Lines Matching refs:lDir
12 real*8 lDir(npro,nshl,3), lStsVec(npro,nshl,nResDims), local
15 call e3StsDir( xl, lDir )
18 lStsVec(:,i,1) = lDir(:,i,1) * lDir(:,i,1)
19 lStsVec(:,i,2) = lDir(:,i,2) * lDir(:,i,2)
20 lStsVec(:,i,3) = lDir(:,i,3) * lDir(:,i,3)
22 lStsVec(:,i,4) = lDir(:,i,1) * lDir(:,i,2)
23 lStsVec(:,i,5) = lDir(:,i,2) * lDir(:,i,3)
24 lStsVec(:,i,6) = lDir(:,i,3) * lDir(:,i,1)
47 real*8 lDir(npro,nshl,3), lStsVec(npro,nshl,nResDims) local
49 call e3StsDir( xl, lDir )
52 lStsVec(:,i,1) = lDir(:,i,1) * rl(:,i,4)
53 lStsVec(:,i,2) = lDir(:,i,2) * rl(:,i,4)
54 lStsVec(:,i,3) = lDir(:,i,3) * rl(:,i,4)
56 lStsVec(:,i,4) = lDir(:,i,1) * rl(:,i,1)
57 lStsVec(:,i,5) = lDir(:,i,2) * rl(:,i,2)
58 lStsVec(:,i,6) = lDir(:,i,3) * rl(:,i,3)
60 lStsVec(:,i,7) = lDir(:,i,1) * rl(:,i,2)
61 & + lDir(:,i,2) * rl(:,i,1)
62 lStsVec(:,i,8) = lDir(:,i,2) * rl(:,i,3)
63 & + lDir(:,i,3) * rl(:,i,2)
64 lStsVec(:,i,9) = lDir(:,i,3) * rl(:,i,1)
65 & + lDir(:,i,1) * rl(:,i,3)
74 subroutine e3StsDir( xl, lDir ) argument
82 real*8 xl(npro,nenl,3), lDir(npro,nshl,3) local
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
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)
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)
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