Mercurial > repos > public > sbplib

annotate 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
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
830
e15a667ffde2 Add quadrature stencils and operator assembly. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
1 function [H, Hinv] = assemble_H(quadrature, h, m)
e15a667ffde2 Add quadrature stencils and operator assembly. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
2 quadrature_variables = read_stencil(quadrature);
e15a667ffde2 Add quadrature stencils and operator assembly. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
3 %TBD: Do we want to reuse assempbe_op here?
e15a667ffde2 Add quadrature stencils and operator assembly. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
4 H = assemble_op(quadrature_variables.inner_stencil, quadrature_variables.closure, m, 1)*h;
e15a667ffde2 Add quadrature stencils and operator assembly. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
5 Hinv = inv(H);
e15a667ffde2 Add quadrature stencils and operator assembly. Add tests
Vidar Stiernström <vidar.stiernstrom@it.uu.se>
parents:
diff changeset
6 end