Lines Matching refs:e
18 integer e, founde local
144 do e = 1, npro
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(:),
152 founde=e
164 do e = 1, npro
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))
171 detaij(e) = 1./detaij(e)
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)))
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(:),
205 founde=e