Mercurial > repos > public > sbplib

comparison +grid/Ti3D.m @ 350:5d5652fe826a feature/hypsyst

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
A commit before I try resolving the performance issues
author Ylva Rydin <ylva.rydin@telia.com>
date Wed, 02 Nov 2016 00:02:01 +0100
parents cd6a29ab3746
children 69b078cf8072
comparison
equal deleted inserted replaced
349:cd6a29ab3746 350:5d5652fe826a
23 ZETA=ZETA'; 23 ZETA=ZETA';
24 24
25 one=0*ETA+1; 25 one=0*ETA+1;
26 zero=0*ETA; 26 zero=0*ETA;
27 27
28 Sw = gw((1-ETA),(1-ZETA)); 28 Sw = gw(ETA,(1-ZETA));
29 Se = ge(ETA,ZETA); 29 Se = ge((1-ETA),(1-ZETA));
30 Ss = gs(XI,(1-ZETA)); 30 Ss = gs(XI,ZETA);
31 Sn = gn((1-XI),ZETA); 31 Sn = gn((1-XI),(1-ZETA));
32 Sb = gb(XI,ETA); 32 Sb = gb((1-XI),ETA);
33 St = gt((1-XI),(1-ETA)); 33 St = gt(XI,ETA);
34 34
35 Ewt = gw(1-ETA,zero); 35 Ewt = gw(ETA,zero);
36 Ewb = gw(1-ETA,one); 36 Ewb = gw(ETA,one);
37 Ews = gw(one,1-ZETA); 37 Ews = gw(zero,1-ZETA);
38 Ewn = gw(zero,1-ZETA); 38 Ewn = gw(one,1-ZETA);
39 Eet = ge(ETA,one); 39 Eet = ge(1-ETA,zero);
40 Eeb = ge(ETA,zero); 40 Eeb = ge(1-ETA,one);
41 Ees = ge(0*one,ZETA); 41 Ees = ge(one,1-ZETA);
42 Een = ge(one,ZETA); 42 Een = ge(zero,1-ZETA);
43 Enb = gn(1-XI,zero); 43 Enb = gn(1-XI,one);
44 Ent = gn(1-XI,one); 44 Ent = gn(1-XI,zero);
45 Est = gs(XI,zero); 45 Est = gs(XI,one);
46 Esb = gs(XI,one); 46 Esb = gs(XI,zero);
47 47
48 Cwbs = gw(one,one); 48 Cwbs = gw(zero,one);
49 Cwbn = gw(zero,one); 49 Cwbn = gw(one,one);
50 Cwts = gw(one,zero); 50 Cwts = gw(zero,zero);
51 Cwtn = gw(zero,zero); 51 Cwtn = gw(one,zero);
52 Cebs = ge(zero,zero); 52 Cebs = ge(one,one);
53 Cebn = ge(one,zero); 53 Cebn = ge(zero,one);
54 Cets = ge(zero,one); 54 Cets = ge(one,zero);
55 Cetn = ge(one,one); 55 Cetn = ge(zero,zero);
56 56
57 57
58 X1 = (1-XI).*Sw(1,:,:) + XI.*Se(1,:,:); 58 X1 = (1-XI).*Sw(1,:,:) + XI.*Se(1,:,:);
59 X2 = (1-ETA).*Ss(1,:,:) + ETA.*Sn(1,:,:); 59 X2 = (1-ETA).*Ss(1,:,:) + ETA.*Sn(1,:,:);
60 X3 = (1-ZETA).*Sb(1,:,:) + ZETA.*St(1,:,:); 60 X3 = (1-ZETA).*Sb(1,:,:) + ZETA.*St(1,:,:);
102 end 102 end
103 103
104 obj.V = @V_fun; 104 obj.V = @V_fun;
105 end 105 end
106 106
107 107 %Should be rewritten so that the input is xi eta zeta
108 function [X,Y,Z] = map(obj,XI,ETA,ZETA) 108 function [X,Y,Z] = map(obj,XI,ETA,ZETA)
109 109
110 V = obj.V; 110 V = obj.V;
111 111
112 p = V(XI,ETA,ZETA); 112 p = V(XI,ETA,ZETA);
245 % grid.place_label(pw,'w'); 245 % grid.place_label(pw,'w');
246 % grid.place_label(pe,'e'); 246 % grid.place_label(pe,'e');
247 % grid.place_label(ps,'s'); 247 % grid.place_label(ps,'s');
248 % grid.place_label(pn,'n'); 248 % grid.place_label(pn,'n');
249 % end 249 % end
250 % end 250 % end
251 end 251 end
252 end 252 end