Properly compute the Hessian quadratic product in BQPD when using a Hessian operator #300
Description
BQPDEvaluationSpace::compute_hessian_quadratic_product forms the quadratic products using the current sparse Hessian matrix, which is empty if we use a Hessian operator.