function [H, Hinv] = assemble_H(quadrature, h, m)
    quadrature_variables = read_stencil(quadrature);
    %TBD: Do we want to reuse assempbe_op here?
    H = assemble_op(quadrature_variables.inner_stencil, quadrature_variables.closure, m, 1)*h;
    Hinv = inv(H);
end