SquareMatrix< T > tensorproduct(const SquareMatrix< T > &op1, const SquareMatrix< T > &op2)
Form the tensor product between two operators.
Definition operator_tensor.h:25
dg::SparseMatrix< int, T, thrust::host_vector > sandwich(const SquareMatrix< T > &left, const dg::SparseMatrix< int, T, thrust::host_vector > &m, const SquareMatrix< T > &right)
Multiply 1d matrices by diagonal block matrices from left and right.
Definition operator_tensor.h:93