changeset 824:ed8d34894589 feature/operator_files

fix bug in assemble_D1
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 10 Sep 2018 17:31:42 +0200
parents 5c8b1a3bd0e6
children 32c360bb480e
files operator_def/assemble_D1.m operator_def/assemble_boundary_op.m
diffstat 2 files changed, 2 insertions(+), 2 deletions(-) [+]
line wrap: on
line diff
--- a/operator_def/assemble_D1.m	Mon Sep 10 17:28:15 2018 +0200
+++ b/operator_def/assemble_D1.m	Mon Sep 10 17:31:42 2018 +0200
@@ -2,5 +2,5 @@
   stencil_variables = read_stencil(stencil);
   
   D1 = assemble_op(stencil_variables.inner_stencil,stencil_variables.boundary_stencils,m)/h;
-  [e_l, e_r] = assemble_boundary_op(stencil_variables);
+  [e_l, e_r] = assemble_boundary_op(stencil_variables.e,m);
   end
\ No newline at end of file
--- a/operator_def/assemble_boundary_op.m	Mon Sep 10 17:28:15 2018 +0200
+++ b/operator_def/assemble_boundary_op.m	Mon Sep 10 17:31:42 2018 +0200
@@ -5,5 +5,5 @@
 op_l(1:l) = op;
 
 op_r = sparse(m,1);
-op_r(m-l+l:m) = flipr(op);
+op_r(m-l+1:m) = fliplr(op);
 end