Lines Matching refs:zi0
11 & zi0(npro,nsd), detaij(npro), dzi0(npro,nsd), local
52 zi0(:,1) = detaij(:)*((al(:,4,2)*al(:,3,3)
60 zi0(:,2) = detaij(:)*((al(:,2,2)*al(:,4,3)
67 zi0(:,3) = detaij(:)*((al(:,3,2)*al(:,2,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)*
109 & zi0(:,2)*zi0(:,3) - xts1
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)*
115 & zi0(:,2)*zi0(:,3) - xts2
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)*
121 & zi0(:,2)*zi0(:,3) - xts3
140 zi0(:,:) = zi0(:,:) + dzi0(:,:)
145 if ((abs(zi0(e,1)).lt.(one+tolpt)).and.
146 & (abs(zi0(e,2)).lt.(one+tolpt)).and.
147 & (abs(zi0(e,3)).lt.(one+tolpt))) then ! got the element
149 call shphex (ipord, zi0(e,:),shape(:),
173 zi0(e,1) = detaij(e)*((al(e,4,2)*al(e,3,3)
180 zi0(e,2) = detaij(e)*((al(e,2,2)*al(e,4,3)
187 zi0(e,3) = detaij(e)*((al(e,3,2)*al(e,2,3)
194 if ((zi0(e,1)+zi0(e,2)+zi0(e,3)).lt.(one+tolpt).and.
195 …& zi0(e,1).lt.(one +tolpt).and. !should not be necessary; the limit in the tet …
196 & zi0(e,1).gt.(zero-tolpt).and.
197 & zi0(e,2).lt.(one +tolpt).and. !should not be necessary
198 & zi0(e,2).gt.(zero-tolpt).and.
199 & zi0(e,3).lt.(one +tolpt).and. !should not be necessary
200 & zi0(e,3).gt.(zero-tolpt)) then
202 call shptet (ipord, zi0(e,:), shape(:),
214 parptts(jj,:)=zi0(founde,:)