OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
upxfem_tagxp.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!|| upxfem_tagxp ../engine/source/elements/xfem/upxfem_tagxp.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| crk_tagxp3 ../engine/source/elements/xfem/crk_tagxp3.F
29!|| crk_tagxp4 ../engine/source/elements/xfem/crk_tagxp4.F
30!|| startimeg ../engine/source/system/timer.F
31!|| stoptimeg ../engine/source/system/timer.F
32!|| upenric3_n3 ../engine/source/elements/xfem/upenric3_nx.F
33!|| upenric3_n4 ../engine/source/elements/xfem/upenric3_nx.F
34!||--- uses -----------------------------------------------------
35!|| crackxfem_mod ../engine/share/modules/crackxfem_mod.F
36!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.f90
37!|| element_mod ../common_source/modules/elements/element_mod.F90
38!||====================================================================
39 SUBROUTINE upxfem_tagxp(XFEM_TAB,
40 . IPARG ,IXC ,NGROUC ,IGROUC ,IXTG ,
41 . IADC_CRK ,IEL_CRK ,INOD_CRK,ELCUTC ,NODEDGE ,
42 . ENRTAG ,CRKEDGE ,XEDGE4N ,XEDGE3N ,ITAB )
43C-----------------------------------------------
44C M o d u l e s
45C-----------------------------------------------
46 USE crackxfem_mod
47 USE elbufdef_mod
48 use element_mod , only : nixc,nixtg
49C-----------------------------------------------
50C I m p l i c i t T y p e s
51C-----------------------------------------------
52#include "implicit_f.inc"
53C-----------------------------------------------
54C C o m m o n B l o c k s
55C-----------------------------------------------
56#include "param_c.inc"
57#include "com01_c.inc"
58#include "com04_c.inc"
59#include "com_xfem1.inc"
60#include "task_c.inc"
61#include "vect01_c.inc"
62C-----------------------------------------------
63C D u m m y A r g u m e n t s
64C-----------------------------------------------
65 INTEGER IPARG(NPARG,*),IXC(NIXC,*),NGROUC,IGROUC(*),IEL_CRK(*),
66 . ELCUTC(2,*),IADC_CRK(*),INOD_CRK(*),IXTG(NIXTG,*),ITAB(*),
67 . NODEDGE(2,*),ENRTAG(NUMNOD,*),XEDGE4N(4,*),XEDGE3N(3,*)
68C
69 TYPE(elbuf_struct_), TARGET ,DIMENSION(NGROUP,NXEL) :: XFEM_TAB
70 TYPE (XFEM_EDGE_) , DIMENSION(*) :: CRKEDGE
71C-----------------------------------------------
72C L o c a l V a r i a b l e s
73C-----------------------------------------------
74 INTEGER IG,NG,JFT,JLT,NEL,NXLAY,IXFEM,OFF,ITG1,ITG2
75C=======================================================================
76 ITG1 = 1+numelc
77 itg2 = 1+4*ecrkxfec
78C-----------------------------------------------
79C SMP dynamic parallel loop
80!$OMP DO SCHEDULE(DYNAMIC,1)
81 DO ig = 1, ngrouc
82 ng = igrouc(ig)
83 off = iparg(8,ng)
84 ixfem = iparg(54,ng) * iparg(70,ng) ! IXFEM * ACTIVEXFEM
85 IF (off == 0 .and. ixfem > 0) THEN
86 IF (iddw > 0) CALL startimeg(ng)
87C---
88 ity = iparg(5,ng)
89 nel = iparg(2,ng)
90 nft = iparg(3,ng)
91 nxlay = iparg(59,ng)
92 lft = 1
93 llt = min(nvsiz,nel)
94 jft = lft
95 jlt = llt
96C---
97 IF (ity == 3) THEN
98 CALL crk_tagxp4(iparg ,ixc ,nft ,jft ,jlt ,
99 . elcutc ,iadc_crk,iel_crk,inod_crk,enrtag,
100 . nxlay ,crkedge ,xedge4n,itab )
101 ELSEIF (ity == 7) THEN
102 CALL crk_tagxp3(
103 . iparg ,ixtg ,nft ,jft ,jlt ,
104 . elcutc(1,itg1),iadc_crk(itg2),iel_crk(itg1),inod_crk,enrtag,
105 . nxlay ,crkedge ,xedge3n ,itab )
106 ENDIF
107C---
108 IF (iddw > 0) CALL stoptimeg(ng)
109 ENDIF
110 ENDDO
111!$OMP END DO
112C-----------------------------------------------
113C SMP dynamic parallel loop
114!$OMP DO SCHEDULE(DYNAMIC,1)
115 DO ig = 1, ngrouc
116 ng = igrouc(ig)
117 off = iparg(8,ng)
118 ixfem = iparg(54,ng) * iparg(70,ng) ! IXFEM * ACTIVEXFEM
119 IF (off == 0 .and. ixfem > 0) THEN
120 IF (iddw > 0) CALL startimeg(ng)
121C---
122 ity = iparg(5,ng)
123 nel = iparg(2,ng)
124 nft = iparg(3,ng)
125 lft = 1
126 llt = min(nvsiz,nel)
127 jft=lft
128 jlt=llt
129C---
130 IF (ity == 3) THEN
131 CALL upenric3_n4(xfem_tab(ng,1:nxel),
132 . iparg ,ixc ,nft ,jft ,jlt ,
133 . elcutc ,iadc_crk,iel_crk,inod_crk,ixfem,
134 . crkedge ,xedge4n )
135 ELSEIF (ity == 7) THEN
136 CALL upenric3_n3(xfem_tab(ng,1:nxel),
137 . iparg ,ixtg ,nft ,jft ,jlt ,
138 . elcutc(1,itg1),iadc_crk(itg2),iel_crk(itg1),inod_crk,ixfem,
139 . crkedge ,xedge3n )
140 ENDIF
141C---
142 IF (iddw > 0) CALL stoptimeg(ng)
143 ENDIF
144 ENDDO
145!$OMP END DO
146c-------------
147 RETURN
148 END
crk_tagxp3
subroutine crk_tagxp3(iparg, ixtg, nft, jft, jlt, elcutc, iad_crktg, iel_crktg, inod_crk, enrtag, nxlay, crkedge, xedge3n, itab)
Definition crk_tagxp3.F:35
crk_tagxp4
subroutine crk_tagxp4(iparg, ixc, nft, jft, jlt, elcutc, iadc_crk, iel_crk, inod_crk, enrtag, nxlay, crkedge, xedge4n, itab)
Definition crk_tagxp4.F:35
startimeg
subroutine startimeg(ng)
Definition timer.F:1371
stoptimeg
subroutine stoptimeg(ng)
Definition timer.F:1419
min
#define min(a, b)
Definition macros.h:20
crackxfem_mod
Definition crackxfem_mod.F:121
upenric3_n4
subroutine upenric3_n4(xfem_tab, iparg, ixc, nft, jft, jlt, elcutc, iadc_crk, iel_crk, inod_crk, ixfem, crkedge, xedge4n)
Definition upenric3_nx.F:36
upenric3_n3
subroutine upenric3_n3(xfem_tab, iparg, ixtg, nft, jft, jlt, elcutc, iad_crktg, iel_crktg, inod_crk, ixfem, crkedge, xedge3n)
Definition upenric3_nx.F:140
upxfem_tagxp
subroutine upxfem_tagxp(xfem_tab, iparg, ixc, ngrouc, igrouc, ixtg, iadc_crk, iel_crk, inod_crk, elcutc, nodedge, enrtag, crkedge, xedge4n, xedge3n, itab)
Definition upxfem_tagxp.F:43