BQPDEvaluationSpace::compute_hessian_quadratic_product forms the quadratic products using the current sparse Hessian matrix, which is empty if we use a Hessian operator.