Mercurial > repos > public > sbplib

diff +rv/ResidualViscosity.m @ 1164:fc2631ba4da5 feature/rv

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Remove unnecessary variable
author Vidar Stiernström <vidar.stiernstrom@it.uu.se>
date Thu, 27 Jun 2019 17:13:28 +0200
parents 856bd6291d17
children beecb580c5bf
line wrap: on
line diff
--- a/+rv/ResidualViscosity.m	Thu Jun 27 11:04:34 2019 +0200
+++ b/+rv/ResidualViscosity.m	Thu Jun 27 17:13:28 2019 +0200
@@ -43,8 +43,7 @@
         end
 
         function viscosity = evaluateViscosity(obj, v, dvdt)
-            firstOrderViscosity = obj.Mfirst*abs(obj.waveSpeed(v));
-            viscosity = min(firstOrderViscosity, obj.fRes(v,dvdt));
+            viscosity = min(obj.Mfirst*abs(obj.waveSpeed(v)), obj.fRes(v,dvdt));
         end
 
         function [viscosity, Df, firstOrderViscosity, residualViscosity] = evaluate(obj, v, dvdt)