Lines Matching refs:al
10 real*8 al(npro,nenl,nsd), local
40 call get_a_not_hex(xl,al) ! get mapping poly. coeff.
44 detaij(:) = -al(:,2,1)*al(:,3,2)*al(:,4,3) +
45 & al(:,2,1)*al(:,4,2)*al(:,3,3) + al(:,2,2)*
46 & al(:,3,1)*al(:,4,3) - al(:,2,2)*al(:,4,1)*
47 & al(:,3,3) - al(:,2,3)*al(:,3,1)*al(:,4,2)+
48 & al(:,2,3)*al(:,4,1)*al(:,3,2)
52 zi0(:,1) = detaij(:)*((al(:,4,2)*al(:,3,3)
53 & - al(:,3,2)*al(:,4,3))*(xts1-al(:,1,1)) +
54 & (al(:,3,1)*al(:,4,3)
55 & - al(:,4,1)*al(:,3,3))*(xts2-al(:,1,2)) +
56 & (al(:,4,1)*al(:,3,2)
57 & - al(:,3,1)*al(:,4,2))*(xts3-al(:,1,3)))
60 zi0(:,2) = detaij(:)*((al(:,2,2)*al(:,4,3)
61 & - al(:,4,2)*al(:,2,3))*(xts1-al(:,1,1)) +
62 & (al(:,4,1)*al(:,2,3)
63 & - al(:,2,1)*al(:,4,3))*(xts2-al(:,1,2)) +
64 & (al(:,2,1)*al(:,4,2)
65 & - al(:,4,1)*al(:,2,2))*(xts3-al(:,1,3)))
67 zi0(:,3) = detaij(:)*((al(:,3,2)*al(:,2,3)
68 & - al(:,2,2)*al(:,3,3))*(xts1-al(:,1,1)) +
69 & (al(:,2,1)*al(:,3,3)
70 & - al(:,3,1)*al(:,2,3))*(xts2-al(:,1,2)) +
71 & (al(:,3,1)*al(:,2,2)
72 & - al(:,2,1)*al(:,3,2))*(xts3-al(:,1,3)))
81 m11(:)=al(:,2,1)+al(:,5,1)*zi0(:,2)+al(:,7,1)*zi0(:,3)
82 & +al(:,8,1)*zi0(:,2)*zi0(:,3)
83 m12(:)=al(:,3,1)+al(:,5,1)*zi0(:,1)+al(:,6,1)*zi0(:,3)
84 & +al(:,8,1)*zi0(:,1)*zi0(:,3)
85 m13(:)=al(:,4,1)+al(:,6,1)*zi0(:,2)+al(:,7,1)*zi0(:,1)
86 & +al(:,8,1)*zi0(:,1)*zi0(:,2)
88 m21(:)=al(:,2,2)+al(:,5,2)*zi0(:,2)+al(:,7,2)*zi0(:,3)
89 & +al(:,8,2)*zi0(:,2)*zi0(:,3)
90 m22(:)=al(:,3,2)+al(:,5,2)*zi0(:,1)+al(:,6,2)*zi0(:,3)
91 & +al(:,8,2)*zi0(:,1)*zi0(:,3)
92 m23(:)=al(:,4,2)+al(:,6,2)*zi0(:,2)+al(:,7,2)*zi0(:,1)
93 & +al(:,8,2)*zi0(:,1)*zi0(:,2)
95 m31(:)=al(:,2,3)+al(:,5,3)*zi0(:,2)+al(:,7,3)*zi0(:,3)
96 & +al(:,8,3)*zi0(:,2)*zi0(:,3)
97 m32(:)=al(:,3,3)+al(:,5,3)*zi0(:,1)+al(:,6,3)*zi0(:,3)
98 & +al(:,8,3)*zi0(:,1)*zi0(:,3)
99 m33(:)=al(:,4,3)+al(:,6,3)*zi0(:,2)+al(:,7,3)*zi0(:,1)
100 & +al(:,8,3)*zi0(:,1)*zi0(:,2)
105 r1(:)=al(:,1,1)+al(:,2,1)*zi0(:,1)+al(:,3,1)*zi0(:,2)+
106 & al(:,4,1)*zi0(:,3)+al(:,5,1)*zi0(:,1)*zi0(:,2)+
107 & al(:,6,1)*zi0(:,2)*zi0(:,3)+al(:,7,1)*
108 & zi0(:,1)*zi0(:,3)+al(:,8,1)*zi0(:,1)*
111 r2(:)=al(:,1,2)+al(:,2,2)*zi0(:,1)+al(:,3,2)*zi0(:,2)+
112 & al(:,4,2)*zi0(:,3)+al(:,5,2)*zi0(:,1)*zi0(:,2)+
113 & al(:,6,2)*zi0(:,2)*zi0(:,3)+al(:,7,2)*
114 & zi0(:,1)*zi0(:,3)+al(:,8,2)*zi0(:,1)*
117 r3(:)=al(:,1,3)+al(:,2,3)*zi0(:,1)+al(:,3,3)*zi0(:,2)+
118 & al(:,4,3)*zi0(:,3)+al(:,5,3)*zi0(:,1)*zi0(:,2)+
119 & al(:,6,3)*zi0(:,2)*zi0(:,3)+al(:,7,3)*
120 & zi0(:,1)*zi0(:,3)+al(:,8,3)*zi0(:,1)*
159 call get_a_not_tet(xl,al)
165 detaij(e) = al(e,2,1)*(-al(e,3,2)*al(e,4,3) +
166 & al(e,4,2)*al(e,3,3)) + al(e,2,2)*
167 & (al(e,3,1)*al(e,4,3) - al(e,4,1)*
168 & al(e,3,3)) + al(e,2,3)*(-al(e,3,1)*al(e,4,2)+
169 & al(e,4,1)*al(e,3,2))
173 zi0(e,1) = detaij(e)*((al(e,4,2)*al(e,3,3)
174 & - al(e,3,2)*al(e,4,3))*(xts1-al(e,1,1)) +
175 & (al(e,3,1)*al(e,4,3)
176 & - al(e,4,1)*al(e,3,3))*(xts2-al(e,1,2)) +
177 & (al(e,4,1)*al(e,3,2)
178 & - al(e,3,1)*al(e,4,2))*(xts3-al(e,1,3)))
180 zi0(e,2) = detaij(e)*((al(e,2,2)*al(e,4,3)
181 & - al(e,4,2)*al(e,2,3))*(xts1-al(e,1,1)) +
182 & (al(e,4,1)*al(e,2,3)
183 & - al(e,2,1)*al(e,4,3))*(xts2-al(e,1,2)) +
184 & (al(e,2,1)*al(e,4,2)
185 & - al(e,4,1)*al(e,2,2))*(xts3-al(e,1,3)))
187 zi0(e,3) = detaij(e)*((al(e,3,2)*al(e,2,3)
188 & - al(e,2,2)*al(e,3,3))*(xts1-al(e,1,1)) +
189 & (al(e,2,1)*al(e,3,3)
190 & - al(e,3,1)*al(e,2,3))*(xts2-al(e,1,2)) +
191 & (al(e,3,1)*al(e,2,2)
192 & - al(e,2,1)*al(e,3,2))*(xts3-al(e,1,3)))