Mercurial > repos > public > sbplib

comparison +rv/+diffops/constructSymmetricD2.m @ 1197:433c89bf19e0 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Merge with default
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Wed, 07 Aug 2019 15:23:42 +0200
parents 65a577db5ca0
children 68ee061639a1
comparison
equal deleted inserted replaced
1196:f6c571d8f22f 1197:433c89bf19e0
18 Dm = ops{1}.Dm; 18 Dm = ops{1}.Dm;
19 Dp = ops{1}.Dp; 19 Dp = ops{1}.Dp;
20 M = Dm - Hi*B; 20 M = Dm - Hi*B;
21 D2 = @(Viscosity) M*Viscosity*Dp; 21 D2 = @(Viscosity) M*Viscosity*Dp;
22 else 22 else
23 % TODO: Fix closure for D2Variable
24 % TODO: Fix Viscosity not being vector 23 % TODO: Fix Viscosity not being vector
25 D2 = @(Viscosity)ops{1}.D2(diag(Viscosity)); 24 d1_r = ops{1}.d1_r';
25 d1_l = ops{1}.d1_l';
26 D2 = @(Viscosity)ops{1}.D2(diag(Viscosity)) + Hi*(Viscosity(1,1)*e_l*d1_l - e_r*Viscosity(end,end)*d1_r);
26 end 27 end
27 case 2 28 case 2
28 % TODO: 29 % TODO:
29 % Currently only implemented for upwind operators. 30 % Currently only implemented for upwind operators.
30 % Remove this part once the time-dependent D2 operator is implemented for other opSets 31 % Remove this part once the time-dependent D2 operator is implemented for other opSets