Mercurial > repos > public > sbplib

comparison +rv/ResidualViscosity.m @ 1161:856bd6291d17 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Use same order for both time derivative and flux term when computing the residual
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Tue, 25 Jun 2019 17:00:27 +0200
parents 1ad7da049b50
children fc2631ba4da5
comparison
equal deleted inserted replaced
1160:76e3bb7836cf 1161:856bd6291d17
16 methods 16 methods
17 % TBD: Decide on how to treat waveSpeed. It would be nice to just pass a constant value without 17 % TBD: Decide on how to treat waveSpeed. It would be nice to just pass a constant value without
18 % wrapping it in a function. 18 % wrapping it in a function.
19 function obj = ResidualViscosity(grid, Df, waveSpeed, Cmax, Cres, h, normalization, postProcess) 19 function obj = ResidualViscosity(grid, Df, waveSpeed, Cmax, Cres, h, normalization, postProcess)
20 default_arg('normalization',@(v)abs(obj.minmaxDiffNeighborhood1d(v)-norm(v-mean(v),inf))); 20 default_arg('normalization',@(v)abs(obj.minmaxDiffNeighborhood1d(v)-norm(v-mean(v),inf)));
21 default_arg('postProcess','none'); 21 default_arg('postProcess','maximum neighbors');
22 obj.Df = Df; 22 obj.Df = Df;
23 obj.waveSpeed = waveSpeed; 23 obj.waveSpeed = waveSpeed;
24 obj.h = h; 24 obj.h = h;
25 obj.Cmax = Cmax; 25 obj.Cmax = Cmax;
26 26