41
42
43
44#include "implicit_f.inc"
45
46
47
48#include "mvsiz_p.inc"
49
50
51
52
53
54
55 INTEGER NEL
56
57
58
59
60 DOUBLE PRECISION
61 . XD1(MVSIZ), XD2(MVSIZ), XD3(MVSIZ), XD4(MVSIZ),
62 . XD5(MVSIZ), XD6(MVSIZ), XD7(MVSIZ), XD8(MVSIZ),
63 . YD1(MVSIZ), YD2(MVSIZ), YD3(MVSIZ), YD4(MVSIZ),
64 . YD5(MVSIZ), YD6(MVSIZ), YD7(MVSIZ), YD8(MVSIZ),
65 . ZD1(MVSIZ), ZD2(MVSIZ), ZD3(MVSIZ), ZD4(MVSIZ),
66 . ZD5(MVSIZ), ZD6(MVSIZ), ZD7(MVSIZ), ZD8(MVSIZ),
67 . SAV(NEL,21),INVJ(MVSIZ,3,3),
68 . V1(MVSIZ,3,3), V2(MVSIZ,3,3), V3(MVSIZ,3,3), V4(MVSIZ,3,3),
69 . V5(MVSIZ,3,3), V6(MVSIZ,3,3), V7(MVSIZ
70
71
72
73
74 INTEGER I,J,JJ
75
76 DOUBLE PRECISION DX_DR,DX_DS,DX_DT,DY_DR,DY_DS,DY_DT,DZ_DR,DZ_DS,DZ_DT
77 DOUBLE PRECISION FMAT(3,3),UM(3,3),IC,I2C,I3C,IU,I2U,I3U
78 DOUBLE PRECISION
79 . SA11,SA21,SA31,SA12,SA22,SA32,SA13,SA23,SA33,
80 . SB11,SB21,SB31,SB12,SB22,SB32,SB13,SB23,SB33,
81 . SC11,SC21,SC31,SC12,SC22,SC32,SC13,SC23,SC33,
82 . SD11,SD21,SD31,SD12,SD22,SD32,SD13,SD23,SD33,
83 . SE11,SE21,SE31,SE12,SE22,SE32,SE13,SE23,SE33,
84 . SF11,SF21,SF31,SF12,SF22,SF32,SF13,SF23,SF33,
85 . SG11,SG21,SG31,SG12,SG22,SG32,SG13,SG23,SG33,
86 . S11,S21,S31,S12,S22,S32,S13,S23,S33
87 DOUBLE PRECISION
88 . B11,B21,B31,B12,,B32,B13,B23,B33,,
89 . BI11,BI21,BI31,BI12,BI22,BI32,BI13,BI23,BI33
90 DOUBLE PRECISION AA
91 DOUBLE PRECISION
92 . DN1_DX(MVSIZ), DN1_DY(MVSIZ), DN1_DZ(MVSIZ),
93 . DN2_DX(MVSIZ), DN2_DY(MVSIZ), DN2_DZ(MVSIZ),
94 . DN3_DX(MVSIZ), DN3_DY(MVSIZ), DN3_DZ(MVSIZ),
95 . DN4_DX(MVSIZ), DN4_DY(MVSIZ), DN4_DZ(MVSIZ),
96 . DN5_DX(MVSIZ), DN5_DY(MVSIZ), DN5_DZ(MVSIZ),
97 . DN6_DX(MVSIZ), DN6_DY(MVSIZ), DN6_DZ(MVSIZ),
98 . DN7_DX(MVSIZ), DN7_DY(MVSIZ), DN7_DZ(MVSIZ),
99 . DN8_DX(MVSIZ), DN8_DY(MVSIZ), DN8_DZ(MVSIZ),
100 . XI1(MVSIZ), XI2(MVSIZ), XI3(MVSIZ), XI4(MVSIZ),
101 . XI5(MVSIZ), XI6(MVSIZ), XI7(MVSIZ), (MVSIZ),
102 . YI1(MVSIZ), YI2(MVSIZ), YI3(MVSIZ), YI4(MVSIZ),
103 . YI5(MVSIZ), YI6(MVSIZ), YI7(MVSIZ), YI8(MVSIZ),
104 . ZI1(MVSIZ), ZI2(MVSIZ), ZI3(MVSIZ), ZI4(MVSIZ),
105 . ZI5(MVSIZ), ZI6(MVSIZ), ZI7(MVSIZ), ZI8(MVSIZ)
106
107 aa = one/eight
108
109 DO i=1,nel
110 xi1(i) = zero
111 yi1(i) = zero
112 zi1(i) = zero
113 xi2(i) = xd2(i)-xd1(i)
114 yi2(i) = yd2(i)-yd1(i)
115 zi2(i) = zd2(i)-zd1(i)
116 xi3(i) = xd3(i)-xd1(i)
117 yi3(i) = yd3(i)-yd1(i)
118 zi3(i) = zd3(i)-zd1(i)
119 xi4(i) = xd4(i)-xd1(i)
120 yi4(i) = yd4(i)-yd1(i)
121 zi4(i) = zd4(i)-zd1(i)
122 xi5(i) = xd5(i)-xd1(i)
123 yi5(i) = yd5(i)-yd1(i)
124 zi5(i) = zd5(i)-zd1(i)
125 xi6(i) = xd6(i)-xd1(i)
126 yi6(i) = yd6(i)-yd1(i)
127 zi6(i) = zd6(i)-zd1(i)
128 xi7(i) = xd7(i)-xd1(i)
129 yi7(i) = yd7(i)-yd1(i)
130 zi7(i) = zd7(i)-zd1(i)
131 xi8(i) = xd8(i)-xd1(i)
132 yi8(i) = yd8(i)-yd1(i)
133 zi8(i) = zd8(i)-zd1(i)
134 ENDDO
135
136 DO i=1,nel
137 dn1_dx(i) = (-invj(i,1,1)-invj(i,2,1)-invj(i,3,1))*aa
138 dn2_dx(i) = (-invj(i,1,1)-invj(i,2,1)+invj(i,3,1))*aa
139 dn3_dx(i) = ( invj(i,1,1)-invj(i,2,1)+invj(i,3,1))*aa
140 dn4_dx(i) = ( invj(i,1,1)-invj(i,2,1)-invj(i,3,1))*aa
141 dn5_dx(i) = (-invj(i,1,1)+invj(i,2,1)-invj(i,3,1))*aa
142 dn6_dx(i) = (-invj(i,1,1)+invj(i,2,1)+invj(i,3,1))*aa
143 dn7_dx(i) = ( invj(i,1,1)+invj(i,2,1)+invj(i,3,1))*aa
144 dn8_dx(i) = ( invj(i,1,1)+invj(i,2,1)-invj(i,3,1))*aa
145
146 dn1_dy(i) = (-invj(i,1,2)-invj(i,2,2)-invj(i,3,2))*aa
147 dn2_dy(i) = (-invj(i,1,2)-invj(i,2,2)+invj(i,3,2))*aa
148 dn3_dy(i) = ( invj(i,1,2)-invj(i,2,2)+invj(i,3,2))*aa
149 dn4_dy(i) = ( invj(i,1,2)-invj(i,2,2)-invj(i,3,2))*aa
150 dn5_dy(i) = (-invj(i,1,2)+invj(i,2,2)-invj(i,3,2))*aa
151 dn6_dy(i) = (-invj(i,1,2)+invj(i,2,2)+invj(i,3,2))*aa
152 dn7_dy(i) = ( invj(i,1,2)+invj(i,2,2)+invj(i,3,2))*aa
153 dn8_dy(i) = ( invj(i,1,2)+invj(i,2,2)-invj(i,3,2))*aa
154
155 dn1_dz(i) = (-invj(i,1,3)-invj(i,2,3)-invj(i,3,3))*aa
156 dn2_dz(i) = (-invj(i,1,3)-invj(i,2,3)+invj(i,3,3))*aa
157 dn3_dz(i) = ( invj(i,1,3)-invj(i,2,3)+invj(i,3,3))*aa
158 dn4_dz(i) = ( invj(i,1,3)-invj(i,2,3)-invj(i,3,3))*aa
159 dn5_dz(i) = (-invj(i,1,3)+invj(i,2,3)-invj(i,3,3))*aa
160 dn6_dz(i) = (-invj(i,1,3)+invj(i,2,3)+invj(i,3,3))*aa
161 dn7_dz(i) = ( invj(i,1,3)+invj(i,2,3)+invj(i,3,3))*aa
162 dn8_dz(i) = ( invj(i,1,3)+invj(i,2,3)-invj(i,3,3))*aa
163
164 sa11 = r(2,1,i)*zi2(i)-r(3,1,i)*yi2(i)
165 sa12 =-r(1,1,i)*zi2(i)+r(3,1,i)*xi2(i)
166 sa13 = r(1,1,i)*yi2(i)-r(2,1,i)*xi2(i)
167 sa21 = r(2,2,i)*zi2(i)-r(3,2,i)*yi2(i)
168 sa22 =-r(1,2,i)*zi2(i)+r(3,2,i)*xi2(i)
169 sa23 = r(1,2,i)*yi2(i)-r(2,2,i)*xi2(i)
170 sa31 = r(2,3,i)*zi2(i)-r(3,3,i)*yi2(i)
171 sa32 =-r(1,3,i)*zi2(i)+r(3,3,i)*xi2(i)
172 sa33 = r(1,3,i)*yi2(i)-r(2,3,i)*xi2(i)
173 sb11 = r(2,1,i)*zi3(i)-r(3,1,i)*yi3(i)
174 sb12 =-r(1,1,i)*zi3(i)+r(3,1,i)*xi3(i)
175 sb13 = r(1,1,i)*yi3(i)-r(2,1,i)*xi3(i)
176 sb21 = r(2,2,i)*zi3(i)-r(3,2,i)*yi3(i)
177 sb22 =-r(1,2,i)*zi3(i)+r(3,2,i)*xi3(i)
178 sb23 = r(1,2,i)*yi3(i)-r(2,2,i)*xi3(i)
179 sb31 = r(2,3,i)*zi3(i)-r(3,3,i)*yi3(i)
180 sb32 =-r(1,3,i)*zi3(i)+r(3,3,i)*xi3(i)
181 sb33 = r(1,3,i)*yi3(i)-r(2,3,i)*xi3(i)
182 sc11 = r(2,1,i)*zi4(i)-r(3,1,i)*yi4(i)
183 sc12 =-r(1,1,i)*zi4(i)+r(3,1,i)*xi4(i)
184 sc13 = r(1,1,i)*yi4(i)-r(2,1,i)*xi4(i)
185 sc21 = r(2,2,i)*zi4(i)-r(3,2,i)*yi4(i)
186 sc22 =-r(1,2,i)*zi4(i)+r(3,2,i)*xi4(i)
187 sc23 = r(1,2,i)*yi4(i)-r(2,2,i)*xi4(i)
188 sc31 = r(2,3,i)*zi4(i)-r(3,3,i)*yi4(i)
189 sc32 =-r(1,3,i)*zi4(i)+r(3,3,i)*xi4(i)
190 sc33 = r(1,3,i)*yi4(i)-r(2,3,i)*xi4(i)
191 sd11 = r(2,1,i)*zi5(i)-r(3,1,i)*yi5(i)
192 sd12 =-r(1,1,i)*zi5(i)+r(3,1,i)*xi5(i)
193 sd13 = r(1,1,i)*yi5(i)-r(2,1,i)*xi5(i)
194 sd21 = r(2,2,i)*zi5(i)-r(3,2,i)*yi5(i)
195 sd22 =-r(1,2,i)*zi5(i)+r(3,2,i)*xi5(i)
196 sd23 = r(1,2,i)*yi5(i)-r(2,2,i)*xi5(i)
197 sd31 = r(2,3,i)*zi5(i)-r(3,3,i)*yi5(i)
198 sd32 =-r(1,3,i)*zi5(i)+r(3,3,i)*xi5(i)
199 sd33 = r(1,3,i)*yi5(i)-r(2,3,i)*xi5(i)
200 se11 = r(2,1,i)*zi6(i)-r(3,1,i)*yi6(i)
201 se12 =-r(1,1,i)*zi6(i)+r(3,1,i)*xi6(i)
202 se13 = r(1,1,i)*yi6(i)-r(2,1,i)*xi6(i)
203 se21 = r(2,2,i)*zi6(i)-r(3,2,i)*yi6(i)
204 se22 =-r(1,2,i)*zi6(i)+r(3,2,i)*xi6(i)
205 se23 = r(1,2,i)*yi6(i)-r(2,2,i)*xi6(i)
206 se31 = r(2,3,i)*zi6(i)-r(3,3,i)*yi6(i)
207 se32 =-r(1,3,i)*zi6(i)+r(3,3,i)*xi6(i)
208 se33 = r(1,3,i)*yi6(i)-r(2,3,i)*xi6(i)
209 sf11 = r(2,1,i)*zi7(i)-r(3,1,i)*yi7(i)
210 sf12 =-r(1,1,i)*zi7(i)+r(3,1,i)*xi7(i)
211 sf13 = r(1,1,i)*yi7(i)-r(2,1,i)*xi7(i)
212 sf21 = r(2,2,i)*zi7(i)-r(3,2,i)*yi7(i)
213 sf22 =-r(1,2,i)*zi7(i)+r(3,2,i)*xi7(i)
214 sf23 = r(1,2,i)*yi7(i)-r(2,2,i)*xi7(i)
215 sf31 = r(2,3,i)*zi7(i)-r(3,3,i)*yi7(i)
216 sf32 =-r(1,3,i)*zi7(i)+r(3,3,i)*xi7(i)
217 sf33 = r(1,3,i)*yi7(i)-r(2,3,i)*xi7(i)
218 sg11 = r(2,1,i)*zi8(i)-r(3,1,i)*yi8(i)
219 sg12 =-r(1,1,i)*zi8(i)+r(3,1,i)*xi8(i)
220 sg13 = r(1,1,i)*yi8(i)-r(2,1,i)*xi8(i)
221 sg21 = r(2,2,i)*zi8(i)-r(3,2,i)*yi8(i)
222 sg22 =-r(1,2,i)*zi8(i)+r(3,2,i)*xi8(i)
223 sg23 = r(1,2,i)*yi8(i)-r(2,2,i)*xi8(i)
224 sg31 = r(2,3,i)*zi8(i)-r(3,3,i)*yi8(i)
225 sg32 =-r(1,3,i)*zi8(i)+r(3,3,i)*xi8(i)
226 sg33 = r(1,3,i)*yi8(i)-r(2,3,i)*xi8(i)
227
228 b11 = -dn2_dz(i)*sa21+dn2_dy(i)*sa31
229 . -dn3_dz(i)*sb21+dn3_dy(i)*sb31
230 . -dn4_dz(i)*sc21+dn4_dy(i)*sc31
231 . -dn5_dz(i)*sd21+dn5_dy(i)*sd31
232 . -dn6_dz(i)*se21+dn6_dy(i)*se31
233 . -dn7_dz(i)*sf21+dn7_dy(i)*sf31
234 . -dn8_dz(i)*sg21+dn8_dy(i)*sg31
235 b12 = -dn2_dz(i)*sa22+dn2_dy(i)*sa32
236 . -dn3_dz(i)*sb22+dn3_dy(i)*sb32
237 . -dn4_dz(i)*sc22+dn4_dy(i)*sc32
238 . -dn5_dz(i)*sd22+dn5_dy(i)*sd32
239 . -dn6_dz(i)*se22+dn6_dy(i)*se32
240 . -dn7_dz(i)*sf22+dn7_dy(i)*sf32
241 . -dn8_dz(i)*sg22+dn8_dy(i)*sg32
242 b13 = -dn2_dz(i)*sa23+dn2_dy(i)*sa33
243 . -dn3_dz(i)*sb23+dn3_dy(i)*sb33
244 . -dn4_dz(i)*sc23+dn4_dy(i)*sc33
245 . -dn5_dz(i)*sd23+dn5_dy(i)*sd33
246 . -dn6_dz(i)*se23+dn6_dy(i)*se33
247 . -dn7_dz(i)*sf23+dn7_dy(i)*sf33
248 . -dn8_dz(i)*sg23+dn8_dy(i)*sg33
249 b21 = dn2_dz(i)*sa11-dn2_dx(i)*sa31
250 . +dn3_dz(i)*sb11-dn3_dx(i)*sb31
251 . +dn4_dz(i)*sc11-dn4_dx(i)*sc31
252 . +dn5_dz(i)*sd11-dn5_dx(i)*sd31
253 . +dn6_dz(i)*se11-dn6_dx(i)*se31
254 . +dn7_dz(i)*sf11-dn7_dx(i)*sf31
255 . +dn8_dz(i)*sg11-dn8_dx(i)*sg31
256 b22 = dn2_dz(i)*sa12-dn2_dx(i)*sa32
257 . +dn3_dz(i)*sb12-dn3_dx(i)*sb32
258 . +dn4_dz(i)*sc12-dn4_dx(i)*sc32
259 . +dn5_dz(i)*sd12-dn5_dx(i)*sd32
260 . +dn6_dz(i)*se12-dn6_dx(i)*se32
261 . +dn7_dz(i)*sf12-dn7_dx(i)*sf32
262 . +dn8_dz(i)*sg12-dn8_dx(i)*sg32
263 b23 = dn2_dz(i)*sa13-dn2_dx(i)*sa33
264 . +dn3_dz(i)*sb13-dn3_dx(i)*sb33
265 . +dn4_dz(i)*sc13-dn4_dx(i)*sc33
266 . +dn5_dz(i)*sd13-dn5_dx(i)*sd33
267 . +dn6_dz(i)*se13-dn6_dx(i)*se33
268 . +dn7_dz(i)*sf13-dn7_dx(i)*sf33
269 . +dn8_dz(i)*sg13-dn8_dx(i)*sg33
270 b31 = -dn2_dy(i)*sa11+dn2_dx(i)*sa21
271 . -dn3_dy(i)*sb11+dn3_dx(i)*sb21
272 . -dn4_dy(i)*sc11+dn4_dx(i)*sc21
273 . -dn5_dy(i)*sd11+dn5_dx(i)*sd21
274 . -dn6_dy(i)*se11+dn6_dx(i)*se21
275 . -dn7_dy(i)*sf11+dn7_dx(i)*sf21
276 . -dn8_dy(i)*sg11+dn8_dx(i)*sg21
277 b32 = -dn2_dy(i)*sa12+dn2_dx(i)*sa22
278 . -dn3_dy(i)*sb12+dn3_dx(i)*sb22
279 . -dn4_dy(i)*sc12+dn4_dx(i)*sc22
280 . -dn5_dy(i)*sd12+dn5_dx(i)*sd22
281 . -dn6_dy(i)*se12+dn6_dx(i)*se22
282 . -dn7_dy(i)*sf12+dn7_dx(i)*sf22
283 . -dn8_dy(i)*sg12+dn8_dx(i)*sg22
284 b33 = -dn2_dy(i)*sa13+dn2_dx(i)*sa23
285 . -dn3_dy(i)*sb13+dn3_dx(i)*sb23
286 . -dn4_dy(i)*sc13+dn4_dx(i)*sc23
287 . -dn5_dy(i)*sd13+dn5_dx(i)*sd23
288 . -dn6_dy(i)*se13+dn6_dx(i)*se23
289 . -dn7_dy(i)*sf13+dn7_dx(i)*sf23
290 . -dn8_dy(i)*sg13+dn8_dx(i)*sg23
291 bb = b11*(b22*b33-b23*b32)
292 . -b12*(b21*b33-b31*b23)
293 . +b13*(b21*b32-b31*b22)
294 IF (bb /= zero) bb = one/bb
295 bi11 = (b22*b33-b32*b23)*bb
296 bi21 = (b31*b23-b21*b33)*bb
297 bi31 = (b21*b32-b31*b22)*bb
298 bi12 = (b13*b32-b12*b33)*bb
299 bi22 = (b11*b33-b13*b31)*bb
300 bi32 = (b12*b31-b11*b32)*bb
301 bi13 = (b12*b23-b13*b22)*bb
302 bi23 = (b21*b13-b11*b23)*bb
303 bi33 = (b11*b22-b12*b21)*bb
304
305 s11 = -dn1_dz(i)*r(1,2,i)+dn1_dy(i)*r(1,3,i)
306 s21 = dn1_dz(i)*r(1,1,i)-dn1_dx(i)*r(1,3,i)
307 s31 = -dn1_dy(i)*r(1,1,i)+dn1_dx(i)*r(1,2,i)
308 s12 = -dn1_dz(i)*r(2,2,i)+dn1_dy(i)*r(2,3,i)
309 s22 = dn1_dz(i)*r(2,1,i)-dn1_dx(i)*r(2,3,i)
310 s32 = -dn1_dy(i)*r(2,1,i)+dn1_dx(i)*r(2,2,i)
311 s13 = -dn1_dz(i)*r(3,2,i)+dn1_dy(i)*r(3,3,i)
312 s23 = dn1_dz(i)*r(3,1,i)-dn1_dx(i)*r(3,3,i)
313 s33 = -dn1_dy(i)*r(3,1,i)+dn1_dx(i)*r(3,2,i)
314
315 v1(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
316 v1(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
317 v1(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
318 v1(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
319 v1(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
320 v1(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
321 v1(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
322 v1(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
323 v1(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
324
325 s11 = -dn2_dz(i)*r(1,2,i)+dn2_dy(i)*r(1,3,i)
326 s21 = dn2_dz(i)*r(1,1,i)-dn2_dx(i)*r(1,3,i)
327 s31 = -dn2_dy(i)*r(1,1,i)+dn2_dx(i)*r(1,2,i)
328 s12 = -dn2_dz(i)*r(2,2,i)+dn2_dy(i)*r(2,3,i)
329 s22 = dn2_dz(i)*r(2,1,i)-dn2_dx(i)*r(2,3,i)
330 s32 = -dn2_dy(i)*r(2,1,i)+dn2_dx(i)*r(2,2,i)
331 s13 = -dn2_dz(i)*r(3,2,i)+dn2_dy(i)*r(3,3,i)
332 s23 = dn2_dz(i)*r(3,1,i)-dn2_dx(i)*r(3,3,i)
333 s33 = -dn2_dy(i)*r(3,1,i)+dn2_dx(i)*r(3,2,i)
334
335 v2(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
336 v2(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
337 v2(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
338 v2(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
339 v2(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
340 v2(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
341 v2(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
342 v2(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
343 v2(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
344
345 s11 = -dn3_dz(i)*r(1,2,i)+dn3_dy(i)*r(1,3,i)
346 s21 = dn3_dz(i)*r(1,1,i)-dn3_dx(i)*r(1,3,i)
347 s31 = -dn3_dy(i)*r(1,1,i)+dn3_dx(i)*r(1,2,i)
348 s12 = -dn3_dz(i)*r(2,2,i)+dn3_dy(i)*r(2,3,i)
349 s22 = dn3_dz(i)*r(2,1,i)-dn3_dx(i)*r(2,3,i)
350 s32 = -dn3_dy(i)*r(2,1,i)+dn3_dx(i)*r(2,2,i)
351 s13 = -dn3_dz(i)*r(3,2,i)+dn3_dy(i)*r(3,3,i)
352 s23 = dn3_dz(i)*r(3,1,i)-dn3_dx(i)*r(3,3,i)
353 s33 = -dn3_dy(i)*r(3,1,i)+dn3_dx(i)*r(3,2,i)
354
355 v3(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
356 v3(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
357 v3(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
358 v3(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
359 v3(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
360 v3(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
361 v3(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
362 v3(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
363 v3(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
364
365 s11 = -dn4_dz(i)*r(1,2,i)+dn4_dy(i)*r(1,3,i)
366 s21 = dn4_dz(i)*r(1,1,i)-dn4_dx(i)*r(1,3,i)
367 s31 = -dn4_dy(i)*r(1,1,i)+dn4_dx(i)*r(1,2,i)
368 s12 = -dn4_dz(i)*r(2,2,i)+dn4_dy(i)*r(2,3,i)
369 s22 = dn4_dz(i)*r(2,1,i)-dn4_dx(i)*r(2,3,i)
370 s32 = -dn4_dy(i)*r(2,1,i)+dn4_dx(i)*r(2,2,i)
371 s13 = -dn4_dz(i)*r(3,2,i)+dn4_dy(i)*r(3,3,i)
372 s23 = dn4_dz(i)*r(3,1,i)-dn4_dx(i)*r(3,3,i)
373 s33 = -dn4_dy(i)*r(3,1,i)+dn4_dx(i)*r(3,2,i)
374
375 v4(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
376 v4(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
377 v4(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
378 v4(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
379 v4(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
380 v4(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
381 v4(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
382 v4(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
383 v4(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
384
385 s11 = -dn5_dz(i)*r(1,2,i)+dn5_dy(i)*r(1,3,i)
386 s21 = dn5_dz(i)*r(1,1,i)-dn5_dx(i)*r(1,3,i)
387 s31 = -dn5_dy(i)*r(1,1,i)+dn5_dx(i)*r(1,2,i)
388 s12 = -dn5_dz(i)*r(2,2,i)+dn5_dy(i)*r(2,3,i)
389 s22 = dn5_dz(i)*r(2,1,i)-dn5_dx(i)*r(2,3,i)
390 s32 = -dn5_dy(i)*r(2,1,i)+dn5_dx(i)*r(2,2,i)
391 s13 = -dn5_dz(i)*r(3,2,i)+dn5_dy(i)*r(3,3,i)
392 s23 = dn5_dz(i)*r(3,1,i)-dn5_dx(i)*r(3,3,i)
393 s33 = -dn5_dy(i)*r(3,1,i)+dn5_dx(i)*r(3,2,i)
394
395 v5(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
396 v5(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
397 v5(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
398 v5(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
399 v5(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
400 v5(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
401 v5(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
402 v5(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
403 v5(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
404
405 s11 = -dn6_dz(i)*r(1,2,i)+dn6_dy(i)*r(1,3,i)
406 s21 = dn6_dz(i)*r(1,1,i)-dn6_dx(i)*r(1,3,i)
407 s31 = -dn6_dy(i)*r(1,1,i)+dn6_dx(i)*r(1,2,i)
408 s12 = -dn6_dz(i)*r(2,2,i)+dn6_dy(i)*r(2,3,i)
409 s22 = dn6_dz(i)*r(2,1,i)-dn6_dx(i)*r(2,3,i)
410 s32 = -dn6_dy(i)*r(2,1,i)+dn6_dx(i)*r(2,2,i)
411 s13 = -dn6_dz(i)*r(3,2,i)+dn6_dy(i)*r(3,3,i)
412 s23 = dn6_dz(i)*r(3,1,i)-dn6_dx(i)*r(3,3,i)
413 s33 = -dn6_dy(i)*r(3,1,i)+dn6_dx(i)*r(3,2,i)
414
415 v6(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
416 v6(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
417 v6(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
418 v6(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
419 v6(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
420 v6(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
421 v6(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
422 v6(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
423 v6(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
424
425 s11 = -dn7_dz(i)*r(1,2,i)+dn7_dy(i)*r(1,3,i)
426 s21 = dn7_dz(i)*r(1,1,i)-dn7_dx(i)*r(1,3,i)
427 s31 = -dn7_dy(i)*r(1,1,i)+dn7_dx(i)*r(1,2,i)
428 s12 = -dn7_dz(i)*r(2,2,i)+dn7_dy(i)*r(2,3,i)
429 s22 = dn7_dz(i)*r(2,1,i)-dn7_dx(i)*r(2,3,i)
430 s32 = -dn7_dy(i)*r(2,1,i)+dn7_dx(i)*r(2,2,i)
431 s13 = -dn7_dz(i)*r(3,2,i)+dn7_dy(i)*r(3,3,i)
432 s23 = dn7_dz(i)*r(3,1,i)-dn7_dx(i)*r(3,3,i)
433 s33 = -dn7_dy(i)*r(3,1,i)+dn7_dx(i)*r(3,2,i)
434
435 v7(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
436 v7(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
437 v7(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
438 v7(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
439 v7(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
440 v7(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
441 v7(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
442 v7(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
443 v7(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
444
445 s11 = -dn8_dz(i)*r(1,2,i)+dn8_dy(i)*r(1,3,i)
446 s21 = dn8_dz(i)*r(1,1,i)-dn8_dx(i)*r(1,3,i)
447 s31 = -dn8_dy(i)*r(1,1,i)+dn8_dx(i)*r(1,2,i)
448 s12 = -dn8_dz(i)*r(2,2,i)+dn8_dy(i)*r(2,3,i)
449 s22 = dn8_dz(i)*r(2,1,i)-dn8_dx(i)*r(2,3,i)
450 s32 = -dn8_dy(i)*r(2,1,i)+dn8_dx(i)*r(2,2,i)
451 s13 = -dn8_dz(i)*r(3,2,i)+dn8_dy(i)*r(3,3,i)
452 s23 = dn8_dz(i)*r(3,1,i)-dn8_dx(i)*r(3,3,i)
453 s33 = -dn8_dy(i)*r(3,1,i)+dn8_dx(i)*r(3,2,i)
454
455 v8(i,1,1) = -(bi11*s11+bi12*s21+bi13*s31)
456 v8(i,2,1) = -(bi21*s11+bi22*s21+bi23*s31)
457 v8(i,3,1) = -(bi31*s11+bi32*s21+bi33*s31)
458 v8(i,1,2) = -(bi11*s12+bi12*s22+bi13*s32)
459 v8(i,2,2) = -(bi21*s12+bi22*s22+bi23*s32)
460 v8(i,3,2) = -(bi31*s12+bi32*s22+bi33*s32)
461 v8(i,1,3) = -(bi11*s13+bi12*s23+bi13*s33)
462 v8(i,2,3) = -(bi21*s13+bi22*s23+bi23*s33)
463 v8(i,3,3) = -(bi31*s13+bi32*s23+bi33*s33)
464
465 trm(i,1,1) = zero
466 trm(i,2,1) = zero
467 trm(i,3,1) = zero
468 trm(i,1,2) = zero
469 trm(i,2,2) = zero
470 trm(i,3,2) = zero
471 trm(i,1,3) = zero
472 trm(i,2,3) = zero
473 trm(i,3,3) = zero
474
475 trm(i,4,1) = sa11*v1(i,1,1)+sa12*v1(i,2,1)+sa13*v1(i,3,1)
476 trm(i,5,1) = sa21*v1(i,1,1)+sa22*v1(i,2,1)+sa23*v1(i,3,1)
477 trm(i,6,1) = sa31*v1(i,1,1)+sa32*v1(i,2,1)+sa33*v1(i,3,1)
478 trm(i,4,2) = sa11*v1(i,1,2)+sa12*v1(i,2,2)+sa13*v1(i,3,2)
479 trm(i,5,2) = sa21*v1(i,1,2)+sa22*v1(i,2,2)+sa23*v1(i,3,2)
480 trm(i,6,2) = sa31*v1(i,1,2)+sa32*v1(i,2,2)+sa33*v1(i,3,2)
481 trm(i,4,3) = sa11*v1(i,1,3)+sa12*v1(i,2,3)+sa13*v1(i,3,3)
482 trm(i,5,3) = sa21*v1(i,1,3)+sa22*v1(i,2,3)+sa23*v1(i,3,3)
483 trm(i,6,3) = sa31*v1(i,1,3)+sa32*v1(i,2,3)+sa33*v1(i,3,3)
484
485 trm(i,7,1) = sb11*v1(i,1,1)+sb12*v1(i,2,1)+sb13*v1(i,3,1)
486 trm(i,8,1) = sb21*v1(i,1,1)+sb22*v1(i,2,1)+sb23*v1(i,3,1)
487 trm(i,9,1) = sb31*v1(i,1,1)+sb32*v1(i,2,1)+sb33*v1(i,3,1)
488 trm(i,7,2) = sb11*v1(i,1,2)+sb12*v1(i,2,2)+sb13*v1(i,3,2)
489 trm(i,8,2) = sb21*v1(i,1,2)+sb22*v1(i,2,2)+sb23*v1(i,3,2)
490 trm(i,9,2) = sb31*v1(i,1,2)+sb32*v1(i,2,2)+sb33*v1(i,3,2)
491 trm(i,7,3) = sb11*v1(i,1,3)+sb12*v1(i,2,3)+sb13*v1(i,3,3)
492 trm(i,8,3) = sb21*v1(i,1,3)+sb22*v1(i,2,3)+sb23*v1(i,3,3)
493 trm(i,9,3) = sb31*v1(i,1,3)+sb32*v1(i,2,3)+sb33*v1(i,3,3)
494
495 trm(i,10,1) = sc11*v1(i,1,1)+sc12*v1(i,2,1)+sc13*v1(i,3,1)
496 trm(i,11,1) = sc21*v1(i,1,1)+sc22*v1(i,2,1)+sc23*v1(i,3,1)
497 trm(i,12,1) = sc31*v1(i,1,1)+sc32*v1(i,2,1)+sc33*v1(i,3,1)
498 trm(i,10,2) = sc11*v1(i,1,2)+sc12*v1(i,2,2)+sc13*v1(i,3,2)
499 trm(i,11,2) = sc21*v1(i,1,2)+sc22*v1(i,2,2)+sc23*v1(i,3,2)
500 trm(i,12,2) = sc31*v1(i,1,2)+sc32*v1(i,2,2)+sc33*v1(i,3,2)
501 trm(i,10,3) = sc11*v1(i,1,3)+sc12*v1(i,2,3)+sc13*v1(i,3,3)
502 trm(i,11,3) = sc21*v1(i,1,3)+sc22*v1(i,2,3)+sc23*v1(i,3,3)
503 trm(i,12,3) = sc31*v1(i,1,3)+sc32*v1(i,2,3)+sc33*v1(i,3,3)
504
505 trm(i,13,1) = sd11*v1(i,1,1)+sd12*v1(i,2,1)+sd13*v1(i,3,1)
506 trm(i,14,1) = sd21*v1(i,1,1)+sd22*v1(i,2,1)+sd23*v1(i,3,1)
507 trm(i,15,1) = sd31*v1(i,1,1)+sd32*v1(i,2,1)+sd33*v1(i,3,1)
508 trm(i,13,2) = sd11*v1(i,1,2)+sd12*v1(i,2,2)+sd13*v1(i,3,2)
509 trm(i,14,2) = sd21*v1(i,1,2)+sd22*v1(i,2,2)+sd23*v1(i,3,2)
510 trm(i,15,2) = sd31*v1(i,1,2)+sd32*v1(i,2,2)+sd33*v1(i,3,2)
511 trm(i,13,3) = sd11*v1(i,1,3)+sd12*v1(i,2,3)+sd13*v1(i,3,3)
512 trm(i,14,3) = sd21*v1(i,1,3)+sd22*v1(i,2,3)+sd23*v1(i,3,3)
513 trm(i,15,3) = sd31*v1(i,1,3)+sd32*v1(i,2,3)+sd33*v1(i,3,3)
514
515 trm(i,16,1) = se11*v1(i,1,1)+se12*v1(i,2,1)+se13*v1(i,3,1)
516 trm(i,17,1) = se21*v1(i,1,1)+se22*v1(i,2,1)+se23*v1(i,3,1)
517 trm(i,18,1) = se31*v1(i,1,1)+se32*v1(i,2,1)+se33*v1(i,3,1)
518 trm(i,16,2) = se11*v1(i,1,2)+se12*v1(i,2,2)+se13*v1(i,3,2)
519 trm(i,17,2) = se21*v1(i,1,2)+se22*v1(i,2,2)+se23*v1(i,3,2)
520 trm(i,18,2) = se31*v1(i,1,2)+se32*v1(i,2,2)+se33*v1(i,3,2)
521 trm(i,16,3) = se11*v1(i,1,3)+se12*v1(i,2,3)+se13*v1(i,3,3)
522 trm(i,17,3) = se21*v1(i,1,3)+se22*v1(i,2,3)+se23*v1(i,3,3)
523 trm(i,18,3) = se31*v1(i,1,3)+se32*v1(i,2,3)+se33*v1(i,3,3)
524
525 trm(i,19,1) = sf11*v1(i,1,1)+sf12*v1(i,2,1)+sf13*v1(i,3,1)
526 trm(i,20,1) = sf21*v1(i,1,1)+sf22*v1(i,2,1)+sf23*v1(i,3,1)
527 trm(i,21,1) = sf31*v1(i,1,1)+sf32*v1(i,2,1)+sf33*v1(i,3,1)
528 trm(i,19,2) = sf11*v1(i,1,2)+sf12*v1(i,2,2)+sf13*v1(i,3,2)
529 trm(i,20,2) = sf21*v1(i,1,2)+sf22*v1(i,2,2)+sf23*v1(i,3,2)
530 trm(i,21,2) = sf31*v1(i,1,2)+sf32*v1(i,2,2)+sf33*v1(i,3,2)
531 trm(i,19,3) = sf11*v1(i,1,3)+sf12*v1(i,2,3)+sf13*v1(i,3,3)
532 trm(i,20,3) = sf21*v1(i,1,3)+sf22*v1(i,2,3)+sf23*v1(i,3,3)
533 trm(i,21,3) = sf31*v1(i,1,3)+sf32*v1(i,2,3)+sf33*v1(i,3,3)
534
535 trm(i,22,1) = sg11*v1(i,1,1)+sg12*v1(i,2,1)+sg13*v1(i,3,1)
536 trm(i,23,1) = sg21*v1(i,1,1)+sg22*v1(i,2,1)+sg23*v1(i,3,1)
537 trm(i,24,1) = sg31*v1(i,1,1)+sg32*v1(i,2,1)+sg33*v1(i,3,1)
538 trm(i,22,2) = sg11*v1(i,1,2)+sg12*v1(i,2,2)+sg13*v1(i,3,2)
539 trm(i,23,2) = sg21*v1(i,1,2)+sg22*v1(i,2,2)+sg23*v1(i,3,2)
540 trm(i,24,2) = sg31*v1(i,1,2)+sg32*v1(i,2,2)+sg33*v1(i,3,2)
541 trm(i,22,3) = sg11*v1(i,1,3)+sg12*v1(i,2,3)+sg13*v1(i,3,3)
542 trm(i,23,3) = sg21*v1(i,1,3)+sg22*v1(i,2,3)+sg23*v1(i,3,3)
543 trm(i,24,3) = sg31*v1(i,1,3)+sg32*v1(i,2,3)+sg33*v1(i,3,3)
544
545 trm(i,1,4) = zero
546 trm(i,2,4) = zero
547 trm(i,3,4) = zero
548 trm(i,1,5) = zero
549 trm(i,2,5) = zero
550 trm(i,3,5) = zero
551 trm(i,1,6) = zero
552 trm(i,2,6) = zero
553 trm(i,3,6) = zero
554
555 trm(i,4,4) = sa11*v2(i,1,1)+sa12*v2(i,2,1)+sa13*v2(i,3,1)
556 trm(i,5,4) = sa21*v2(i,1,1)+sa22*v2(i,2,1)+sa23*v2(i,3,1)
557 trm(i,6,4) = sa31*v2(i,1,1)+sa32*v2(i,2,1)+sa33*v2(i,3,1)
558 trm(i,4,5) = sa11*v2(i,1,2)+sa12*v2(i,2,2)+sa13*v2(i,3,2)
559 trm(i,5,5) = sa21*v2(i,1,2)+sa22*v2(i,2,2)+sa23*v2(i,3,2)
560 trm(i,6,5) = sa31*v2(i,1,2)+sa32*v2(i,2,2)+sa33*v2(i,3,2)
561 trm(i,4,6) = sa11*v2(i,1,3)+sa12*v2(i,2,3)+sa13*v2(i,3,3)
562 trm(i,5,6) = sa21*v2(i,1,3)+sa22*v2(i,2,3)+sa23*v2(i,3,3)
563 trm(i,6,6) = sa31*v2(i,1,3)+sa32*v2(i,2,3)+sa33*v2(i,3,3)
564
565 trm(i,7,4) = sb11*v2(i,1,1)+sb12*v2(i,2,1)+sb13*v2(i,3,1)
566 trm(i,8,4) = sb21*v2(i,1,1)+sb22*v2(i,2,1)+sb23*v2(i,3,1)
567 trm(i,9,4) = sb31*v2(i,1,1)+sb32*v2(i,2,1)+sb33*v2(i,3,1)
568 trm(i,7,5) = sb11*v2(i,1,2)+sb12*v2(i,2,2)+sb13*v2(i,3,2)
569 trm(i,8,5) = sb21*v2(i,1,2)+sb22*v2(i,2,2)+sb23*v2(i,3,2)
570 trm(i,9,5) = sb31*v2(i,1,2)+sb32*v2(i,2,2)+sb33*v2(i,3,2)
571 trm(i,7,6) = sb11*v2(i,1,3)+sb12*v2(i,2,3)+sb13*v2(i,3,3)
572 trm(i,8,6) = sb21*v2(i,1,3)+sb22*v2(i,2,3)+sb23*v2(i,3,3)
573 trm(i,9,6) = sb31*v2(i,1,3)+sb32*v2(i,2,3)+sb33*v2(i,3,3)
574
575 trm(i,10,4) = sc11*v2(i,1,1)+sc12*v2(i,2,1)+sc13*v2(i,3,1)
576 trm(i,11,4) = sc21*v2(i,1,1)+sc22*v2(i,2,1)+sc23*v2(i,3,1)
577 trm(i,12,4) = sc31*v2(i,1,1)+sc32*v2(i,2,1)+sc33*v2(i,3,1)
578 trm(i,10,5) = sc11*v2(i,1,2)+sc12*v2(i,2,2)+sc13*v2(i,3,2)
579 trm(i,11,5) = sc21*v2(i,1,2)+sc22*v2(i,2,2)+sc23*v2(i,3,2)
580 trm(i,12,5) = sc31*v2(i,1,2)+sc32*v2(i,2,2)+sc33*v2(i,3,2)
581 trm(i,10,6) = sc11*v2(i,1,3)+sc12*v2(i,2,3)+sc13*v2(i,3,3)
582 trm(i,11,6) = sc21*v2(i,1,3)+sc22*v2(i,2,3)+sc23*v2(i,3,3)
583 trm(i,12,6) = sc31*v2(i,1,3)+sc32*v2(i,2,3)+sc33*v2(i,3,3)
584
585 trm(i,13,4) = sd11*v2(i,1,1)+sd12*v2(i,2,1)+sd13*v2(i,3,1)
586 trm(i,14,4) = sd21*v2(i,1,1)+sd22*v2(i,2,1)+sd23*v2(i,3,1)
587 trm(i,15,4) = sd31*v2(i,1,1)+sd32*v2(i,2,1)+sd33*v2(i,3,1)
588 trm(i,13,5) = sd11*v2(i,1,2)+sd12*v2(i,2,2)+sd13*v2(i,3,2)
589 trm(i,14,5) = sd21*v2(i,1,2)+sd22*v2(i,2,2)+sd23*v2(i,3,2)
590 trm(i,15,5) = sd31*v2(i,1,2)+sd32*v2(i,2,2)+sd33*v2(i,3,2)
591 trm(i,13,6) = sd11*v2(i,1,3)+sd12*v2(i,2,3)+sd13*v2(i,3,3)
592 trm(i,14,6) = sd21*v2(i,1,3)+sd22*v2(i,2,3)+sd23*v2(i,3,3)
593 trm(i,15,6) = sd31*v2(i,1,3)+sd32*v2(i,2,3)+sd33*v2(i,3,3)
594
595 trm(i,16,4) = se11*v2(i,1,1)+se12*v2(i,2,1)+se13*v2(i,3,1)
596 trm(i,17,4) = se21*v2(i,1,1)+se22*v2(i,2,1)+se23*v2(i,3,1)
597 trm(i,18,4) = se31*v2(i,1,1)+se32*v2(i,2,1)+se33*v2(i,3,1)
598 trm(i,16,5) = se11*v2(i,1,2)+se12*v2(i,2,2)+se13*v2(i,3,2)
599 trm(i,17,5) = se21*v2(i,1,2)+se22*v2(i,2,2)+se23*v2(i,3,2)
600 trm(i,18,5) = se31*v2(i,1,2)+se32*v2(i,2,2)+se33*v2(i,3,2)
601 trm(i,16,6) = se11*v2(i,1,3)+se12*v2(i,2,3)+se13*v2(i,3,3)
602 trm(i,17,6) = se21*v2(i,1,3)+se22*v2(i,2,3)+se23*v2(i,3,3)
603 trm(i,18,6) = se31*v2(i,1,3)+se32*v2(i,2,3)+se33*v2(i,3,3)
604
605 trm(i,19,4) = sf11*v2(i,1,1)+sf12*v2(i,2,1)+sf13*v2(i,3,1)
606 trm(i,20,4) = sf21*v2(i,1,1)+sf22*v2(i,2,1)+sf23*v2(i,3,1)
607 trm(i,21,4) = sf31*v2(i,1,1)+sf32*v2(i,2,1)+sf33*v2(i,3,1)
608 trm(i,19,5) = sf11*v2(i,1,2)+sf12*v2(i,2,2)+sf13*v2(i,3,2)
609 trm(i,20,5) = sf21*v2(i,1,2)+sf22*v2(i,2,2)+sf23*v2(i,3,2)
610 trm(i,21,5) = sf31*v2(i,1,2)+sf32*v2(i,2,2)+sf33*v2(i,3,2)
611 trm(i,19,6) = sf11*v2(i,1,3)+sf12*v2(i,2,3)+sf13*v2(i,3,3)
612 trm(i,20,6) = sf21*v2(i,1,3)+sf22*v2(i,2,3)+sf23*v2(i,3,3)
613 trm(i,21,6) = sf31*v2(i,1,3)+sf32*v2(i,2,3)+sf33*v2(i,3,3)
614
615 trm(i,22,4) = sg11*v2(i,1,1)+sg12*v2(i,2,1)+sg13*v2(i,3,1)
616 trm(i,23,4) = sg21*v2(i,1,1)+sg22*v2(i,2,1)+sg23*v2(i,3,1)
617 trm(i,24,4) = sg31*v2(i,1,1)+sg32*v2(i,2,1)+sg33*v2(i,3,1)
618 trm(i,22,5) = sg11*v2(i,1,2)+sg12*v2(i,2,2)+sg13*v2(i,3,2)
619 trm(i,23,5) = sg21*v2(i,1,2)+sg22*v2(i,2,2)+sg23*v2(i,3,2)
620 trm(i,24,5) = sg31*v2(i,1,2)+sg32*v2(i,2,2)+sg33*v2(i,3,2)
621 trm(i,22,6) = sg11*v2(i,1,3)+sg12*v2(i,2,3)+sg13*v2(i,3,3)
622 trm(i,23,6) = sg21*v2(i,1,3)+sg22*v2(i,2,3)+sg23*v2(i,3,3)
623 trm(i,24,6) = sg31*v2(i,1,3)+sg32*v2(i,2,3)+sg33*v2(i,3,3)
624
625 trm(i,1,7) = zero
626 trm(i,2,7) = zero
627 trm(i,3,7) = zero
628 trm(i,1,8) = zero
629 trm(i,2,8) = zero
630 trm(i,3,8) = zero
631 trm(i,1,9) = zero
632 trm(i,2,9) = zero
633 trm(i,3,9) = zero
634
635 trm(i,4,7) = sa11*v3(i,1,1)+sa12*v3(i,2,1)+sa13*v3(i,3,1)
636 trm(i,5,7) = sa21*v3(i,1,1)+sa22*v3(i,2,1)+sa23*v3(i,3,1)
637 trm(i,6,7) = sa31*v3(i,1,1)+sa32*v3(i,2,1)+sa33*v3(i,3,1)
638 trm(i,4,8) = sa11*v3(i,1,2)+sa12*v3(i,2,2)+sa13*v3(i,3,2)
639 trm(i,5,8) = sa21*v3(i,1,2)+sa22*v3(i,2,2)+sa23*v3(i,3,2)
640 trm(i,6,8) = sa31*v3(i,1,2)+sa32*v3(i,2,2)+sa33*v3(i,3,2)
641 trm(i,4,9) = sa11*v3(i,1,3)+sa12*v3(i,2,3)+sa13*v3(i,3,3)
642 trm(i,5,9) = sa21*v3(i,1,3)+sa22*v3(i,2,3)+sa23*v3(i,3,3)
643 trm(i,6,9) = sa31*v3(i,1,3)+sa32*v3(i,2,3)+sa33*v3(i,3,3)
644
645 trm(i,7,7) = sb11*v3(i,1,1)+sb12*v3(i,2,1)+sb13*v3(i,3,1)
646 trm(i,8,7) = sb21*v3(i,1,1)+sb22*v3(i,2,1)+sb23*v3(i,3,1)
647 trm(i,9,7) = sb31*v3(i,1,1)+sb32*v3(i,2,1)+sb33*v3(i,3,1)
648 trm(i,7,8) = sb11*v3(i,1,2)+sb12*v3(i,2,2)+sb13*v3(i,3,2)
649 trm(i,8,8) = sb21*v3(i,1,2)+sb22*v3(i,2,2)+sb23*v3(i,3,2)
650 trm(i,9,8) = sb31*v3(i,1,2)+sb32*v3(i,2,2)+sb33*v3(i,3,2)
651 trm(i,7,9) = sb11*v3(i,1,3)+sb12*v3(i,2,3)+sb13*v3(i,3,3)
652 trm(i,8,9) = sb21*v3(i,1,3)+sb22*v3(i,2,3)+sb23*v3(i,3,3)
653 trm(i,9,9) = sb31*v3(i,1,3)+sb32*v3(i,2,3)+sb33*v3(i,3,3)
654
655 trm(i,10,7) = sc11*v3(i,1,1)+sc12*v3(i,2,1)+sc13*v3(i,3,1)
656 trm(i,11,7) = sc21*v3(i,1,1)+sc22*v3(i,2,1)+sc23*v3(i,3,1)
657 trm(i,12,7) = sc31*v3(i,1,1)+sc32*v3(i,2,1)+sc33*v3(i,3,1)
658 trm(i,10,8) = sc11*v3(i,1,2)+sc12*v3(i,2,2)+sc13*v3(i,3,2)
659 trm(i,11,8) = sc21*v3(i,1,2)+sc22*v3(i,2,2)+sc23*v3(i,3,2)
660 trm(i,12,8) = sc31*v3(i,1,2)+sc32*v3(i,2,2)+sc33*v3(i,3,2)
661 trm(i,10,9) = sc11*v3(i,1,3)+sc12*v3(i,2,3)+sc13*v3(i,3,3)
662 trm(i,11,9) = sc21*v3(i,1,3)+sc22*v3(i,2,3)+sc23*v3(i,3,3)
663 trm(i,12,9) = sc31*v3(i,1,3)+sc32*v3(i,2,3)+sc33*v3(i,3,3)
664
665 trm(i,13,7) = sd11*v3(i,1,1)+sd12*v3(i,2,1)+sd13*v3(i,3,1)
666 trm(i,14,7) = sd21*v3(i,1,1)+sd22*v3(i,2,1)+sd23*v3(i,3,1)
667 trm(i,15,7) = sd31*v3(i,1,1)+sd32*v3(i,2,1)+sd33*v3(i,3,1)
668 trm(i,13,8) = sd11*v3(i,1,2)+sd12*v3(i,2,2)+sd13*v3(i,3,2)
669 trm(i,14,8) = sd21*v3(i,1,2)+sd22*v3(i,2,2)+sd23*v3(i,3,2)
670 trm(i,15,8) = sd31*v3(i,1,2)+sd32*v3(i,2,2)+sd33*v3(i,3,2)
671 trm(i,13,9) = sd11*v3(i,1,3)+sd12*v3(i,2,3)+sd13*v3(i,3,3)
672 trm(i,14,9) = sd21*v3(i,1,3)+sd22*v3(i,2,3)+sd23*v3(i,3,3)
673 trm(i,15,9) = sd31*v3(i,1,3)+sd32*v3(i,2,3)+sd33*v3(i,3,3)
674
675 trm(i,16,7) = se11*v3(i,1,1)+se12*v3(i,2,1)+se13*v3(i,3,1)
676 trm(i,17,7) = se21*v3(i,1,1)+se22*v3(i,2,1)+se23*v3(i,3,1)
677 trm(i,18,7) = se31*v3(i,1,1)+se32*v3(i,2,1)+se33*v3(i,3,1)
678 trm(i,16,8) = se11*v3(i,1,2)+se12*v3(i,2,2)+se13*v3(i,3,2)
679 trm(i,17,8) = se21*v3(i,1,2)+se22*v3(i,2,2)+se23*v3(i,3,2)
680 trm(i,18,8) = se31*v3(i,1,2)+se32*v3(i,2,2)+se33*v3(i,3,2)
681 trm(i,16,9) = se11*v3(i,1,3)+se12*v3(i,2,3)+se13*v3(i,3,3)
682 trm(i,17,9) = se21*v3(i,1,3)+se22*v3(i,2,3)+se23*v3(i,3,3)
683 trm(i,18,9) = se31*v3(i,1,3)+se32*v3(i,2,3)+se33*v3(i,3,3)
684
685 trm(i,19,7) = sf11*v3(i,1,1)+sf12*v3(i,2,1)+sf13*v3(i,3,1)
686 trm(i,20,7) = sf21*v3(i,1,1)+sf22*v3(i,2,1)+sf23*v3(i,3,1)
687 trm(i,21,7) = sf31*v3(i,1,1)+sf32*v3(i,2,1)+sf33*v3(i,3,1)
688 trm(i,19,8) = sf11*v3(i,1,2)+sf12*v3(i,2,2)+sf13*v3(i,3,2)
689 trm(i,20,8) = sf21*v3(i,1,2)+sf22*v3(i,2,2)+sf23*v3(i,3,2)
690 trm(i,21,8) = sf31*v3(i,1,2)+sf32*v3(i,2,2)+sf33*v3(i,3,2)
691 trm(i,19,9) = sf11*v3(i,1,3)+sf12*v3(i,2,3)+sf13*v3(i,3,3)
692 trm(i,20,9) = sf21*v3(i,1,3)+sf22*v3(i,2,3)+sf23*v3(i,3,3)
693 trm(i,21,9) = sf31*v3(i,1,3)+sf32*v3(i,2,3)+sf33*v3(i,3,3)
694
695 trm(i,22,7) = sg11*v3(i,1,1)+sg12*v3(i,2,1)+sg13*v3(i,3,1)
696 trm(i,23,7) = sg21*v3(i,1,1)+sg22*v3(i,2,1)+sg23*v3(i,3,1)
697 trm(i,24,7) = sg31*v3(i,1,1)+sg32*v3(i,2,1)+sg33*v3(i,3,1)
698 trm(i,22,8) = sg11*v3(i,1,2)+sg12*v3(i,2,2)+sg13*v3(i,3,2)
699 trm(i,23,8) = sg21*v3(i,1,2)+sg22*v3(i,2,2)+sg23*v3(i,3,2)
700 trm(i,24,8) = sg31*v3(i,1,2)+sg32*v3(i,2,2)+sg33*v3(i,3,2)
701 trm(i,22,9) = sg11*v3(i,1,3)+sg12*v3(i,2,3)+sg13*v3(i,3,3)
702 trm(i,23,9) = sg21*v3(i,1,3)+sg22*v3(i,2,3)+sg23*v3(i,3,3)
703 trm(i,24,9) = sg31*v3(i,1,3)+sg32*v3(i,2,3)+sg33*v3(i,3,3)
704
705 trm(i,1,10) = zero
706 trm(i,2,10) = zero
707 trm(i,3,10) = zero
708 trm(i,1,11) = zero
709 trm(i,2,11) = zero
710 trm(i,3,11) = zero
711 trm(i,1,12) = zero
712 trm(i,2,12) = zero
713 trm(i,3,12) = zero
714
715 trm(i,4,10) = sa11*v4(i,1,1)+sa12*v4(i,2,1)+sa13*v4(i,3,1)
716 trm(i,5,10) = sa21*v4(i,1,1)+sa22*v4(i,2,1)+sa23*v4(i,3,1)
717 trm(i,6,10) = sa31*v4(i,1,1)+sa32*v4(i,2,1)+sa33*v4(i,3,1)
718 trm(i,4,11) = sa11*v4(i,1,2)+sa12*v4(i,2,2)+sa13*v4(i,3,2)
719 trm(i,5,11) = sa21*v4(i,1,2)+sa22*v4(i,2,2)+sa23*v4(i,3,2)
720 trm(i,6,11) = sa31*v4(i,1,2)+sa32*v4(i,2,2)+sa33*v4(i,3,2)
721 trm(i,4,12) = sa11*v4(i,1,3)+sa12*v4(i,2,3)+sa13*v4(i,3,3)
722 trm(i,5,12) = sa21*v4(i,1,3)+sa22*v4(i,2,3)+sa23*v4(i,3,3)
723 trm(i,6,12) = sa31*v4(i,1,3)+sa32*v4(i,2,3)+sa33*v4(i,3,3)
724
725 trm(i,7,10) = sb11*v4(i,1,1)+sb12*v4(i,2,1)+sb13*v4(i,3,1)
726 trm(i,8,10) = sb21*v4(i,1,1)+sb22*v4(i,2,1)+sb23*v4(i,3,1)
727 trm(i,9,10) = sb31*v4(i,1,1)+sb32*v4(i,2,1)+sb33*v4(i,3,1)
728 trm(i,7,11) = sb11*v4(i,1,2)+sb12*v4(i,2,2)+sb13*v4(i,3,2)
729 trm(i,8,11) = sb21*v4(i,1,2)+sb22*v4(i,2,2)+sb23*v4(i,3,2)
730 trm(i,9,11) = sb31*v4(i,1,2)+sb32*v4(i,2,2)+sb33*v4(i,3,2)
731 trm(i,7,12) = sb11*v4(i,1,3)+sb12*v4(i,2,3)+sb13*v4(i,3,3)
732 trm(i,8,12) = sb21*v4(i,1,3)+sb22*v4(i,2,3)+sb23*v4(i,3,3)
733 trm(i,9,12) = sb31*v4(i,1,3)+sb32*v4(i,2,3)+sb33*v4(i,3,3)
734
735 trm(i,10,10) = sc11*v4(i,1,1)+sc12*v4(i,2,1)+sc13*v4(i,3,1)
736 trm(i,11,10) = sc21*v4(i,1,1)+sc22*v4(i,2,1)+sc23*v4(i,3,1)
737 trm(i,12,10) = sc31*v4(i,1,1)+sc32*v4(i,2,1)+sc33*v4(i,3,1)
738 trm(i,10,11) = sc11*v4(i,1,2)+sc12*v4(i,2,2)+sc13*v4(i,3,2)
739 trm(i,11,11) = sc21*v4(i,1,2)+sc22*v4(i,2,2)+sc23*v4(i,3,2)
740 trm(i,12,11) = sc31*v4(i,1,2)+sc32*v4(i,2,2)+sc33*v4(i,3,2)
741 trm(i,10,12) = sc11*v4(i,1,3)+sc12*v4(i,2,3)+sc13*v4(i,3,3)
742 trm(i,11,12) = sc21*v4(i,1,3)+sc22*v4(i,2,3)+sc23*v4(i,3,3)
743 trm(i,12,12) = sc31*v4(i,1,3)+sc32*v4(i,2,3)+sc33*v4(i,3,3)
744
745 trm(i,13,10) = sd11*v4(i,1,1)+sd12*v4(i
746 trm(i,14,10) = sd21*v4(i,1,1)+sd22*v4(i,2,1)+sd23*v4(i,3,1)
747 trm(i,15,10) = sd31*v4(i,1,1)+sd32*v4(i,2,1)+sd33*v4(i,3,1)
748 trm(i,13,11) = sd11*v4(i,1,2)+sd12*v4(i,2,2)+sd13*v4(i,3,2)
749 trm(i,14,11) = sd21*v4(i,1,2)+sd22*v4(i,2,2)+sd23*v4(i,3,2)
750 trm(i,15,11) = sd31*v4(i,1,2)+sd32*v4(i,2,2)+sd33*v4(i,3,2)
751 trm(i,13,12) = sd11*v4(i,1,3)+sd12*v4(i,2,3)+sd13*v4(i,3,3)
752 trm(i,14,12) = sd21*v4(i,1,3)+sd22*v4(i,2,3)+sd23*v4(i,3,3)
753 trm(i,15,12) = sd31*v4(i,1,3)+sd32*v4(i,2,3)+sd33*v4(i,3,3)
754
755 trm(i,16,10) = se11*v4(i,1,1)+se12*v4(i,2,1)+se13*v4(i,3,1)
756 trm(i,17,10) = se21*v4(i,1,1)+se22*v4(i,2,1)+se23*v4(i,3,1)
757 trm(i,18,10) = se31*v4(i,1,1)+se32*v4(i,2,1)+se33*v4(i,3,1)
758 trm(i,16,11) = se11*v4(i,1,2)+se12*v4(i,2,2)+se13*v4(i,3,2)
759 trm(i,17,11) = se21*v4(i,1,2)+se22*v4(i,2,2)+se23*v4(i,3,2)
760 trm(i,18,11) = se31*v4(i,1,2)+se32*v4(i,2,2)+se33*v4(i,3,2)
761 trm(i,16,12) = se11*v4(i,1,3)+se12*v4(i,2,3)+se13*v4(i,3,3)
762 trm(i,17,12) = se21*v4(i,1,3)+se22*v4(i,2,3)+se23*v4(i,3,3)
763 trm(i,18,12) = se31*v4(i,1,3)+se32*v4(i,2,3)+se33*v4(i,3,3)
764
765 trm(i,19,10) = sf11*v4(i,1,1)+sf12*v4(i,2,1)+sf13*v4(i,3,1)
766 trm(i,20,10) = sf21*v4(i,1,1)+sf22*v4(i,2,1)+sf23*v4(i,3,1)
767 trm(i,21,10) = sf31*v4(i,1,1)+sf32*v4(i,2,1)+sf33*v4(i,3,1)
768 trm(i,19,11) = sf11*v4(i,1,2)+sf12*v4(i,2,2)+sf13*v4(i,3,2)
769 trm(i,20,11) = sf21*v4(i,1,2)+sf22*v4(i,2,2)+sf23*v4(i,3,2)
770 trm(i,21,11) = sf31*v4(i,1,2)+sf32*v4(i,2,2)+sf33*v4(i,3,2)
771 trm(i,19,12) = sf11*v4(i,1,3)+sf12*v4(i,2,3)+sf13*v4
772 trm(i,20,12) = sf21*v4(i,1,3)+sf22*v4(i,2,3)+sf23*v4(i,3,3)
773 trm(i,21,12) = sf31*v4(i,1,3)+sf32*v4(i,2,3)+sf33*v4(i,3,3)
774
775 trm(i,22,10) = sg11*v4(i,1,1)+sg12*v4(i,2,1)+sg13*v4(i,3,1)
776 trm(i,23,10) = sg21*v4(i,1,1)+sg22*v4(i,2,1)+sg23*v4(i,3,1)
777 trm(i,24,10) = sg31*v4(i,1,1)+sg32*v4(i,2,1)+sg33*v4(i,3,1)
778 trm(i,22,11) = sg11*v4(i,1,2)+sg12*v4(i,2,2)+sg13*v4(i,3,2)
779 trm(i,23,11) = sg21*v4(i,1,2)+sg22*v4(i,2,2)+sg23*v4(i,3,2)
780 trm(i,24,11) = sg31*v4(i,1,2)+sg32*v4(i,2,2)+sg33*v4(i,3,2)
781 trm(i,22,12) = sg11*v4(i,1,3)+sg12*v4(i,2,3)+sg13*v4(i,3,3)
782 trm(i,23,12) = sg21*v4(i,1,3)+sg22*v4(i,2,3)+sg23*v4(i,3,3)
783 trm(i,24,12) = sg31*v4(i,1,3)+sg32*v4(i,2,3)+sg33*v4(i,3,3)
784
785 trm(i,1,13) = zero
786 trm(i,2,13) = zero
787 trm(i,3,13) = zero
788 trm(i,1,14) = zero
789 trm(i,2,14) = zero
790 trm(i,3,14) = zero
791 trm(i,1,15) = zero
792 trm(i,2,15) = zero
793 trm(i,3,15) = zero
794
795 trm(i,4,13) = sa11*v5(i,1,1)+sa12*v5(i,2,1)+sa13*v5(i,3,1)
796 trm(i,5,13) = sa21*v5(i,1,1)+sa22*v5(i,2,1)+sa23*v5(i,3,1)
797 trm(i,6,13) = sa31*v5(i,1,1)+sa32*v5(i,2,1)+sa33*v5(i,3,1)
798 trm(i,4,14) = sa11*v5(i,1,2)+sa12*v5(i,2,2)+sa13*v5(i,3,2)
799 trm(i,5,14) = sa21*v5(i,1,2)+sa22*v5(i,2,2)+sa23*v5(i,3,2)
800 trm(i,6,14) = sa31*v5(i,1,2)+sa32*v5(i,2,2)+sa33*v5(i,3,2)
801 trm(i,4,15) = sa11*v5(i,1,3)+sa12*v5(i,2,3)+sa13*v5(i,3,3)
802 trm(i,5,15) = sa21*v5(i,1,3)+sa22*v5(i,2,3)+sa23*v5(i,3,3)
803 trm(i,6,15) = sa31*v5(i,1,3)+sa32*v5(i,2,3)+sa33*v5(i,3,3)
804
805 trm(i,7,13) = sb11*v5(i,1,1)+sb12*v5(i,2,1)+sb13*v5(i,3,1)
806
807 trm(i,9,13) = sb31*v5(i,1,1)+sb32*v5(i,2,1)+sb33
808 trm(i,7,14) = sb11*v5(i,1,2)+sb12*v5
809 trm(i,8,14) = sb21*v5(i,1,2)+sb22*v5(i,2,2)+sb23*v5(i,3,2)
810 trm(i,9,14) = sb31*v5(i,1,2)+sb32*v5(i,2,2)+sb33*v5(i,3,2)
811 trm(i,7,15) = sb11*v5(i,1,3)+sb12*v5(i,2,3)+sb13*v5(i,3,3)
812 trm(i,8,15) = sb21*v5(i,1,3)+sb22*v5(i,2,3)+sb23*v5(i,3,3)
813 trm(i,9,15) = sb31*v5(i,1,3)+sb32*v5(i,2,3)+sb33*v5(i,3,3)
814
815 trm(i,10,13) = sc11*v5(i,1,1)+sc12*v5(i,2,1)+sc13*v5(i,3,1)
816 trm(i,11,13) = sc21*v5(i,1,1)+sc22*v5(i,2,1)+sc23*v5(i,3,1)
817 trm(i,12,13) = sc31*v5(i,1,1)+sc32*v5(i,2,1)+sc33*v5(i,3,1)
818 trm(i,10,14) = sc11*v5(i,1,2)+sc12*v5(i,2,2)+sc13*v5(i
819 trm(i,11,14) = sc21*v5(i,1,2)+sc22*v5(i,2,2)+sc23*v5(i,3,2)
820 trm(i,12,14) = sc31*v5(i,1,2)+sc32*v5(i,2,2)+sc33*v5(i,3,2)
821 trm(i,10,15) = sc11*v5(i,1,3)+sc12*v5(i,2,3)+sc13*v5(i,3,3)
822 trm(i,11,15) = sc21*v5(i,1,3)+sc22*v5(i,2,3)+sc23*v5(i,3,3)
823 trm(i,12,15) = sc31*v5(i,1,3)+sc32*v5(i,2,3)+sc33*v5(i,3,3)
824
825 trm(i,13,13) = sd11*v5(i,1,1)+sd12*v5(i,2,1)+sd13*v5(i,3,1)
826 trm(i,14,13) = sd21*v5(i,1,1)+sd22*v5(i,2,1)+sd23*v5(i,3,1)
827 trm(i,15,13) = sd31*v5(i,1,1)+sd32*v5(i,2,1)+sd33*v5(i,3,1)
828 trm(i,13,14) = sd11*v5(i,1,2)+sd12*v5(i,2,2)+sd13*v5(i,
829 trm(i,14,14) = sd21*v5(i,1,2)+sd22*v5(i,2,2)+sd23*v5(i,3,2)
830 trm(i,15,14) = sd31*v5(i,1,2)+sd32*v5(i,2,2)+sd33*v5
831 trm(i,13,15) = sd11*v5(i,1,3)+sd12*v5(i,2,3)+sd13*v5(i,3,3)
832 trm(i,14,15) = sd21*v5(i,1,3)+sd22*v5(i,2,3)+sd23*v5(i,3,3)
833 trm(i,15,15) = sd31*v5(i,1,3)+sd32*v5(i,2,3)+sd33*v5(i,3,3)
834
835 trm(i,16,13) = se11*v5(i,1,1)+se12*v5(i,2,1)+se13*v5(i,3,1)
836 trm(i,17,13) = se21*v5(i,1,1)+se22*v5(i,2,1)+se23*v5(i
837 trm(i,18,13) = se31*v5(i,1,1)+se32*v5(i,2,1)+se33*v5(i
838 trm(i,16,14) = se11*v5(i,1,2)+se12*v5(i,2,2)+se13*v5(i,3,2)
839 trm(i,17,14) = se21*v5(i,1,2)+se22*v5(i,2,2)+se23*v5(i,3,2)
840 trm(i,18,14) = se31*v5(i,1,2)+se32*v5(i,2,2)+se33*v5(i,3,2)
841 trm(i,16,15) = se11*v5(i,1,3)+se12*v5(i,2,3)+se13*v5(i,3,3)
842 trm(i,17,15) = se21*v5(i,1,3)+se22*v5(i,2,3)+se23*v5(i,3,3)
843 trm(i,18,15) = se31*v5(i,1,3)+se32*v5(i,2,3)+se33*v5(i,3,3)
844
845 trm(i,19,13) = sf11*v5(i,1,1)+sf12*v5(i,2,1)+sf13*v5(i,3,1)
846 trm(i,20,13) = sf21*v5(i,1,1)+sf22*v5(i,2,1)+sf23*v5(i,3,1)
847 trm(i,21,13) = sf31*v5(i,1,1)+sf32*v5(i,2,1)+sf33*v5(i,3,1)
848 trm(i,19,14) = sf11*v5(i,1,2)+sf12*v5(i,2,2)+sf13*v5(i,3,2)
849 trm(i,20,14) = sf21*v5(i,1,2)+sf22*v5(i,2,2)+sf23*v5(i,3,2)
850 trm(i,21,14) = sf31*v5(i,1,2)+sf32*v5(i,2,2)+sf33*v5(i
851 trm(i,19,15) = sf11*v5(i,1,3)+sf12*v5(i,2,3)+sf13*v5(i,3,3)
852 trm(i,20,15) = sf21*v5(i,1,3)+sf22*v5(i,2,3)+sf23*v5(i,3,3)
853 trm(i,21,15) = sf31*v5(i,1,3)+sf32*v5(i,2,3)+sf33*v5(i,3,3)
854
855 trm(i,22,13) = sg11*v5(i,1,1)+sg12*v5(i,2,1)+sg13*v5(i,3,1)
856 trm(i,23,13) = sg21*v5(i,1,1)+sg22*v5(i,2,1)+sg23*v5(i,3,1)
857 trm(i,24,13) = sg31*v5(i,1,1)+sg32*v5(i,2,1)+sg33*v5(i,3,1)
858 trm(i,22,14) = sg11*v5(i,1,2)+sg12*v5(i,2,2)+sg13*v5(i,3,2)
859 trm(i,23,14) = sg21*v5(i,1,2)+sg22*v5(i,2,2)+sg23*v5(i,3,2)
860 trm(i,24,14) = sg31*v5(i,1,2)+sg32*v5(i,2,2)+sg33*v5(i,3,2)
861 trm(i,22,15) = sg11*v5(i,1,3)+sg12*v5(i,2,3)+sg13*v5(i,3,3)
862 trm(i,23,15) = sg21*v5(i,1,3)+sg22*v5(i,2,3)+sg23*v5(i,3,3)
863 trm(i,24,15) = sg31*v5(i,1,3)+sg32*v5(i,2,3)+sg33*v5(i,3,3)
864
865 trm(i,1,16) = zero
866 trm(i,2,16) = zero
867 trm(i,3,16) = zero
868 trm(i,1,17) = zero
869 trm(i,2,17) = zero
870 trm(i,3,17) = zero
871 trm(i,1,18) = zero
872 trm(i,2,18) = zero
873 trm(i,3,18) = zero
874
875 trm(i,4,16) = sa11*v6(i,1,1)+sa12*v6(i,2,1)+sa13*v6(i,3,1)
876 trm(i,5,16) = sa21*v6(i,1,1)+sa22*v6(i,2,1)+sa23*v6(i,3,1)
877 trm(i,6,16) = sa31*v6(i,1,1)+sa32*v6(i,2,1)+sa33*v6(i,3,1)
878 trm(i,4,17) = sa11*v6(i,1,2)+sa12*v6(i,2,2)+sa13*v6(i,3,2)
879 trm(i,5,17) = sa21*v6(i,1,2)+sa22*v6(i,2,2)+sa23*v6(i,3,2)
880 trm(i,6,17) = sa31*v6(i,1,2)+sa32*v6(i,2,2)+sa33*v6(i,3,2)
881 trm(i,4,18) = sa11*v6(i,1,3)+sa12*v6(i,2,3)+sa13*v6(i,3,3)
882 trm(i,5,18) = sa21*v6(i,1,3)+sa22*v6(i,2,3)+sa23*v6(i,3,3)
883 trm(i,6,18) = sa31*v6(i,1,3)+sa32*v6(i,2,3)+sa33*v6(i,3,3)
884
885 trm(i,7,16) = sb11*v6(i,1,1)+sb12*v6(i,2,1)+sb13*v6(i,3,1)
886 trm(i,8,16) = sb21*v6(i,1,1)+sb22*v6(i,2,1)+sb23*v6(i,3,1)
887 trm(i,9,16) = sb31*v6(i,1,1)+sb32*v6(i,2,1)+sb33*v6(i,3,1)
888 trm(i,7,17) = sb11*v6(i,1,2)+sb12*v6(i,2,2)+sb13*v6(i,3,2)
889 trm(i,8,17) = sb21*v6(i,1,2)+sb22*v6(i,2,2)+sb23*v6(i,3,2)
890 trm(i,9,17) = sb31*v6(i,1,2)+sb32*v6(i,2,2)+sb33*v6(i,3,2)
891 trm(i,7,18) = sb11*v6(i,1,3)+sb12*v6(i,2,3)+sb13*v6(i,3,3)
892 trm(i,8,18) = sb21*v6(i,1,3)+sb22*v6(i,2,3)+sb23*v6(i,3,3)
893 trm(i,9,18) = sb31*v6(i,1,3)+sb32*v6(i,2,3)+sb33*v6(i,3,3)
894
895 trm(i,10,16) = sc11*v6(i,1,1)+sc12*v6(i,2,1)+sc13*v6(i,3,1)
896 trm(i,11,16) = sc21*v6(i,1,1)+sc22*v6(i,2,1)+sc23*v6(i,3,1)
897 trm(i,12,16) = sc31*v6(i,1,1)+sc32*v6(i,2,1)+sc33*v6(i,3,1)
898 trm(i,10,17) = sc11*v6(i,1,2)+sc12*v6(i,2,2)+sc13*v6(i,3,2)
899 trm(i,11,17) = sc21*v6(i,1,2)+sc22*v6(i,2,2)+sc23*v6(i,3,2)
900 trm(i,12,17) = sc31*v6(i,1,2)+sc32*v6(i,2,2)+sc33*v6(i,3,2)
901 trm(i,10,18) = sc11*v6(i,1,3)+sc12*v6(i,2,3)+sc13*v6(i,3,3)
902 trm(i,11,18) = sc21*v6(i,1,3)+sc22*v6(i,2,3)+sc23*v6(i,3,3)
903 trm(i,12,18) = sc31*v6(i,1,3)+sc32*v6(i,2,3)+sc33*v6(i,3,3)
904
905 trm(i,13,16) = sd11*v6(i,1,1)+sd12*v6(i,2,1)+sd13*v6(i,3,1)
906 trm(i,14,16) = sd21*v6(i,1,1)+sd22*v6(i,2,1)+sd23*v6(i,3,1)
907 trm(i,15,16) = sd31*v6(i,1,1)+sd32*v6(i,2,1)+sd33*v6(i,3,1)
908 trm(i,13,17) = sd11*v6(i,1,2)+sd12*v6(i,2,2)+sd13*v6(i,3,2)
909 trm(i,14,17) = sd21*v6(i,1,2)+sd22*v6(i,2,2)+sd23*v6(i,3,2)
910 trm(i,15,17) = sd31*v6(i,1,2)+sd32*v6(i,2,2)+sd33*v6(i,3,2)
911 trm(i,13,18) = sd11*v6(i,1,3)+sd12*v6(i,2,3)+sd13*v6(i,3,3)
912 trm(i,14,18) = sd21*v6(i,1,3)+sd22*v6(i,2,3)+sd23*v6(i,3,3)
913 trm(i,15,18) = sd31*v6(i,1,3)+sd32*v6(i,2,3)+sd33*v6(i,3,3)
914
915 trm(i,16,16) = se11*v6(i,1,1)+se12*v6(i,2,1)+se13*v6(i,3,1)
916 trm(i,17,16) = se21*v6(i,1,1)+se22*v6(i,2,1)+se23*v6(i,3,1)
917 trm(i,18,16) = se31*v6(i,1,1)+se32*v6(i,2,1)+se33*v6(i,3,1)
918 trm(i,16,17) = se11*v6(i,1,2)+se12*v6(i,2,2)+se13*v6(i,3,2)
919 trm(i,17,17) = se21*v6(i,1,2)+se22*v6(i,2,2)+se23*v6(i,3,2)
920 trm(i,18,17) = se31*v6(i,1,2)+se32*v6(i,2,2)+se33*v6(i,3,2)
921 trm(i,16,18) = se11*v6(i,1,3)+se12*v6(i,2,3)+se13*v6(i,3,3)
922 trm(i,17,18) = se21*v6(i,1,3)+se22*v6(i,2,3)+se23*v6(i,3,3)
923 trm(i,18,18) = se31*v6(i,1,3)+se32*v6(i,2,3)+se33*v6(i,3,3)
924
925 trm(i,19,16) = sf11*v6(i,1,1)+sf12*v6(i,2,1)+sf13*v6(i,3,1)
926 trm(i,20,16) = sf21*v6(i,1,1)+sf22*v6(i,2,1)+sf23*v6(i,3,1)
927 trm(i,21,16) = sf31*v6(i,1,1)+sf32*v6(i,2,1)+sf33*v6(i,3,1)
928 trm(i,19,17) = sf11*v6(i,1,2)+sf12*v6(i,2,2)+sf13*v6(i,3,2)
929 trm(i,20,17) = sf21*v6(i,1,2)+sf22*v6(i,2,2)+sf23*v6(i,3,2)
930 trm(i,21,17) = sf31*v6(i,1,2)+sf32*v6(i,2,2)+sf33*v6(i,3,2)
931 trm(i,19,18) = sf11*v6(i,1,3)+sf12*v6(i,2,3)+sf13*v6(i,3,3)
932 trm(i,20,18) = sf21*v6(i,1,3)+sf22*v6(i,2,3)+sf23*v6(i,3,3)
933 trm(i,21,18) = sf31*v6(i,1,3)+sf32*v6(i,2,3)+sf33*v6(i,3,3)
934
935 trm(i,22,16) = sg11*v6(i,1,1)+sg12*v6(i,2,1)+sg13*v6(i,3,1)
936 trm(i,23,16) = sg21*v6(i,1,1)+sg22*v6(i,2,1)+sg23*v6(i,3,1)
937 trm(i,24,16) = sg31*v6(i,1,1)+sg32*v6(i,2,1)+sg33*v6(i,3,1)
938 trm(i,22,17) = sg11*v6(i,1,2)+sg12*v6(i,2,2)+sg13*v6(i,3,2)
939 trm(i,23,17) = sg21*v6(i,1,2)+sg22*v6(i,2,2)+sg23*v6(i,3,2)
940 trm(i,24,17) = sg31*v6(i,1,2)+sg32*v6(i,2,2)+sg33*v6(i,3,2)
941 trm(i,22,18) = sg11*v6(i,1,3)+sg12*v6(i,2,3)+sg13*v6(i,3,3)
942 trm(i,23,18) = sg21*v6(i,1,3)+sg22*v6(i,2,3)+sg23*v6(i,3,3)
943 trm(i,24,18) = sg31*v6(i,1,3)+sg32*v6(i,2,3)+sg33*v6(i,3,3)
944
945 trm(i,1,19) = zero
946 trm(i,2,19) = zero
947 trm(i,3,19) = zero
948 trm(i,1,20) = zero
949 trm(i,2,20) = zero
950 trm(i,3,20) = zero
951 trm(i,1,21) = zero
952 trm(i,2,21) = zero
953 trm(i,3,21) = zero
954
955 trm(i,4,19) = sa11*v7(i,1,1)+sa12*v7(i,2,1)+sa13*v7(i,3,1)
956 trm(i,5,19) = sa21*v7(i,1,1)+sa22*v7(i,2,1)+sa23*v7(i,3,1)
957 trm(i,6,19) = sa31*v7(i,1,1)+sa32*v7(i,2,1)+sa33*v7(i,3,1)
958 trm(i,4,20) = sa11*v7(i,1,2)+sa12*v7(i,2,2)+sa13*v7(i,3,2)
959 trm(i,5,20) = sa21*v7(i,1,2)+sa22*v7(i,2,2)+sa23*v7(i,3,2)
960 trm(i,6,20) = sa31*v7(i,1,2)+sa32*v7(i,2,2)+sa33*v7(i,3,2)
961 trm(i,4,21) = sa11*v7(i,1,3)+sa12*v7(i,2,3)+sa13*v7(i,3,3)
962 trm(i,5,21) = sa21*v7(i,1,3)+sa22*v7(i,2,3)+sa23*v7(i,3,3)
963 trm(i,6,21) = sa31*v7(i,1,3)+sa32*v7(i,2,3)+sa33*v7(i,3,3)
964
965 trm(i,7,19) = sb11*v7(i,1,1)+sb12*v7(i,2,1)+sb13*v7(i,3,1)
966 trm(i,8,19) = sb21*v7(i,1,1)+sb22*v7(i,2,1)+sb23*v7(i,3,1)
967 trm(i,9,19) = sb31*v7(i,1,1)+sb32*v7(i,2,1)+sb33*v7(i,3,1)
968 trm(i,7,20) = sb11*v7(i,1,2)+sb12*v7(i,2,2)+sb13*v7(i,3,2)
969 trm(i,8,20) = sb21*v7(i,1,2)+sb22*v7(i,2,2)+sb23*v7(i,3,2)
970 trm(i,9,20) = sb31*v7(i,1,2)+sb32*v7(i,2,2)+sb33*v7(i,3,2)
971 trm(i,7,21) = sb11*v7(i,1,3)+sb12*v7(i,2,3)+sb13*v7(i,3,3)
972 trm(i,8,21) = sb21*v7(i,1,3)+sb22*v7(i,2,3)+sb23*v7(i,3,3)
973 trm(i,9,21) = sb31*v7(i,1,3)+sb32*v7(i,2,3)+sb33*v7(i,3,3)
974
975 trm(i,10,19) = sc11*v7(i,1,1)+sc12*v7(i,2,1)+sc13*v7(i,3,1)
976 trm(i,11,19) = sc21*v7(i,1,1)+sc22*v7(i,2,1)+sc23*v7(i,3,1)
977 trm(i,12,19) = sc31*v7(i,1,1)+sc32*v7(i,2,1)+sc33*v7(i,3,1)
978 trm(i,10,20) = sc11*v7(i,1,2)+sc12*v7(i,2,2)+sc13*v7(i,3,2)
979 trm(i,11,20) = sc21*v7(i,1,2)+sc22*v7(i,2,2)+sc23*v7(i,3,2)
980 trm(i,12,20) = sc31*v7(i,1,2)+sc32*v7(i,2,2)+sc33*v7(i,3,2)
981 trm(i,10,21) = sc11*v7(i,1,3)+sc12*v7(i,2,3)+sc13*v7(i,3,3)
982 trm(i,11,21) = sc21*v7(i,1,3)+sc22*v7(i,2,3)+sc23*v7(i,3,3)
983 trm(i,12,21) = sc31*v7(i,1,3)+sc32*v7(i,2,3)+sc33*v7(i,3,3)
984
985 trm(i,13,19) = sd11*v7(i,1,1)+sd12*v7(i,2,1)+sd13*v7(i,3,1)
986 trm(i,14,19) = sd21*v7(i,1,1)+sd22*v7(i,2,1)+sd23*v7(i,3,1)
987 trm(i,15,19) = sd31*v7(i,1,1)+sd32*v7(i,2,1)+sd33*v7(i,3,1)
988 trm(i,13,20) = sd11*v7(i,1,2)+sd12*v7(i,2,2)+sd13*v7(i,3,2)
989 trm(i,14,20) = sd21*v7(i,1,2)+sd22*v7(i,2,2)+sd23*v7(i,3,2)
990 trm(i,15,20) = sd31*v7(i,1,2)+sd32*v7(i,2,2)+sd33*v7(i,3,2)
991 trm(i,13,21) = sd11*v7(i,1,3)+sd12*v7(i,2,3)+sd13*v7(i,3,3)
992 trm(i,14,21) = sd21*v7(i,1,3)+sd22*v7(i,2,3)+sd23*v7(i,3,3)
993 trm(i,15,21) = sd31*v7(i,1,3)+sd32*v7(i,2,3)+sd33*v7(i,3,3)
994
995 trm(i,16,19) = se11*v7(i,1,1)+se12*v7(i,2,1)+se13*v7(i,3,1)
996 trm(i,17,19) = se21*v7(i,1,1)+se22*v7(i,2,1)+se23*v7(i,3,1)
997 trm(i,18,19) = se31*v7(i,1,1)+se32*v7(i,2,1)+se33*v7(i,3,1)
998 trm(i,16,20) = se11*v7(i,1,2)+se12*v7(i,2,2)+se13*v7(i,3,2)
999 trm(i,17,20) = se21*v7(i,1,2)+se22*v7(i,2,2)+se23*v7(i,3,2)
1000 trm(i,18,20) = se31*v7(i,1,2)+se32*v7(i,2,2)+se33*v7(i,3,2)
1001 trm(i,16,21) = se11*v7(i,1,3)+se12*v7(i,2,3)+se13*v7(i,3,3)
1002 trm(i,17,21) = se21*v7(i,1,3)+se22*v7(i,2,3)+se23*v7(i,3,3)
1003 trm(i,18,21) = se31*v7(i,1,3)+se32*v7(i,2,3)+se33*v7(i,3,3)
1004
1005 trm(i,19,19) = sf11*v7(i,1,1)+sf12*v7(i,2,1)+sf13*v7(i,3,1)
1006 trm(i,20,19) = sf21*v7(i,1,1)+sf22*v7(i,2,1)+sf23*v7(i,3,1)
1007 trm(i,21,19) = sf31*v7(i,1,1)+sf32*v7(i,2,1)+sf33*v7(i,3,1)
1008 trm(i,19,20) = sf11*v7(i,1,2)+sf12*v7(i,2,2)+sf13*v7(i,3,2)
1009 trm(i,20,20) = sf21*v7(i,1,2)+sf22*v7(i,2,2)+sf23*v7(i,3,2)
1010 trm(i,21,20) = sf31*v7(i,1,2)+sf32*v7(i,2,2)+sf33*v7(i,3,2)
1011 trm(i,19,21) = sf11*v7(i,1,3)+sf12*v7(i,2,3)+sf13*v7(i,3,3)
1012 trm(i,20,21) = sf21*v7(i,1,3)+sf22*v7(i,2,3)+sf23*v7(i,3,3)
1013 trm(i,21,21) = sf31*v7(i,1,3)+sf32*v7(i,2,3)+sf33*v7(i,3,3)
1014
1015 trm(i,22,19) = sg11*v7(i,1,1)+sg12*v7(i,2,1)+sg13*v7(i,3,1)
1016 trm(i,23,19) = sg21*v7(i,1,1)+sg22*v7(i,2,1)+sg23*v7(i,3,1)
1017 trm(i,24,19) = sg31*v7(i,1,1)+sg32*v7(i,2,1)+sg33*v7(i,3,1)
1018 trm(i,22,20) = sg11*v7(i,1,2)+sg12*v7(i,2,2)+sg13*v7(i,3,2)
1019 trm(i,23,20) = sg21*v7(i,1,2)+sg22*v7(i,2,2)+sg23*v7(i,3,2)
1020 trm(i,24,20) = sg31*v7(i,1,2)+sg32*v7(i,2,2)+sg33*v7(i,3,2)
1021 trm(i,22,21) = sg11*v7(i,1,3)+sg12*v7(i,2,3)+sg13*v7(i,3,3)
1022 trm(i,23,21) = sg21*v7(i,1,3)+sg22*v7(i,2,3)+sg23*v7(i,3,3)
1023 trm(i,24,21) = sg31*v7(i,1,3)+sg32*v7(i,2,3)+sg33*v7(i,3,3)
1024
1025 trm(i,1,22) = zero
1026 trm(i,2,22) = zero
1027 trm(i,3,22) = zero
1028 trm(i,1,23) = zero
1029 trm(i,2,23) = zero
1030 trm(i,3,23) = zero
1031 trm(i,1,24) = zero
1032 trm(i,2,24) = zero
1033 trm(i,3,24) = zero
1034
1035 trm(i,4,22) = sa11*v8(i,1,1)+sa12*v8(i,2,1)+sa13*v8(i,3,1)
1036 trm(i,5,22) = sa21*v8(i,1,1)+sa22*v8(i,2,1)+sa23*v8(i,3,1)
1037 trm(i,6,22) = sa31*v8(i,1,1)+sa32*v8(i,2,1)+sa33*v8(i,3,1)
1038 trm(i,4,23) = sa11*v8(i,1,2)+sa12*v8(i,2,2)+sa13*v8(i,3,2)
1039 trm(i,5,23) = sa21*v8(i,1,2)+sa22*v8(i,2,2)+sa23*v8(i,3,2)
1040 trm(i,6,23) = sa31*v8(i,1,2)+sa32*v8(i,2,2)+sa33*v8(i,3,2)
1041 trm(i,4,24) = sa11*v8(i,1,3)+sa12*v8(i,2,3)+sa13*v8(i,3,3)
1042 trm(i,5,24) = sa21*v8(i,1,3)+sa22*v8(i,2,3)+sa23*v8(i,3,3)
1043 trm(i,6,24) = sa31*v8(i,1,3)+sa32*v8(i,2,3)+sa33*v8(i,3,3)
1044
1045 trm(i,7,22) = sb11*v8(i,1,1)+sb12*v8(i,2,1)+sb13*v8(i,3,1)
1046 trm(i,8,22) = sb21*v8(i,1,1)+sb22*v8(i,2,1)+sb23*v8(i,3,1)
1047 trm(i,9,22) = sb31*v8(i,1,1)+sb32*v8(i,2,1)+sb33*v8(i,3,1)
1048 trm(i,7,23) = sb11*v8(i,1,2)+sb12*v8(i,2,2)+sb13*v8(i,3,2)
1049 trm(i,8,23) = sb21*v8(i,1,2)+sb22*v8(i,2,2)+sb23*v8(i,3,2)
1050 trm(i,9,23) = sb31*v8(i,1,2)+sb32*v8(i,2,2)+sb33*v8(i,3,2)
1051 trm(i,7,24) = sb11*v8(i,1,3)+sb12*v8(i,2,3)+sb13*v8(i,3,3)
1052 trm(i,8,24) = sb21*v8(i,1,3)+sb22*v8(i,2,3)+sb23*v8(i,3,3)
1053 trm(i,9,24) = sb31*v8(i,1,3)+sb32*v8(i,2,3)+sb33*v8(i,3,3)
1054
1055 trm(i,10,22) = sc11*v8(i,1,1)+sc12*v8(i,2,1)+sc13*v8(i,3,1)
1056 trm(i,11,22) = sc21*v8(i,1,1)+sc22*v8(i,2,1)+sc23*v8(i,3,1)
1057 trm(i,12,22) = sc31*v8(i,1,1)+sc32*v8(i,2,1)+sc33*v8(i,3,1)
1058 trm(i,10,23) = sc11*v8(i,1,2)+sc12*v8(i,2,2)+sc13*v8(i,3,2)
1059 trm(i,11,23) = sc21*v8(i,1,2)+sc22*v8(i,2,2)+sc23*v8(i,3,2)
1060 trm(i,12,23) = sc31*v8(i,1,2)+sc32*v8(i,2,2)+sc33*v8(i,3,2)
1061 trm(i,10,24) = sc11*v8(i,1,3)+sc12*v8(i,2,3)+sc13*v8(i,3,3)
1062 trm(i,11,24) = sc21*v8(i,1,3)+sc22*v8(i,2,3)+sc23*v8(i,3,3)
1063 trm(i,12,24) = sc31*v8(i,1,3)+sc32*v8(i,2,3)+sc33*v8(i,3,3)
1064
1065 trm(i,13,22) = sd11*v8(i,1,1)+sd12*v8(i,2,1)+sd13*v8(i,3,1)
1066 trm(i,14,22) = sd21*v8(i,1,1)+sd22*v8(i,2,1)+sd23*v8(i,3,1)
1067 trm(i,15,22) = sd31*v8(i,1,1)+sd32*v8(i,2,1)+sd33*v8(i,3,1)
1068 trm(i,13,23) = sd11*v8(i,1,2)+sd12*v8(i,2,2)+sd13*v8(i,3,2)
1069 trm(i,14,23) = sd21*v8(i,1,2)+sd22*v8(i,2,2)+sd23*v8(i,3,2)
1070 trm(i,15,23) = sd31*v8(i,1,2)+sd32*v8(i,2,2)+sd33*v8(i,3,2)
1071 trm(i,13,24) = sd11*v8(i,1,3)+sd12*v8(i,2,3)+sd13*v8(i,3,3)
1072 trm(i,14,24) = sd21*v8(i,1,3)+sd22*v8(i,2,3)+sd23*v8(i,3,3)
1073 trm(i,15,24) = sd31*v8(i,1,3)+sd32*v8(i,2,3)+sd33*v8(i,3,3)
1074
1075 trm(i,16,22) = se11*v8(i,1,1)+se12*v8(i,2,1)+se13*v8(i,3,1)
1076 trm(i,17,22) = se21*v8(i,1,1)+se22*v8(i,2,1)+se23*v8(i,3,1)
1077 trm(i,18,22) = se31*v8(i,1,1)+se32*v8(i,2,1)+se33*v8(i,3,1)
1078 trm(i,16,23) = se11*v8(i,1,2)+se12*v8(i,2,2)+se13*v8(i,3,2)
1079 trm(i,17,23) = se21*v8(i,1,2)+se22*v8(i,2,2)+se23*v8(i,3,2)
1080 trm(i,18,23) = se31*v8(i,1,2)+se32*v8(i,2,2)
1081 trm(i,16,24) = se11*v8(i,1,3)+se12*v8(i,2,3)+se13*v8(i,3,3)
1082 trm(i,17,24) = se21*v8(i,1,3)+se22*v8(i,2,3)+se23*v8(i,3,3)
1083 trm(i,18,24) = se31*v8(i,1,3)+se32*v8(i,2,3)+se33*v8(i,3,3)
1084
1085 trm(i,19,22) = sf11*v8(i,1,1)+sf12*v8(i,2,1)+sf13*v8(i,3,1)
1086 trm(i,20,22) = sf21*v8(i,1,1)+sf22*v8(i,2,1)+sf23*v8(i,3,1)
1087 trm(i,21,22) = sf31*v8(i,1,1)+sf32*v8(i,2,1)+sf33*v8(i,3,1)
1088 trm(i,19,23) = sf11*v8(i,1,2)+sf12*v8(i,2,2)+sf13*v8(i,3,2)
1089 trm(i,20,23) = sf21*v8(i,1,2)+sf22*v8(i,2,2)+sf23*v8(i,3,2)
1090 trm(i,21,23) = sf31*v8(i,1,2)+sf32*v8(i,2,2)+sf33*v8(i,3,2)
1091 trm(i,19,24) = sf11*v8(i,1,3)+sf12*v8(i,2,3)+sf13*v8(i,3,3)
1092 trm(i,20,24) = sf21*v8(i,1,3)+sf22*v8(i,2,3)+sf23*v8(i,3,3)
1093 trm(i,21,24) = sf31*v8(i,1,3)+sf32*v8(i,2,3)+sf33*v8(i,3,3)
1094
1095 trm(i,22,22) = sg11*v8(i,1,1)+sg12*v8(i,2,1)+sg13*v8(i,3,1)
1096 trm(i,23,22) = sg21*v8(i,1,1)+sg22*v8(i,2,1)+sg23*v8(i,3,1)
1097 trm(i,24,22) = sg31*v8(i,1,1)+sg32*v8(i,2,1)+sg33*v8(i,3,1)
1098 trm(i,22,23) = sg11*v8(i,1,2)+sg12*v8(i,2,2)+sg13*v8(i,3,2)
1099 trm(i,23,23) = sg21*v8(i,1,2)+sg22*v8(i,2,2)+sg23*v8(i,3,2)
1100 trm(i,24,23) = sg31*v8(i,1,2)+sg32*v8(i,2,2)+sg33*v8(i,3,2)
1101 trm(i,22,24) = sg11*v8(i,1,3)+sg12*v8(i,2,3)+sg13*v8(i,3,3)
1102 trm(i,23,24) = sg21*v8(i,1,3)+sg22*v8(i,2,3)+sg23*v8(i,3,3)
1103 trm(i,24,24) = sg31*v8(i,1,3)+sg32*v8(i,2,3)+sg33*v8(i,3,3)
1104 enddo
1105
1106 DO j=1,8
1107 jj = 3*(j-1)
1108 DO i=1,nel
1109 trm(i,jj+1,jj+1) = r(1,1,i)+trm(i,jj+1,jj+1)
1110 trm(i,jj+2,jj+1) = r(1,2,i)+trm(i,jj+2,jj+1)
1111 trm(i,jj+3,jj+1) = r(1,3,i)+trm(i,jj+3,jj+1)
1112 trm(i,jj+1,jj+2) = r(2,1,i)+trm(i,jj+1,jj+2)
1113 trm(i,jj+2,jj+2) = r(2,2,i)+trm(i,jj+2,jj+2)
1114 trm(i,jj+3,jj+2) = r(2,3,i)+trm(i,jj+3,jj+2)
1115 trm(i,jj+1,jj+3) = r(3,1,i)+trm(i,jj+1,jj+3)
1116 trm(i,jj+2,jj+3) = r(3,2,i)+trm(i,jj+2,jj+3)
1117 trm(i,jj+3,jj+3) = r(3,3,i)+trm(i,jj+3,jj+3)
1118 ENDDO
1119 ENDDO
1120
1121 RETURN
1.15.0