void finish_tensor_trace( SpinBlock& b, SpinBlock* sysdot, SparseMatrix& sysdot_op, SparseMatrix& op, std::string& build_pattern ) { // Build and store new operator assert( ! op.get_built() ); op.set_built() = true; op.set_build_pattern() = build_pattern; //FIXME magic number 2 op.set_deltaQuantum(1, op.get_quantum_ladder().at( build_pattern ).at(2) ); op.allocate(b.get_stateInfo()); SpinAdapted::operatorfunctions::TensorTrace(sysdot, sysdot_op, &b, &(b.get_stateInfo()), op); }
void finish_tensor_product( SpinBlock& b, SpinBlock* sysdot, const SparseMatrix& sysdot_op1, const SparseMatrix& sysdot_op2, SparseMatrix& op, bool include_parity, std::string& build_pattern ) { // Build and store new operator assert( ! op.get_built() ); op.set_built() = true; op.set_build_pattern() = build_pattern; op.set_deltaQuantum(1, op.get_quantum_ladder().at( build_pattern ).at(2) ); op.allocate(b.get_stateInfo()); // Do tensor product double parity = 1.0; if ( include_parity ) parity = getCommuteParity( sysdot_op1.get_deltaQuantum(0), sysdot_op2.get_deltaQuantum(0), op.get_deltaQuantum(0) ); SpinAdapted::operatorfunctions::TensorProduct(sysdot, sysdot_op1, sysdot_op2, &b, &(b.get_stateInfo()), op, parity); }