Mercurial > repos > public > sbplib

changeset 822:4808c4bd844e feature/operator_files

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add assemble_boundary_op
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 10 Sep 2018 17:16:22 +0200
parents 95c26000c0ba
children 5c8b1a3bd0e6
files operator_def/assemble_D1.m operator_def/assemble_boundary_op.m
diffstat 2 files changed, 11 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/operator_def/assemble_D1.m	Mon Sep 10 16:57:17 2018 +0200
+++ b/operator_def/assemble_D1.m	Mon Sep 10 17:16:22 2018 +0200
@@ -1,6 +1,6 @@
 function [D1,e_l,e_r] = assemble_D1(stencil,h,m)
   stencil_variables = read_stencil(stencil);
   
-  D1 = assemble_op(stencil_variables.inner_stencil,stencil_variables.boundary_block)/h;
-  
+  D1 = assemble_op(stencil_variables.inner_stencil,stencil_variables.boundary_block,m)/h;
+  [e_l, e_r] = assemble_boundary_op(stencil_variables);
   end
\ No newline at end of file
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/operator_def/assemble_boundary_op.m	Mon Sep 10 17:16:22 2018 +0200
@@ -0,0 +1,9 @@
+function [op_l, op_r] = assemble_boundary_op(op,m)
+l = length(op);
+
+op_l = sparse(m,1);
+op_l(1:l) = op;
+
+op_r = sparse(m,1);
+op_r(m-l+l:m) = flipr(op);
+end