Mercurial > repos > public > sbplib

changeset 799:8c65ef13df89 feature/poroelastic

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Bort med martins hemska whitespace
author Jonatan Werpers <jonatan@werpers.com>
date Thu, 26 Jul 2018 17:51:48 -0700
parents 5cf9fdf4c98f
children 87ea9cac3287
files spdiagsVariablePeriodic.m
diffstat 1 files changed, 5 insertions(+), 5 deletions(-) [+]
line wrap: on
line diff
diff -r 5cf9fdf4c98f -r 8c65ef13df89 spdiagsVariablePeriodic.m
--- a/spdiagsVariablePeriodic.m	Thu Jul 26 10:53:05 2018 -0700
+++ b/spdiagsVariablePeriodic.m	Thu Jul 26 17:51:48 2018 -0700
@@ -2,17 +2,17 @@
     % Creates an m x m periodic discretization matrix.
     % vals - m x ndiags matrix of values
     % diags - 1 x ndiags vector of the 'center diagonals' that vals end up on
-    % vals that are not on main diagonal are going to spill over to 
+    % vals that are not on main diagonal are going to spill over to
     % off-diagonal corners.
 
     default_arg('diags',0);
 
-    [m, ~] = size(vals); 
+    [m, ~] = size(vals);
 
     A = sparse(m,m);
 
     for i = 1:length(diags)
-        
+
         d = diags(i);
         a = vals(:,i);
 
@@ -21,7 +21,7 @@
             a_bulk = a(1+abs(d):end);
             a_corner = a(1:1+abs(d)-1);
             corner_diag = m-abs(d);
-            A = A + spdiagVariable(a_bulk, d); 
+            A = A + spdiagVariable(a_bulk, d);
             A = A + spdiagVariable(a_corner, corner_diag);
 
         % Super-diagonals
@@ -29,7 +29,7 @@
             a_bulk = a(1:end-d);
             a_corner = a(end-d+1:end);
             corner_diag = -m + d;
-            A = A + spdiagVariable(a_bulk, d); 
+            A = A + spdiagVariable(a_bulk, d);
             A = A + spdiagVariable(a_corner, corner_diag);
 
         % Main diagonal