Mercurial > repos > public > sbplib

comparison +rv/ResidualViscosity.m @ 1194:bd5383809917 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Comment out filtering option in ResidualViscosity
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Mon, 05 Aug 2019 10:45:08 +0200
parents 921595039ab8
children
comparison
equal deleted inserted replaced
1193:921595039ab8 1194:bd5383809917
30 obj.Mres = obj.Cres*obj.h^2; 30 obj.Mres = obj.Cres*obj.h^2;
31 obj.Mfirst = obj.Cmax*obj.h; 31 obj.Mfirst = obj.Cmax*obj.h;
32 switch postProcess 32 switch postProcess
33 case {'', 'none'} 33 case {'', 'none'}
34 obj.fRes = @(v,dvdt) obj.Mres*abs(dvdt + obj.Df(v))./obj.normalization(v); 34 obj.fRes = @(v,dvdt) obj.Mres*abs(dvdt + obj.Df(v))./obj.normalization(v);
35 case {'filt', 'filter'} 35 % TBD: Keep?
36 % TBD: Keep? 36 % case {'filt', 'filter'}
37 order = 37 % order = 4;
38 F = obj.shapiroFilter(obj.g, order); 38 % F = obj.shapiroFilter(obj.g, order);
39 obj.Mres = F*obj.Mres; 39 % obj.Mres = F*obj.Mres;
40 obj.fRes = @(v,dvdt) obj.Mres*abs(dvdt + obj.Df(v))./obj.normalization(v); 40 % obj.fRes = @(v,dvdt) obj.Mres*abs(dvdt + obj.Df(v))./obj.normalization(v);
41 case {'max', 'maximum neighbors'} 41 case {'max', 'maximum neighbors'}
42 switch g.D() 42 switch g.D()
43 case 1 43 case 1
44 obj.fRes = @(v,dvdt) movmax(obj.Mres*abs(dvdt + obj.Df(v))./obj.normalization(v),3); 44 obj.fRes = @(v,dvdt) movmax(obj.Mres*abs(dvdt + obj.Df(v))./obj.normalization(v),3);
45 case 2 45 case 2