Mercurial > repos > public > sbplib

diff +rv/+diffops/constructSymmetricD2.m @ 1225:68ee061639a1 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Make sure matrices are sparse.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 06 Nov 2019 14:52:07 +0100
parents 433c89bf19e0
children
line wrap: on
line diff
--- a/+rv/+diffops/constructSymmetricD2.m	Wed Nov 06 14:51:33 2019 +0100
+++ b/+rv/+diffops/constructSymmetricD2.m	Wed Nov 06 14:52:07 2019 +0100
@@ -12,12 +12,12 @@
         case 1
             e_r = ops{1}.e_r;
             e_l = ops{1}.e_l;
-            Hi = ops{1}.HI;
-            B = e_r*e_r' - e_l*e_l';
+            Hi = ops{1}.HI;            
             if isequal(opSet,@sbp.D1Upwind)
-                Dm = ops{1}.Dm;
-                Dp = ops{1}.Dp;
-                M =  Dm - Hi*B;
+                B = e_r*e_r' - e_l*e_l';
+                Dm = sparse(ops{1}.Dm);
+                Dp = sparse(ops{1}.Dp);
+                M =  sparse(Dm - Hi*B);
                 D2 = @(Viscosity) M*Viscosity*Dp;
             else
                 % TODO: Fix Viscosity not being vector