OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i21main_tri.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "com04_c.inc"
#include "com08_c.inc"
#include "param_c.inc"
#include "task_c.inc"
#include "timeri_c.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i21main_tri (timers, ipari, x, nin, itask, weight, retri, num_imp, ind_imp, intstamp, mwag, intbuf_tab, nspmd)

Function/Subroutine Documentation

◆ i21main_tri()

subroutine i21main_tri ( type(timer_) timers,
integer, dimension(npari,ninter) ipari,
x,
integer nin,
integer itask,
integer, dimension(*) weight,
integer retri,
integer num_imp,
integer, dimension(*) ind_imp,
type(intstamp_data) intstamp,
integer, dimension(*) mwag,
type(intbuf_struct_) intbuf_tab,
integer, intent(in) nspmd )
Parameters
[in]nspmdnumber of processor

Definition at line 43 of file i21main_tri.F.

47C-----------------------------------------------
48C M o d u l e s
49C-----------------------------------------------
50 USE timer_mod
51 USE intstamp_mod
52 USE intbufdef_mod
53 use check_sorting_criteria_mod , only : check_sorting_criteria
54C-----------------------------------------------
55C I m p l i c i t T y p e s
56C-----------------------------------------------
57#include "implicit_f.inc"
58#include "comlock.inc"
59C-----------------------------------------------
60C C o m m o n B l o c k s
61C-----------------------------------------------
62#include "com04_c.inc"
63#include "com08_c.inc"
64#include "param_c.inc"
65#include "task_c.inc"
66#include "timeri_c.inc"
67 COMMON /i21mainc/result21,i_memg
68 INTEGER RESULT21,I_MEMG
69C-----------------------------------------------
70C D u m m y A r g u m e n t s
71C-----------------------------------------------
72 TYPE(TIMER_) :: TIMERS
73 integer, intent(in) :: nspmd !< number of processor
74 INTEGER NIN ,ITASK, RETRI, NUM_IMP ,IND_IMP(*),
75 . IPARI(NPARI,NINTER),
76 . WEIGHT(*), MWAG(*)
77 TYPE(INTSTAMP_DATA) INTSTAMP
78C REAL
79 my_real
80 . x(*)
81
82 TYPE(INTBUF_STRUCT_) INTBUF_TAB
83C-----------------------------------------------
84C L o c a l V a r i a b l e s
85C-----------------------------------------------
86 INTEGER J17_T,
87 . I, IP0, IP1, IP2, IP21, K11_T, I_SK_OLD,
88 . ADD1, NB_N_B, INIT, NOINT, INACTI, MULTIMP, IGAP, I_STOK,
89 . INTTH , ADXLOC
90 INTEGER
91 . NRTM_T,ESHIFT,ILD, NCONTACT, NRTM_L,
92 . I_MEM,CAND_N_OLD, IROT
93C REAL
94 my_real
95 . gap,maxbox,minbox,tzinf,
96 . xmax, ymax, zmax, xmin, ymin, zmin, gapmin, gapmax,
97 . dd(3)
98 INTEGER :: NRTM,NSN,NTY,NMN
99 logical :: need_computation
100C-----------------------------------------------
101 ! --------------
102 ! check if the current interface needs to be sorted
103 call check_sorting_criteria( need_computation,nin,npari,nspmd,
104 . itask,ipari(1,nin),tt,intbuf_tab )
105 if( .not.need_computation ) return
106 ! --------------
107 i_mem = 0
108 i_memg = 0
109C
110 nrtm =ipari(4,nin)
111 nsn =ipari(5,nin)
112 nty =ipari(7,nin)
113 noint =ipari(15,nin)
114 nmn =ipari(8,nin)
115 intth = ipari(47,nin)
116
117 inacti =ipari(22,nin)
118 multimp=ipari(23,nin)
119
120 ncontact=multimp*nsn
121C
122 gap =intbuf_tab%VARIABLES(2)
123 gapmin=intbuf_tab%VARIABLES(13)
124 gapmax=intbuf_tab%VARIABLES(16)
125C
126 IF(nsn > 0) THEN
127 adxloc = 1+3*(nsn+2)
128 ELSE
129 adxloc = 1
130 ENDIF
131C
132C -------------------------------------------------------------
133C
134 retri=1
135C
136 dd(1)=intstamp%D(1)
137 dd(2)=intstamp%D(2)
138 dd(3)=intstamp%D(3)
139C sauvegarde de XSAV (tableau BUFIN(JD(19)))
140 irot=intstamp%IROT
141 CALL i21xsave(
142 1 x ,intbuf_tab%NSV,nsn ,itask ,intbuf_tab%XSAV,
143 2 dd ,irot ,intstamp%ROT,intstamp%BRACKET)
144C -------------------------------------------------------------
145C
146 maxbox = intbuf_tab%VARIABLES(9)
147 minbox = intbuf_tab%VARIABLES(12)
148 tzinf = intbuf_tab%VARIABLES(8)
149C
150C -------------------------------------------------------------
151 i_sk_old=0
152 intbuf_tab%I_STOK(1)=zero
153C BARRIER II_STOK et RESULT
154 result21 = 0
155C
156 cand_n_old = intbuf_tab%I_STOK(1)
157 40 CONTINUE
158 CALL my_barrier
159C
160 inacti = ipari(22,nin)
161C -------------------------------------------------------------
162 nrtm_t=nrtm/nthread
163C eshift : decalage sur cand_e
164 eshift = itask*nrtm_t
165C il faut conserver K11 global
166
167 k11_t = 1+4*itask*nrtm_t
168
169 j17_t = 1+itask*nrtm_t
170 IF(itask==nthread-1)nrtm_t=nrtm-(nthread-1)*(nrtm/nthread)
171 init = 1
172 ild = 0
173 50 CONTINUE
174C
175 IF (imonm > 0) CALL startime(timers,30)
176 IF(nrtm_t/=0) THEN
177 CALL i21buce(
178 1 x ,intbuf_tab%IRECTM(k11_t),intbuf_tab%NSV ,ipari(22,nin),
179 2 nrtm_t ,nsn ,intbuf_tab%CAND_E ,intbuf_tab%CAND_N ,gap,
180 3 noint ,intbuf_tab%I_STOK(1) ,tzinf ,maxbox ,minbox ,
181 4 ncontact ,xmin ,xmax ,ymin ,
182 5 ymax ,zmin ,zmax ,nb_n_b ,eshift ,
183 6 ild ,init ,weight ,intbuf_tab%STFNS ,nin ,
184 7 intbuf_tab%STF(j17_t) ,ipari(21,nin),intbuf_tab%GAP_S,gapmin,gapmax ,
185 8 ipari(39,nin),num_imp ,intbuf_tab%XM0 ,intbuf_tab%NOD_NORMAL ,
186 9 intbuf_tab%VARIABLES(23) ,intbuf_tab%VARIABLES(22) ,intbuf_tab%VARIABLES(27) ,
187 . intbuf_tab%VARIABLES(28) ,intbuf_tab%VARIABLES(29) ,
188 a nrtm_l ,intbuf_tab%XSAV(adxloc),i_mem,intbuf_tab%VARIABLES(32) ,
189 . nmn ,
190 b intth ,intbuf_tab%MNDD ,intbuf_tab%MSR_L ,itask,intbuf_tab%IRECTM,
191 c ipari(48,nin) ,intbuf_tab%VARIABLES(46))
192
193 ENDIF
194
195C Upgrade MultiMP
196 IF (i_mem == 2)THEN
197#include "lockon.inc"
198 i_memg = i_mem
199#include "lockoff.inc"
200 ENDIF
201C New barrier needed for Dynamic MultiMP
202 CALL my_barrier
203
204 IF(i_memg /=0)THEN
205!$OMP SINGLE
206 multimp = ipari(23,nin) + 4
207 CALL upgrade_multimp(nin,multimp,intbuf_tab)
208!$OMP END SINGLE
209 i_memg = 0
210 i_mem = 0
211 intbuf_tab%I_STOK(1) =cand_n_old
212 multimp=ipari(23,nin)
213 ncontact=multimp*nsn
214 GOTO 40
215 ENDIF
216
217 CALL update_struct_int21(
218 1 intbuf_tab%CAND_E, intbuf_tab%I_STOK(1), ipari(48,nin) , nin , nmn,
219 2 intth, intbuf_tab%MNDD, intbuf_tab%MSR_L , intbuf_tab%IRECTM )
220
221 IF (imonm > 0) CALL stoptime(timers,30)
222C
223#include "lockon.inc"
224 intbuf_tab%VARIABLES(9) = min(maxbox,intbuf_tab%VARIABLES(9))
225 intbuf_tab%VARIABLES(12) = min(minbox,intbuf_tab%VARIABLES(12))
226 intbuf_tab%VARIABLES(8) = min(tzinf,intbuf_tab%VARIABLES(8))
227c +1 en 21 :
228 intbuf_tab%VARIABLES(5) = ep30
229 result21 = result21 + ild
230#include "lockoff.inc"
231C--------------------------------------------------------------
232C--------------------------------------------------------------
233 CALL my_barrier
234 IF (result21/=0) THEN
235 CALL my_barrier
236 IF (itask==0) THEN
237 intbuf_tab%I_STOK(1) = i_sk_old
238 result21 = 0
239 ENDIF
240 CALL my_barrier
241 ild = 0
242 init = 0
243 maxbox = intbuf_tab%VARIABLES(9)
244 minbox = intbuf_tab%VARIABLES(12)
245 tzinf = intbuf_tab%VARIABLES(8)
246
247 CALL my_barrier
248 GOTO 50
249 ENDIF
250C-----------------------------------------------------------
251C Remise a 0 de IFPEN pour les nds qui ne sont plus candidats
252C--------------------------------------------------------------
253 IF(itask==0)THEN
254 ip0 = 1
255 ip1 = ip0 + nsn + 3
256C MWA = MWAG SUR TASK 0
257 i_stok = intbuf_tab%I_STOK(1)
258 CALL i21fpen(
259 1 nsn ,i_stok ,intbuf_tab%CAND_N,intbuf_tab%CAND_E,
260 . intbuf_tab%PENIS,
261 2 intbuf_tab%FTSAVX,intbuf_tab%FTSAVY,intbuf_tab%FTSAVZ,mwag(ip0) ,
262 . intbuf_tab%IFPEN,
263 3 inacti)
264 ENDIF
265C--------------------------------------------------------------
266
267C
268 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
i21buce
subroutine i21buce(x, irect, nsv, inacti, nrtm, nsn, cand_e, cand_n, gap, noint, ii_stok, tzinf, maxbox, minbox, ncontact, xmin, xmax, ymin, ymax, zmin, zmax, nb_n_b, eshift, ild, init, weight, stfn, nin, stf, igap, gap_s, gapmin, gapmax, icurv, num_imp, xm0, nod_normal, depth, margeref, lxm, lym, lzm, nrtm_l, xloc, i_mem, drad, nmn, intth, mndd, msr_l, itask, irectt, iform, dgapload)
Definition i21buce.F:49
i21fpen
subroutine i21fpen(nsn, i_stok, cand_n, cand_e, peni, ftxsav, ftysav, ftzsav, cand_tag, ifpen, inacti)
Definition i21fpen.F:32
i21xsave
subroutine i21xsave(x, nsv, nsn, itask, xsav, dd, irot, rot, bracket)
Definition i21xsave.F:31
ymax
subroutine ymax(idn, fac, npc, pld, stiffmin, stiffmax, stiffini, stiffavg)
Definition law100_upd.F:272
min
#define min(a, b)
Definition macros.h:20
intstamp_mod
Definition intstamp_mod.F:38
my_barrier
subroutine my_barrier
Definition machine.F:31
startime
subroutine startime(event, itask)
Definition timer.F:93
stoptime
subroutine stoptime(event, itask)
Definition timer.F:135
update_struct_int21
subroutine update_struct_int21(cand_e, ii_stok, iform, nin, nmn, intth, mndd, msr_l, irectt)
Definition update_struct_int21.F:37
upgrade_multimp
subroutine upgrade_multimp(ni, multimp_parameter, intbuf_tab)
Definition upgrade_multimp.F:48