30
31
32
33#include "implicit_f.inc"
34
35
36
37#include "mvsiz_p.inc"
38
39
40
41 INTEGER JFT, JLT
42
44 . sig(mvsiz,6), g1x(*),g1y(*),g1z(*),
45 . g2x(*),g2y(*),g2z(*),g3x(*),g3y(*),g3z(*),off(*)
46
47
48
49 INTEGER I
50
52 . sx,sy,sz,sign(mvsiz,6)
53
54 DO i=jft,jlt
55 IF (off(i) <=one) cycle
56 sx = sig(i,1)*g1x(i)+sig(i,4)*g1y(i)+sig(i,6)*g1z(i)
57 sy = sig(i,4)*g1x(i)+sig(i,2)*g1y(i)+sig(i,5)*g1z(i)
58 sz = sig(i,6)*g1x(i)+sig(i,5)*g1y(i)+sig(i,3)*g1z(i)
59 sign(i,1) = sx*g1x(i)+sy*g1y(i)+sz*g1z(i)
60 sign(i,4) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
61 sign(i,6) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
62 sx = sig(i,1)*g2x(i)+sig(i,4)*g2y(i)+sig(i,6)*g2z(i)
63 sy = sig(i,4)*g2x(i)+sig(i,2)*g2y(i)+sig(i,5)*g2z(i)
64 sz = sig(i,6)*g2x(i)+sig(i,5)*g2y(i)+sig(i,3)*g2z(i)
65 sign(i,2) = sx*g2x(i)+sy*g2y(i)+sz*g2z(i)
66 sign(i,5) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
67 sx = sig(i,1)*g3x(i)+sig(i,4)*g3y(i)+sig(i,6)*g3z(i)
68 sy = sig(i,4)*g3x(i)+sig(i,2)*g3y(i)+sig(i,5)*g3z(i)
69 sz = sig(i,6)*g3x(i)+sig(i,5)*g3y(i)+sig(i,3)*g3z(i)
70 sign(i,3) = sx*g3x(i)+sy*g3y(i)+sz*g3z(i)
71 sig(i,1) = sign(i,1)
72 sig(i,2) = sign(i,2)
73 sig(i,3) = sign(i,3)
74 sig(i,4) = sign(i,4)
75 sig(i,5) = sign(i,5)
76 sig(i,6) = sign(i,6)
77 ENDDO
78
79 RETURN