Mercurial > repos > public > sbplib

view stripeMatrixPeriodic.m @ 741:5a9acf282b34 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add scheme Heat2Dcurvilinear. Neumann and Dirichlet seem to work. Only tested for stretched Cartesian grids though.
author Martin Almquist <malmquist@stanford.edu>
date Wed, 09 May 2018 19:29:12 -0700
parents 7368affc8f78
children bbf303c1f0cf
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 = spdiagsVariablePeriodic(D,diags);
end