Mercurial > repos > public > sbplib

diff +scheme/Burgers1d.m @ 1196:f6c571d8f22f feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Add option to discretize Burgers1d using quasi-linear form
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 07 Aug 2019 13:28:21 +0200
parents 635386c073b9
children 433c89bf19e0
line wrap: on
line diff
--- a/+scheme/Burgers1d.m	Wed Aug 07 13:27:36 2019 +0200
+++ b/+scheme/Burgers1d.m	Wed Aug 07 13:28:21 2019 +0200
@@ -31,6 +31,8 @@
                 obj.D1 = (ops.Dp + ops.Dm)/2;
                 DissOp = (ops.Dm - ops.Dp)/2;
                 switch pde_form
+                    case 'quasi-linear'
+                        obj.D = @(v) -((spdiag(v)*obj.D1  + fluxSplitting(v)*DissOp)*v);
                     case 'skew-symmetric'
                         obj.D = @(v) -(1/3*obj.D1*(v.*v) + (1/3*spdiag(v)*obj.D1  + fluxSplitting(v)*DissOp)*v);
                     case 'conservative'
@@ -39,6 +41,8 @@
             else 
                 obj.D1 = ops.D1;
                 switch pde_form
+                    case 'quasi-linear'
+                        obj.D = @(v) -(spdiag(v)*obj.D1*v);
                     case 'skew-symmetric'
                         obj.D = @(v) -(1/3*obj.D1*(v.*v) + 1/3*spdiag(v)*obj.D1*v);
                     case 'conservative'