Mercurial > repos > public > sbplib

diff +scheme/Burgers1d.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/+scheme/Burgers1d.m	Wed Nov 06 14:51:33 2019 +0100
+++ b/+scheme/Burgers1d.m	Wed Nov 06 14:52:07 2019 +0100
@@ -28,8 +28,8 @@
             ops = opSet(m, xlim, order);
 
             if (isequal(opSet, @sbp.D1Upwind))
-                obj.D1 = (ops.Dp + ops.Dm)/2;
-                DissOp = (ops.Dm - ops.Dp)/2;
+                obj.D1 = sparse((ops.Dp + ops.Dm)/2);
+                DissOp = sparse((ops.Dm - ops.Dp)/2);
                 switch pde_form
                     case 'quasi-linear'
                         obj.D = @(v) -((spdiag(v)*obj.D1  + fluxSplitting(v)*DissOp)*v);
@@ -39,7 +39,7 @@
                         obj.D = @(v) -(1/2*obj.D1*(v.*v) + fluxSplitting(v)*DissOp*v);
                 end
             else 
-                obj.D1 = ops.D1;
+                obj.D1 = sparse(ops.D1);
                 switch pde_form
                     case 'quasi-linear'
                         obj.D = @(v) -(spdiag(v)*obj.D1*v);