Mercurial > repos > public > sbplib

changeset 1222:55463e3c1e4a rv_diffOp_test

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge bug fix from feature/rv
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 05 Mar 2019 10:59:39 +0100
parents 0c906f7ab8bf (current diff) 635386c073b9 (diff)
children 9fddc8749445
files
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
diff -r 0c906f7ab8bf -r 55463e3c1e4a +scheme/Burgers1d.m
--- a/+scheme/Burgers1d.m	Tue Mar 05 10:56:29 2019 +0100
+++ b/+scheme/Burgers1d.m	Tue Mar 05 10:59:39 2019 +0100
@@ -32,17 +32,17 @@
                 DissOp = (ops.Dm - ops.Dp)/2;
                 switch pde_form
                     case 'skew-symmetric'
-                        obj.D = @(v) -(1/3*obj.D1*v.*v + (1/3*spdiag(v)*obj.D1  + fluxSplitting(v)*DissOp)*v);
+                        obj.D = @(v) -(1/3*obj.D1*(v.*v) + (1/3*spdiag(v)*obj.D1  + fluxSplitting(v)*DissOp)*v);
                     case 'conservative'
-                        obj.D = @(v) -(1/2*obj.D1*v.*v + fluxSplitting(v)*DissOp*v);
+                        obj.D = @(v) -(1/2*obj.D1*(v.*v) + fluxSplitting(v)*DissOp*v);
                 end
             else 
                 obj.D1 = ops.D1;
                 switch pde_form
                     case 'skew-symmetric'
-                        obj.D = @(v) -(1/3*obj.D1*v.*v + 1/3*spdiag(v)*obj.D1*v);
+                        obj.D = @(v) -(1/3*obj.D1*(v.*v) + 1/3*spdiag(v)*obj.D1*v);
                     case 'conservative'
-                        obj.D = @(v) -1/2*obj.D1*v.*v;
+                        obj.D = @(v) -1/2*obj.D1*(v.*v);
                 end
             end
             obj.grid = g;