Mercurial > repos > public > sbplib
view spdiagVariable.m @ 724:a9e8c9d71307 feature/poroelastic
Modify penalty for Neumann in Head2d so that data is for normal derivative and not u_x or u_y.
author | Martin Almquist <malmquist@stanford.edu> |
---|---|
date | Wed, 14 Mar 2018 13:04:25 -0700 |
parents | 7368affc8f78 |
children |
line wrap: on
line source
function A = spdiagVariable(a,i) default_arg('i',0); if isrow(a) a = a'; end n = length(a)+abs(i); if i > 0 a = [sparse(i,1); a]; elseif i < 0 a = [a; sparse(abs(i),1)]; end A = spdiags(a,i,n,n); end