Mercurial > repos > public > sbplib

annotate stripeMatrixPeriodic.m @ 684:fcf004066ea9 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add property lim to grid/Cartesian to make periodic discr easier.
author Martin Almquist <malmquist@stanford.edu>
date Thu, 08 Feb 2018 16:40:48 -0800
parents 7368affc8f78
children bbf303c1f0cf
Ignore whitespace changes - Everywhere: Within whitespace: At end of lines:
rev   line source
681
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
1 % Creates a periodic discretization matrix of size n x n
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
2 % with the values of val on the diagonals diag.
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
3 % A = stripeMatrix(val,diags,n)
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
4 function A = stripeMatrixPeriodic(val,diags,n)
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
5
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
6 D = ones(n,1)*val;
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
7 A = spdiagsVariablePeriodic(D,diags);
7368affc8f78 Add D2 variable periodic for second order.
Martin Almquist <malmquist@stanford.edu>
parents:
diff changeset
8 end