OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
stat_brick_spmd.F
Go to the documentation of this file.
1Copyright> OpenRadioss
2Copyright> Copyright (C) 1986-2025 Altair Engineering Inc.
3Copyright>
4Copyright> This program is free software: you can redistribute it and/or modify
5Copyright> it under the terms of the GNU Affero General Public License as published by
6Copyright> the Free Software Foundation, either version 3 of the License, or
7Copyright> (at your option) any later version.
8Copyright>
9Copyright> This program is distributed in the hope that it will be useful,
10Copyright> but WITHOUT ANY WARRANTY; without even the implied warranty of
11Copyright> MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
12Copyright> GNU Affero General Public License for more details.
13Copyright>
14Copyright> You should have received a copy of the GNU Affero General Public License
15Copyright> along with this program. If not, see <https://www.gnu.org/licenses/>.
16Copyright>
17Copyright>
18Copyright> Commercial Alternative: Altair Radioss Software
19Copyright>
20Copyright> As an alternative to this open-source version, Altair also offers Altair Radioss
21Copyright> software under a commercial license. Contact Altair to discuss further if the
22Copyright> commercial version may interest you: https://www.altair.com/radioss/.
23!||====================================================================
24!|| stat_brick_spmd ../engine/source/output/sta/stat_brick_spmd.F
25!||--- called by ------------------------------------------------------
26!|| genstat ../engine/source/output/sta/genstat.F
27!||--- calls -----------------------------------------------------
28!|| my_orders ../common_source/tools/sort/my_orders.c
29!|| spmd_iget_partn_sta ../engine/source/mpi/output/spmd_stat.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!|| element_mod ../common_source/modules/elements/element_mod.F90
33!|| my_alloc_mod ../common_source/tools/memory/my_alloc.F90
34!||====================================================================
35 SUBROUTINE stat_brick_spmd(ITAB,IPART,IGEO,IXS,IPARTS,
36 . IPART_STATE,NODTAG,STAT_INDXS,
37 . IPARG, LENGS,IXS10,IXS16,IXS20,
38 . ELBUF_TAB,IDEL )
39C-----------------------------------------------
40C M o d u l e s
41C-----------------------------------------------
42 USE elbufdef_mod
43 USE my_alloc_mod
44 use element_mod , only : nixs
45C-----------------------------------------------
46C I m p l i c i t T y p e s
47C-----------------------------------------------
48#include "implicit_f.inc"
49C-----------------------------------------------
50C C o m m o n B l o c k s
51C-----------------------------------------------
52#include "com01_c.inc"
53#include "com04_c.inc"
54#include "param_c.inc"
55#include "scr16_c.inc"
56#include "scr17_c.inc"
57#include "spmd_c.inc"
58#include "task_c.inc"
59#include "units_c.inc"
60C-----------------------------------------------
61C D u m m y A r g u m e n t s
62C-----------------------------------------------
63 INTEGER ITAB(*), IPART(LIPART1,*),
64 . IGEO(NPROPGI,*), IXS(NIXS,*),
65 . IPARTS(*), IPART_STATE(*),
66 . nodtag(*), stat_indxs(*),
67 . iparg(nparg,*),lengs,ixs10(6,*),ixs16(8,*),
68 . ixs20(12,*),idel
69 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
70C-----------------------------------------------
71C L o c a l V a r i a b l e s
72C-----------------------------------------------
73 INTEGER I, N, JJ, IPRT, IPRT0, K
74 INTEGER NG, NEL, NFT, LFT, LLT, ITY, ISOLNOD, IOFF
75 INTEGER,DIMENSION(:,:),ALLOCATABLE :: IADG
76 INTEGER,DIMENSION(:),ALLOCATABLE :: IADD
77 INTEGER,DIMENSION(:),ALLOCATABLE :: NP
78 INTEGER,DIMENSION(:),ALLOCATABLE :: NPGLOB
79 INTEGER,DIMENSION(:,:),ALLOCATABLE :: CLEF
80 INTEGER WORK(70000)
81 TYPE(g_bufel_) ,POINTER :: GBUF
82C-----------------------------------------------
83C 8 NODES BRICK
84C-----------------------------------------------
85 CALL my_alloc(iadg,nspmd,npart)
86 CALL my_alloc(iadd,npart+1)
87 CALL my_alloc(np,24*numelsg)
88 CALL my_alloc(npglob,24*numelsg)
89 CALL my_alloc(clef,2,numelsg)
90C-----------------------------------------------
91 iadd = 0
92 npglob(1:24*numelsg) = 0
93C
94C In SPMD you must send the info even qd we have 0 local elem
95 jj = 0
96 DO ng=1,ngroup
97 ity =iparg(5,ng)
98 isolnod = iparg(28,ng)
99 nel =iparg(2,ng)
100 nft =iparg(3,ng)
101 gbuf => elbuf_tab(ng)%GBUF
102 lft=1
103 llt=nel
104 IF(ity == 1) THEN
105 DO i=lft,llt
106 n = i + nft
107
108 iprt=iparts(n)
109 IF(ipart_state(iprt)==0)cycle
110
111 np(jj+1) = ixs(nixs,n)
112 IF (isolnod == 10)THEN
113 np(jj+2) = itab(ixs(2,n))
114 np(jj+3) = itab(ixs(4,n))
115 np(jj+4) = itab(ixs(7,n))
116 np(jj+5) = itab(ixs(6,n))
117 IF(ixs10(1,n - numels8) /= 0)THEN
118 np(jj+6) = itab(ixs10(1,n - numels8))
119 ELSE
120 np(jj+6) = 0
121 ENDIF
122 IF(ixs10(2,n - numels8) /= 0)THEN
123 np(jj+7) = itab(ixs10(2,n - numels8))
124 ELSE
125 np(jj+7) = 0
126 ENDIF
127 IF(ixs10(3,n - numels8) /= 0)THEN
128 np(jj+8) = itab(ixs10(3,n - numels8))
129 ELSE
130 np(jj+8) = 0
131 ENDIF
132 IF(ixs10(4,n - numels8) /= 0)THEN
133 np(jj+9) = itab(ixs10(4,n - numels8))
134 ELSE
135 np(jj+9) = 0
136 ENDIF
137 IF(ixs10(5,n - numels8) /= 0)THEN
138 np(jj+10) = itab(ixs10(5,n - numels8))
139 ELSE
140 np(jj+10) = 0
141 ENDIF
142 IF(ixs10(6,n - numels8) /= 0)THEN
143 np(jj+11) = itab(ixs10(6,n - numels8))
144 ELSE
145 np(jj+11) = 0
146 ENDIF
147 np(jj+12) = 0
148 np(jj+13) = 0
149 np(jj+14) = 0
150 np(jj+15) = 0
151 np(jj+16) = 0
152 np(jj+17) = 0
153 np(jj+18) = 0
154 np(jj+19) = 0
155 np(jj+20) = 0
156 np(jj+21) = 0
157 ELSEIF (isolnod == 16)THEN
158 np(jj+2) = itab(ixs(2,n))
159 np(jj+3) = itab(ixs(3,n))
160 np(jj+4) = itab(ixs(4,n))
161 np(jj+5) = itab(ixs(5,n))
162 np(jj+6) = itab(ixs(6,n))
163 np(jj+7) = itab(ixs(7,n))
164 np(jj+8) = itab(ixs(8,n))
165 np(jj+9) = itab(ixs(9,n))
166 IF(ixs16(1,n - (numels8+numels10+numels20)) /= 0)THEN
167 np(jj+10) = itab(ixs16(1,n - (numels8+numels10+numels20)))
168 ELSE
169 np(jj+10) = 0
170 ENDIF
171 IF(ixs16(2,n - (numels8+numels10+numels20)) /= 0)THEN
172 np(jj+11) = itab(ixs16(2,n - (numels8+numels10+numels20)))
173 ELSE
174 np(jj+11) = 0
175 ENDIF
176 IF(ixs16(3,n - (numels8+numels10+numels20)) /= 0)THEN
177 np(jj+12) = itab(ixs16(3,n - (numels8+numels10+numels20)))
178 ELSE
179 np(jj+12) = 0
180 ENDIF
181 IF(ixs16(4,n - (numels8+numels10+numels20)) /= 0)THEN
182 np(jj+13) = itab(ixs16(4,n - (numels8+numels10+numels20)))
183 ELSE
184 np(jj+13) = 0
185 ENDIF
186 IF(ixs16(5,n - (numels8+numels10+numels20)) /= 0)THEN
187 np(jj+14) = itab(ixs16(5,n - (numels8+numels10+numels20)))
188 ELSE
189 np(jj+14) = 0
190 ENDIF
191 IF(ixs16(6,n - (numels8+numels10+numels20)) /= 0)THEN
192 np(jj+15) = itab(ixs16(6,n - (numels8+numels10+numels20)))
193 ELSE
194 np(jj+15) = 0
195 ENDIF
196 IF(ixs16(7,n - (numels8+numels10+numels20)) /= 0)THEN
197 np(jj+16) = itab(ixs16(7,n - (numels8+numels10+numels20)))
198 ELSE
199 np(jj+16) = 0
200 ENDIF
201 IF(ixs16(8,n - (numels8+numels10+numels20)) /= 0)THEN
202 np(jj+17) = itab(ixs16(8,n - (numels8+numels10+numels20)))
203 ELSE
204 np(jj+17) = 0
205 ENDIF
206 np(jj+18) = 0
207 np(jj+19) = 0
208 np(jj+20) = 0
209 np(jj+21) = 0
210 ELSEIF (isolnod == 20)THEN
211 np(jj+2) = itab(ixs(2,n))
212 np(jj+3) = itab(ixs(3,n))
213 np(jj+4) = itab(ixs(4,n))
214 np(jj+5) = itab(ixs(5,n))
215 np(jj+6) = itab(ixs(6,n))
216 np(jj+7) = itab(ixs(7,n))
217 np(jj+8) = itab(ixs(8,n))
218 np(jj+9) = itab(ixs(9,n))
219 IF(ixs20(1,n - (numels8+numels10)) /= 0)THEN
220 np(jj+10) = itab(ixs20(1,n - (numels8+numels10)))
221 ELSE
222 np(jj+10) = 0
223 ENDIF
224 IF(ixs20(2,n - (numels8+numels10)) /= 0)THEN
225 np(jj+11) = itab(ixs20(2,n - (numels8+numels10)))
226 ELSE
227 np(jj+11) = 0
228 ENDIF
229 IF(ixs20(3,n - (numels8+numels10)) /= 0)THEN
230 np(jj+12) = itab(ixs20(3,n - (numels8+numels10)))
231 ELSE
232 np(jj+12) = 0
233 ENDIF
234 IF(ixs20(4,n - (numels8+numels10)) /= 0)THEN
235 np(jj+13) = itab(ixs20(4,n - (numels8+numels10)))
236 ELSE
237 np(jj+13) = 0
238 ENDIF
239 IF(ixs20(5,n - (numels8+numels10)) /= 0)THEN
240 np(jj+14) = itab(ixs20(5,n - (numels8+numels10)))
241 ELSE
242 np(jj+14) = 0
243 ENDIF
244 IF(ixs20(6,n - (numels8+numels10)) /= 0)THEN
245 np(jj+15) = itab(ixs20(6,n - (numels8+numels10)))
246 ELSE
247 np(jj+15) = 0
248 ENDIF
249 IF(ixs20(7,n - (numels8+numels10)) /= 0)THEN
250 np(jj+16) = itab(ixs20(7,n - (numels8+numels10)))
251 ELSE
252 np(jj+16) = 0
253 ENDIF
254 IF(ixs20(8,n - (numels8+numels10)) /= 0)THEN
255 np(jj+17) = itab(ixs20(8,n - (numels8+numels10)))
256 ELSE
257 np(jj+17) = 0
258 ENDIF
259 IF(ixs20(9,n - (numels8+numels10)) /= 0)THEN
260 np(jj+18) = itab(ixs20(9,n - (numels8+numels10)))
261 ELSE
262 np(jj+18) = 0
263 ENDIF
264 IF(ixs20(10,n - (numels8+numels10)) /= 0)THEN
265 np(jj+19) = itab(ixs20(10,n - (numels8+numels10)))
266 ELSE
267 np(jj+19) = 0
268 ENDIF
269 IF(ixs20(11,n - (numels8+numels10)) /= 0)THEN
270 np(jj+20) = itab(ixs20(11,n - (numels8+numels10)))
271 ELSE
272 np(jj+20) = 0
273 ENDIF
274 IF(ixs20(12,n - (numels8+numels10)) /= 0)THEN
275 np(jj+21) = itab(ixs20(12,n - (numels8+numels10)))
276 ELSE
277 np(jj+21) = 0
278 ENDIF
279 ELSE
280 np(jj+2) = itab(ixs(2,n))
281 np(jj+3) = itab(ixs(3,n))
282 np(jj+4) = itab(ixs(4,n))
283 np(jj+5) = itab(ixs(5,n))
284 np(jj+6) = itab(ixs(6,n))
285 np(jj+7) = itab(ixs(7,n))
286 np(jj+8) = itab(ixs(8,n))
287 np(jj+9) = itab(ixs(9,n))
288 np(jj+10) = 0
289 np(jj+11) = 0
290 np(jj+12) = 0
291 np(jj+13) = 0
292 np(jj+14) = 0
293 np(jj+15) = 0
294 np(jj+16) = 0
295 np(jj+17) = 0
296 np(jj+18) = 0
297 np(jj+19) = 0
298 np(jj+20) = 0
299 np(jj+21) = 0
300 ENDIF
301 np(jj+22) = iprt
302 np(jj+23) = isolnod
303 np(jj+24) = iabs(nint(gbuf%OFF(i)))
304 jj = jj + 24
305
306 stat_numels =stat_numels+1
307
308 IF (isolnod == 10)THEN
309 nodtag(ixs(2,n))=1
310 nodtag(ixs(4,n))=1
311 nodtag(ixs(7,n))=1
312 nodtag(ixs(6,n))=1
313 IF (ixs10(1,n - numels8) /= 0) nodtag(ixs10(1,n - numels8))=1
314 IF (ixs10(2,n - numels8) /= 0) nodtag(ixs10(2,n - numels8))=1
315 IF (ixs10(3,n - numels8) /= 0) nodtag(ixs10(3,n - numels8))=1
316 IF (ixs10(4,n - numels8) /= 0) nodtag(ixs10(4,n - numels8))=1
317 IF (ixs10(5,n - numels8) /= 0) nodtag(ixs10(5,n - numels8))=1
318 IF (ixs10(6,n - numels8) /= 0) nodtag(ixs10(6,n - numels8))=1
319 ELSEIF (isolnod == 16)THEN
320 nodtag(ixs(2,n))=1
321 nodtag(ixs(3,n))=1
322 nodtag(ixs(4,n))=1
323 nodtag(ixs(5,n))=1
324 nodtag(ixs(6,n))=1
325 nodtag(ixs(7,n))=1
326 nodtag(ixs(8,n))=1
327 nodtag(ixs(9,n))=1
328 IF (ixs16(1,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(1,n - (numels8+numels10+numels20)))=1
329 IF (ixs16(2,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(2,n - (numels8+numels10+numels20)))=1
330 IF (ixs16(3,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(3,n - (numels8+numels10+numels20)))=1
331 IF (ixs16(4,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(4,n - (numels8+numels10+numels20)))=1
332 IF (ixs16(5,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(5,n - (numels8+numels10+numels20)))=1
333 IF (ixs16(6,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(6,n - (numels8+numels10+numels20)))=1
334 IF (ixs16(7,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(7,n - (numels8+numels10+numels20)))=1
335 IF (ixs16(8,n - (numels8+numels10+numels20)) /= 0) nodtag(ixs16(8,n - (numels8+numels10+numels20)))=1
336 ELSEIF (isolnod == 20)THEN
337 nodtag(ixs(2,n))=1
338 nodtag(ixs(3,n))=1
339 nodtag(ixs(4,n))=1
340 nodtag(ixs(5,n))=1
341 nodtag(ixs(6,n))=1
342 nodtag(ixs(7,n))=1
343 nodtag(ixs(8,n))=1
344 nodtag(ixs(9,n))=1
345 IF (ixs20(1,n - (numels8+numels10)) /= 0) nodtag(ixs20(1,n - (numels8+numels10)))=1
346 IF (ixs20(2,n - (numels8+numels10)) /= 0) nodtag(ixs20(2,n - (numels8+numels10)))=1
347 IF (ixs20(3,n - (numels8+numels10)) /= 0) nodtag(ixs20(3,n - (numels8+numels10)))=1
348 IF (ixs20(4,n - (numels8+numels10)) /= 0) nodtag(ixs20(4,n - (numels8+numels10)))=1
349 IF (ixs20(5,n - (numels8+numels10)) /= 0) nodtag(ixs20(5,n - (numels8+numels10)))=1
350 IF (ixs20(6,n - (numels8+numels10)) /= 0) nodtag(ixs20(6,n - (numels8+numels10)))=1
351 IF (ixs20(7,n - (numels8+numels10)) /= 0) nodtag(ixs20(7,n - (numels8+numels10)))=1
352 IF (ixs20(8,n - (numels8+numels10)) /= 0) nodtag(ixs20(8,n - (numels8+numels10)))=1
353 IF (ixs20(9,n - (numels8+numels10)) /= 0) nodtag(ixs20(9,n - (numels8+numels10)))=1
354 IF (ixs20(10,n - (numels8+numels10)) /= 0) nodtag(ixs20(10,n - (numels8+numels10)))=1
355 IF (ixs20(11,n - (numels8+numels10)) /= 0) nodtag(ixs20(11,n - (numels8+numels10)))=1
356 IF (ixs20(12,n - (numels8+numels10)) /= 0) nodtag(ixs20(12,n - (numels8+numels10)))=1
357 ELSE
358 nodtag(ixs(2,n))=1
359 nodtag(ixs(3,n))=1
360 nodtag(ixs(4,n))=1
361 nodtag(ixs(5,n))=1
362 nodtag(ixs(6,n))=1
363 nodtag(ixs(7,n))=1
364 nodtag(ixs(8,n))=1
365 nodtag(ixs(9,n))=1
366 ENDIF
367 END DO
368 END IF
369 END DO
370C-----
371 stat_numels_g=0
372 CALL spmd_iget_partn_sta(24,stat_numels,stat_numels_g,lengs,np,
373 . iadg,npglob,stat_indxs)
374C-----
375 IF (ispmd==0) THEN
376 DO n=1,stat_numels_g
377 stat_indxs(n)=n
378 clef(1,n)=npglob(19*(n-1)+22)
379 clef(2,n)=npglob(19*(n-1)+1)
380 END DO
381 CALL my_orders(0,work,clef,stat_indxs,stat_numels_g,2)
382
383 iprt0=0
384 DO n=1,stat_numels_g
385 k=stat_indxs(n)
386 jj=24*(k-1)
387 iprt=npglob(jj+22)
388 ioff=npglob(jj+24)
389 IF (npglob(jj+23) == 4) THEN
390 IF(idel==0.OR.(idel==1.AND.ioff >= 1)) THEN
391 IF(iprt /= iprt0)THEN
392 WRITE(iugeo,'(A,I10)')'/TETRA4/',ipart(4,iprt)
393 WRITE(iugeo,'(A)')
394 . '# TETRA4ID NOD1 NOD2 NOD3 NOD4'
395 iprt0=iprt
396 END IF
397 WRITE(iugeo,'(5I10)') npglob(jj+1),npglob(jj+2),npglob(jj+4),
398 . npglob(jj+8),npglob(jj+6)
399 ENDIF
400 ELSEIF (npglob(jj+23) == 6) THEN
401 IF(idel==0.OR.(idel==1.AND.ioff >= 1)) THEN
402 IF(iprt /= iprt0)THEN
403 WRITE(iugeo,'(A,I10)')'/PENTA6/',ipart(4,iprt)
404 WRITE(iugeo,'(A)')
405 . '# PENTA6ID NOD1 NOD2 NOD3 NOD4 NOD5 NOD6'
406 iprt0=iprt
407 END IF
408 WRITE(iugeo,'(7I10)') npglob(jj+1),npglob(jj+2),npglob(jj+3),
409 . npglob(jj+4),npglob(jj+6),npglob(jj+7),
410 . npglob(jj+8)
411 ENDIF
412 ELSEIF (npglob(jj+23) == 8) THEN
413 IF(idel==0.OR.(idel==1.AND.ioff >= 1)) THEN
414 IF(iprt /= iprt0)THEN
415 WRITE(iugeo,'(A,I10)')'/BRICK/',ipart(4,iprt)
416 WRITE(iugeo,'(A)')
417 . '# BRICKID NOD1 NOD2 NOD3 NOD4 NOD5 NOD6 NOD7 NOD8'
418 iprt0=iprt
419 END IF
420 WRITE(iugeo,'(9I10)') npglob(jj+1),npglob(jj+2),npglob(jj+3),
421 . npglob(jj+4),npglob(jj+5),npglob(jj+6),
422 . npglob(jj+7),npglob(jj+8),npglob(jj+9)
423 ENDIF
424 ELSEIF (npglob(jj+23) == 10) THEN
425 IF(idel==0.OR.(idel==1.AND.ioff >= 1)) THEN
426 IF(iprt /= iprt0)THEN
427 WRITE(iugeo,'(A,I10)')'/TETRA10/',ipart(4,iprt)
428 WRITE(iugeo,'(A)')
429 . '#TETRA10ID'
430 WRITE(iugeo,'(A)')
431 . '# NOD1 NOD2 NOD3 NOD4 NOD5 NOD6 NOD7 NOD8 NOD9 NOD10'
432 iprt0=iprt
433 END IF
434 WRITE(iugeo,'(I10)') npglob(jj+1)
435 WRITE(iugeo,'(10I10)') npglob(jj+2),npglob(jj+3),npglob(jj+4),
436 . npglob(jj+5),npglob(jj+6),npglob(jj+7),
437 . npglob(jj+8),npglob(jj+9),npglob(jj+10),
438 . npglob(jj+11)
439 ENDIF
440 ELSEIF (npglob(jj+23) == 16) THEN
441 IF(idel==0.OR.(idel==1.AND.ioff >= 1)) THEN
442 IF(iprt /= iprt0)THEN
443 WRITE(iugeo,'(A,I10)')'/SHEL16/',ipart(4,iprt)
444 WRITE(iugeo,'(A)')
445 . '#TSHEL16ID NOD1 NOD2 NOD3 NOD4 NOD5 NOD6 NOD7 NOD8'
446 WRITE(iugeo,'(A)')
447 . '# NOD9 NOD10 NOD11 NOD12'
448 WRITE(iugeo,'(A)')
449 . '# NOD13 NOD14 NOD15 NOD16'
450 iprt0=iprt
451 ENDIF
452 WRITE(iugeo,'(9I10)') npglob(jj+1),npglob(jj+2),npglob(jj+3),
453 . npglob(jj+4),npglob(jj+5),npglob(jj+6),
454 . npglob(jj+7), npglob(jj+8),npglob(jj+9)
455 WRITE(iugeo,'(4I10)') npglob(jj+10),npglob(jj+11),
456 . npglob(jj+12), npglob(jj+13)
457 WRITE(iugeo,'(4I10)') npglob(jj+14),npglob(jj+15),
458 . npglob(jj+16),npglob(jj+17)
459 ENDIF
460 ELSEIF (npglob(jj+23) == 20) THEN
461 IF(idel==0.OR.(idel==1.AND.ioff >= 1)) THEN
462 IF(iprt /= iprt0)THEN
463 WRITE(iugeo,'(A,I10)')'/BRIC20/',ipart(4,iprt)
464 WRITE(iugeo,'(A)')
465 . '#BRICK20ID NOD1 NOD2 NOD3 NOD4 NOD5 NOD6 NOD7 NOD8'
466 WRITE(iugeo,'(A)')
467 . '# NOD9 NOD10 NOD11 NOD12 NOD13 NOD14 NOD15 NOD16'
468 WRITE(iugeo,'(A)')
469 . '# NOD17 NOD18 NOD19 NOD20'
470 iprt0=iprt
471 END IF
472 WRITE(iugeo,'(9I10)') npglob(jj+1),npglob(jj+2),npglob(jj+3),
473 . npglob(jj+4),npglob(jj+5),npglob(jj+6),
474 . npglob(jj+7), npglob(jj+8),npglob(jj+9)
475 WRITE(iugeo,'(8I10)')npglob(jj+10),npglob(jj+11),
476 . npglob(jj+12),npglob(jj+13),
477 . npglob(jj+14),npglob(jj+15),
478 . npglob(jj+16),npglob(jj+17)
479 WRITE(iugeo,'(4I10)')npglob(jj+18),npglob(jj+19),
480 . npglob(jj+20),npglob(jj+21)
481 ENDIF
482 ENDIF
483 END DO
484 ENDIF
485C-----------------------------------------------
486 DEALLOCATE(iadg)
487 DEALLOCATE(iadd)
488 DEALLOCATE(np)
489 DEALLOCATE(npglob)
490 DEALLOCATE(clef)
491C-----------------------------------------------
492 RETURN
493 END
my_orders
void my_orders(int *mode, int *iwork, int *data, int *index, int *n, int *irecl)
Definition my_orders.c:82
spmd_iget_partn_sta
subroutine spmd_iget_partn_sta(size, stat_numel, stat_lenelg, leng, np, iadg, npglob, stat_indx)
Definition spmd_stat.F:129
stat_brick_spmd
subroutine stat_brick_spmd(itab, ipart, igeo, ixs, iparts, ipart_state, nodtag, stat_indxs, iparg, lengs, ixs10, ixs16, ixs20, elbuf_tab, idel)
Definition stat_brick_spmd.F:39