OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i7mainf.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!|| i7mainf ../engine/source/interfaces/int07/i7mainf.F
25!||--- called by ------------------------------------------------------
26!|| intfop2 ../engine/source/interfaces/interf/intfop2.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| frictionparts_model_isot ../engine/source/interfaces/int07/frictionparts_model.F
31!|| frictionparts_model_ortho ../engine/source/interfaces/int07/frictionparts_model.F
32!|| i17main_pena ../engine/source/interfaces/int17/i17main_pena.F
33!|| i18dst3 ../engine/source/interfaces/int18/i18dst3.F
34!|| i18for3 ../engine/source/interfaces/int18/i18for3.F
35!|| i7ass33 ../engine/source/interfaces/int07/i7ass3.F
36!|| i7cdcor3 ../engine/source/interfaces/int07/i7cdcor3.F
37!|| i7cor3 ../engine/source/interfaces/int07/i7cor3.F
38!|| i7curvsz ../engine/source/interfaces/int07/i7curv.F
39!|| i7dst3 ../engine/source/interfaces/int07/i7dst3.F
40!|| i7for3 ../engine/source/interfaces/int07/i7for3.F
41!|| i7norm ../engine/source/interfaces/int07/i7curv.F
42!|| i7norme ../engine/source/interfaces/int07/i7rcurv.F
43!|| i7normn ../engine/source/interfaces/int07/i7rcurv.F
44!|| i7normnp ../engine/source/interfaces/int07/i7rcurv.F
45!|| i7normp ../engine/source/interfaces/int07/i7curv.F
46!|| i7rcurv ../engine/source/interfaces/int07/i7rcurv.F
47!|| i7therm ../engine/source/interfaces/int07/i7therm.f
48!|| my_barrier ../engine/source/system/machine.f
49!|| spmd_exch_n ../engine/source/mpi/generic/spmd_exch_n.F
50!|| spmd_i7curvsz ../engine/source/mpi/interfaces/spmd_i7curvsz.F
51!|| startime ../engine/source/system/timer_mod.F90
52!|| stoptime ../engine/source/system/timer_mod.F90
53!|| sum_6_float_sens ../engine/source/system/parit.F
54!||--- uses -----------------------------------------------------
55!|| ale_connectivity_mod ../common_source/modules/ale/ale_connectivity_mod.F
56!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
57!|| groupdef_mod ../common_source/modules/groupdef_mod.F
58!|| h3d_mod ../engine/share/modules/h3d_mod.F
59!|| i18dst3_mod ../engine/source/interfaces/int18/i18dst3.F
60!|| i18for3_mod ../engine/source/interfaces/int18/i18for3.f
61!|| intbuf_fric_mod ../common_source/modules/interfaces/intbuf_fric_mod.F90
62!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
63!|| message_mod ../engine/share/message_module/message_mod.F
64!|| multi_fvm_mod ../common_source/modules/ale/multi_fvm_mod.F90
65!|| output_mod ../common_source/modules/output/output_mod.f90
66!|| timer_mod ../engine/source/system/timer_mod.F90
67!||====================================================================
68 SUBROUTINE i7mainf(OUTPUT, TIMERS,
69 1 IPARI ,X ,A ,ALE_CONNECTIVITY ,XCELL ,
70 2 ICODT ,FSAV ,V ,MS ,DT2T ,
71 3 NELTST ,ITYPTST ,ITAB ,STIFN ,FSKYI ,
72 4 ISKY ,FCONT ,NIN ,LINDMAX ,KINET ,
73 5 JTASK ,NB_JLT ,NB_JLT_NEW ,NB_STOK_N ,ELBUF_TAB ,
74 6 NISKYFI ,NEWFRONT ,NSTRF ,SECFCUM ,IGROUPS ,
75 7 ICONTACT ,VISCN ,NUM_IMP ,
76 9 NS_IMP ,NE_IMP ,IND_IMP ,FSAVSUB ,NRTMDIM ,
77 A IGRBRIC ,
78 B EMINX ,IXS ,IXS16 ,IXS20 ,FNCONT ,
79 C FTCONT ,IAD_ELEM ,FR_ELEM ,RCONTACT ,ACONTACT ,
80 D PCONTACT ,TEMP ,FTHE ,FTHESKYI ,
81 E PM ,IPARG ,IAD17 ,MSKYI_SMS ,ISKYI_SMS ,
82 F NODNX_SMS ,MS0 ,QFRICINT ,NPC ,TF ,
83 G CONDN ,CONDNSKYI ,INTBUF_TAB ,NODADT_THERM ,THEACCFACT ,
84 H FBSAV6 ,ISENSINT ,DIMFB ,IXIG3D ,KXIG3D ,
85 I WIGE ,KNOT ,IGEO ,MULTI_FVM ,H3D_DATA ,
86 J INTBUF_FRIC_TAB,KNOTLOCPC ,KNOTLOCEL ,ITASK ,TAGNCONT ,
87 K KLOADPINTER ,LOADPINTER ,LOADP_HYD_INTER,DGAPLOADINT ,S_LOADPINTER,
88 L INTEREFRIC ,S_XCELL_REMOTE,XCELL_REMOTE )
89C======================================================================
90C-----------------------------------------------
91C M o d u l e s
92C-----------------------------------------------
93 USE timer_mod
94 USE elbufdef_mod
95 USE message_mod
96 USE intbufdef_mod
97 USE multi_fvm_mod
98 USE h3d_mod
99 USE intbuf_fric_mod
100 USE groupdef_mod
101 USE ale_connectivity_mod
102 USE output_mod
103 USE i18dst3_mod , ONLY : i18dst3
104 USE i18for3_mod , ONLY : i18for3
105C-----------------------------------------------
106C I m p l i c i t T y p e s
107C-----------------------------------------------
108#include "implicit_f.inc"
109C-----------------------------------------------
110C G l o b a l P a r a m e t e r s
111C-----------------------------------------------
112#include "mvsiz_p.inc"
113C-----------------------------------------------
114C C o m m o n B l o c k s
115C-----------------------------------------------
116#include "com01_c.inc"
117#include "com04_c.inc"
118#include "com08_c.inc"
119#include "param_c.inc"
120#include "warn_c.inc"
121#include "task_c.inc"
122#include "parit_c.inc"
123#include "timeri_c.inc"
124#include "impl1_c.inc"
125#include "macro.inc"
126#include "comlock.inc"
127#include "tabsiz_c.inc"
128C-----------------------------------------------
129C D u m m y A r g u m e n t s
130C-----------------------------------------------
131 TYPE(timer_), INTENT(INOUT) :: TIMERS
132 TYPE(OUTPUT_), INTENT(INOUT) :: OUTPUT
133 TYPE(t_ale_connectivity), INTENT(IN) :: ALE_CONNECTIVITY
134 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
135 INTEGER NELTST,ITYPTST,NIN,NEWFRONT,NSTRF(*),NRTMDIM, IAD17, ISENSINT(*),DIMFB,IGROUPS(NUMELS)
136 INTEGER IPARI(NPARI,NINTER), ICODT(*),ICONTACT(*),ITAB(*), ISKY(*), KINET(*),IPARG(NPARG,*),TAGNCONT(NLOADP_HYD_INTER,NUMNOD)
137 INTEGER NB_JLT,NB_JLT_NEW,NB_STOK_N,JTASK,NISKYFI, LINDMAX,ITASK
138 INTEGER NUM_IMP,NS_IMP(*),NE_IMP(*),IND_IMP(*)
139 INTEGER IXS(*) ,IXS16(*) ,IXS20(*)
140 INTEGER IAD_ELEM(2,*),FR_ELEM(*), ISKYI_SMS(*), NODNX_SMS(*),NPC(*),KXIG3D(NIXIG3D,*),IXIG3D(*),IGEO(NPROPGI,*)
141 INTEGER, INTENT(IN) :: S_LOADPINTER
142 INTEGER, INTENT(IN) :: KLOADPINTER(NINTER+1),LOADPINTER(S_LOADPINTER),LOADP_HYD_INTER(NLOADP_HYD)
143 INTEGER, INTENT(IN) :: INTEREFRIC
144 INTEGER, INTENT(IN) :: NODADT_THERM
145 my_real, intent(in) :: theaccfact
146 my_real, INTENT(IN) :: DGAPLOADINT(S_LOADPINTER)
147 my_real EMINX(*), XCELL(3,SXCELL)
148 my_real DT2T,X(3,*), A(3,*), FSAV(*), V(3,*),MS(*),STIFN(*),FSKYI(LSKYI,4),FCONT(3,*),MS0(*),
149 . SECFCUM(7,NUMNOD,NSECT),VISCN(*), FSAVSUB(*),
150 . FNCONT(3,*), FTCONT(3,*), RCONTACT(*), ACONTACT(*),
151 . pcontact(*),temp(*),fthe(*),ftheskyi(lskyi),pm(npropm,*),
152 . mskyi_sms(*), qfricint(*),tf(*),condn(*),condnskyi(lskyi),wige(*),knot(*)
153 my_real knotlocpc(*),knotlocel(*)
154 DOUBLE PRECISION FBSAV6(12,6,DIMFB)
155 INTEGER, INTENT(in) :: S_XCELL_REMOTE
156 my_real, DIMENSION(S_XCELL_REMOTE), INTENT(in) :: xcell_remote
157 TYPE(intbuf_struct_) INTBUF_TAB
158 TYPE(multi_fvm_struct), INTENT(INOUT) :: MULTI_FVM
159 TYPE(h3d_database) :: H3D_DATA
160 TYPE(intbuf_fric_struct_), TARGET, DIMENSION(NINTERFRIC) :: INTBUF_FRIC_TAB
161 TYPE (GROUP_) , DIMENSION(NGRBRIC) :: IGRBRIC
162C-----------------------------------------------
163C L o c a l V a r i a b l e s
164C-----------------------------------------------
165 INTEGER I, J, I_STOK, JLT_NEW, JLT , NFT, IVIS2,
166 . ibc, isecin, ibag, iadm,
167 . igap, inacti, ifq, mfrot, igsti, nisub, itied,
168 . jlt_tied, nb_loc, i_stok_loc,debut,
169 . ilagm, lenr, lent, maxcc,intth,iform, h, ierror,
170 . sym_flag_type19, intfric ,nsetprts ,npartfric ,iorthfric,
171 . nforth ,nfisot
172 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
173 . nsvg(mvsiz), cn_loc(mvsiz),ce_loc(mvsiz),
174 . cand_n_n(mvsiz),cand_e_n(mvsiz),kini(mvsiz),
175 . index2(lindmax),
176 . isdsiz(nspmd+1),ircsiz(nspmd+1),itag(numnod),
177 . ieleci(mvsiz), nsms(mvsiz), ipartfricsi(mvsiz),
178 . ipartfricmi(mvsiz),indexisot(mvsiz),indexorth(mvsiz),
179 . irep_fricmi(mvsiz)
180 my_real startt, fric, gap, stopt,visc,viscf,stiglo,gapmin, kmin, kmax, gapmax,rstif,fheats,fheatm,tint,dtmini
181C-----------------------------------------------
182 my_real nx1(mvsiz), nx2(mvsiz), nx3(mvsiz), nx4(mvsiz),
183 . ny1(mvsiz), ny2(mvsiz), ny3(mvsiz), ny4(mvsiz),
184 . nz1(mvsiz), nz2(mvsiz), nz3(mvsiz), nz4(mvsiz),
185 . lb1(mvsiz), lb2(mvsiz), lb3(mvsiz), lb4(mvsiz),
186 . lc1(mvsiz), lc2(mvsiz), lc3(mvsiz), lc4(mvsiz),
187 . p1(mvsiz), p2(mvsiz), p3(mvsiz), p4(mvsiz),
188 . x1(mvsiz), x2(mvsiz), x3(mvsiz), x4(mvsiz),
189 . y1(mvsiz), y2(mvsiz), y3(mvsiz), y4(mvsiz),
190 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
191 . xi(mvsiz), yi(mvsiz), zi(mvsiz), stif(mvsiz),
192 . h1(mvsiz), h2(mvsiz), h3(mvsiz), h4(mvsiz),
193 . gapv(mvsiz),vxi(mvsiz),vyi(mvsiz),vzi(mvsiz),msi(mvsiz),
194 . tempi(mvsiz),phi(mvsiz),areasi(mvsiz),
195 . efrict(mvsiz),fni(mvsiz),condint(mvsiz),msi_append(2*mvsiz)
196 my_real rcurvi(mvsiz), anglmi(mvsiz), anglt, padm
197 INTEGER NRTMFT, NRTMLT, NMNFT, NMNLT, NRADM,IRSTH,IFRIC
198 my_real NNX1(MVSIZ), NNX2(MVSIZ), NNX3(MVSIZ), NNX4(MVSIZ),
199 . NNY1(MVSIZ), NNY2(MVSIZ), NNY3(MVSIZ), NNY4(MVSIZ),
200 . nnz1(mvsiz), nnz2(mvsiz), nnz3(mvsiz), nnz4(mvsiz),
201 . cmaj(mvsiz)
202 INTEGER ICURV,SFSAVPARIT,NINLOADP
203 my_real DRAD, DRAD2, FRAD,XTHE,XFRIC,DGAPLOADP
204 my_real, DIMENSION(:,:,:), ALLOCATABLE :: FSAVPARIT
205 my_real KS(MVSIZ),K1(MVSIZ),K2(MVSIZ),K3(MVSIZ),K4(MVSIZ),
206 . CS(MVSIZ),C1(MVSIZ),C2(MVSIZ),C3(MVSIZ),C4(MVSIZ),
207 . KT(MVSIZ),C(MVSIZ),CF(MVSIZ),
208 . FXI(MVSIZ), FYI(MVSIZ), FZI(MVSIZ),
209 . FX1(MVSIZ), FX2(MVSIZ), FX3(MVSIZ), FX4(MVSIZ),
210 . FY1(MVSIZ), FY2(MVSIZ), FY3(MVSIZ), FY4(MVSIZ),
211 . FZ1(MVSIZ), FZ2(MVSIZ), FZ3(MVSIZ), FZ4(MVSIZ),
212 . PHI1(MVSIZ),PHI2(MVSIZ),PHI3(MVSIZ),PHI4(MVSIZ)
213 my_real XFILTR_FRIC,FRIC_COEFS(MVSIZ,10),VISCFFRIC(MVSIZ),FRICC(MVSIZ),
214 . FRIC_COEFS2(MVSIZ,10),VISCFFRIC2(MVSIZ),FRICC2(MVSIZ),
215 . DIR1(MVSIZ,3),DIR2(MVSIZ,3),DIR_FRICMI(MVSIZ,2)
216 INTEGER, DIMENSION(:) ,POINTER :: TABCOUPLEPARTS_FRIC
217 INTEGER, DIMENSION(:) ,POINTER :: TABPARTS_FRIC
218 INTEGER, DIMENSION(:) ,POINTER :: ADPARTS_FRIC
219 INTEGER, DIMENSION(:) ,POINTER :: IFRICORTH
220 my_real, DIMENSION(:) ,POINTER :: TABCOEF_FRIC
221 INTEGER,TARGET, DIMENSION(1):: TABCOUPLEPARTS_FRIC_BID
222 INTEGER,TARGET, DIMENSION(1):: TABPARTS_FRIC_BID
223 INTEGER,TARGET, DIMENSION(1):: ADPARTS_FRIC_BID
224 INTEGER,TARGET, DIMENSION(1):: IFRICORTH_BID
225 my_real,TARGET, DIMENSION(1):: tabcoef_fric_bid
226 INTEGER :: NSN, NTY, NOINT
227 INTEGER :: NRTM
228 LOGICAL TYPE18
229 INTEGER :: IDIR
230C-----------------------------------------------
231C S o u r c e L i n e s
232C-----------------------------------------------
233 type18 = .false.
234 nrtm = ipari(4,nin)
235 nsn = ipari(5,nin)
236 nty = ipari(7,nin)
237 ibc = ipari(11,nin)
238 ivis2 = ipari(14,nin)
239 IF(ipari(33,nin) == 1) RETURN !precondition
240 noint = ipari(15,nin)
241 igap = ipari(21,nin)
242 inacti = ipari(22,nin) !inact=7 <=> /INTER/TYPE18
243 isecin = ipari(28,nin)
244 mfrot = ipari(30,nin)
245 ifq = ipari(31,nin)
246 ibag = ipari(32,nin)
247 igsti = ipari(34,nin)
248 nisub = ipari(36,nin)
249 icurv = ipari(39,nin)
250 itied = ipari(85,nin)
251 sym_flag_type19 = ipari(71,nin)
252C heat interface
253 intth = ipari(47,nin)
254 iform = ipari(48,nin)
255 irsth = ipari(42,nin)
256C adaptive meshing
257 iadm = ipari(44,nin)
258 nradm = ipari(49,nin)
259 padm = intbuf_tab%VARIABLES(24)
260 anglt = intbuf_tab%VARIABLES(25)
261 stiglo = -intbuf_tab%STFAC(1)
262 startt = intbuf_tab%VARIABLES(3)
263 stopt = intbuf_tab%VARIABLES(11)
264 IF(startt > tt) RETURN !precondition
265 IF(tt > stopt) RETURN !precondition
266 fric = intbuf_tab%VARIABLES(1)
267 gap = intbuf_tab%VARIABLES(2)
268 gapmin = intbuf_tab%VARIABLES(13)
269 visc = intbuf_tab%VARIABLES(14)
270 viscf= intbuf_tab%VARIABLES(15)
271 ifric = 0
272 IF(intth > 0) ifric =ipari(50,nin)
273 xfric = intbuf_tab%VARIABLES(34)
274 gapmax = intbuf_tab%VARIABLES(16)
275 kmin = intbuf_tab%VARIABLES(17)
276 kmax = intbuf_tab%VARIABLES(18)
277 rstif = intbuf_tab%VARIABLES(20)
278 fheats = intbuf_tab%VARIABLES(21)
279 tint = intbuf_tab%VARIABLES(22)
280 dtmini = intbuf_tab%VARIABLES(41)
281 fheatm = intbuf_tab%VARIABLES(23)
282 xthe = intbuf_tab%VARIABLES(33)
283 frad = zero
284 drad = zero
285 drad2 = zero
286 IF(nty == 7)THEN
287 frad = intbuf_tab%VARIABLES(31)
288 drad = intbuf_tab%VARIABLES(32)
289 drad2 = drad*drad
290 ENDIF
291
292C--- Corresponding Friction model
293 intfric=ipari(72,nin)
294 iorthfric = 0
295 nsetprts = 0
296 npartfric = 0
297 xfiltr_fric = zero
298
299 IF(intfric /= 0) THEN
300 tabcoupleparts_fric => intbuf_fric_tab(intfric)%TABCOUPLEPARTS_FRIC
301 tabcoef_fric => intbuf_fric_tab(intfric)%TABCOEF_FRIC
302 tabparts_fric => intbuf_fric_tab(intfric)%TABPARTS_FRIC
303 adparts_fric => intbuf_fric_tab(intfric)%ADPARTS_FRIC
304 xfiltr_fric = intbuf_fric_tab(intfric)%XFILTR_FRIC
305 nsetprts = intbuf_fric_tab(intfric)%NSETPRTS
306 npartfric = intbuf_fric_tab(intfric)%S_TABPARTS_FRIC
307 iorthfric = intbuf_fric_tab(intfric)%IORTHFRIC
308 ifricorth => intbuf_fric_tab(intfric)%IFRICORTH
309 ELSE
310 tabcoupleparts_fric => tabcoupleparts_fric_bid
311 tabparts_fric => tabparts_fric_bid
312 tabcoef_fric => tabcoef_fric_bid
313 adparts_fric => adparts_fric_bid
314 ifricorth => ifricorth_bid
315 IF (ifq/=0) xfiltr_fric = intbuf_tab%XFILTR(1)
316 ENDIF
317
318 ninloadp = ipari(95,nin) ! load pressure related to inter
319 dgaploadp = intbuf_tab%VARIABLES(46)
320
321 IF(nty == 7)THEN
322 IF(inacti == 7) type18 = .true.
323 IF(type18 .AND. ipari(34,nin) == -2) RETURN
324!----------------------------------------------------
325! Compute nodal normals using quadratic curvature
326!----------------------------------------------------
327 IF(icurv == 3)THEN
328 CALL my_barrier()
329 IF(jtask == 1)THEN
330 ALLOCATE(intbuf_tab%NODNORM_NORMAL (3,numnod))
331 IF(iparit == 0)THEN
332 CALL i7norm(ipari(4,nin),intbuf_tab%IRECTM,numnod,x,intbuf_tab%NODNORM_NORMAL,ipari(6,nin),intbuf_tab%MSR)
333 IF(nspmd > 1)THEN
334 lenr = iad_elem(1,nspmd+1)-iad_elem(1,1)
335 CALL spmd_exch_n(intbuf_tab%NODNORM_NORMAL,iad_elem,fr_elem,lenr)
336 END IF
337 ELSE
338! parallel assembly(parith/on using spmd): optimize if necessary
339 lenr = iad_elem(1,nspmd+1)-iad_elem(1,1)
340 IF(nspmd > 1)THEN
341! memory space estimation for parith/on(spmd communication)
342 CALL spmd_i7curvsz(
343 1 ipari(4,nin),intbuf_tab%IRECTM,numnod,iad_elem,fr_elem,
344 2 isdsiz ,ircsiz ,itag ,lenr ,lent ,
345 3 maxcc )
346 ELSE
347! Memory space estimation for parith/on (without spmd)
348 CALL i7curvsz(
349 1 ipari(4,nin),intbuf_tab%IRECTM,numnod,itag,lent,maxcc)
350 END IF
351 CALL i7normp(
352 1 ipari(4,nin),intbuf_tab%IRECTM,numnod ,x ,intbuf_tab%NODNORM_NORMAL,
353 2 ipari(6,nin),intbuf_tab%MSR ,lent ,maxcc,isdsiz ,
354 3 ircsiz ,iad_elem ,fr_elem,itag )
355 END IF
356 END IF
357 CALL my_barrier()
358 ENDIF
359!----------------------------------------------------
360! Curvature radius: compute normalized nodal normals
361! Warning / Precondition : IADM /= 0 and ICURV /= 0 cannot be used together (Starter error).
362!----------------------------------------------------
363 IF(iadm /= 0)THEN
364 CALL my_barrier()
365 IF(jtask==1)THEN
366 ALLOCATE(intbuf_tab%MODRCURV(nrtmdim),intbuf_tab%MODANGLM(nrtmdim))
367 ALLOCATE(intbuf_tab%NODNORM_NORMAL (3,numnod))
368 IF(iparit == 0)THEN
369 CALL i7normn(
370 . ipari(4,nin),intbuf_tab%IRECTM,numnod,x ,intbuf_tab%NODNORM_NORMAL,
371 . ipari(6,nin),intbuf_tab%MSR)
372 IF(nspmd > 1)THEN
373 lenr = iad_elem(1,nspmd+1)-iad_elem(1,1)
374 CALL spmd_exch_n(intbuf_tab%NODNORM_NORMAL,iad_elem,fr_elem,lenr)
375 END IF
376 ELSE
377C Traitement d'assemblage parith/on spmd a optimiser si besoin
378 lenr = iad_elem(1,nspmd+1)-iad_elem(1,1)
379 IF(nspmd > 1)THEN
380C Sizing memory space for parith/on and spmd
381 CALL spmd_i7curvsz(
382 1 ipari(4,nin),intbuf_tab%IRECTM,numnod,iad_elem,fr_elem,
383 2 isdsiz ,ircsiz ,itag ,lenr ,lent ,
384 3 maxcc )
385 ELSE
386C Sizing memory space for parith/on without spmd
387 CALL i7curvsz(
388 1 ipari(4,nin),intbuf_tab%IRECTM,numnod,itag,lent,maxcc)
389 END IF
390 CALL i7normnp(
391 1 ipari(4,nin),intbuf_tab%IRECTM,numnod ,x ,intbuf_tab%NODNORM_NORMAL,
392 2 ipari(6,nin),intbuf_tab%MSR,lent ,maxcc,isdsiz ,
393 3 ircsiz ,iad_elem ,fr_elem,itag )
394 END IF
395 END IF
396 CALL my_barrier()
397 nmnft=1+(jtask-1)*ipari(6,nin)/nthread
398 nmnlt=jtask*ipari(6,nin)/nthread
399 CALL i7norme(nmnft,nmnlt,intbuf_tab%NODNORM_NORMAL,intbuf_tab%MSR)
400 CALL my_barrier()
401 nrtmft=1+(jtask-1)*ipari(4,nin)/nthread
402 nrtmlt=jtask*ipari(4,nin)/nthread
403 CALL i7rcurv(nrtmft, nrtmlt ,x ,intbuf_tab%NODNORM_NORMAL ,intbuf_tab%IRECTM ,
404 . intbuf_tab%MODRCURV , nradm ,intbuf_tab%MODANGLM ,anglt )
405 CALL my_barrier()
406 END IF
407C----------------------------------------------------
408C
409 i_stok = intbuf_tab%I_STOK(1)
410 debut = 0
411 i_stok_loc = 0
412
413C static decoupage
414 IF(inacti /= 7 .OR. ibag == 0)THEN
415 nb_loc = i_stok / nthread
416 IF (jtask == nthread) THEN
417 i_stok_loc = i_stok-nb_loc*(nthread-1)
418 ELSE
419 i_stok_loc = nb_loc
420 ENDIF
421 debut = (jtask-1)*nb_loc
422 END IF
423 i_stok = 0
424
425 IF (impl_s == 1) THEN
426 num_imp = 0
427 visc =zero
428 viscf =zero
429 ENDIF
430C
431C updating ISTOK
432C
433 IF(itied == 0)THEN
434 IF (inacti == 5 .OR. inacti == 6 .OR. type18)THEN
435 DO i = debut+1, debut+i_stok_loc
436 IF(intbuf_tab%CAND_N(i)<0) THEN
437 i_stok = i_stok + 1
438 index2(i_stok) = i
439C inbuf == cand_n
440 intbuf_tab%CAND_N(i) = -intbuf_tab%CAND_N(i)
441 ELSE
442C Discount at 0 of Cand_P
443 intbuf_tab%CAND_P(i) = zero
444 ENDIF
445 ENDDO
446 ELSE
447 DO i = debut+1, debut+i_stok_loc
448 IF(intbuf_tab%CAND_N(i)<0) THEN
449 i_stok = i_stok + 1
450 index2(i_stok) = i
451C inbuf == cand_n
452 intbuf_tab%CAND_N(i) = -intbuf_tab%CAND_N(i)
453 ENDIF
454 ENDDO
455 ENDIF
456 ELSE ! IF(ITIED == 0)THEN
457 IF (inacti == 5 .OR. inacti == 6)THEN
458 DO i = debut+1, debut+i_stok_loc
459 IF(intbuf_tab%CAND_F(8*(i-1)+1)/=zero) THEN
460 i_stok = i_stok + 1
461 index2(i_stok) = i
462 intbuf_tab%CAND_N(i) = iabs(intbuf_tab%CAND_N(i))
463 ELSEIF(intbuf_tab%CAND_N(i)<0) THEN
464 i_stok = i_stok + 1
465 index2(i_stok) = i
466C inbuf == cand_n
467 intbuf_tab%CAND_N(i) = -intbuf_tab%CAND_N(i)
468 ELSE
469C Discount at 0 of Cand_P
470 intbuf_tab%CAND_P(i) = zero
471 intbuf_tab%CAND_F(8*(i-1)+1:8*(i-1)+8) = zero
472 ENDIF
473 ENDDO
474 ELSE
475 DO i = debut+1, debut+i_stok_loc
476 IF(intbuf_tab%CAND_F(8*(i-1)+1)/=zero) THEN
477 i_stok = i_stok + 1
478 index2(i_stok) = i
479 intbuf_tab%CAND_N(i) = iabs(intbuf_tab%CAND_N(i))
480 ELSEIF(intbuf_tab%CAND_N(i)<0) THEN
481 i_stok = i_stok + 1
482 index2(i_stok) = i
483C inbuf == cand_n
484 intbuf_tab%CAND_N(i) = -intbuf_tab%CAND_N(i)
485 ELSE
486 intbuf_tab%CAND_F(8*(i-1)+1:8*(i-1)+8) = zero
487 ENDIF
488 ENDDO
489 END IF
490 ENDIF
491
492c------------------------------------------------
493 IF (debug(3) >= 1) THEN
494 nb_jlt = nb_jlt + i_stok_loc
495 nb_stok_n = nb_stok_n + i_stok
496 ENDIF
497
498 sfsavparit = 0
499 DO i=1,nisub+1
500 IF(isensint(i)/=0) THEN
501 sfsavparit = sfsavparit + 1
502 ENDIF
503 ENDDO
504 IF (sfsavparit /= 0) THEN
505 ALLOCATE(fsavparit(nisub+1,11,i_stok),stat=ierror)
506 IF(ierror/=0) THEN
507 CALL ancmsg(msgid=19,anmode=aninfo,c1='(/INTER/TYPE7)')
508 CALL arret(2)
509 ENDIF
510 DO j=1,i_stok
511 DO i=1,11
512 DO h=1,nisub+1
513 fsavparit(h,i,j) = zero
514 ENDDO
515 ENDDO
516 ENDDO
517 ELSE
518 ALLOCATE(fsavparit(0,0,0),stat=ierror)
519 IF(ierror/=0) THEN
520 CALL ancmsg(msgid=19,anmode=aninfo,c1='(/INTER/TYPE7)')
521 CALL arret(2)
522 ENDIF
523 ENDIF
524
525 !--------------------------------------------------------
526 ! FORCE COMPUTATION :
527 ! inter 7 & inter 18
528 !--------------------------------------------------------
529C!$OMP ORDERED
530 DO nft = 0 , i_stok - 1 , nvsiz
531 jlt = min( nvsiz, i_stok - nft )
532 ! preparing retained candidates
533 CALL i7cdcor3(jlt,index2(nft+1),intbuf_tab%CAND_E,intbuf_tab%CAND_N,cand_e_n,cand_n_n)
534 ! cand_n & cand_e replaced with cand_n_n & cand_e_n
535 CALL i7cor3(
536 1 jlt ,x ,intbuf_tab%IRECTM,intbuf_tab%NSV,cand_e_n,
537 2 cand_n_n ,intbuf_tab%STFM,intbuf_tab%STFNS,x1 ,x2 ,
538 3 x3 ,x4 ,y1 ,y2 ,y3 ,
539 4 y4 ,z1 ,z2 ,z3 ,z4 ,
540 5 xi ,yi ,zi ,stif ,ix1 ,
541 6 ix2 ,ix3 ,ix4 ,nsvg ,igap ,
542 7 gap ,intbuf_tab%GAP_S,intbuf_tab%GAP_M,gapv ,
543 9 ms ,vxi ,vyi ,
544 a vzi ,msi ,nsn ,v ,kinet ,
545 b kini ,nty ,nin ,igsti ,kmin ,
546 c kmax ,gapmax ,gapmin ,iadm ,intbuf_tab%MODRCURV ,
547 d rcurvi ,intbuf_tab%MODANGLM ,anglmi ,intth ,temp ,
548 e tempi ,phi ,intbuf_tab%AREAS,intbuf_tab%IELEC,areasi ,
549 f ieleci ,nodnx_sms ,nsms ,intbuf_tab%GAP_SL,intbuf_tab%GAP_ML,
550 g intfric ,intbuf_tab%IPARTFRICS,ipartfricsi,intbuf_tab%IPARTFRICM,ipartfricmi,
551 h iorthfric,intbuf_tab%IREP_FRICM,intbuf_tab%DIR_FRICM,irep_fricmi,dir_fricmi )
552 ! -- PIT
553
554 jlt_new = 0
555 jlt_tied= 0
556
557 IF(type18)THEN
558 CALL i18dst3(
559 1 jlt ,cand_n_n ,cand_e_n , cn_loc ,ce_loc ,
560 2 x1 ,x2 ,x3 , x4 ,y1 ,
561 3 y2 ,y3 ,y4 , z1 ,z2 ,
562 4 z3 ,z4 ,xi , yi ,zi ,
563 5 nx1 ,nx2 ,nx3 , nx4 ,ny1 ,
564 6 ny2 ,ny3 ,ny4 , nz1 ,nz2 ,
565 7 nz3 ,nz4 ,lb1 , lb2 ,lb3 ,
566 8 lb4 ,lc1 ,lc2 , lc3 ,lc4 ,
567 9 p1 ,p2 ,p3 , p4 ,ix1 ,
568 a ix2 ,ix3 ,ix4 , nsvg ,stif ,
569 b jlt_new ,gapv ,intbuf_tab%CAND_P ,ale_connectivity%NE_CONNECT,
570 c index2(nft+1) ,vxi ,vyi ,itab ,xcell ,
571 d vzi ,msi ,kini ,
572 e igap ,multi_fvm ,s_xcell_remote ,xcell_remote)
573
574 IF( multi_fvm%IS_INT18_LAW151 ) THEN
575 msi_append(1:mvsiz) = zero
576 msi_append(mvsiz+1:2*mvsiz) = msi(1:mvsiz)
577 ELSE
578 msi_append(1:mvsiz) = msi(1:mvsiz)
579 msi_append(mvsiz+1:2*mvsiz) = zero
580 ENDIF
581
582 jlt = jlt_new
583 IF (imonm > 0 .AND. jtask == 1) CALL startime(timers,20)
584 IF(jlt_new /= 0) THEN
585 ipari(29,nin) = 1
586 idir = ipari(34,nin)
587 IF (debug(3) >= 1)nb_jlt_new = nb_jlt_new + jlt_new
588 CALL i18for3(output,
589 1 jlt ,a ,v ,ibc ,icodt ,
590 2 fsav ,gap ,fric ,ms ,visc ,
591 3 viscf ,noint ,intbuf_tab%STFNS ,itab ,cn_loc ,
592 4 stiglo ,stifn ,stif ,fskyi ,isky ,
593 5 nx1 ,nx2 ,nx3 ,nx4 ,ny1 ,
594 6 ny2 ,ny3 ,ny4 ,nz1 ,nz2 ,
595 7 nz3 ,nz4 ,lb1 ,lb2 ,lb3 ,
596 8 lb4 ,lc1 ,lc2 ,lc3 ,lc4 ,
597 9 p1 ,p2 ,p3 ,p4 ,fcont ,
598 b ix1 ,ix2 ,ix3 ,ix4 ,nsvg ,
599 c ivis2 ,neltst ,ityptst ,dt2t ,ixs ,
600 d gapv ,intbuf_tab%CAND_P ,index2(nft+1) ,niskyfi ,
601 e kinet ,newfront ,isecin ,nstrf ,secfcum ,
602 f x ,intbuf_tab%IRECTM ,ce_loc ,mfrot ,ifq ,
603 g intbuf_tab%FRIC_P ,intbuf_tab%FTSAVX ,intbuf_tab%FTSAVY ,intbuf_tab%FTSAVZ ,
604 h intbuf_tab%IFPEN , icontact ,igroups ,iparg ,
605 j viscn ,vxi ,vyi ,vzi ,msi_append ,
606 k kini ,nin ,nisub ,intbuf_tab%LISUB ,intbuf_tab%ADDSUBS,
607 l intbuf_tab%ADDSUBM,intbuf_tab%LISUBS ,intbuf_tab%LISUBM ,fsavsub ,intbuf_tab%CAND_N ,
608 m ipari(33,nin) ,ipari(39,nin) ,fncont ,ms0 ,
609 n jtask ,isensint ,fsavparit ,nft ,multi_fvm ,
610 o h3d_data ,elbuf_tab ,idir )
611 ENDIF
612
613 ELSE
614 CALL i7dst3(
615 1 jlt ,cand_n_n,cand_e_n,cn_loc ,ce_loc ,
616 2 x1 ,x2 ,x3 ,x4 ,y1 ,
617 3 y2 ,y3 ,y4 ,z1 ,z2 ,
618 4 z3 ,z4 ,xi ,yi ,zi ,
619 5 nx1 ,nx2 ,nx3 ,nx4 ,ny1 ,
620 6 ny2 ,ny3 ,ny4 ,nz1 ,nz2 ,
621 7 nz3 ,nz4 ,lb1 ,lb2 ,lb3 ,
622 8 lb4 ,lc1 ,lc2 ,lc3 ,lc4 ,
623 9 p1 ,p2 ,p3 ,p4 ,ix1 ,
624 a ix2 ,ix3 ,ix4 ,nsvg ,stif ,
625 b jlt_new ,gapv ,inacti ,intbuf_tab%CAND_P,
626 c index2(nft+1),vxi ,vyi ,
627 d vzi ,msi ,kini ,icurv ,intbuf_tab%IRECTM,
628 e nnx1 ,nnx2 ,nnx3 ,nnx4 ,nny1 ,
629 f nny2 ,nny3 ,nny4 ,nnz1 ,nnz2 ,
630 g nnz3 ,nnz4 ,intbuf_tab%NODNORM_NORMAL ,iadm,rcurvi ,
631 h anglmi ,intth ,tempi ,phi ,areasi ,
632 i ieleci ,nsms ,cmaj ,drad2 ,
633 j intfric ,ipartfricsi,
634 k ipartfricmi ,itied ,jlt_tied,intbuf_tab%CAND_F,iorthfric,
635 l irep_fricmi ,dir_fricmi,dgaploadp)
636
637 jlt = jlt_new
638 IF (imonm > 0 .AND. jtask == 1) CALL startime(timers,20)
639 IF(jlt_new /= 0) THEN
640 ipari(29,nin) = 1
641 IF (debug(3)>=1) nb_jlt_new = nb_jlt_new + jlt_new
642
643 !-------------------------------------------------------------------------------
644 ! Friction model : computation of friction coefficients based on Material of connected Parts
645 !-------------------------------------------------------------------------------
646 IF(itask==1) CALL startime(timers,macro_timer_fric)
647 IF(iorthfric > 0) THEN
648 CALL frictionparts_model_ortho(
649 1 intfric ,jlt ,ipartfricsi ,ipartfricmi ,adparts_fric ,
650 2 nsetprts ,tabcoupleparts_fric,npartfric,tabparts_fric,tabcoef_fric ,
651 3 fric ,viscf ,intbuf_tab%FRIC_P,fric_coefs , fricc ,
652 4 viscffric ,nty ,mfrot ,iorthfric , fric_coefs2,
653 5 fricc2 ,viscffric2 ,ifricorth ,nforth , nfisot ,
654 6 indexorth ,indexisot ,jlt_tied ,irep_fricmi , dir_fricmi,
655 7 ix3 ,ix4 ,x1 ,y1 , z1 ,
656 8 x2 ,y2 ,z2 ,x3 , y3 ,
657 9 z3 ,x4 ,y4 ,z4 ,ce_loc ,
658 a dir1 ,dir2 )
659 ELSE
660 nforth = 0
661 nfisot = 0
662 CALL frictionparts_model_isot(
663 1 intfric ,jlt ,ipartfricsi ,ipartfricmi ,adparts_fric ,
664 2 nsetprts ,tabcoupleparts_fric,npartfric,tabparts_fric,tabcoef_fric ,
665 3 fric ,viscf ,intbuf_tab%FRIC_P,fric_coefs , fricc ,
666 4 viscffric ,nty ,mfrot ,iorthfric , ifric ,
667 5 jlt_tied , tint ,tempi ,npc ,tf ,
668 6 temp , h1 ,h2 ,h3 ,h4 ,
669 7 ix1 , ix2 ,ix3 ,ix4 ,iform )
670 ENDIF
671 IF(itask==1) CALL stoptime(timers,macro_timer_fric)
672
673 CALL i7for3(output,
674 1 jlt ,a ,v ,ibc ,icodt ,
675 2 fsav ,gap ,fric ,ms ,visc ,
676 3 viscf ,noint ,intbuf_tab%STFNS,itab ,cn_loc ,
677 4 stiglo ,stifn ,stif ,fskyi ,isky ,
678 5 nx1 ,nx2 ,nx3 ,nx4 ,ny1 ,
679 6 ny2 ,ny3 ,ny4 ,nz1 ,nz2 ,
680 7 nz3 ,nz4 ,lb1 ,lb2 ,lb3 ,
681 8 lb4 ,lc1 ,lc2 ,lc3 ,lc4 ,
682 9 p1 ,p2 ,p3 ,p4 ,fcont ,
683 a ix1 ,ix2 ,ix3 ,ix4 ,nsvg ,
684 b ivis2 ,neltst ,ityptst ,dt2t ,gapv ,
685 c inacti,intbuf_tab%CAND_P,index2(nft+1),kinet ,newfront ,
686 d isecin ,nstrf ,x ,intbuf_tab%IRECTM,ce_loc ,
687 e mfrot ,ifq ,intbuf_tab%FTSAVX ,intbuf_tab%FTSAVY,
688 f intbuf_tab%FTSAVZ,xfiltr_fric,intbuf_tab%IFPEN ,ibag ,
689 h icontact ,viscn ,vxi ,vyi ,vzi ,
690 i msi ,kini ,nin ,nisub ,intbuf_tab%LISUB,
691 j intbuf_tab%ADDSUBS,intbuf_tab%ADDSUBM,intbuf_tab%LISUBS ,
692 . intbuf_tab%LISUBM,fsavsub,
693 k intbuf_tab%CAND_N,ipari(33,nin),ipari(39,nin),intbuf_tab%NODNORM_NORMAL ,fncont,
694 l ftcont ,x1 ,x2 ,x3 ,x4 ,
695 m y1 , y2 ,y3 ,y4 ,z1 ,
696 n z2 ,z3 ,z4 ,xi ,yi ,
697 o zi ,iadm ,rcurvi ,rcontact ,acontact ,
698 p pcontact ,anglmi ,padm ,intth ,temp ,
699 q tempi ,iform ,npc ,tf ,cmaj ,
700 r dtmini ,drad ,fheats ,fheatm ,efrict ,
701 s qfricint(nin),fni ,ifric ,jtask ,h1 ,
702 u h2 ,h3 ,h4 ,ks ,kt ,
703 v k1 ,k2 ,k3 ,k4 ,c1 ,
704 w c2 ,c3 ,c4 ,cs ,c ,
705 x cf ,tint ,xfric ,fxi ,fyi ,
706 y fzi ,fx1 ,fy1 ,fz1 ,fx2 ,
707 z fy2 ,fz2 ,fx3 ,fy3 ,fz3 ,
708 1 fx4 ,fy4 ,fz4 ,isensint ,fsavparit ,
709 5 nft ,sym_flag_type19,h3d_data,fricc ,viscffric ,
710 6 fric_coefs ,itied ,jlt_tied ,intbuf_tab%CAND_F,iorthfric,
711 7 fric_coefs2 ,fricc2 ,viscffric2 ,nforth ,nfisot ,
712 8 indexorth ,indexisot ,dir1 ,dir2 ,tagncont ,
713 9 kloadpinter ,loadpinter ,loadp_hyd_inter,intbuf_tab%TYPSUB ,
714 . intbuf_tab%INFLG_SUBS,
715 a intbuf_tab%INFLG_SUBM,ninloadp,dgaploadint,s_loadpinter,dgaploadp,
716 b interefric )
717
718 IF( intth > 0 ) THEN
719 CALL i7therm(jlt ,iparg ,pm ,ixs ,iform ,x ,
720 1 xi ,yi ,zi ,x1 ,y1 ,z1 ,
721 2 x2 ,y2 ,z2 ,x3 ,y3 ,z3 ,
722 3 x4 ,y4 ,z4 ,ix1 ,ix2 ,ix3 ,
723 4 ix4 ,rstif ,tempi, intbuf_tab%IELEC ,
724 5 phi ,tint , areasi,ieleci ,frad,drad ,
725 6 gapv ,fni ,irsth,xthe ,npc ,tf ,
726 6 condint,phi1,phi2 ,phi3 ,phi4 ,fheats,
727 7 fheatm,efrict,temp ,h1 ,h2 ,h3 ,
728 8 h4 ,theaccfact)
729 ENDIF
730
731 CALL i7ass33(
732 1 jlt ,a ,noint ,itab ,stifn ,
733 2 stif ,fskyi ,isky ,fcont ,ix1 ,
734 3 ix2 ,ix3 ,ix4 ,nsvg ,neltst ,
735 4 ityptst ,dt2t ,niskyfi ,isecin ,nstrf ,
736 5 secfcum ,viscn ,nin ,fxi ,fyi ,
737 6 fzi ,fx1 ,fy1 ,fz1 ,fx2 ,
738 7 fy2 ,fz2 ,fx3 ,fy3 ,fz3 ,
739 8 fx4 ,fy4 ,fz4 ,h1 ,h2 ,
740 9 h3 ,h4 ,ks ,kt ,k1 ,
741 a k2 ,k3 ,k4 ,cs ,cf ,
742 b c1 ,c2 ,c3 ,c4 ,c ,
743 c intth ,phi ,phi1 ,phi2 ,phi3 ,
744 d phi4 ,fthe ,ftheskyi ,mskyi_sms ,iskyi_sms ,
745 e nsms ,dtmini ,jtask ,
746 f condn ,condint ,condnskyi ,ixig3d ,kxig3d ,
747 j wige ,knot ,igeo ,intbuf_tab%NIGE ,
748 k intbuf_tab%RIGE ,x ,h3d_data ,knotlocpc ,
749 l knotlocel ,iform ,nodadt_therm)
750
751 ENDIF
752 ENDIF
753
754C-----------------------------------------------
755C
756 IF (imonm > 0 .AND. jtask == 1) CALL stoptime(timers,20)
757 IF(impl_s = =1) THEN
758 DO i = 1 ,jlt_new
759 ns_imp(i+num_imp)=cn_loc(i)
760 ne_imp(i+num_imp)=ce_loc(i)
761 ind_imp(i+num_imp)=index2(i+nft)
762 ENDDO
763 num_imp=num_imp+jlt_new
764 ENDIF
765 ENDDO
766C!$OMP END ORDERED
767 IF (sfsavparit /= 0)THEN
768 CALL sum_6_float_sens(fsavparit, nisub+1, 11, i_stok,1,i_stok, fbsav6, 12, 6, dimfb, isensint )
769 ENDIF
770 IF (ALLOCATED(fsavparit)) DEALLOCATE (fsavparit)
771
772 IF(icurv == 3 .OR. iadm /= 0)THEN
773 CALL my_barrier
774 IF(jtask==1)DEALLOCATE(intbuf_tab%NODNORM_NORMAL)
775 END IF
776 IF(iadm/=0)THEN
777 CALL my_barrier
778 IF(jtask==1)DEALLOCATE(intbuf_tab%MODRCURV,intbuf_tab%MODANGLM)
779 END IF
780 ELSEIF(nty==17)THEN
781 ilagm = ipari(33,nin)
782C I7KGLO = 1 => mis ds inttri
783 IF(ilagm == 0)THEN
784 CALL i17main_pena(output,
785 1 nin ,ipari ,intbuf_tab ,x ,v ,
786 2 a ,jtask-1,igrbric ,eminx(iad17),ms ,
787 3 ixs ,ixs16 ,ixs20 ,stifn ,fskyi ,isky ,
788 4 fsav ,fcont ,niskyfi ,h3d_data)
789 ENDIF
790 ENDIF
791C
792 RETURN
793 END
794
frictionparts_model_ortho
subroutine frictionparts_model_ortho(intfric, jlt, ipartfricsi, ipartfricmi, adparts_fric, nset, tabcoupleparts_fric, npartfric, tabparts_fric, tabcoef_fric, fric, viscf, frot_p, fric_coefs, fricc, viscffric, nty, mfrot, iorthfric, fric_coefs2, fricc2, viscffric2, ifricorth, nforth, nfisot, indexorth, indexisot, jlt_tied, irep_fricmi, dir_fricmi, ix3, ix4, x1, y1, z1, x2, y2, z2, x3, y3, z3, x4, y4, z4, ce_loc, dir1, dir2)
Definition frictionparts_model.F:43
frictionparts_model_isot
subroutine frictionparts_model_isot(intfric, jlt, ipartfricsi, ipartfricmi, adparts_fric, nset, tabcoupleparts_fric, npartfric, tabparts_fric, tabcoef_fric, fric, viscf, frot_p, fric_coefs, fricc, viscffric, nty, mfrot, iorthfric, ifric, jlt_tied, tint, tempi, npc, tf, temp, h1, h2, h3, h4, ix1, ix2, ix3, ix4, iform)
Definition frictionparts_model.F:350
i17main_pena
subroutine i17main_pena(output, nin, ipari, intbuf_tab, x, v, a, itask, igrbric, eminx, ms, ixs, ixs16, ixs20, stifn, fskyi, isky, fsav, fcont, niskyfi, h3d_data)
Definition i17main_pena.F:381
i7ass33
subroutine i7ass33(jlt, a, noint, itab, stifn, stif, fskyi, isky, fcont, ix1, ix2, ix3, ix4, nsvg, neltst, ityptst, dt2t, niskyfi, isecin, nstrf, secfcum, viscn, nin, fxi, fyi, fzi, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, h1, h2, h3, h4, ks, kt, k1, k2, k3, k4, cs, cf, c1, c2, c3, c4, c, intth, phi, phi1, phi2, phi3, phi4, fthe, ftheskyi, mskyi_sms, iskyi_sms, nsms, dtmini, jtask, condn, condint, condnskyi, ixig3d, kxig3d, wige, knot, igeo, nige, rige, x, h3d_data, knotlocpc, knotlocel, iform, nodadt_therm)
Definition i7ass3.F:58
i7cdcor3
subroutine i7cdcor3(jlt, index, cand_e, cand_n, cand_e_n, cand_n_n)
Definition i7cdcor3.F:38
i7curvsz
subroutine i7curvsz(nrtm, irect, numnod, itag, lent, maxcc)
Definition i7curv.F:219
i7normp
subroutine i7normp(nrtm, irect, numnod, x, nod_normal, nmn, msr, lent, maxcc, isdsiz, ircsiz, iad_elem, fr_elem, itag)
Definition i7curv.F:102
i7norm
subroutine i7norm(nrtm, irect, numnod, x, nod_normal, nmn, msr)
Definition i7curv.F:30
i7for3
subroutine i7for3(output, jlt, a, v, ibcc, icodt, fsav, gap, fric, ms, visc, viscf, noint, stfn, itab, cn_loc, stiglo, stifn, stif, fskyi, isky, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, p1, p2, p3, p4, fcont, ix1, ix2, ix3, ix4, nsvg, ivis2, neltst, ityptst, dt2t, gapv, inacti, cand_p, index, kinet, newfront, isecin, nstrf, x, irect, ce_loc, mfrot, ifq, cand_fx, cand_fy, cand_fz, alpha0, ifpen, ibag, icontact, viscn, vxi, vyi, vzi, msi, kini, nin, nisub, lisub, addsubs, addsubm, lisubs, lisubm, fsavsub, cand_n, ilagm, icurv, nod_normal, fncont, ftcont, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, iadm, rcurvi, rcontact, acontact, pcontact, anglmi, padm, intth, temp, tempi, iform, npc, tf, cmaj, dtmini, drad, fheats, fheatm, efrict, qfric, fni, ifric, jtask, h1, h2, h3, h4, ks, kt, k1, k2, k3, k4, c1, c2, c3, c4, cs, c, cf, tint, xfric, fxi, fyi, fzi, fx1, fy1, fz1, fx2, fy2, fz2, fx3, fy3, fz3, fx4, fy4, fz4, isensint, fsavparit, nft, sym_flag_type19, h3d_data, fricc, viscffric, fric_coefs, itied, jlt_tied, cand_f, iorthfric, fric_coefs2, fricc2, viscffric2, nforth, nfisot, indexorth, indexisot, dir1, dir2, tagncont, kloadpinter, loadpinter, loadp_hyd_inter, typsub, inflg_subs, inflg_subm, ninloadp, dgaploadint, s_loadpinter, dgaploadp, interefric)
Definition i7for3.F:78
i7mainf
subroutine i7mainf(output, timers, ipari, x, a, ale_connectivity, xcell, icodt, fsav, v, ms, dt2t, neltst, ityptst, itab, stifn, fskyi, isky, fcont, nin, lindmax, kinet, jtask, nb_jlt, nb_jlt_new, nb_stok_n, elbuf_tab, niskyfi, newfront, nstrf, secfcum, igroups, icontact, viscn, num_imp, ns_imp, ne_imp, ind_imp, fsavsub, nrtmdim, igrbric, eminx, ixs, ixs16, ixs20, fncont, ftcont, iad_elem, fr_elem, rcontact, acontact, pcontact, temp, fthe, ftheskyi, pm, iparg, iad17, mskyi_sms, iskyi_sms, nodnx_sms, ms0, qfricint, npc, tf, condn, condnskyi, intbuf_tab, nodadt_therm, theaccfact, fbsav6, isensint, dimfb, ixig3d, kxig3d, wige, knot, igeo, multi_fvm, h3d_data, intbuf_fric_tab, knotlocpc, knotlocel, itask, tagncont, kloadpinter, loadpinter, loadp_hyd_inter, dgaploadint, s_loadpinter, interefric, s_xcell_remote, xcell_remote)
Definition i7mainf.F:89
i7rcurv
subroutine i7rcurv(nrtmft, nrtmlt, x, nod_normal, irect, rcurv, nradm, anglm, anglt)
Definition i7rcurv.F:274
i7normnp
subroutine i7normnp(nrtm, irect, numnod, x, nod_normal, nmn, msr, lent, maxcc, isdsiz, ircsiz, iad_elem, fr_elem, itag)
Definition i7rcurv.F:107
i7normn
subroutine i7normn(nrtm, irect, numnod, x, nod_normal, nmn, msr)
Definition i7rcurv.F:30
i7norme
subroutine i7norme(nmnft, nmnlt, nod_normal, msr)
Definition i7rcurv.F:229
i7therm
subroutine i7therm(jlt, iparg, pm, ixs, iform, x, xi, yi, zi, x1, y1, z1, x2, y2, z2, x3, y3, z3, x4, y4, z4, ix1, ix2, ix3, ix4, rstif, tempi, ieles, phi, tint, areas, ieleci, frad, drad, gapv, fni, ifunctk, xthe, npc, tf, condint, phi1, phi2, phi3, phi4, fheats, fheatm, efrict, temp, h1, h2, h3, h4, theaccfact)
Definition i7therm.F:44
min
#define min(a, b)
Definition macros.h:20
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
ale_connectivity_mod
Definition ale_connectivity_mod.F:225
groupdef_mod
Definition groupdef_mod.F:668
h3d_mod
Definition h3d_mod.F:311
i18dst3_mod
Definition i18dst3.F:29
i18dst3_mod::i18dst3
subroutine i18dst3(jlt, cand_n, cand_e, cn_loc, ce_loc, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, p1, p2, p3, p4, ix1, ix2, ix3, ix4, nsvg, stif, jlt_new, gapv, cand_p, ale_ne_connect, index, vxi, vyi, itab, xcell, vzi, msi, kini, igap, multi_fvm, s_xcell_remote, xcell_remote)
Definition i18dst3.F:56
i18for3_mod
Definition i18for3.F:28
i18for3_mod::i18for3
subroutine i18for3(output, jlt, a, v, ibcc, icodt, fsav, gap, fric, ms, visc, viscf, noint, stfn, itab, cn_loc, stfval, stifn, stif, fskyi, isky, nx1, nx2, nx3, nx4, ny1, ny2, ny3, ny4, nz1, nz2, nz3, nz4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, p1, p2, p3, p4, fcont, ix1, ix2, ix3, ix4, nsvg, ivis2, neltst, ityptst, dt2t, ixs, gapv, cand_p, index, niskyfi, kinet, newfront, isecin, nstrf, secfcum, x, irect, ce_loc, mfrot, ifq, frot_p, cand_fx, cand_fy, cand_fz, ifpen, icontact, igroups, iparg, viscn, vxi, vyi, vzi, msi, kini, nin, nisub, lisub, addsubs, addsubm, lisubs, lisubm, fsavsub, cand_n, ilagm, icurv, fncont, ms0, jtask, isensint, fsavparit, nft, multi_fvm, h3d_data, elbuf_tab, idir)
Definition i18for3.F:70
message_mod
Definition message_mod.F:1257
sum_6_float_sens
subroutine sum_6_float_sens(f, a, b, c, jft, jlt, f6, d, e, g, isensint)
Definition parit.F:540
spmd_exch_n
subroutine spmd_exch_n(xnorm, iad_elem, fr_elem, lenr)
Definition spmd_exch_n.F:37
spmd_i7curvsz
subroutine spmd_i7curvsz(nrtm, irect, numnod, iad_elem, fr_elem, isdsiz, ircsiz, itag, lenr, lent, maxcc)
Definition spmd_i7curvsz.F:34
i7cor3
subroutine i7cor3(x, irect, nsv, cand_e, cand_n, stf, stfn, gapv, igap, gap, gap_s, gap_m, istf, gapmin, gapmax, gap_s_l, gap_m_l, drad, ix1, ix2, ix3, ix4, nsvg, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, dgapload, last)
Definition i7cor3.F:43
i7dst3
subroutine i7dst3(ix3, ix4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, x0, y0, z0, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, last)
Definition i7dst3.F:46
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:895
arret
subroutine arret(nn)
Definition arret.F:86
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
h3d_mod::h3d_database
Definition h3d_mod.F:438