Mercurial > repos > public > sbplib

comparison +scheme/Burgers1D.m @ 831:d0934d1143b7 feature/burgers1d

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Fix bug in initialization of differential operators
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 11 Sep 2018 13:24:08 +0200
parents fae41958af4f
children 5573913a0949
comparison
equal deleted inserted replaced
815:fae41958af4f 831:d0934d1143b7
31 error('Other operator types not yet supported', operator_type); 31 error('Other operator types not yet supported', operator_type);
32 end 32 end
33 33
34 switch pde_form 34 switch pde_form
35 case 'skew-symmetric' 35 case 'skew-symmetric'
36 D = @(v, viscosity) -1/3*v.*D1*v - 1/3*D1*v.^2 + D2(obj.params.eps + viscosity)*v; 36 D = @(v, viscosity) -1/3*v.*D1*v - 1/3*D1*v.^2 + D2(params.eps + viscosity)*v;
37 case 'conservative' 37 case 'conservative'
38 D = @(v, viscosity) -1/2*D1*v.^2 + D2(obj.params.eps + viscosity)*v; 38 D = @(v, viscosity) -1/2*D1*v.^2 + D2(params.eps + viscosity)*v;
39 end 39 end
40 40
41 obj.m = m; 41 obj.m = m;
42 obj.h = h; 42 obj.h = h;
43 obj.order = order; 43 obj.order = order;