Mercurial > repos > public > sbplib

comparison spdiagsVariablePeriodic.m @ 683:50e77b15d841 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bugfix in spdiagsVariablePeriodic
author Martin Almquist <malmquist@stanford.edu>
date Thu, 08 Feb 2018 16:02:20 -0800
parents 7368affc8f78
children 8c65ef13df89
comparison
equal deleted inserted replaced
681:7368affc8f78 683:50e77b15d841
25 A = A + spdiagVariable(a_corner, corner_diag); 25 A = A + spdiagVariable(a_corner, corner_diag);
26 26
27 % Super-diagonals 27 % Super-diagonals
28 elseif d > 0 28 elseif d > 0
29 a_bulk = a(1:end-d); 29 a_bulk = a(1:end-d);
30 a_corner = a(end-d+1,end); 30 a_corner = a(end-d+1:end);
31 corner_diag = -m + d; 31 corner_diag = -m + d;
32 A = A + spdiagVariable(a_bulk, d); 32 A = A + spdiagVariable(a_bulk, d);
33 A = A + spdiagVariable(a_corner, corner_diag); 33 A = A + spdiagVariable(a_corner, corner_diag);
34 34
35 % Main diagonal 35 % Main diagonal