Mercurial > repos > public > sbplib

view operator_def/assemble_op.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 d1e5143d67ed
children
line wrap: on
line source

function [D] = assemble_op(inner_stencil,boundary_block,m,sign)  
  [height_b, width_b] = size(boundary_block);
  width_i = length(inner_stencil);
  
  max_offset = floor(width_i/2);
  diags = -max_offset:max_offset;
   
  D = stripeMatrix(inner_stencil, diags, m);
  
  D(1:height_b,1:width_b) = boundary_block;
  D(m-height_b+1:m,m-width_b+1:m) = rot90( sign*boundary_block ,2);
end