Mercurial > repos > public > sbplib

comparison +scheme/Hypsyst3dCurve.m @ 512:4ef2d2a493f1 feature/quantumTriangles

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
add H_xyz to hypsyst to compute H-norm
author Ylva Rydin <ylva.rydin@telia.com>
date Mon, 26 Jun 2017 19:23:19 +0200
parents feebfca90080
children
comparison
equal deleted inserted replaced
511:57f3493f851b 512:4ef2d2a493f1
214 end 214 end
215 215
216 obj.Hxii = kr(I_n, ops_xi.HI, I_eta,I_zeta); 216 obj.Hxii = kr(I_n, ops_xi.HI, I_eta,I_zeta);
217 obj.Hetai = kr(I_n, I_xi, ops_eta.HI,I_zeta); 217 obj.Hetai = kr(I_n, I_xi, ops_eta.HI,I_zeta);
218 obj.Hzetai = kr(I_n, I_xi,I_eta, ops_zeta.HI); 218 obj.Hzetai = kr(I_n, I_xi,I_eta, ops_zeta.HI);
219 Hi = obj.Hxii*obj.Hetai*obj.Hzetai;
220 obj.H = inv(Hi);
219 221
220 obj.index_w = (kr(ops_xi.e_l, O_eta,O_zeta)==1); 222 obj.index_w = (kr(ops_xi.e_l, O_eta,O_zeta)==1);
221 obj.index_e = (kr(ops_xi.e_r, O_eta,O_zeta)==1); 223 obj.index_e = (kr(ops_xi.e_r, O_eta,O_zeta)==1);
222 obj.index_s = (kr(O_xi, ops_eta.e_l,O_zeta)==1); 224 obj.index_s = (kr(O_xi, ops_eta.e_l,O_zeta)==1);
223 obj.index_n = (kr(O_xi, ops_eta.e_r,O_zeta)==1); 225 obj.index_n = (kr(O_xi, ops_eta.e_r,O_zeta)==1);