Mercurial > repos > public > sbplib
diff +scheme/Hypsyst3d.m @ 368:53abf04f5e4e feature/hypsyst
Made some small changes.
| author | Ylva Rydin <ylva.rydin@telia.com> |
|---|---|
| date | Wed, 25 Jan 2017 15:37:12 +0100 |
| parents | 9567a9dd220d |
| children | 9d1fc984f40d |
line wrap: on
line diff
--- a/+scheme/Hypsyst3d.m Tue Dec 27 09:11:29 2016 +0100 +++ b/+scheme/Hypsyst3d.m Wed Jan 25 15:37:12 2017 +0100 @@ -107,9 +107,9 @@ switch operator case 'upwind' - alphaA = max(eig(A(params,obj.x(end),obj.y(end),obj.z(end)))); - alphaB = max(eig(B(params,obj.x(end),obj.y(end),obj.z(end)))); - alphaC = max(eig(C(params,obj.x(end),obj.y(end),obj.z(end)))); + alphaA = max(abs(eig(A(params,obj.x(end),obj.y(end),obj.z(end))))); + alphaB = max(abs(eig(B(params,obj.x(end),obj.y(end),obj.z(end))))); + alphaC = max(abs(eig(C(params,obj.x(end),obj.y(end),obj.z(end))))); Ap = (obj.Aevaluated+alphaA*I_N)/2; Am = (obj.Aevaluated-alphaA*I_N)/2;