29
30
31
32#include "implicit_f.inc"
33
34
35
36#include "mvsiz_p.inc"
37
38
39
40 INTEGER JFT, JLT, ISYM
42 . vq(2,2,*), kk(2,2,*)
43
44
45
46 INTEGER I, , EP
48 . k(2,2,mvsiz)
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73 DO ep=jft,jlt
74 DO i=1,2
75 DO j=1,2
76 k(i,j,ep) =
77 1 vq(i,1,ep)*(kk(1,1,ep)*vq(j,1,ep)+kk(1,2,ep)*vq(j,2,ep)) +
78 2 vq(i,2,ep)*(kk(2,1,ep)*vq(j,1,ep)+kk(2,2,ep)*vq(j,2,ep))
79 ENDDO
80 ENDDO
81 ENDDO
82 DO ep=jft,jlt
83 DO i=1,2
84 DO j=1,2
85 kk(i,j,ep) = k(i,j,ep)
86 ENDDO
87 ENDDO
88 ENDDO
89
90
91 RETURN