OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
s8zlke3.F File Reference
#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine s8zlke3 (px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, px5, px6, px7, px8, py5, py6, py7, py8, pz5, pz6, pz7, pz8, pxy1, pxy2, pxy3, pxy4, pxy5, pxy6, pxy7, pxy8, pyx1, pyx2, pyx3, pyx4, pyx5, pyx6, pyx7, pyx8, pxz1, pxz2, pxz3, pxz4, pxz5, pxz6, pxz7, pxz8, pzx1, pzx2, pzx3, pzx4, pzx5, pzx6, pzx7, pzx8, pyz1, pyz2, pyz3, pyz4, pyz5, pyz6, pyz7, pyz8, pzy1, pzy2, pzy3, pzy4, pzy5, pzy6, pzy7, pzy8, bxy1, bxy2, bxy3, bxy4, bxy5, bxy6, bxy7, bxy8, byx1, byx2, byx3, byx4, byx5, byx6, byx7, byx8, bxz1, bxz2, bxz3, bxz4, bxz5, bxz6, bxz7, bxz8, bzx1, bzx2, bzx3, bzx4, bzx5, bzx6, bzx7, bzx8, byz1, byz2, byz3, byz4, byz5, byz6, byz7, byz8, bzy1, bzy2, bzy3, bzy4, bzy5, bzy6, bzy7, bzy8, k11, k12, k13, k14, k15, k16, k17, k18, k22, k23, k24, k25, k26, k27, k28, k33, k34, k35, k36, k37, k38, k44, k45, k46, k47, k48, k55, k56, k57, k58, k66, k67, k68, k77, k78, k88, dd, gg, icp, off, iksup, dg, g33, nel)

Function/Subroutine Documentation

◆ s8zlke3()

subroutine s8zlke3 ( px1,
px2,
px3,
px4,
py1,
py2,
py3,
py4,
pz1,
pz2,
pz3,
pz4,
px5,
px6,
px7,
px8,
py5,
py6,
py7,
py8,
pz5,
pz6,
pz7,
pz8,
pxy1,
pxy2,
pxy3,
pxy4,
pxy5,
pxy6,
pxy7,
pxy8,
pyx1,
pyx2,
pyx3,
pyx4,
pyx5,
pyx6,
pyx7,
pyx8,
pxz1,
pxz2,
pxz3,
pxz4,
pxz5,
pxz6,
pxz7,
pxz8,
pzx1,
pzx2,
pzx3,
pzx4,
pzx5,
pzx6,
pzx7,
pzx8,
pyz1,
pyz2,
pyz3,
pyz4,
pyz5,
pyz6,
pyz7,
pyz8,
pzy1,
pzy2,
pzy3,
pzy4,
pzy5,
pzy6,
pzy7,
pzy8,
bxy1,
bxy2,
bxy3,
bxy4,
bxy5,
bxy6,
bxy7,
bxy8,
byx1,
byx2,
byx3,
byx4,
byx5,
byx6,
byx7,
byx8,
bxz1,
bxz2,
bxz3,
bxz4,
bxz5,
bxz6,
bxz7,
bxz8,
bzx1,
bzx2,
bzx3,
bzx4,
bzx5,
bzx6,
bzx7,
bzx8,
byz1,
byz2,
byz3,
byz4,
byz5,
byz6,
byz7,
byz8,
bzy1,
bzy2,
bzy3,
bzy4,
bzy5,
bzy6,
bzy7,
bzy8,
k11,
k12,
k13,
k14,
k15,
k16,
k17,
k18,
k22,
k23,
k24,
k25,
k26,
k27,
k28,
k33,
k34,
k35,
k36,
k37,
k38,
k44,
k45,
k46,
k47,
k48,
k55,
k56,
k57,
k58,
k66,
k67,
k68,
k77,
k78,
k88,
dd,
gg,
integer icp,
off,
integer iksup,
dg,
g33,
integer, intent(in) nel )

Definition at line 35 of file s8zlke3.F.

77C-----------------------------------------------
78C I m p l i c i t T y p e s
79C-----------------------------------------------
80#include "implicit_f.inc"
81C-----------------------------------------------
82C C o m m o n B l o c k s
83C-----------------------------------------------
84C-----------------------------------------------
85C D u m m y A r g u m e n t s
86C-----------------------------------------------
87 INTEGER, INTENT(IN) :: NEL
88 INTEGER ICP,IKSUP
89C REAL
90 my_real
91 . px1(*), px2(*), px3(*), px4(*),
92 . py1(*), py2(*), py3(*), py4(*),
93 . pz1(*), pz2(*), pz3(*), pz4(*),
94 . px5(*), px6(*), px7(*), px8(*),
95 . py5(*), py6(*), py7(*), py8(*),
96 . pz5(*), pz6(*), pz7(*), pz8(*),
97 . pxy1(*),pxy2(*),pxy3(*),pxy4(*),
98 . pxy5(*),pxy6(*),pxy7(*),pxy8(*),
99 . pyx1(*),pyx2(*),pyx3(*),pyx4(*),
100 . pyx5(*),pyx6(*),pyx7(*),pyx8(*),
101 . pxz1(*),pxz2(*),pxz3(*),pxz4(*),
102 . pxz5(*),pxz6(*),pxz7(*),pxz8(*),
103 . pzx1(*),pzx2(*),pzx3(*),pzx4(*),
104 . pzx5(*),pzx6(*),pzx7(*),pzx8(*),
105 . pyz1(*),pyz2(*),pyz3(*),pyz4(*),
106 . pyz5(*),pyz6(*),pyz7(*),pyz8(*),
107 . pzy1(*),pzy2(*),pzy3(*),pzy4(*),
108 . pzy5(*),pzy6(*),pzy7(*),pzy8(*),
109 . bxy1(*),bxy2(*),bxy3(*),bxy4(*),
110 . bxy5(*),bxy6(*),bxy7(*),bxy8(*),
111 . byx1(*),byx2(*),byx3(*),byx4(*),
112 . byx5(*),byx6(*),byx7(*),byx8(*),
113 . bxz1(*),bxz2(*),bxz3(*),bxz4(*),
114 . bxz5(*),bxz6(*),bxz7(*),bxz8(*),
115 . bzx1(*),bzx2(*),bzx3(*),bzx4(*),
116 . bzx5(*),bzx6(*),bzx7(*),bzx8(*),
117 . byz1(*),byz2(*),byz3(*),byz4(*),
118 . byz5(*),byz6(*),byz7(*),byz8(*),
119 . bzy1(*),bzy2(*),bzy3(*),bzy4(*),
120 . bzy5(*),bzy6(*),bzy7(*),bzy8(*),
121 . k11(9,*) ,k12(9,*) ,k13(9,*) ,k14(9,*) ,k15(9,*) ,
122 . k16(9,*) ,k17(9,*) ,k18(9,*) ,k22(9,*) ,k23(9,*) ,
123 . k24(9,*) ,k25(9,*) ,k26(9,*) ,k27(9,*) ,k28(9,*) ,
124 . k33(9,*) ,k34(9,*) ,k35(9,*) ,k36(9,*) ,k37(9,*) ,
125 . k38(9,*) ,k44(9,*) ,k45(9,*) ,k46(9,*) ,k47(9,*) ,
126 . k48(9,*) ,k55(9,*) ,k56(9,*) ,k57(9,*) ,k58(9,*) ,
127 . k66(9,*) ,k67(9,*) ,k68(9,*) ,k77(9,*) ,k78(9,*) ,
128 . k88(9,*) ,dd(3,3,*),off(*),gg(*),dg(9,*) ,g33(9,*)
129C-----------------------------------------------
130C L o c a l V a r i a b l e s
131C-----------------------------------------------
132 INTEGER I, J,IS
133C REAL
134C-----------------------------------------------
135C -------symm(diag)------
136 is =1
137 IF (icp==1) THEN
138 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px1, py1, pz1, dd, k11,is)
139 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px2, py2, pz2, dd, k22,is)
140 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px3, py3, pz3, dd, k33,is)
141 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px4, py4, pz4, dd, k44,is)
142 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px5, py5, pz5, dd, k55,is)
143 CALL s8zked3(1,nel,px6, py6 ,pz6 ,px6, py6, pz6, dd, k66,is)
144 CALL s8zked3(1,nel,px7, py7 ,pz7 ,px7, py7, pz7, dd, k77,is)
145 CALL s8zked3(1,nel,px8, py8 ,pz8 ,px8, py8, pz8, dd, k88,is)
146 ELSE
147 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px1, py1, pz1,
148 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
149 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,dd,k11,is)
150 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px2, py2, pz2,
151 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
152 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,dd,k22,is)
153 CALL s8zkeb3(1,nel,px3, py3 ,pz3 ,px3, py3, pz3,
154 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
155 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,dd,k33,is)
156 CALL s8zkeb3(1,nel,px4, py4 ,pz4 ,px4, py4, pz4,
157 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
158 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,dd,k44,is)
159 CALL s8zkeb3(1,nel,px5, py5 ,pz5 ,px5, py5, pz5,
160 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
161 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,dd,k55,is)
162 CALL s8zkeb3(1,nel,px6, py6 ,pz6 ,px6, py6, pz6,
163 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
164 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,dd,k66,is)
165 CALL s8zkeb3(1,nel,px7, py7 ,pz7 ,px7, py7, pz7,
166 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
167 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k77,is)
168 CALL s8zkeb3(1,nel,px8, py8 ,pz8 ,px8, py8, pz8,
169 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
170 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k88,is)
171 ENDIF
172 IF (iksup==0) THEN
173 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
174 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,gg,k11,is)
175 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
176 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,gg,k22,is)
177 CALL s8zkeg3(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
178 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,gg,k33,is)
179 CALL s8zkeg3(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
180 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,gg,k44,is)
181 CALL s8zkeg3(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
182 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,gg,k55,is)
183 CALL s8zkeg3(1,nel,pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
184 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,gg,k66,is)
185 CALL s8zkeg3(1,nel,pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,
186 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k77,is)
187 CALL s8zkeg3(1,nel,pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,
188 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k88,is)
189 ENDIF
190C -------non symm------
191 is =0
192 IF (icp==1) THEN
193 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px2, py2, pz2, dd, k12,is)
194 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px3, py3, pz3, dd, k13,is)
195 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px4, py4, pz4, dd, k14,is)
196 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px5, py5, pz5, dd, k15,is)
197 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px6, py6, pz6, dd, k16,is)
198 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px7, py7, pz7, dd, k17,is)
199 CALL s8zked3(1,nel,px1, py1 ,pz1 ,px8, py8, pz8, dd, k18,is)
200 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px3, py3, pz3, dd, k23,is)
201 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px4, py4, pz4, dd, k24,is)
202 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px5, py5, pz5, dd, k25,is)
203 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px6, py6, pz6, dd, k26,is)
204 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px7, py7, pz7, dd, k27,is)
205 CALL s8zked3(1,nel,px2, py2 ,pz2 ,px8, py8, pz8, dd, k28,is)
206 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px4, py4, pz4, dd, k34,is)
207 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px5, py5, pz5, dd, k35,is)
208 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px6, py6, pz6, dd, k36,is)
209 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px7, py7, pz7, dd, k37,is)
210 CALL s8zked3(1,nel,px3, py3 ,pz3 ,px8, py8, pz8, dd, k38,is)
211 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px5, py5, pz5, dd, k45,is)
212 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px6, py6, pz6, dd, k46,is)
213 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px7, py7, pz7, dd, k47,is)
214 CALL s8zked3(1,nel,px4, py4 ,pz4 ,px8, py8, pz8, dd, k48,is)
215 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px6, py6, pz6, dd, k56,is)
216 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px7, py7, pz7, dd, k57,is)
217 CALL s8zked3(1,nel,px5, py5 ,pz5 ,px8, py8, pz8, dd, k58,is)
218 CALL s8zked3(1,nel,px6, py6 ,pz6 ,px7, py7, pz7, dd, k67,is)
219 CALL s8zked3(1,nel,px6, py6 ,pz6 ,px8, py8, pz8, dd, k68,is)
220 CALL s8zked3(1,nel,px7, py7 ,pz7 ,px8, py8, pz8, dd, k78,is)
221 ELSE
222 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px2, py2, pz2,
223 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
224 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,dd,k12,is)
225 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px3, py3, pz3,
226 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
227 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,dd,k13,is)
228 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px4, py4, pz4,
229 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
230 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,dd,k14,is)
231 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px5, py5, pz5,
232 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
233 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,dd,k15,is)
234 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px6, py6, pz6,
235 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
236 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,dd,k16,is)
237 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px7, py7, pz7,
238 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
239 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k17,is)
240 CALL s8zkeb3(1,nel,px1, py1 ,pz1 ,px8, py8, pz8,
241 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
242 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k18,is)
243 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px3, py3, pz3,
244 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
245 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,dd,k23,is)
246 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px4, py4, pz4,
247 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
248 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,dd,k24,is)
249 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px5, py5, pz5,
250 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
251 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,dd,k25,is)
252 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px6, py6, pz6,
253 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
254 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,dd,k26,is)
255 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px7, py7, pz7,
256 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
257 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k27,is)
258 CALL s8zkeb3(1,nel,px2, py2 ,pz2 ,px8, py8, pz8,
259 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
260 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k28,is)
261 CALL s8zkeb3(1,nel,px3, py3 ,pz3 ,px4, py4, pz4,
262 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
263 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,dd,k34,is)
264 CALL s8zkeb3(1,nel,px3, py3 ,pz3 ,px5, py5, pz5,
265 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
266 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,dd,k35,is)
267 CALL s8zkeb3(1,nel,px3, py3 ,pz3 ,px6, py6, pz6,
268 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
269 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,dd,k36,is)
270 CALL s8zkeb3(1,nel,px3, py3 ,pz3 ,px7, py7, pz7,
271 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
272 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k37,is)
273 CALL s8zkeb3(1,nel,px3, py3 ,pz3 ,px8, py8, pz8,
274 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
275 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k38,is)
276 CALL s8zkeb3(1,nel,px4, py4 ,pz4 ,px5, py5, pz5,
277 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
278 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,dd,k45,is)
279 CALL s8zkeb3(1,nel,px4, py4 ,pz4 ,px6, py6, pz6,
280 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
281 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,dd,k46,is)
282 CALL s8zkeb3(1,nel,px4, py4 ,pz4 ,px7, py7, pz7,
283 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
284 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k47,is)
285 CALL s8zkeb3(1,nel,px4, py4 ,pz4 ,px8, py8, pz8,
286 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
287 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k48,is)
288 CALL s8zkeb3(1,nel,px5, py5 ,pz5 ,px6, py6, pz6,
289 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
290 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,dd,k56,is)
291 CALL s8zkeb3(1,nel,px5, py5 ,pz5 ,px7, py7, pz7,
292 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
293 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k57,is)
294 CALL s8zkeb3(1,nel,px5, py5 ,pz5 ,px8, py8, pz8,
295 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
296 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k58,is)
297 CALL s8zkeb3(1,nel,px6, py6 ,pz6 ,px7, py7, pz7,
298 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
299 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,dd,k67,is)
300 CALL s8zkeb3(1,nel,px6, py6 ,pz6 ,px8, py8, pz8,
301 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
302 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k68,is)
303 CALL s8zkeb3(1,nel,px7, py7 ,pz7 ,px8, py8, pz8,
304 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
305 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,dd,k78,is)
306 ENDIF
307C
308C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
309 IF (iksup==0) THEN
310 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
311 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,gg,k12,is)
312 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
313 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,gg,k13,is)
314 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
315 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,gg,k14,is)
316 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
317 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,gg,k15,is)
318 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
319 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,gg,k16,is)
320 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
321 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k17,is)
322 CALL s8zkeg3(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
323 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k18,is)
324 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
325 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,gg,k23,is)
326 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
327 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,gg,k24,is)
328 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
329 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,gg,k25,is)
330 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
331 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,gg,k26,is)
332 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
333 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k27,is)
334 CALL s8zkeg3(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
335 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k28,is)
336 CALL s8zkeg3(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
337 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,gg,k34,is)
338 CALL s8zkeg3(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
339 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,gg,k35,is)
340 CALL s8zkeg3(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
341 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,gg,k36,is)
342 CALL s8zkeg3(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
343 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k37,is)
344 CALL s8zkeg3(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
345 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k38,is)
346 CALL s8zkeg3(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
347 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,gg,k45,is)
348 CALL s8zkeg3(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
349 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,gg,k46,is)
350 CALL s8zkeg3(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
351 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k47,is)
352 CALL s8zkeg3(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
353 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k48,is)
354 CALL s8zkeg3(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
355 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,gg,k56,is)
356 CALL s8zkeg3(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
357 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k57,is)
358 CALL s8zkeg3(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
359 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k58,is)
360 CALL s8zkeg3(1,nel,pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
361 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,gg,k67,is)
362 CALL s8zkeg3(1,nel,pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
363 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k68,is)
364 CALL s8zkeg3(1,nel,pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,
365 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,gg,k78,is)
366 ELSE
367 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
368 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,g33,k12,is)
369 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
370 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,g33,k13,is)
371 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
372 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,g33,k14,is)
373 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
374 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,g33,k15,is)
375 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
376 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,g33,k16,is)
377 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
378 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k17,is)
379 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
380 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k18,is)
381 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
382 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,g33,k23,is)
383 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
384 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,g33,k24,is)
385 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
386 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,g33,k25,is)
387 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
388 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,g33,k26,is)
389 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
390 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k27,is)
391 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
392 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k28,is)
393 CALL s8zkeg33(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
394 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,g33,k34,is)
395 CALL s8zkeg33(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
396 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,g33,k35,is)
397 CALL s8zkeg33(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
398 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,g33,k36,is)
399 CALL s8zkeg33(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
400 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k37,is)
401 CALL s8zkeg33(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
402 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k38,is)
403 CALL s8zkeg33(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
404 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,g33,k45,is)
405 CALL s8zkeg33(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
406 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,g33,k46,is)
407 CALL s8zkeg33(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
408 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k47,is)
409 CALL s8zkeg33(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
410 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k48,is)
411 CALL s8zkeg33(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
412 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,g33,k56,is)
413 CALL s8zkeg33(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
414 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k57,is)
415 CALL s8zkeg33(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
416 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k58,is)
417 CALL s8zkeg33(1,nel,pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
418 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k67,is)
419 CALL s8zkeg33(1,nel,pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
420 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k68,is)
421 CALL s8zkeg33(1,nel,pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,
422 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k78,is)
423C---+----1----+----2----+----3----+----4----+----5----+----6----+----7--
424 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px2, py2, pz2,
425 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
426 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
427 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
428 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,dg,k12,is,icp)
429 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px3, py3, pz3,
430 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
431 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
432 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
433 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,dg,k13,is,icp)
434 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px4, py4, pz4,
435 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
436 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
437 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
438 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,dg,k14,is,icp)
439 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px5, py5, pz5,
440 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
441 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
442 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
443 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,dg,k15,is,icp)
444 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px6, py6, pz6,
445 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
446 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
447 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
448 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,dg,k16,is,icp)
449 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px7, py7, pz7,
450 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
451 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
452 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
453 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k17,is,icp)
454 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px8, py8, pz8,
455 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
456 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
457 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
458 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k18,is,icp)
459 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px3, py3, pz3,
460 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
461 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
462 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
463 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,dg,k23,is,icp)
464 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px4, py4, pz4,
465 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
466 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
467 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
468 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,dg,k24,is,icp)
469 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px5, py5, pz5,
470 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
471 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
472 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
473 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,dg,k25,is,icp)
474 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px6, py6, pz6,
475 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
476 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
477 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
478 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,dg,k26,is,icp)
479 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px7, py7, pz7,
480 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
481 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
482 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
483 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k27,is,icp)
484 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px8, py8, pz8,
485 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
486 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
487 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
488 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k28,is,icp)
489 CALL s8zkebg3(1,nel,px3, py3 ,pz3 ,px4, py4, pz4,
490 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
491 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
492 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
493 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,dg,k34,is,icp)
494 CALL s8zkebg3(1,nel,px3, py3 ,pz3 ,px5, py5, pz5,
495 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
496 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
497 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
498 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,dg,k35,is,icp)
499 CALL s8zkebg3(1,nel,px3, py3 ,pz3 ,px6, py6, pz6,
500 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
501 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
502 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
503 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,dg,k36,is,icp)
504 CALL s8zkebg3(1,nel,px3, py3 ,pz3 ,px7, py7, pz7,
505 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
506 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
507 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
508 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k37,is,icp)
509 CALL s8zkebg3(1,nel,px3, py3 ,pz3 ,px8, py8, pz8,
510 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
511 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
512 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
513 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k38,is,icp)
514 CALL s8zkebg3(1,nel,px4, py4 ,pz4 ,px5, py5, pz5,
515 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
516 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
517 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
518 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,dg,k45,is,icp)
519 CALL s8zkebg3(1,nel,px4, py4 ,pz4 ,px6, py6, pz6,
520 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
521 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
522 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
523 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,dg,k46,is,icp)
524 CALL s8zkebg3(1,nel,px4, py4 ,pz4 ,px7, py7, pz7,
525 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
526 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
527 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
528 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k47,is,icp)
529 CALL s8zkebg3(1,nel,px4, py4 ,pz4 ,px8, py8, pz8,
530 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
531 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
532 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
533 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k48,is,icp)
534 CALL s8zkebg3(1,nel,px5, py5 ,pz5 ,px6, py6, pz6,
535 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
536 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
537 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
538 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,dg,k56,is,icp)
539 CALL s8zkebg3(1,nel,px5, py5 ,pz5 ,px7, py7, pz7,
540 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
541 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
542 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
543 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k57,is,icp)
544 CALL s8zkebg3(1,nel,px5, py5 ,pz5 ,px8, py8, pz8,
545 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
546 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
547 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
548 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k58,is,icp)
549 CALL s8zkebg3(1,nel,px6, py6 ,pz6 ,px7, py7, pz7,
550 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
551 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
552 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
553 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k67,is,icp)
554 CALL s8zkebg3(1,nel,px6, py6 ,pz6 ,px8, py8, pz8,
555 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
556 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
557 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
558 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k68,is,icp)
559 CALL s8zkebg3(1,nel,px7, py7 ,pz7 ,px8, py8, pz8,
560 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
561 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
562 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,
563 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k78,is,icp)
564C -------symm(diag)------
565 is =1
566c IS =0!..DG is non-symmetric
567 CALL s8zkebg3(1,nel,px1, py1 ,pz1 ,px1, py1, pz1,
568 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
569 . bxy1, byx1 ,bxz1,bzx1 ,byz1,bzy1,
570 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
571 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,dg,k11,is,icp)
572 CALL s8zkebg3(1,nel,px2, py2 ,pz2 ,px2, py2, pz2,
573 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
574 . bxy2, byx2 ,bxz2,bzx2 ,byz2,bzy2,
575 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
576 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,dg,k22,is,icp)
577 CALL s8zkebg3(1,nel,px3, py3 ,pz3 ,px3, py3, pz3,
578 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
579 . bxy3, byx3 ,bxz3,bzx3 ,byz3,bzy3,
580 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
581 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,dg,k33,is,icp)
582 CALL s8zkebg3(1,nel,px4, py4 ,pz4 ,px4, py4, pz4,
583 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
584 . bxy4, byx4 ,bxz4,bzx4 ,byz4,bzy4,
585 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
586 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,dg,k44,is,icp)
587 CALL s8zkebg3(1,nel,px5, py5 ,pz5 ,px5, py5, pz5,
588 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
589 . bxy5, byx5 ,bxz5,bzx5 ,byz5,bzy5,
590 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
591 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,dg,k55,is,icp)
592 CALL s8zkebg3(1,nel,px6, py6 ,pz6 ,px6, py6, pz6,
593 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
594 . bxy6, byx6 ,bxz6,bzx6 ,byz6,bzy6,
595 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
596 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,dg,k66,is,icp)
597 CALL s8zkebg3(1,nel,px7, py7 ,pz7 ,px7, py7, pz7,
598 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
599 . bxy7, byx7 ,bxz7,bzx7 ,byz7,bzy7,
600 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,
601 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,dg,k77,is,icp)
602 CALL s8zkebg3(1,nel,px8, py8 ,pz8 ,px8, py8, pz8,
603 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
604 . bxy8, byx8 ,bxz8,bzx8 ,byz8,bzy8,
605 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,
606 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,dg,k88,is,icp)
607 CALL s8zkeg33(1,nel,pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,
608 . pxy1, pyx1 ,pxz1,pzx1 ,pyz1,pzy1,g33,k11,is)
609 CALL s8zkeg33(1,nel,pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,
610 . pxy2, pyx2 ,pxz2,pzx2 ,pyz2,pzy2,g33,k22,is)
611 CALL s8zkeg33(1,nel,pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,
612 . pxy3, pyx3 ,pxz3,pzx3 ,pyz3,pzy3,g33,k33,is)
613 CALL s8zkeg33(1,nel,pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,
614 . pxy4, pyx4 ,pxz4,pzx4 ,pyz4,pzy4,g33,k44,is)
615 CALL s8zkeg33(1,nel,pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,
616 . pxy5, pyx5 ,pxz5,pzx5 ,pyz5,pzy5,g33,k55,is)
617 CALL s8zkeg33(1,nel,pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,
618 . pxy6, pyx6 ,pxz6,pzx6 ,pyz6,pzy6,g33,k66,is)
619 CALL s8zkeg33(1,nel,pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,
620 . pxy7, pyx7 ,pxz7,pzx7 ,pyz7,pzy7,g33,k77,is)
621 CALL s8zkeg33(1,nel,pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,
622 . pxy8, pyx8 ,pxz8,pzx8 ,pyz8,pzy8,g33,k88,is)
623 ENDIF
624C ---------------
625 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
s8zkeb3
subroutine s8zkeb3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, dd, kij, is)
Definition s8zkeb3.F:33
s8zkebg3
subroutine s8zkebg3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, bxyi, byxi, bxzi, bzxi, byzi, bzyi, bxyj, byxj, bxzj, bzxj, byzj, bzyj, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, dg, kij, is, icp)
Definition s8zkebg3.F:34
s8zked3
subroutine s8zked3(lft, llt, pxi, pyi, pzi, pxj, pyj, pzj, dd, kij, is)
Definition s8zked3.F:34
s8zkeg33
subroutine s8zkeg33(lft, llt, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, g33, kij, is)
Definition s8zkeg33.F:34
s8zkeg3
subroutine s8zkeg3(lft, llt, pxyi, pyxi, pxzi, pzxi, pyzi, pzyi, pxyj, pyxj, pxzj, pzxj, pyzj, pzyj, g, kij, is)
Definition s8zkeg3.F:34