Mercurial > repos > public > sbplib

view stripeMatrixPeriodic.m @ 1275:8ff3a95ad7cc feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Change metric computation in CurvilinearUpwind so that one can use odd orders for the diffOp
author Martin Almquist <malmquist@stanford.edu>
date Wed, 03 Jun 2020 10:29:53 -0700
parents bbf303c1f0cf
children
line wrap: on
line source

% Creates a periodic discretization matrix of size n x n
%  with the values of val on the diagonals diag.
%   A = stripeMatrix(val,diags,n)
function A = stripeMatrixPeriodic(val,diags,n)

    D = ones(n,1)*val;
    A = spdiagsPeriodic(D,diags);
end