Mercurial > repos > public > sbplib

comparison +scheme/Hypsyst3d.m @ 368:53abf04f5e4e feature/hypsyst

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
Made some small changes.
author Ylva Rydin <ylva.rydin@telia.com>
date Wed, 25 Jan 2017 15:37:12 +0100
parents 9567a9dd220d
children 9d1fc984f40d
comparison
equal deleted inserted replaced
367:05947fc2505c 368:53abf04f5e4e
105 obj.h = [ops_x.h ops_y.h ops_x.h]; 105 obj.h = [ops_x.h ops_y.h ops_x.h];
106 obj.order = order; 106 obj.order = order;
107 107
108 switch operator 108 switch operator
109 case 'upwind' 109 case 'upwind'
110 alphaA = max(eig(A(params,obj.x(end),obj.y(end),obj.z(end)))); 110 alphaA = max(abs(eig(A(params,obj.x(end),obj.y(end),obj.z(end)))));
111 alphaB = max(eig(B(params,obj.x(end),obj.y(end),obj.z(end)))); 111 alphaB = max(abs(eig(B(params,obj.x(end),obj.y(end),obj.z(end)))));
112 alphaC = max(eig(C(params,obj.x(end),obj.y(end),obj.z(end)))); 112 alphaC = max(abs(eig(C(params,obj.x(end),obj.y(end),obj.z(end)))));
113 113
114 Ap = (obj.Aevaluated+alphaA*I_N)/2; 114 Ap = (obj.Aevaluated+alphaA*I_N)/2;
115 Am = (obj.Aevaluated-alphaA*I_N)/2; 115 Am = (obj.Aevaluated-alphaA*I_N)/2;
116 Dpx = kr(I_n, ops_x.Dp, I_y,I_z); 116 Dpx = kr(I_n, ops_x.Dp, I_y,I_z);
117 Dmx = kr(I_n, ops_x.Dm, I_y,I_z); 117 Dmx = kr(I_n, ops_x.Dm, I_y,I_z);