Mercurial > repos > public > sbplib

changeset 1153:635386c073b9 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in Burgers schemes resulting in incorrect wave speed.
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 05 Mar 2019 10:57:26 +0100
parents 010bb2677230
children 3108963cc42c 55463e3c1e4a
files +scheme/Burgers1d.m
diffstat 1 files changed, 4 insertions(+), 4 deletions(-) [+]
line wrap: on
line diff
--- a/+scheme/Burgers1d.m	Tue Mar 05 10:53:34 2019 +0100
+++ b/+scheme/Burgers1d.m	Tue Mar 05 10:57:26 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;