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_routtines.c
33!||--- uses -----------------------------------------------------
34!|| my_alloc_mod ../common_source/tools/memory/my_alloc.F90
35!|| plyxfem_mod ../engine/share/modules/plyxfem_mod.F
36!||====================================================================
37 SUBROUTINE parsor_ply(NEL_PLY, X, D, XNORM, CDG,
38 . IPARG, IXC, IXTG, INVERT, EL2FA,
39 . MATER, IPARTC, NODGLOB, IDPLY, IADPLY,
40 . IADPLYG, PLYNUMC, NBF_PXFEMG)
41 USE plyxfem_mod
42 USE my_alloc_mod
43C-----------------------------------------------
44C I m p l i c i t T y p e s
45C-----------------------------------------------
46#include "implicit_f.inc"
47C-----------------------------------------------
48C C o m m o n B l o c k s
49C-----------------------------------------------
50#include "com01_c.inc"
51#include "com04_c.inc"
52#include "param_c.inc"
53#include "task_c.inc"
54C-----------------------------------------------
55C D u m m y A r g u m e n t s
56C-----------------------------------------------
57C REAL
58 my_real
59 . x(*),d(*),xnorm(3,*),cdg(*)
60 INTEGER IXC(NIXC,*),IXTG(NIXTG,*),IPARG(NPARG,*),
61 . invert(*), el2fa(*),mater(*),
62 . ipartc(*),nodglob(*),
63 . nel_ply,idply,iadply(*),iadplyg(*),plynumc,
64 . nbf_pxfemg
65C-----------------------------------------------
66 INTEGER GET_PLY_NODID
67 EXTERNAL GET_PLY_NODID
68 INTEGER GET_PLY_NOD
69 EXTERNAL get_ply_nod
70C-----------------------------------------------
71C REAL
72 my_real off
73 INTEGER II(4),IE,NG, ITY, LFT, LLT, N, I, J,
74 . IPRT, NEL, IAD, NPAR, NFT, IMID,IALEL,MTN,
75 . nn1,nn2,nn3,nn4,nn5,nn6,nn7,nn8,nn9,nn10,
76 . jj, k, sh_ih,buf,ishplyxfem,
77 . ipt,elc,plys,iply,ie_el2fa,
78 . npcc,idplyn
79 INTEGER EMPL
80 INTEGER,DIMENSION(:),ALLOCATABLE::NP
81 INTEGER,DIMENSION(:),ALLOCATABLE::PLYELEMS
82C-----------------------------------------------
83C PLY ELEMS INIT
84C-----------------------------------------------
85 CALL my_alloc(plyelems, numelc)
86 CALL my_alloc(np, plynumc*4)
87 npcc = 0
88 DO plys = 1,nplypart
89 iply = indx_ply(plys)
90 npcc = npcc + plyshell(iply)%PLYNUMSHELL
91 ENDDO
92
93 jj = 0
94 nel_ply = 0
95 ie=0
96 DO plys = 1,nplypart
97 iply = indx_ply(plys)
98 plyelems=0
99 DO i=1,plyshell(iply)%PLYNUMSHELL
100 ipt = plyshell(iply)%SHELLIPT(i)
101 elc = plyshell(iply)%SHID(i)
102 plyelems(elc)=ipt
103 ENDDO
104C-----------------------------------------------
105 ie_el2fa=0
106 npar = 0
107C
108C en spmd il faut envoyer l'info meme qd on a 0 elem local
109 DO 490 ng=1,ngroup
110 mtn =iparg(1,ng)
111 nel =iparg(2,ng)
112 nft =iparg(3,ng)
113 iad =iparg(4,ng)
114 ity =iparg(5,ng)
115 ishplyxfem = iparg(50,ng)
116 lft=1
117 llt=nel
118cccc KPT =IPARG(6,NG)
119C-----------------------------------------------
120C COQUES
121C-----------------------------------------------
122 IF(ishplyxfem > 0) THEN
123 IF(ity==3)THEN
124c KPT =IPARG(6,NG)
125c IHBE = IPARG(23,NG)
126c ISTRAIN=IPARG(44,NG)
127c IEXPAN=IPARG(49,NG)
128c CALL CBUFAD (IAD-1,NEL,ISTRAIN,KPT,IHBE,0,0,0,IEXPAN,
129c . IPARG(1,NG))
130 DO 130 i=lft,llt
131 n = i + nft
132 IF(plyelems(n) == 0)GOTO 130
133cc OFF = BUFEL(IAD+I+SH_IH*NEL-1)
134 ii(1) = ixc(2,n)
135 ii(2) = ixc(3,n)
136 ii(3) = ixc(4,n)
137 ii(4) = ixc(5,n)
138 ie = ie + 1
139 ie_el2fa = ie_el2fa + 1
140c CALL FACNOR(X,D,II,XNORM,CDG,INVERT(IE))
141 IF (nspmd == 1) THEN
142 ii(1) = get_ply_nodid(iply,ii(1),idply)- 1
143 ii(2) = get_ply_nodid(iply,ii(2),idply)- 1
144 ii(3) = get_ply_nodid(iply,ii(3),idply)- 1
145 ii(4) = get_ply_nodid(iply,ii(4),idply)- 1
146 CALL write_i_c(ii,4)
147 ELSE
148 empl = get_ply_nod(iply,ii(1))
149 np(jj+1) = plynod(iply)%PLYNODID(empl)
150
151 empl = get_ply_nod(iply,ii(2))
152 np(jj+2) = plynod(iply)%PLYNODID(empl)
153
154 empl = get_ply_nod(iply,ii(3))
155 np(jj+3) = plynod(iply)%PLYNODID(empl)
156
157 empl = get_ply_nod(iply,ii(4))
158 np(jj+4) = plynod(iply)%PLYNODID(empl)
159
160 END IF
161 el2fa(nel_ply+ie_el2fa) = ie
162 jj = jj + 4
163 130 CONTINUE
164C-----------------------------------------------
165C COQUES 3 NOEUDS
166C-----------------------------------------------
167 ELSEIF(ity==7)THEN
168 ENDIF
169 ENDIF
170 490 CONTINUE
171C
172
173 iadply(plys) = ie
174 nel_ply = nel_ply + plyshell(iply)%PLYNUMSHELL
175 ENDDO
176
177 IF (nspmd > 1) THEN
178C construction tableau global des parts sur p0
179
180 IF (ispmd==0) THEN
181 CALL spmd_iglob_partn(iadply,nplypart,iadplyg,nplypart)
182 buf = nbf_pxfemg*4
183 idplyn = idply - 1
184 CALL spmd_iget_partn_ply(4,jj,np,nplypart,iadplyg,buf,1,idplyn)
185 ELSE
186 CALL spmd_iglob_partn(iadply,nplypart,iadplyg,nplypart)
187 buf=1
188 idplyn = idply - 1
189 CALL spmd_iget_partn_ply(4,jj,np,nplypart,iadplyg,buf,1,idplyn)
190 ENDIF
191
192 ENDIF
193 DEALLOCATE(plyelems)
194 DEALLOCATE(np)
195C-----------------------------------------------
196 RETURN
197 END
my_real
#define my_real
Definition cppsort.cpp:32
invert
subroutine invert(matrix, inverse, n, errorflag)
Definition hm_read_rbe3.F:1354
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:41
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_routtines.c:645