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);
}