Lines Matching refs:e
83 integer e local
90 do e = 1, npro
92 x12 = xl(e,2,1) - xl(e,1,1)
93 x13 = xl(e,3,1) - xl(e,1,1)
94 x14 = xl(e,4,1) - xl(e,1,1)
95 x23 = xl(e,3,1) - xl(e,2,1)
96 x24 = xl(e,4,1) - xl(e,2,1)
97 x34 = xl(e,4,1) - xl(e,3,1)
99 y12 = xl(e,2,2) - xl(e,1,2)
100 y13 = xl(e,3,2) - xl(e,1,2)
101 y14 = xl(e,4,2) - xl(e,1,2)
102 y23 = xl(e,3,2) - xl(e,2,2)
103 y24 = xl(e,4,2) - xl(e,2,2)
104 y34 = xl(e,4,2) - xl(e,3,2)
106 z12 = xl(e,2,3) - xl(e,1,3)
107 z13 = xl(e,3,3) - xl(e,1,3)
108 z14 = xl(e,4,3) - xl(e,1,3)
109 z23 = xl(e,3,3) - xl(e,2,3)
110 z24 = xl(e,4,3) - xl(e,2,3)
111 z34 = xl(e,4,3) - xl(e,3,3)
118 lDir(e,1,1) = fct * (y14*(z12 - z13) + y12*(z13 - z14)
120 lDir(e,1,2) = fct * ( x14*(-z12 + z13) + x13*(z12 - z14)
122 lDir(e,1,3) = fct * ( x14*(y12 - y13)
125 lDir(e,2,1) = fct * (-(y13*z12) + y14*z12 + y12*z13
127 lDir(e,2,2) = fct * (x13*z12 - x14*z12-x12*z13 + x12*z14
129 lDir(e,2,3) = fct * (-(x13*y12) + x14*y12 + x12*y13
132 lDir(e,3,1) = fct * (y12*z13 - y14*z13 + y13*(-z12 + z14)
134 lDir(e,3,2) = fct * (-(x12*z13) + x14*z13 + x13*(z12 - z14)
136 lDir(e,3,3) = fct * (x12*y13 - x14*y13 + x13*(-y12 + y14)
139 lDir(e,4,1) = fct * (y14*(z12 - z13) - y12*z14 + y13*z14
141 lDir(e,4,2) = fct * (x14*(-z12 + z13) + x12*z14 - x13*z14
143 lDir(e,4,3) = fct * (x14*(y12 - y13) - x12*y14 + x13*y14
152 do e = 1, npro
154 x12 = xl(e,2,1) - xl(e,1,1)
155 x13 = xl(e,3,1) - xl(e,1,1)
156 x14 = xl(e,4,1) - xl(e,1,1)
157 x23 = xl(e,3,1) - xl(e,2,1)
158 x24 = xl(e,4,1) - xl(e,2,1)
159 x34 = xl(e,4,1) - xl(e,3,1)
161 y12 = xl(e,2,2) - xl(e,1,2)
162 y13 = xl(e,3,2) - xl(e,1,2)
163 y14 = xl(e,4,2) - xl(e,1,2)
164 y23 = xl(e,3,2) - xl(e,2,2)
165 y24 = xl(e,4,2) - xl(e,2,2)
166 y34 = xl(e,4,2) - xl(e,3,2)
168 z12 = xl(e,2,3) - xl(e,1,3)
169 z13 = xl(e,3,3) - xl(e,1,3)
170 z14 = xl(e,4,3) - xl(e,1,3)
171 z23 = xl(e,3,3) - xl(e,2,3)
172 z24 = xl(e,4,3) - xl(e,2,3)
173 z34 = xl(e,4,3) - xl(e,3,3)
178 lDir(e,1,1) = fct * (y14*(z12 - z13) + y12*(z13 - z14)
180 lDir(e,1,2) = fct * ( x14*(-z12 + z13) + x13*(z12 - z14)
182 lDir(e,1,3) = fct * ( x14*(y12 - y13)
185 lDir(e,2,1) = fct * (-(y13*z12) + y14*z12 + y12*z13
187 lDir(e,2,2) = fct * (x13*z12 - x14*z12-x12*z13 + x12*z14
189 lDir(e,2,3) = fct * (-(x13*y12) + x14*y12 + x12*y13
192 lDir(e,3,1) = fct * (y12*z13 - y14*z13 + y13*(-z12 + z14)
194 lDir(e,3,2) = fct * (-(x12*z13) + x14*z13 + x13*(z12 - z14)
196 lDir(e,3,3) = fct * (x12*y13 - x14*y13 + x13*(-y12 + y14)
199 lDir(e,4,1) = fct * (y14*(z12 - z13) - y12*z14 + y13*z14
201 lDir(e,4,2) = fct * (x14*(-z12 + z13) + x12*z14 - x13*z14
203 lDir(e,4,3) = fct * (x14*(y12 - y13) - x12*y14 + x13*y14
208 lDir(e,5,1) = pt25*(-(y13*z12) + y14*z12 + y12*(z13-z14))
209 lDir(e,5,2) = pt25*(x13*z12 - x14*z12 + x12*(-z13 + z14))
210 lDir(e,5,3) = pt25*(-(x13*y12) + x14*y12 + x12*(y13 - y14))
212 lDir(e,6,1) = pt25*(-(y13*z12) + y12*z13 + y24*z23-y23*z24)
213 lDir(e,6,2) = pt25*(x13*z12 - x12*z13 - x24*z23 + x23*z24)
214 lDir(e,6,3) = pt25*(-(x13*y12) + x12*y13 + x24*y23-x23*y24)
216 lDir(e,7,1) = pt25*((y12 - y14)*z13 + y13*(-z12 + z14))
217 lDir(e,7,2) = pt25*((-x12 + x14)*z13 + x13*(z12 - z14))
218 lDir(e,7,3) = pt25*((x12 - x14)*y13 + x13*(-y12 + y14))
220 lDir(e,8,1) = pt25*(y14*(z12 - z13) + (-y12 + y13)*z14)
221 lDir(e,8,2) = pt25*(x14*(-z12 + z13) + (x12 - x13)*z14)
222 lDir(e,8,3) = pt25*(x14*(y12 - y13) + (-x12 + x13)*y14)
224 lDir(e,9,1) = pt25*(y14*z12 - y12*z14 + y24*z23 - y23*z24)
225 lDir(e,9,2) = pt25*(-(x14*z12) + x12*z14 - x24*z23+x23*z24)
226 lDir(e,9,3) = pt25*(x14*y12 - x12*y14 + x24*y23 - x23*y24)
228 lDir(e,10,1) = pt25*(-(y14*z13) + y13*z14+y24*z23-y23*z24)
229 lDir(e,10,2) = pt25*(x14*z13 - x13*z14-x24*z23 + x23*z24)
230 lDir(e,10,3) = pt25*(-(x14*y13) + x13*y14+x24*y23-x23*y24)
239 do e = 1, npro
241 x12 = xl(e,2,1) - xl(e,1,1)
242 x13 = xl(e,3,1) - xl(e,1,1)
243 x14 = xl(e,4,1) - xl(e,1,1)
244 x23 = xl(e,3,1) - xl(e,2,1)
245 x24 = xl(e,4,1) - xl(e,2,1)
248 y12 = xl(e,2,2) - xl(e,1,2)
249 y13 = xl(e,3,2) - xl(e,1,2)
250 y14 = xl(e,4,2) - xl(e,1,2)
251 y23 = xl(e,3,2) - xl(e,2,2)
252 y24 = xl(e,4,2) - xl(e,2,2)
255 z12 = xl(e,2,3) - xl(e,1,3)
256 z13 = xl(e,3,3) - xl(e,1,3)
257 z14 = xl(e,4,3) - xl(e,1,3)
258 z23 = xl(e,3,3) - xl(e,2,3)
259 z24 = xl(e,4,3) - xl(e,2,3)
265 lDir(e,1,1) = fct * (y14*(z12 - z13) + y12*(z13 - z14)
267 lDir(e,1,2) = fct * ( x14*(-z12 + z13) + x13*(z12 - z14)
269 lDir(e,1,3) = fct * ( x14*(y12 - y13)
272 lDir(e,2,1) = fct * (-(y13*z12) + y14*z12 + y12*z13
274 lDir(e,2,2) = fct * (x13*z12 - x14*z12-x12*z13 + x12*z14
276 lDir(e,2,3) = fct * (-(x13*y12) + x14*y12 + x12*y13
279 lDir(e,3,1) = fct * (y12*z13 - y14*z13 + y13*(-z12 + z14)
281 lDir(e,3,2) = fct * (-(x12*z13) + x14*z13 + x13*(z12 - z14)
283 lDir(e,3,3) = fct * (x12*y13 - x14*y13 + x13*(-y12 + y14)
286 lDir(e,4,1) = fct * (y14*(z12 - z13) - y12*z14 + y13*z14
288 lDir(e,4,2) = fct * (x14*(-z12 + z13) + x12*z14 - x13*z14
290 lDir(e,4,3) = fct * (x14*(y12 - y13) - x12*y14 + x13*y14
295 lDir(e,5,1) = pt25*(-(y13*z12) + y14*z12 + y12*(z13-z14))
296 lDir(e,5,2) = pt25*(x13*z12 - x14*z12 + x12*(-z13 + z14))
297 lDir(e,5,3) = pt25*(-(x13*y12) + x14*y12 + x12*(y13 - y14))
298 lDir(e,6,1) = lDir(e,5,1)
299 lDir(e,6,2) = lDir(e,5,2)
300 lDir(e,6,3) = lDir(e,5,3)
302 lDir(e,7,1) = pt25*(-(y13*z12) + y12*z13 + y24*z23-y23*z24)
303 lDir(e,7,2) = pt25*(x13*z12 - x12*z13 - x24*z23 + x23*z24)
304 lDir(e,7,3) = pt25*(-(x13*y12) + x12*y13 + x24*y23-x23*y24)
305 lDir(e,8,1) = lDir(e,7,1)
306 lDir(e,8,2) = lDir(e,7,2)
307 lDir(e,8,3) = lDir(e,7,3)
309 lDir(e,9,1) = pt25*((y12 - y14)*z13 + y13*(-z12 + z14))
310 lDir(e,9,2) = pt25*((-x12 + x14)*z13 + x13*(z12 - z14))
311 lDir(e,9,3) = pt25*((x12 - x14)*y13 + x13*(-y12 + y14))
312 lDir(e,10,1) = lDir(e,9,1)
313 lDir(e,10,2) = lDir(e,9,2)
314 lDir(e,10,3) = lDir(e,9,3)
316 lDir(e,11,1) = pt25*(y14*(z12 - z13) + (-y12 + y13)*z14)
317 lDir(e,11,2) = pt25*(x14*(-z12 + z13) + (x12 - x13)*z14)
318 lDir(e,11,3) = pt25*(x14*(y12 - y13) + (-x12 + x13)*y14)
319 lDir(e,12,1) = lDir(e,11,1)
320 lDir(e,12,2) = lDir(e,11,2)
321 lDir(e,12,3) = lDir(e,11,3)
324 lDir(e,13,1) = pt25*(y14*z12 - y12*z14 + y24*z23 - y23*z24)
325 lDir(e,13,2) = pt25*(-(x14*z12) + x12*z14-x24*z23+x23*z24)
326 lDir(e,13,3) = pt25*(x14*y12 - x12*y14 + x24*y23 - x23*y24)
327 lDir(e,14,1) = lDir(e,13,1)
328 lDir(e,14,2) = lDir(e,13,2)
329 lDir(e,14,3) = lDir(e,13,3)
331 lDir(e,15,1) = pt25*(-(y14*z13) + y13*z14+y24*z23-y23*z24)
332 lDir(e,15,2) = pt25*(x14*z13 - x13*z14-x24*z23 + x23*z24)
333 lDir(e,15,3) = pt25*(-(x14*y13) + x13*y14+x24*y23-x23*y24)
334 lDir(e,16,1) = lDir(e,15,1)
335 lDir(e,16,2) = lDir(e,15,2)
336 lDir(e,16,3) = lDir(e,15,3)
340 lDir(e,17,1) = pt5*(-(y13*z12) + y12*z13)
341 lDir(e,17,2) = pt5*(x13*z12 - x12*z13)
342 lDir(e,17,3) = pt5*(-(x13*y12) + x12*y13)
344 lDir(e,18,1) = pt5*(y14*z12 - y12*z14)
345 lDir(e,18,2) = pt5*(-(x14*z12) + x12*z14)
346 lDir(e,18,3) = pt5*(x14*y12 - x12*y14)
348 lDir(e,19,1) = pt5*(y24*z23 - y23*z24)
349 lDir(e,19,2) = pt5*(-(x24*z23) + x23*z24)
350 lDir(e,19,3) = pt5*(x24*y23 - x23*y24)
352 lDir(e,20,1) = pt5*(-(y14*z13) + y13*z14)
353 lDir(e,20,2) = pt5*(x14*z13 - x13*z14)
354 lDir(e,20,3) = pt5*(-(x14*y13) + x13*y14)
362 do e = 1, npro
363 x13 = xl(e,1,1) - xl(e,3,1)
364 x16 = xl(e,1,1) - xl(e,6,1)
365 x18 = xl(e,1,1) - xl(e,8,1)
366 x24 = xl(e,2,1) - xl(e,4,1)
367 x25 = xl(e,2,1) - xl(e,5,1)
368 x27 = xl(e,2,1) - xl(e,7,1)
369 x36 = xl(e,3,1) - xl(e,6,1)
370 x38 = xl(e,3,1) - xl(e,8,1)
371 x45 = xl(e,4,1) - xl(e,5,1)
372 x47 = xl(e,4,1) - xl(e,7,1)
373 x57 = xl(e,5,1) - xl(e,7,1)
374 x68 = xl(e,6,1) - xl(e,8,1)
376 y13 = xl(e,1,2) - xl(e,3,2)
377 y16 = xl(e,1,2) - xl(e,6,2)
378 y18 = xl(e,1,2) - xl(e,8,2)
379 y24 = xl(e,2,2) - xl(e,4,2)
380 y25 = xl(e,2,2) - xl(e,5,2)
381 y27 = xl(e,2,2) - xl(e,7,2)
382 y36 = xl(e,3,2) - xl(e,6,2)
383 y38 = xl(e,3,2) - xl(e,8,2)
384 y45 = xl(e,4,2) - xl(e,5,2)
385 y47 = xl(e,4,2) - xl(e,7,2)
386 y57 = xl(e,5,2) - xl(e,7,2)
387 y68 = xl(e,6,2) - xl(e,8,2)
389 z13 = xl(e,1,3) - xl(e,3,3)
390 z16 = xl(e,1,3) - xl(e,6,3)
391 z18 = xl(e,1,3) - xl(e,8,3)
392 z24 = xl(e,2,3) - xl(e,4,3)
393 z25 = xl(e,2,3) - xl(e,5,3)
394 z27 = xl(e,2,3) - xl(e,7,3)
395 z36 = xl(e,3,3) - xl(e,6,3)
396 z38 = xl(e,3,3) - xl(e,8,3)
397 z45 = xl(e,4,3) - xl(e,5,3)
398 z47 = xl(e,4,3) - xl(e,7,3)
399 z57 = xl(e,5,3) - xl(e,7,3)
400 z68 = xl(e,6,3) - xl(e,8,3)
439 lDir(e,1,1) = fct * (-y24 * z45 + y36 * z24 - y68 * z45
441 lDir(e,2,1) = fct * (-y16 * z63 + y54 * z16 - y47 * z63
443 lDir(e,3,1) = fct * (-y42 * z27 + y18 * z42 - y86 * z27
445 lDir(e,4,1) = fct * (-y38 * z81 + y72 * z38 - y25 * z81
447 lDir(e,5,1) = fct * (-y61 * z18 + y27 * z61 - y74 * z18
449 lDir(e,6,1) = fct * (-y57 * z72 + y81 * z57 - y13 * z72
451 lDir(e,7,1) = fct * (-y83 * z36 + y45 * z83 - y52 * z36
453 lDir(e,8,1) = fct * (-y75 * z54 + y63 * z75 - y31 * z54
456 lDir(e,1,2) = fct * (-z24 * x45 + z36 * x24 - z68 * x45
458 lDir(e,2,2) = fct * (-z16 * x63 + z54 * x16 - z47 * x63
460 lDir(e,3,2) = fct * (-z42 * x27 + z18 * x42 - z86 * x27
462 lDir(e,4,2) = fct * (-z38 * x81 + z72 * x38 - z25 * x81
464 lDir(e,5,2) = fct * (-z61 * x18 + z27 * x61 - z74 * x18
466 lDir(e,6,2) = fct * (-z57 * x72 + z81 * x57 - z13 * x72
468 lDir(e,7,2) = fct * (-z83 * x36 + z45 * x83 - z52 * x36
470 lDir(e,8,2) = fct * (-z75 * x54 + z63 * x75 - z31 * x54
473 lDir(e,1,3) = fct * (-x24 * y45 + x36 * y24 - x68 * y45
475 lDir(e,2,3) = fct * (-x16 * y63 + x54 * y16 - x47 * y63
477 lDir(e,3,3) = fct * (-x42 * y27 + x18 * y42 - x86 * y27
479 lDir(e,4,3) = fct * (-x38 * y81 + x72 * y38 - x25 * y81
481 lDir(e,5,3) = fct * (-x61 * y18 + x27 * y61 - x74 * y18
483 lDir(e,6,3) = fct * (-x57 * y72 + x81 * y57 - x13 * y72
485 lDir(e,7,3) = fct * (-x83 * y36 + x45 * y83 - x52 * y36
487 lDir(e,8,3) = fct * (-x75 * y54 + x63 * y75 - x31 * y54