OperatorSpace_test.cpp File Reference

Namespaces
namespace	Eigen

Functions
template<typename charT , typename traits >
std::basic_ostream< charT, traits > &	Eigen::operator<< (std::basic_ostream< charT, traits > &strm, const std::complex< double > &c)
Eigen::IOFormat	OctaveFmt (Eigen::StreamPrecision, 0, ", ", ";\n", "", "", "[", "]")
int	main (void)

Function Documentation

main()

int main

(

void

)

               {
    int L = 7, dimLoc = 2;
    HilbertSpace<> baseSpace(dimLoc);
    std::cout << baseSpace.dim()     << std::endl;
    std::cout << &baseSpace << std::endl;
 
    {
        std::cout << "Testing class OperatorSpace." << std::endl;
        OperatorSpace<> opSpace(baseSpace);
        std::cout << opSpace.dim()     << "\n"
                << &opSpace.baseSpace()  << std::endl;
        std::cout << &baseSpace << std::endl;
 
        long opNum = 3;
        Eigen::SparseMatrix<std::complex<double>> op(opSpace.baseSpace().dim(), opSpace.baseSpace().dim());
        for(int j = 0;j < opSpace.baseSpace().dim(); ++j) {
            std::pair<long, std::complex<double> > p = opSpace.action(opNum,j);
            if( std::abs(p.second) != 0 ) op.insert(p.first, j) = p.second;
        }
        std::cout << op.pruned(1,0.5) << std::endl;
    }
 
    {
        std::cout << "Testing class ManyBodyOperatorSpace." << std::endl;
        int L = 4;
        ManyBodySpinSpace mbHSpace(L, baseSpace);
        std::cout << std::endl;
        ManyBodyOperatorSpace<> mbOpSpace(mbHSpace);
        std::cout << std::endl;
        std::cout << mbOpSpace.dim() << "\n"
                << mbOpSpace.baseSpace().dim()    << " " << mbHSpace.dim() << "\n"
                << mbOpSpace.dimLoc() << "\n"
                << mbOpSpace.baseSpace().dimLoc() << " " << mbHSpace.dimLoc() << "\n"
                << mbOpSpace.sysSize() << " "
                << mbOpSpace.baseSpace().sysSize() << " " << mbHSpace.sysSize() << "\n" << std::endl;
 
        std::cout << "&baseSpace                        = " << &baseSpace << "\n"
                  << "&mbHSpace                         = " << &mbHSpace << "\n"
                  << "&mbHSpace.locSpace()              = " << &mbHSpace.locSpace() << "\n"
                  << "&mbOpSpace.baseSpace()            = " << &mbOpSpace.baseSpace() << "\n"
                  << "&mbOpSpace.locSpace()             = " << &mbOpSpace.locSpace() << "\n"
                  << "&mbOpSpace.locSpace().baseSpace() = " << &mbOpSpace.locSpace().baseSpace() << "\n"
                  << "&mbOpSpace                        = " << &mbOpSpace << "\n";
 
        std::cout << mbOpSpace.ordinalToConfig(1) << std::endl;
    }
 
    return EXIT_SUCCESS;
}

OctaveFmt()

Eigen::IOFormat OctaveFmt	(	Eigen::StreamPrecision	,
		0	,
		"	,
		"	,
		";\n"	,
		""	,
		""	,
		""	[", "]
	)