OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
parsor_ply.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!|| parsor_ply ../engine/source/output/anim/generate/parsor_ply.F
25!||--- called by ------------------------------------------------------
26!|| genani ../engine/source/output/anim/generate/genani.F
27!||--- calls -----------------------------------------------------
28!|| get_ply_nod ../engine/source/output/anim/generate/anin_ply_init.F
29!|| get_ply_nodid ../engine/source/output/anim/generate/anin_ply_init.F
30!|| spmd_iget_partn_ply ../engine/source/mpi/anim/spmd_iget_partn_ply.F
31!|| spmd_iglob_partn ../engine/source/mpi/anim/spmd_iglob_partn.F
32!|| write_i_c ../common_source/tools/input_output/write_routines.c
33!||--- uses -----------------------------------------------------
34!|| element_mod ../common_source/modules/elements/element_mod.F90
35!|| my_alloc_mod ../common_source/tools/memory/my_alloc.F90
36!|| plyxfem_mod ../engine/share/modules/plyxfem_mod.F
37!||====================================================================
38 SUBROUTINE parsor_ply(NEL_PLY, X, D, XNORM, CDG,
39 . IPARG, IXC, IXTG, INVERT, EL2FA,
40 . MATER, IPARTC, NODGLOB, IDPLY, IADPLY,
41 . IADPLYG, PLYNUMC, NBF_PXFEMG)
42 USE plyxfem_mod
43 USE my_alloc_mod
44 use element_mod , only : nixc,nixtg
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 "task_c.inc"
56C-----------------------------------------------
57C D u m m y A r g u m e n t s
58C-----------------------------------------------
59C REAL
60 my_real
61 . x(*),d(*),xnorm(3,*),cdg(*)
62 INTEGER IXC(NIXC,*),IXTG(NIXTG,*),IPARG(NPARG,*),
63 . invert(*), el2fa(*),mater(*),
64 . ipartc(*),nodglob(*),
65 . nel_ply,idply,iadply(*),iadplyg(*),plynumc,
66 . nbf_pxfemg
67C-----------------------------------------------
68 INTEGER GET_PLY_NODID
69 EXTERNAL GET_PLY_NODID
70 INTEGER GET_PLY_NOD
71 EXTERNAL get_ply_nod
72C-----------------------------------------------
73C REAL
74 INTEGER II(4), IE, NG, ITY, LFT, LLT, N, I,
75 . NEL, IAD, NPAR, NFT,MTN,
76 . NN90,
77 . jj,buf,ishplyxfem,
78 . ipt,elc,plys,iply,ie_el2fa,
79 . npcc,idplyn
80 INTEGER EMPL
81 INTEGER,DIMENSION(:),ALLOCATABLE::NP
82 INTEGER,DIMENSION(:),ALLOCATABLE::PLYELEMS
83C-----------------------------------------------
84C PLY ELEMS INIT
85C-----------------------------------------------
86 CALL my_alloc(plyelems, numelc)
87 CALL my_alloc(np, plynumc*4)
88 npcc = 0
89 DO plys = 1,nplypart
90 iply = indx_ply(plys)
91 npcc = npcc + plyshell(iply)%PLYNUMSHELL
92 ENDDO
93
94 jj = 0
95 nel_ply = 0
96 ie=0
97 DO plys = 1,nplypart
98 iply = indx_ply(plys)
99 plyelems=0
100 DO i=1,plyshell(iply)%PLYNUMSHELL
101 ipt = plyshell(iply)%SHELLIPT(i)
102 elc = plyshell(iply)%SHID(i)
103 plyelems(elc)=ipt
104 ENDDO
105C-----------------------------------------------
106 ie_el2fa=0
107 npar = 0
108C
109C In SPMD you must send the info even qd we have 0 local elem
110 DO 490 ng=1,ngroup
111 mtn =iparg(1,ng)
112 nel =iparg(2,ng)
113 nft =iparg(3,ng)
114 iad =iparg(4,ng)
115 ity =iparg(5,ng)
116 ishplyxfem = iparg(50,ng)
117 lft=1
118 llt=nel
119cccc KPT =IPARG(6,NG)
120C-----------------------------------------------
121C COQUES
122C-----------------------------------------------
123 IF(ishplyxfem > 0) THEN
124 IF(ity==3)THEN
125c KPT =IPARG(6,NG)
126c IHBE = IPARG(23,NG)
127c ISTRAIN=IPARG(44,NG)
128c IEXPAN=IPARG(49,NG)
129c CALL CBUFAD (IAD-1,NEL,ISTRAIN,KPT,IHBE,0,0,0,IEXPAN,
130c . IPARG(1,NG))
131 DO 130 i=lft,llt
132 n = i + nft
133 IF(plyelems(n) == 0)GOTO 130
134cc OFF = BUFEL(IAD+I+SH_IH*NEL-1)
135 ii(1) = ixc(2,n)
136 ii(2) = ixc(3,n)
137 ii(3) = ixc(4,n)
138 ii(4) = ixc(5,n)
139 ie = ie + 1
140 ie_el2fa = ie_el2fa + 1
141c CALL FACNOR(X,D,II,XNORM,CDG,INVERT(IE))
142 IF (nspmd == 1) THEN
143 ii(1) = get_ply_nodid(iply,ii(1),idply)- 1
144 ii(2) = get_ply_nodid(iply,ii(2),idply)- 1
145 ii(3) = get_ply_nodid(iply,ii(3),idply)- 1
146 ii(4) = get_ply_nodid(iply,ii(4),idply)- 1
147 CALL write_i_c(ii,4)
148 ELSE
149 empl = get_ply_nod(iply,ii(1))
150 np(jj+1) = plynod(iply)%PLYNODID(empl)
151
152 empl = get_ply_nod(iply,ii(2))
153 np(jj+2) = plynod(iply)%PLYNODID(empl)
154
155 empl = get_ply_nod(iply,ii(3))
156 np(jj+3) = plynod(iply)%PLYNODID(empl)
157
158 empl = get_ply_nod(iply,ii(4))
159 np(jj+4) = plynod(iply)%PLYNODID(empl)
160
161 END IF
162 el2fa(nel_ply+ie_el2fa) = ie
163 jj = jj + 4
164 130 CONTINUE
165C-----------------------------------------------
166C COQUES 3 NOEUDS
167C-----------------------------------------------
168 ELSEIF(ity==7)THEN
169 ENDIF
170 ENDIF
171 490 CONTINUE
172C
173
174 iadply(plys) = ie
175 nel_ply = nel_ply + plyshell(iply)%PLYNUMSHELL
176 ENDDO
177
178 IF (nspmd > 1) THEN
179C build global part array on proc 0
180
181 IF (ispmd==0) THEN
182 CALL spmd_iglob_partn(iadply,nplypart,iadplyg,nplypart)
183 buf = nbf_pxfemg*4
184 idplyn = idply - 1
185 CALL spmd_iget_partn_ply(4,jj,np,nplypart,iadplyg,buf,1,idplyn)
186 ELSE
187 CALL spmd_iglob_partn(iadply,nplypart,iadplyg,nplypart)
188 buf=1
189 idplyn = idply - 1
190 CALL spmd_iget_partn_ply(4,jj,np,nplypart,iadplyg,buf,1,idplyn)
191 ENDIF
192
193 ENDIF
194 DEALLOCATE(plyelems)
195 DEALLOCATE(np)
196C-----------------------------------------------
197 RETURN
198 END
my_real
#define my_real
Definition cppsort.cpp:32
invert
subroutine invert(matrix, inverse, n, errorflag)
Definition hm_read_rbe3.F:1368
plyxfem_mod
Definition plyxfem_mod.F:31
plyxfem_mod::plynod
type(plynods), dimension(:), allocatable plynod
Definition plyxfem_mod.F:44
plyxfem_mod::indx_ply
integer, dimension(:), allocatable indx_ply
Definition plyxfem_mod.F:60
plyxfem_mod::nplypart
integer nplypart
Definition plyxfem_mod.F:59
plyxfem_mod::plyshell
type(plyshells), dimension(:), allocatable plyshell
Definition plyxfem_mod.F:56
parsor_ply
subroutine parsor_ply(nel_ply, x, d, xnorm, cdg, iparg, ixc, ixtg, invert, el2fa, mater, ipartc, nodglob, idply, iadply, iadplyg, plynumc, nbf_pxfemg)
Definition parsor_ply.F:42
spmd_iget_partn_ply
subroutine spmd_iget_partn_ply(size, nbf_l, np, nbpart, iadg, srbuf, iflag, idply)
Definition spmd_iget_partn_ply.F:38
spmd_iglob_partn
subroutine spmd_iglob_partn(iad, nbpart, iadg, sbuf)
Definition spmd_iglob_partn.F:39
write_i_c
void write_i_c(int *w, int *len)
Definition write_routines.c:645