Mercurial > repos > public > sbplib

comparison +sbp/+implementations/d2_sparsity_pattern_inds.m @ 1344:b4e5e45bd239 feature/D2_boundary_opt

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove round off zeros from D2Nonequidistant operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Sat, 15 Oct 2022 15:48:20 +0200
parents
children
comparison
equal deleted inserted replaced
1343:09a5783a3d37 1344:b4e5e45bd239
1 function inds = d2_sparsity_pattern_inds(m, order, BP, interior_offset, closure_offset)
2 sparsity_pattern = sparse(m,m);
3 inner_stencil_inds = -order/2-interior_offset:order/2+interior_offset;
4 inner_stencil_pattern = ones(m,length(inner_stencil_inds));
5 sparsity_pattern = spdiags(inner_stencil_pattern,inner_stencil_inds,sparsity_pattern);
6 sparsity_pattern(1:BP,1:BP+closure_offset) = 1;
7 sparsity_pattern(end-BP+1:end,end-(BP+closure_offset)+1:end) = 1;
8 for k = 1:closure_offset
9 sparsity_pattern(BP+k,1:BP+closure_offset+k) = 1;
10 sparsity_pattern(end-(BP+k)+1,end-(BP+k+closure_offset)+1:end) = 1;
11 end
12 inds = find(sparsity_pattern);
13 end