OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
outp_sp_t.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!|| outp_sp_t ../engine/source/output/sty/outp_sp_t.F
25!||--- called by ------------------------------------------------------
26!|| genoutp ../engine/source/output/sty/genoutp.F
27!||--- calls -----------------------------------------------------
28!|| initbuf ../engine/share/resol/initbuf.F
29!|| spmd_rgather9_1comm ../engine/source/mpi/interfaces/spmd_outp.F
30!||--- uses -----------------------------------------------------
31!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
32!|| initbuf_mod ../engine/share/resol/initbuf.F
33!||====================================================================
34 SUBROUTINE outp_sp_t(KEY,TEXT,ELBUF_TAB,IPARG,DD_IAD,SIZLOC,SIZP0,SIZ_WR)
35C-----------------------------------------------
36C M o d u l e s
37C-----------------------------------------------
38 USE initbuf_mod
39 USE elbufdef_mod
40C-----------------------------------------------
41C I m p l i c i t T y p e s
42C-----------------------------------------------
43#include "implicit_f.inc"
44C-----------------------------------------------
45C C o m m o n B l o c k s
46C-----------------------------------------------
47#include "vect01_c.inc"
48#include "com01_c.inc"
49#include "param_c.inc"
50#include "units_c.inc"
51#include "task_c.inc"
52#include "scr16_c.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56 CHARACTER*10 KEY
57 CHARACTER*40 TEXT
58 INTEGER IPARG(NPARG,*), DD_IAD(NSPMD+1,*),
59 . SIZLOC,SIZP0,SIZ_WR
60 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
61C-----------------------------------------------
62C L o c a l V a r i a b l e s
63C-----------------------------------------------
64 INTEGER I, J, NG, NEL, II(6), JJ, JJ_OLD, NGF, NGL, NN,
65 . LEN,RESP0,WRTLEN,RES,COMPTEUR,L,K
66 INTEGER, DIMENSION(NSPGROUP) :: JJ_LOC
67 INTEGER, DIMENSION(NSPGROUP+1,NSPMD) :: ADRESS
68 my_real
69 . wa(sizloc),wap0(siz_wr),wap0_loc(sizp0)
70 my_real
71 . func(6)
72 TYPE(g_bufel_) ,POINTER :: GBUF
73C=======================================================================
74 IF (ispmd == 0) THEN
75 WRITE(iugeo,'(2A)')'/SPHCEL /TENSOR /',key
76 WRITE(iugeo,'(A)')text
77 IF (outyy_fmt == 2) THEN
78 WRITE(iugeo,'(2A)') '#FORMAT: (1P6E12.5) ',
79 . '(TX(I),TY(I),TZ(I),TXY(I),TYZ(I),TZX(I),I=1,NUMSPH)'
80 ELSE
81 WRITE(iugeo,'(2A)') '#FORMAT: (1P6E20.13) ',
82 . '(TX(I),TY(I),TZ(I),TXY(I),TYZ(I),TZX(I),I=1,NUMSPH)'
83 END IF
84 END IF
85C
86 jj_old = 1
87 ngf = 1
88 ngl = 0
89 resp0=0
90 jj = 0
91 compteur = 0
92 DO nn=1,nspgroup
93 ngl = ngl + dd_iad(ispmd+1,nn)
94 DO ng = ngf, ngl
95 ity =iparg(5,ng)
96 IF(ity == 51) THEN
97 CALL initbuf(iparg ,ng ,
98 2 mtn ,nel ,nft ,iad ,ity ,
99 3 npt ,jale ,ismstr ,jeul ,jtur ,
100 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
101 5 nvaux ,jpor ,jcvt ,jclose ,jplasol ,
102 6 irep ,iint ,igtyp ,israt ,isrot ,
103 7 icsen ,isorth ,isorthg ,ifailure,jsms )
104 gbuf => elbuf_tab(ng)%GBUF
105 lft=1
106 llt=nel
107!
108 DO i=1,6
109 ii(i) = nel*(i-1)
110 ENDDO
111!
112 DO i=lft,llt
113 wa(jj+(i-1)*6+1) = gbuf%SIG(ii(1)+i)
114 wa(jj+(i-1)*6+2) = gbuf%SIG(ii(2)+i)
115 wa(jj+(i-1)*6+3) = gbuf%SIG(ii(3)+i)
116 wa(jj+(i-1)*6+4) = gbuf%SIG(ii(4)+i)
117 wa(jj+(i-1)*6+5) = gbuf%SIG(ii(5)+i)
118 wa(jj+(i-1)*6+6) = gbuf%SIG(ii(6)+i)
119 ENDDO
120 jj = jj + 6*llt
121 ENDIF
122 ENDDO
123 ngf = ngl + 1
124 jj_loc(nn) = jj - compteur ! size of each group
125 compteur = jj
126 ENDDO
127! ++++++++++
128 IF( nspmd>1 ) THEN
129 CALL spmd_rgather9_1comm(wa,jj,jj_loc,wap0_loc,sizp0,adress)
130 ELSE
131 wap0_loc(1:jj) = wa(1:jj)
132 adress(1,1) = 1
133 DO nn = 2,nspgroup+1
134 adress(nn,1) = jj_loc(nn-1) + adress(nn-1,1)
135 ENDDO
136 ENDIF
137! ++++++++++
138 IF(ispmd==0) THEN
139 resp0 = 0
140 DO nn=1,nspgroup
141 compteur = 0
142 DO k = 1,nspmd
143 IF((adress(nn+1,k)-1-adress(nn,k))>=0) THEN
144 DO l = adress(nn,k),adress(nn+1,k)-1
145 compteur = compteur + 1
146 wap0(compteur+resp0) = wap0_loc(l)
147 ENDDO ! l=... , ...
148 ENDIF !if(size_loc>0)
149 ENDDO ! k=1,nspmd
150
151 jj_old = compteur+resp0
152 IF(jj_old > 0) THEN
153 res=mod(jj_old,6)
154 wrtlen=jj_old-res
155 IF (wrtlen > 0) THEN
156 IF (outyy_fmt == 2) THEN
157 WRITE(iugeo,'(1P6E12.5)')(wap0(j),j=1,wrtlen)
158 ELSE
159 WRITE(iugeo,'(1p6e20.13)')(WAP0(J),J=1,WRTLEN)
160 END IF
161 ENDIF
162 DO I=1,RES
163 WAP0(I)=WAP0(WRTLEN+I)
164 ENDDO
165 RESP0=RES
166 ENDIF ! jj_old>0
167 ENDDO ! nn = 1, nspgroup
168 IF (RESP0 > 0) THEN
169 IF (OUTYY_FMT == 2) THEN
170 WRITE(IUGEO,'(1p6e12.5)')(WAP0(J),J=1,RESP0)
171 ELSE
172 WRITE(IUGEO,'(1p6e20.13)')(WAP0(J),J=1,RESP0)
173 ENDIF
174 ENDIF
175 END IF ! ispmd = 0
176 RETURN
177 END
178!||====================================================================
179!|| outp_sp_tt ../engine/source/output/sty/outp_sp_t.F
180!||--- called by ------------------------------------------------------
181!|| genoutp ../engine/source/output/sty/genoutp.F
182!||--- calls -----------------------------------------------------
183!|| initbuf ../engine/share/resol/initbuf.F
184!|| spmd_rgather9_1comm ../engine/source/mpi/interfaces/spmd_outp.F
185!||--- uses -----------------------------------------------------
186!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
187!|| initbuf_mod ../engine/share/resol/initbuf.F
188!||====================================================================
189 SUBROUTINE OUTP_SP_TT(KEY,TEXT,ELBUF_TAB,IPARG,DD_IAD,
190 2 IPM ,KXSP ,SPBUF,SIZLOC,SIZP0,SIZ_WR)
191C-----------------------------------------------
192C M o d u l e s
193C-----------------------------------------------
194 USE INITBUF_MOD
195 USE ELBUFDEF_MOD
196C-----------------------------------------------
197C I m p l i c i t T y p e s
198C-----------------------------------------------
199#include "implicit_f.inc"
200C-----------------------------------------------
201C C o m m o n B l o c k s
202C-----------------------------------------------
203#include "vect01_c.inc"
204#include "scr16_c.inc"
205#include "com01_c.inc"
206#include "param_c.inc"
207#include "units_c.inc"
208#include "task_c.inc"
209#include "sphcom.inc"
210C-----------------------------------------------
211C D u m m y A r g u m e n t s
212C-----------------------------------------------
213 CHARACTER*10 KEY
214 CHARACTER*40 TEXT
215 INTEGER KXSP(NISP,*),IPARG(NPARG,*), DD_IAD(NSPMD+1,*),
216 . IPM(NPROPMI,*),SIZ_WR
217 INTEGER SIZLOC,SIZP0
218 my_real
219 . SPBUF(NSPBUF,*)
220 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP), TARGET :: ELBUF_TAB
221C-----------------------------------------------
222C L o c a l V a r i a b l e s
223C-----------------------------------------------
224 INTEGER I, J, K, II(6), JJ, JJ_OLD,
225 . NG, NEL, IADD, N,NGF, NGL, NN, LEN,MLW,COMPTEUR,L
226 INTEGER, DIMENSION(NSPGROUP) :: JJ_LOC
227 INTEGER, DIMENSION(NSPGROUP+1,NSPMD) :: ADRESS
228 my_real
229 . WA(SIZLOC),WAP0(SIZ_WR),WAP0_LOC(SIZP0)
230 my_real
231 . FUNC(6)
232 TYPE(G_BUFEL_) ,POINTER :: GBUF
233C=======================================================================
234 IF (ISPMD == 0) THEN
235 WRITE(IUGEO,'(2a)')'/sphcel /tensor /',KEY
236 WRITE(IUGEO,'(a)')TEXT
237 IF (OUTYY_FMT == 2) THEN
238 WRITE(IUGEO,'(a)')
239 . '#FORMAT:(3E12.5),
240 . eint(i),rho(i),h(i),i=1,numsph '
241 WRITE(IUGEO,'(2a)')'#FORMAT: (1P6E12.5/E12.5) ',
242 . '(TX(I),TY(I),TZ(I),TXY(I),TYZ(I),TZX(I),',
243 . '#EPSP(I),I=1,NUMSPH)'
244 ELSE
245 WRITE(iugeo,'(A)')
246 . '#FORMAT:(3E20.13),
247 . EINT(I),RHO(I),H(I),I=1,NUMSPH '
248 WRITE(iugeo,'(2A)')'#FORMAT: (1P6E20.13/E20.13) ',
249 . '(TX(I),TY(I),TZ(I),TXY(I),TYZ(I),TZX(I),',
250 . '#EPSP(I),I=1,NUMSPH)'
251 END IF
252 END IF
253C
254 jj_old = 1
255 ngf = 1
256 ngl = 0
257 jj = 0
258 compteur = 0
259 DO nn=1,nspgroup
260 ngl = ngl + dd_iad(ispmd+1,nn)
261 DO ng = ngf, ngl
262 ity =iparg(5,ng)
263 IF (ity == 51) THEN
264 CALL initbuf(iparg ,ng ,
265 2 mtn ,nel ,nft ,iad ,ity ,
266 3 npt ,jale ,ismstr ,jeul ,jtur ,
267 4 jthe ,jlag ,jmult ,jhbe ,jivf ,
268 5 nvaux ,jpor ,jcvt ,jclose ,jplasol ,
269 6 irep ,iint ,igtyp ,israt ,isrot ,
270 7 icsen ,isorth ,isorthg ,ifailure,jsms )
271 gbuf => elbuf_tab(ng)%GBUF
272 lft=1
273 llt=nel
274!
275 DO i=1,6
276 ii(i) = nel*(i-1)
277 ENDDO
278!
279 DO i=lft,llt
280 wa(jj+1) = gbuf%EINT(i)
281 wa(jj+2) = gbuf%RHO(i)
282 wa(jj+3) = spbuf(1,nft+i)
283 wa(jj+4) = gbuf%SIG(ii(1)+i)
284 wa(jj+5) = gbuf%SIG(ii(2)+i)
285 wa(jj+6) = gbuf%SIG(ii(3)+i)
286 wa(jj+7) = gbuf%SIG(ii(4)+i)
287 wa(jj+8) = gbuf%SIG(ii(5)+i)
288 wa(jj+9) = gbuf%SIG(ii(6)+i)
289 IF (gbuf%G_PLA > 0) THEN
290 wa(jj+10) = gbuf%PLA(i)
291 wa(jj+11) = one
292 ELSE
293 wa(jj+10) = zero
294 wa(jj+11) = - one
295 ENDIF
296 jj=jj+11
297 ENDDO
298 ENDIF
299 ENDDO
300 ngf = ngl + 1
301 jj_loc(nn) = jj - compteur ! size of each group
302 compteur = jj
303 ENDDO ! nn = 1, nspdgroup
304
305! ++++++++++
306 IF( nspmd>1 ) THEN
307 CALL spmd_rgather9_1comm(wa,jj,jj_loc,wap0_loc,sizp0,adress)
308 ELSE
309 wap0_loc(1:jj) = wa(1:jj)
310 adress(1,1) = 1
311 DO nn = 2,nspgroup+1
312 adress(nn,1) = jj_loc(nn-1) + adress(nn-1,1)
313 ENDDO
314 ENDIF
315! ++++++++++
316 IF(ispmd==0) THEN
317 DO nn=1,nspgroup
318 compteur = 0
319 DO k = 1,nspmd
320 IF((adress(nn+1,k)-1-adress(nn,k))>=0) THEN
321 DO l = adress(nn,k),adress(nn+1,k)-1
322 compteur = compteur + 1
323 wap0(compteur) = wap0_loc(l)
324 ENDDO ! l=... , ...
325 ENDIF !if(size_loc>0)
326 ENDDO ! k=1,nspmd
327
328 jj_old = compteur
329 IF(jj_old>0) THEN
330 IF (outyy_fmt == 2) THEN
331 j = 1
332 DO WHILE (j<jj_old)
333 WRITE(iugeo,'(1P3E12.5)')(wap0(j-1+k),k=1,3) ! j-1+k=1,...,3
334 j=j+3 ! j = 4
335 WRITE(iugeo,'(1P6E12.5)')(wap0(j-1+k),k=1,6)! j-1+k = 4,...,9
336 IF (wap0(j+7) > zero) WRITE(iugeo,'(1P1E12.5)') (wap0(j+6)) ! test on G_PLA (wa(11)) then write wa(10)
337 j=j+8 ! j = 12
338 ENDDO
339 ELSE
340 j = 1
341 DO WHILE (j<=jj_old)
342 WRITE(iugeo,'(1P3E20.13)')(wap0(j-1+k),k=1,3) ! j-1+k=1,...,3
343 j=j+3 ! j = 4
344 WRITE(iugeo,'(1P6E20.13)')(wap0(j-1+k),k=1,6) ! j-1+k = 4,...,9
345 IF (wap0(j+7) > zero) WRITE(iugeo,'(1P1E20.13)') (wap0(j+6)) ! test on G_PLA (wa(11)) then write wa(10)
346 j=j+8 ! j = 12
347 ENDDO
348 END IF
349 END IF
350 ENDDO ! nn = 1, nspgroup
351 END IF ! ispmd=0
352C-----------
353 RETURN
354 END
355!||====================================================================
356!|| count_arsz_spt ../engine/source/output/sty/outp_sp_t.F
357!||--- called by ------------------------------------------------------
358!|| genoutp ../engine/source/output/sty/genoutp.F
359!|| outp_arsz_spt ../engine/source/mpi/interfaces/spmd_outp.F
360!||====================================================================
361 SUBROUTINE count_arsz_spt(IPARG,DD_IAD,WASZ,SIZ_WRITE_LOC)
362C-----------------------------------------------
363C I m p l i c i t T y p e s
364C-----------------------------------------------
365#include "implicit_f.inc"
366C-----------------------------------------------
367C C o m m o n B l o c k s
368C-----------------------------------------------
369#include "param_c.inc"
370#include "com01_c.inc"
371#include "task_c.inc"
372#include "scr16_c.inc"
373C-----------------------------------------------
374C D u m m y A r g u m e n t s
375C-----------------------------------------------
376 INTEGER IPARG(NPARG,*),DD_IAD(NSPMD+1,*),WASZ,SIZ_WRITE_LOC(NSPGROUP+1)
377C-----------------------------------------------
378C L o c a l V a r i a b l e s
379C-----------------------------------------------
380 INTEGER NGF,NGL,NN,ITY,NEL,NG,JJ,
381 . P0ARSZ2,WASZ2
382C-----------------------------------------------
383 wasz = 0
384 IF (outp_spt( 1) == 1 ) THEN
385
386 ngf = 1
387 ngl = 0
388 DO nn=1,nspgroup
389 jj = 0
390 ngl = ngl + dd_iad(ispmd+1,nn)
391 DO ng = ngf, ngl
392 ity =iparg(5,ng)
393 IF(ity == 51) THEN
394 nel = iparg(2,ng)
395 jj = jj + 6*nel
396 ENDIF
397 ENDDO
398 ngf = ngl + 1
399 wasz = wasz + jj
400 siz_write_loc(nn) = jj
401 ENDDO
402 siz_write_loc(nspgroup+1) = wasz
403
404 ENDIF
405
406 RETURN
407 END
408!||====================================================================
409!|| count_arsz_sptt ../engine/source/output/sty/outp_sp_t.F
410!||--- called by ------------------------------------------------------
411!|| genoutp ../engine/source/output/sty/genoutp.F
412!|| outp_arsz_sptt ../engine/source/mpi/interfaces/spmd_outp.F
413!||====================================================================
414 SUBROUTINE count_arsz_sptt(IPARG,DD_IAD,WASZ,SIZ_WRITE_LOC)
415C-----------------------------------------------
416C I m p l i c i t T y p e s
417C-----------------------------------------------
418#include "implicit_f.inc"
419C-----------------------------------------------
420C C o m m o n B l o c k s
421C-----------------------------------------------
422#include "param_c.inc"
423#include "com01_c.inc"
424#include "task_c.inc"
425#include "scr16_c.inc"
426C-----------------------------------------------
427C D u m m y A r g u m e n t s
428C-----------------------------------------------
429 INTEGER IPARG(NPARG,*),DD_IAD(NSPMD+1,*),WASZ,SIZ_WRITE_LOC(NSPGROUP+1)
430C-----------------------------------------------
431C L o c a l V a r i a b l e s
432C-----------------------------------------------
433 INTEGER NGF,NGL,NN,ITY,NEL,NG,JJ,
434 . P0ARSZ2,WASZ2
435C-----------------------------------------------
436 wasz = 0
437 IF (outp_spt( 1) == 1 ) THEN
438
439 ngf = 1
440 ngl = 0
441 DO nn=1,nspgroup
442 jj = 0
443 ngl = ngl + dd_iad(ispmd+1,nn)
444 DO ng = ngf, ngl
445 ity =iparg(5,ng)
446 IF(ity == 51) THEN
447 nel = iparg(2,ng)
448 jj = jj + 11*nel
449 ENDIF
450 ENDDO
451 ngf = ngl + 1
452 wasz = wasz + jj
453 siz_write_loc(nn) = jj
454 ENDDO
455 siz_write_loc(nspgroup+1) = wasz
456 ENDIF
457 RETURN
458 END
my_real
#define my_real
Definition cppsort.cpp:32
initbuf_mod
Definition initbuf.F:160
initbuf_mod::initbuf
subroutine initbuf(iparg, ng, mtn, llt, nft, iad, ity, npt, jale, ismstr, jeul, jtur, jthe, jlag, jmult, jhbe, jivf, mid, jpor, jcvt, jclose, jpla, irep, iint, igtyp, israt, isrot, icsen, isorth, isorthg, ifailure, jsms)
Definition initbuf.F:261
outp_sp_t
subroutine outp_sp_t(key, text, elbuf_tab, iparg, dd_iad, sizloc, sizp0, siz_wr)
Definition outp_sp_t.F:35
count_arsz_sptt
subroutine count_arsz_sptt(iparg, dd_iad, wasz, siz_write_loc)
Definition outp_sp_t.F:415
count_arsz_spt
subroutine count_arsz_spt(iparg, dd_iad, wasz, siz_write_loc)
Definition outp_sp_t.F:362
spmd_rgather9_1comm
subroutine spmd_rgather9_1comm(v, sizv, len, vp0, sizv0, adress)
Definition spmd_outp.F:1181