Mercurial > repos > public > sbplib

diff 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 diff
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/operator_def/assemble_H.m	Mon Sep 10 19:27:27 2018 +0200
@@ -0,0 +1,6 @@
+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
\ No newline at end of file