Mercurial > repos > public > sbplib

view operator_def/assemble_H.m @ 830:e15a667ffde2 feature/operator_files

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add quadrature stencils and operator assembly. Add tests
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 10 Sep 2018 19:27:27 +0200
parents
children
line wrap: on
line source

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