33
34
35
36#include "implicit_f.inc"
37
38
39
40 INTEGER JFT, JLT,NEL
41
43 . sig(nel,6),sign(nel,6), g1x(*),g1y(*),g1z(*),
44 . g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*)
45
46
47
48 INTEGER I
49
51 . sx,sy,sz
52
53 DO i=jft,jlt
54 sx = sig(i,1)*g1x(i)+sig(i,4)*g1y(i)+sig(i,6)*g1z(i)
55 sy = sig(i,4)*g1x(i)+sig(i,2)*g1y(i)+sig(i,5)*g1z(i)
56 sz = sig(i,6)*g1x(i)+sig(i,5)*g1y(i)+sig(i,3)*g1z(i)
57 sign(i,1) = sx*g1x(i)+sy*g1y(i)+sz*g1z(i)
58 sign(i,4) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
59 sign(i,6) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
60 sx = sig(i,1)*g2x(i)+sig(i,4)*g2y(i)+sig(i,6)*g2z(i)
61 sy = sig(i,4)*g2x(i)+sig(i,2)*g2y(i)+sig(i,5)*g2z(i)
62 sz = sig(i,6)*g2x(i)+sig(i,5)*g2y(i)+sig(i,3)*g2z(i)
63 sign(i,2) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
64 sign(i,5) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
65 sx = sig(i,1)*g3x(i)+sig(i,4)*g3y(i)+sig(i,6)*g3z(i)
66 sy = sig(i,4)*g3x(i)+sig(i,2)*g3y(i)+sig(i,5)*g3z(i)
67 sz = sig(i,6)*g3x(i)+sig(i,5)*g3y(i)+sig(i,3)*g3z(i)
68 sign(i,3) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
69 ENDDO
70
71 RETURN