OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
imp_sol_init.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!|| imp_sol_init ../engine/source/implicit/imp_sol_init.F
25!||--- called by ------------------------------------------------------
26!|| resol ../engine/source/engine/resol.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| buf_dim1 ../engine/source/implicit/produt_v.F
31!|| cp_dm ../engine/source/implicit/produt_v.F
32!|| cp_real ../engine/source/implicit/produt_v.F
33!|| dim_glob_k ../engine/source/implicit/ind_glob_k.F
34!|| dim_int7 ../engine/source/implicit/ind_glob_k.F
35!|| dim_span ../engine/source/implicit/ind_glob_k.F
36!|| fil_span0 ../engine/source/implicit/ind_glob_k.F
37!|| imp_frii ../engine/source/mpi/implicit/imp_fri.F
38!|| imp_trans0 ../engine/source/output/restart/wrrest.F
39!|| ind_glob_k ../engine/source/implicit/ind_glob_k.F
40!|| spmd_mstop ../engine/source/mpi/init/spmd_mstop.F
41!|| zero1 ../engine/source/system/zero.F
42!||--- uses -----------------------------------------------------
43!|| elbufdef_mod ../common_source/modules/mat_elem/elbufdef_mod.F90
44!|| element_mod ../common_source/modules/elements/element_mod.F90
45!|| groupdef_mod ../common_source/modules/groupdef_mod.F
46!|| impbufdef_mod ../engine/share/modules/impbufdef_mod.F
47!|| intbufdef_mod ../common_source/modules/interfaces/intbufdef_mod.F90
48!|| message_mod ../engine/share/message_module/message_mod.F
49!||====================================================================
50 SUBROUTINE imp_sol_init(
51 1 GEO ,NPBY ,LPBY ,ITAB ,
52 2 IPARI ,IXS ,IXQ ,IXC ,IXT ,
53 4 IXP ,IXR ,IXTG ,IXTG1 ,IXS10 ,
54 5 IXS20 ,IXS16 ,IPARG ,
55 6 ELBUF ,NINT7 ,NBINTC ,X ,DMCP ,
56 7 FR_ELEM ,IAD_ELEM ,FR_I2M ,IAD_I2M ,
57 8 NPRW ,NUM_IMP1 ,NUM_IMPL ,MONVOL ,IGRSURF ,
58 9 FR_MV ,IPM ,IGEO ,IAD_RBY ,
59 A FR_RBY ,SH4TREE ,SH3TREE ,IRBE3 ,LRBE3 ,
60 B FR_RBE3M ,IAD_RBE3M ,IRBE2 ,LRBE2 ,IBFV ,
61 C VEL ,ELBUF_TAB ,IFRAME ,INTBUF_TAB,
62 D NDDL0 ,NNZK0 ,IMPBUF_TAB )
63C-----------------------------------------------
64C M o d u l e s
65C-----------------------------------------------
66 USE elbufdef_mod
67 USE intbufdef_mod
68 USE groupdef_mod
69 USE message_mod
70 USE impbufdef_mod
71 use element_mod , only : nixs,nixq,nixc,nixt,nixp,nixr,nixtg
72C-----------------------------------------------
73C I m p l i c i t T y p e s
74C-----------------------------------------------
75#include "implicit_f.inc"
76C-----------------------------------------------
77C C o m m o n B l o c k s
78C-----------------------------------------------
79#include "com01_c.inc"
80#include "com04_c.inc"
81#include "com08_c.inc"
82#include "param_c.inc"
83#include "scr05_c.inc"
84#include "impl1_c.inc"
85#include "task_c.inc"
86#include "buckcom.inc"
87#include "units_c.inc"
88C-----------------------------------------------------------------
89C D u m m y A r g u m e n t s
90C-----------------------------------------------
91 INTEGER NDDL0,NNZK0,IPARG(NPARG,*),FR_ELEM(*) ,IAD_ELEM(2,*)
92 INTEGER NPBY(NNPBY,*),LPBY(*),ITAB(*)
93 INTEGER NINT7,NBINTC,IPARI(NPARI,*),
94 . FR_I2M(*),IAD_I2M(*),FR_RBY(*),IAD_RBY(*)
95 INTEGER MONVOL(*),
96 . FR_MV(NSPMD+2,NVOLU),NPRW(*),FR_RBE3M(*),IAD_RBE3M(*)
97 INTEGER IPM(NPROPMI,*),IGEO(NPROPGI,*),IFRAME(LISKN,*)
98 INTEGER
99 . IXS(NIXS,*),IXQ(NIXQ,*),IXC(NIXC,*), IXT(NIXT,*), IXP(NIXP,*),
100 . IXR(NIXR,*), IXTG(NIXTG,*),IXS10(6,*),IXS20(12,*),
101 . IXS16(8,*),IXTG1(4,*),IRBE3(*),LRBE3(*),
102 . SH4TREE(*), SH3TREE(*),
103 . irbe2(*),lrbe2(*),ibfv(*),num_imp1(*),num_impl(ninter,nthread)
104C REAL
105 my_real
106 . geo(npropg,*),elbuf(*),vel(*),x(*),dmcp(*)
107 TYPE (ELBUF_STRUCT_), DIMENSION(NGROUP) :: ELBUF_TAB
108 TYPE (INTBUF_STRUCT_) INTBUF_TAB(*)
109 TYPE (SURF_) , DIMENSION(NSURF) :: IGRSURF
110 TYPE (IMPBUF_STRUCT_) ,TARGET :: IMPBUF_TAB
111C-----------------------------------------------
112C L o c a l V a r i a b l e s
113C-----------------------------------------------
114 INTEGER I,J,K,N,M,L,NDOFI,NDOFJ,NKINE,NMIJ2,IP,NPN,NPP,IER1,
115 . LI1,LI2,LI3,LI4,LI5,LI6,LI7,LI8,LI9,LI10,LI11,LI12,
116 . LIF,LI13,LI14,LI15,LI16,LI17,IER2
117 INTEGER NTMP,L1,NNDLNKINE,NNMAX,NKMAX,NNDL
118 INTEGER, POINTER :: NDDL,NNZK,NRBYAC,NINT2,NMC,NMC2,NMONV
119 INTEGER, DIMENSION(:) ,POINTER :: IADK,JDIK,IDDL,NDOF,INLOC,LSIZE,I_IMP,
120 . irbyac,nsc,iint2,nkud,imonv,ikinw
121 my_real, DIMENSION(:) ,POINTER :: diag_k,lt_k,diag_m,lt_m,lb,
122 . lb0,bkud,d_imp,elbuf_c,bufmat_c,
123 . x_c,dd,ddr
124C=======================================================================
125 IF (n2d>0) THEN
126 IF(ispmd==0)THEN
127 CALL ancmsg(msgid=161,anmode=aninfo)
128 ENDIF
129C CALL ARRET(2)
130 ENDIF
131 IF (iresp==1) THEN
132 IF(ispmd==0)THEN
133 CALL ancmsg(msgid=162,anmode=aninfo)
134 ENDIF
135 CALL arret(2)
136 ENDIF
137#ifndef MPI
138 IF (neig>0) THEN
139 IF(ispmd==0)THEN
140 CALL ancmsg(msgid=294,anmode=aninfo)
141 ENDIF
142 CALL arret(2)
143 END IF
144 IF (nbuck>0) THEN
145 IF(ispmd==0)THEN
146 CALL ancmsg(msgid=295,anmode=aninfo)
147 ENDIF
148 CALL arret(2)
149 END IF
150#endif
151c
152 ALLOCATE(impbuf_tab%IDDL(numnod))
153 ALLOCATE(impbuf_tab%NDOF(numnod))
154 ALLOCATE(impbuf_tab%INLOC(numnod))
155 ALLOCATE(impbuf_tab%IRBYAC(2*nrbykin))
156 ALLOCATE(impbuf_tab%NSC(nrbykin))
157 ALLOCATE(impbuf_tab%IINT2(ninter))
158 ALLOCATE(impbuf_tab%NKUD(nfxvel))
159 ALLOCATE(impbuf_tab%IMONV(nvolu))
160 nddl => impbuf_tab%NDDL
161 nnzk => impbuf_tab%NNZK
162 nrbyac => impbuf_tab%NRBYAC
163 nint2 => impbuf_tab%NINT2
164 nmc => impbuf_tab%NMC
165 nmc2 => impbuf_tab%NMC2
166 nmonv => impbuf_tab%NMONV
167 iddl => impbuf_tab%IDDL
168 ndof => impbuf_tab%NDOF
169 inloc => impbuf_tab%INLOC
170 lsize => impbuf_tab%LSIZE
171 i_imp => impbuf_tab%I_IMP
172 irbyac => impbuf_tab%IRBYAC
173 nsc => impbuf_tab%NSC
174 nsc = 0
175 iint2 => impbuf_tab%IINT2
176 nkud => impbuf_tab%NKUD
177 imonv => impbuf_tab%IMONV
178 CALL dim_glob_k(
179 1 geo ,npby ,lpby ,itab ,nrbyac ,
180 2 irbyac ,nint2 ,iint2 ,ipari ,
181 3 ixs ,ixq ,ixc ,ixt ,
182 4 ixp ,ixr ,ixtg ,ixtg1 ,ixs10 ,
183 5 ixs20 ,ixs16 ,iparg ,ndof ,
184 6 iddl ,nddl ,nnzk ,elbuf ,inloc ,
185 7 lsize ,fr_elem ,iad_elem ,fr_i2m ,iad_i2m ,
186 8 nprw ,nmonv ,imonv ,monvol ,igrsurf ,
187 9 fr_mv ,ipm ,igeo ,iad_rby ,
188 a fr_rby ,sh4tree ,sh3tree ,irbe3 ,lrbe3 ,
189 b fr_rbe3m ,iad_rbe3m ,irbe2 ,lrbe2 ,ibfv ,
190 c vel ,elbuf_tab ,iframe ,intbuf_tab )
191 s_iadk=nddl+1
192 s_jdik=nnzk
193 nddl0 = nddl
194 nnzk0 = nnzk
195 ALLOCATE(impbuf_tab%IADK(s_iadk))
196 ALLOCATE(impbuf_tab%JDIK(s_jdik))
197 iadk => impbuf_tab%IADK
198 jdik => impbuf_tab%JDIK
199c S_ISIJ=IMPBUF_TAB%LSIZE(4)
200c S_IMIJ=IMPBUF_TAB%LSIZE(5)
201c S_NSS=IMPBUF_TAB%LSIZE(1)
202c S_ISS=IMPBUF_TAB%LSIZE(3)
203c S_ISIJ2=IMPBUF_TAB%LSIZE(7)
204c S_NSS2=IMPBUF_TAB%LSIZE(2)
205c S_ISS2=IMPBUF_TAB%LSIZE(6)
206c S_NSC2=NINT2
207c S_NROWK=IMPBUF_TAB%LSIZE(8)
208c S_ICOK=(LSIZE(8)-LCOKM)*LSIZE(9)
209c S_ICOKM=LCOKM*LSIZE(10)
210c S_NMIJ2=4*LSIZE(11)
211c S_NSS3=IMPBUF_TAB%LSIZE(14)
212c S_ISB2=IMPBUF_TAB%LSIZE(15)
213c S_NSRB2=IMPBUF_TAB%LSIZE(16)
214C
215 nkine = lsize(8)
216 nnmax = lsize(9)
217 nkmax = lsize(10)
218 nmij2 = lsize(11)
219 npn = lsize(12)
220 npp = lsize(13)
221C
222c ALLOCATE(IMPBUF_TAB%ISIJ(S_ISIJ))
223c ALLOCATE(IMPBUF_TAB%IMIJ(S_IMIJ))
224c ALLOCATE(IMPBUF_TAB%NSS(S_NSS))
225c ALLOCATE(IMPBUF_TAB%ISIJ2(S_ISIJ2))
226c ALLOCATE(IMPBUF_TAB%NSS2(S_NSS2))
227c ALLOCATE(IMPBUF_TAB%ISS2(S_ISS2))
228c ALLOCATE(IMPBUF_TAB%NSC2(S_NSC2))
229c ALLOCATE(IMPBUF_TAB%ISS2(S_ISS2))
230c ALLOCATE(IMPBUF_TAB%NROWK(S_NROWK))
231c ALLOCATE(IMPBUF_TAB%ICOK(S_ICOK))
232c ALLOCATE(IMPBUF_TAB%ICOKM(S_ICOKM))
233c ALLOCATE(IMPBUF_TAB%NMIJ2(S_NMIJ2))
234c ALLOCATE(IMPBUF_TAB%NSS3(S_NSS3))
235c ALLOCATE(IMPBUF_TAB%ISB2(S_ISB2))
236c ALLOCATE(IMPBUF_TAB%NSRB2(S_NSRB2))
237C----finally not to splite IKINW, otherwise too many changes
238 li1 =1
239 li2 = li1+lsize(4)
240 li3 = li2+lsize(5)
241 li4 = li3+lsize(1)
242 li5 = li4+lsize(3)
243 li6 = li5+lsize(7)
244 li7 = li6+lsize(2)
245 li8 = li7+lsize(6)
246 li9 = li8+nint2
247 li10 = li9+lsize(8)
248 li11 = li10+(lsize(8)-lcokm)*lsize(9)
249 li12 = li11+lcokm*lsize(10)
250 li13 = li12+4*lsize(11)
251 li14 = li13+lsize(14)
252 li15 = li14+lsize(15)
253 lif = li15+lsize(16)
254 s_ikinw = lif
255 ALLOCATE(impbuf_tab%IKINW(s_ikinw))
256 ikinw => impbuf_tab%IKINW
257 CALL ind_glob_k(npby ,lpby ,
258 1 itab ,nrbyac ,irbyac ,nsc ,ikinw(li1),
259 2 nmc ,ikinw(li2),ikinw(li3),ikinw(li4),nint2 ,
260 3 iint2 ,ipari ,intbuf_tab,ikinw(li8),ikinw(li5),
261 4 ikinw(li6),ikinw(li7),iparg ,elbuf ,elbuf_tab ,
262 5 ixs ,ixq ,ixc ,ixt ,ixp ,
263 6 ixr ,ixtg ,ixtg1 ,ixs10 ,ixs20 ,
264 7 ixs16 ,iddl ,ndof ,iadk ,
265 8 jdik ,nddl ,nnzk ,lsize(9) ,lsize(8) ,
266 9 inloc ,lsize(10),ikinw(li9),ikinw(li10),ikinw(li11),
267 a lsize(11) ,ikinw(li12),li1 ,lsize(12) ,lsize(13) ,
268 b fr_elem ,iad_elem ,ipm ,igeo ,irbe3 ,
269 c lrbe3 ,ikinw(li13),fr_i2m ,iad_i2m ,fr_rbe3m ,
270 d iad_rbe3m ,irbe2 ,lrbe2 ,ikinw(li14),ikinw(li15))
271 ntmp=max(nkmax,nnmax)
272 IF (iroddl/=0) THEN
273 maxb = min(6*(ntmp+1),nddl)
274 maxb1 = min(6*(nnmax+1),nddl)
275 ELSE
276 maxb = min(3*(ntmp+1),nddl)
277 maxb1 = min(3*(nnmax+1),nddl)
278 ENDIF
279 maxb0 = maxb
280 nndl=3*numnod
281c IF (IMON>0) CALL STOPTIME(TIMERS,31)
282 IF(ispmd==0)THEN
283 WRITE(istdo,*)
284 WRITE(istdo,*)' **************************'
285 WRITE(istdo,*)' ** IMPLICIT OPTION USED **'
286 WRITE(istdo,*)' **************************'
287 WRITE(istdo,*)
288 ENDIF
289C-----integer : 1:IKC,2:IKUD,3:W_DDL,4:IADM,5:JDIM,6:NDOFI,7:IDDLI,8:INBUF_C
290 IF(nfxvel>0.AND.neig==0)THEN
291 l1=nfxvel*max(3,maxb)
292 ELSE
293 l1=0
294 ENDIF
295 s_ikc=nddl+6
296 s_ikud=l1
297C -----------W_DDL is reserved anyway------
298 s_w_ddl=nddl
299 max_l = 0
300 IF (nmonv>0.AND.isolv>=3.AND.neig==0) THEN
301 IF(ispmd==0)THEN
302 CALL ancmsg(msgid=163,anmode=aninfo)
303 ENDIF
304 CALL arret(2)
305 ENDIF
306C ----- ---------
307 IF ((isolv==1.OR.isolv>4).AND.n_pat>1) THEN
308 CALL fil_span0(nrbyac,irbyac,npby,iddl,ndof,nddl)
309 CALL dim_span(n_pat,nddl,iadk,jdik,max_l,maxb1)
310 s_iadm=nddl+1
311 ELSE
312 s_iadm=0
313 ENDIF
314 s_jdim=max_l
315 ALLOCATE(impbuf_tab%IKC(s_ikc))
316 ALLOCATE(impbuf_tab%IKUD(s_ikud))
317 ALLOCATE(impbuf_tab%W_DDL(s_w_ddl))
318 ALLOCATE(impbuf_tab%IADM(s_iadm))
319 ALLOCATE(impbuf_tab%JDIM(s_jdim))
320 s_ndofi=numnod
321 s_iddli=0
322 nint7 = 0
323 s_cand_n = 0
324 s_cand_e = 0
325 s_indsubt= 0
326 IF (ninter/=0.AND.neig==0) THEN
327 CALL dim_int7(ninter,ipari ,intbuf_tab ,nint7)
328 IF (nint7>0) THEN
329 IF (isolv==4) THEN
330 CALL ancmsg(msgid=214,anmode=aninfo)
331 CALL arret(2)
332 ENDIF
333 s_cand_n = nint7*nthread !NS_IMP
334 s_cand_e = s_cand_n !NE_IMP
335 s_indsubt= s_cand_n !IND_IMP
336 s_iddli=numnod
337 ENDIF
338 ENDIF
339 ALLOCATE(impbuf_tab%CAND_N(s_cand_n))
340 ALLOCATE(impbuf_tab%CAND_E(s_cand_e))
341 ALLOCATE(impbuf_tab%INDSUBT(s_indsubt),stat=ier1)
342 s_inbuf_c =0
343 ALLOCATE(impbuf_tab%NDOFI(s_ndofi))
344 ALLOCATE(impbuf_tab%IDDLI(s_iddli))
345 ALLOCATE(impbuf_tab%INBUF_C(s_inbuf_c))
346 IF (nspmd>1.AND.nbintc>0) CALL imp_frii(ninter)
347c IIF=II9
348C-----reel : 1,2,3,4:DIAG_K,LT_K,DIAG_M,LT_M,5,6:LB,DB,7:BKUD,8,9:D_IMP,DR_IMP
349C---- 10,11,12:ELBUF_C,BUFMAT_C,X_C,13,14:DD,DDR,15,16:X_actual,free
350C-----------17-21:FEXT,DG,DGR,DI,DIR;(Riks arc-length);22:BUFIN_C,23,24:AC,ACR,
351 s_diag_k=nddl
352 s_lt_k=nnzk
353 s_diag_m=nddl
354 s_lt_m=max(nnzk,max_l)
355 s_lb=nddl
356 IF (iline/=1.AND.tt==zero.AND.isprb==0) THEN
357 IF (nmonv>0) isigini=1
358 ELSE
359 isigini=0
360 ENDIF
361 IF ((isprb==1.OR.isigini==1.OR.ilintf>0)
362 . .AND.neig==0) THEN
363 s_lb0=nddl
364 ELSE
365 s_lb0=0
366 ENDIF
367 s_bkud=l1
368 s_d_imp=nndl
369 IF (iroddl/=0.AND.neig==0) THEN
370 s_dr_imp=nndl
371 ELSE
372 s_dr_imp=0
373 ENDIF
374 IF (iline/=1) THEN
375 CALL buf_dim1(s_elbuf_c,s_bufmat_c)
376 s_x_c=nndl
377 s_dd=nndl
378 IF (iroddl/=0) THEN
379 s_ddr=nndl
380 ELSE
381 s_ddr=0
382 ENDIF
383 ELSEIF (ilintf>0) THEN
384 CALL buf_dim1(s_elbuf_c,s_bufmat_c)
385 s_x_c=nndl
386 s_dd=0
387 s_ddr=0
388 ELSE
389 s_elbuf_c=0
390 s_bufmat_c=0
391 s_x_c=0
392 s_dd=0
393 s_ddr=0
394 ENDIF
395 s_x_a=0
396 IF (neig==0) THEN
397 IF (ismdisp>0) s_x_a=nndl
398 ENDIF
399C-----------17-21:FEXT,DG,DGR,DG0,DGR0;(Riks arc-length);22:BUFIN_C,23,24:AC,ACR,
400 s_fext=0
401 s_dg=0
402 s_dgr=0
403 s_dg0=0
404 s_dgr0=0
405 IF (idtc==3) THEN
406 s_fext=nddl
407 IF (iroddl/=0) THEN
408 s_dg=nndl
409 s_dgr=nndl
410 s_dg0=nndl
411 s_dgr0=nndl
412 ELSE
413 s_dg=nndl
414 s_dg0=nndl
415 ENDIF
416 ENDIF
417 s_bufin_c=0
418 s_ac=0
419 s_acr=0
420 IF (neig==0) THEN
421 s_ac=nndl
422 IF (iroddl/=0) s_acr=nndl
423 ENDIF
424C
425 ALLOCATE(impbuf_tab%DIAG_K(s_diag_k))
426 ALLOCATE(impbuf_tab%LT_K(s_lt_k))
427 ALLOCATE(impbuf_tab%DIAG_M(s_diag_m))
428 ALLOCATE(impbuf_tab%LT_M(s_lt_m))
429 ALLOCATE(impbuf_tab%LB(s_lb))
430 ALLOCATE(impbuf_tab%LB0(s_lb0))
431 ALLOCATE(impbuf_tab%BKUD(s_bkud))
432 ALLOCATE(impbuf_tab%D_IMP(s_d_imp))
433 ALLOCATE(impbuf_tab%DR_IMP(s_dr_imp))
434 ALLOCATE(impbuf_tab%ELBUF_C(s_elbuf_c))
435 ALLOCATE(impbuf_tab%BUFMAT_C(s_bufmat_c))
436 ALLOCATE(impbuf_tab%X_C(s_x_c))
437 ALLOCATE(impbuf_tab%DD(s_dd))
438 ALLOCATE(impbuf_tab%DDR(s_ddr))
439 ALLOCATE(impbuf_tab%X_A(s_x_a))
440 ALLOCATE(impbuf_tab%FEXT(s_fext))
441 ALLOCATE(impbuf_tab%DG(s_dg))
442 ALLOCATE(impbuf_tab%DGR(s_dgr))
443 ALLOCATE(impbuf_tab%DG0(s_dg0))
444 ALLOCATE(impbuf_tab%DGR0(s_dgr0))
445 ALLOCATE(impbuf_tab%BUFIN_C(s_bufin_c))
446 ALLOCATE(impbuf_tab%AC(s_ac))
447 ALLOCATE(impbuf_tab%ACR(s_acr),stat=ier2)
448c
449 IF (ier1/=0.OR.ier2/=0) THEN
450 CALL ancmsg(msgid=19,anmode=aninfo,
451 . c1='FOR IMPLICIT')
452 CALL arret(2)
453 ENDIF
454C---------INITIALISATION------
455 impbuf_tab%D_IMP=zero
456 impbuf_tab%IKC=0
457 nddl0 = nddl
458 nnzk0 = nnzk
459 isetk=1
460 idsc=1
461C---------IT(TOTAL),ITC,IWAIT,IDIV,R02,EIMP,E02-,Ichang-----
462 i_imp=0
463 it_bcs = 0
464 it_pcg = 0
465 impbuf_tab%R_IMP(1:25)=zero
466C ---- X_actual,V_zero-----
467 IF (neig==0) THEN
468 IF (s_x_a>0) CALL cp_real(nndl,x,impbuf_tab%X_A)
469 IF (idyna==0) CALL cp_dm(numgeo,geo,igeo,dmcp,1)
470 END IF
471C----i--initialisation for for travers nodes detection-----
472 IF (ninter/=0.AND.neig==0) THEN
473 nt_imp1=0
474 DO i=1,ninter
475 num_imp1(i)=0
476 ENDDO
477 DO j=1,nthread
478 DO i=1,ninter
479 num_impl(i,j)=0
480 ENDDO
481 ENDDO
482 ENDIF
483 CALL imp_trans0(impbuf_tab%R_IMP,impl_s0)
484 IF (iline/=1) CALL zero1(impbuf_tab%DD,nndl)
485c IBUCK= IBUCKL
486 IF (nbuck>0.AND.bisolv==2) THEN
487 WRITE(istdo,'(A)')
488 .' !! BISOLV =2 REQUIRES OPTION /IMPL/GRAPH IN RADIOSS STARTER'
489 WRITE(iout,'(A)')
490 .' !! BISOLV =2 REQUIRES OPTION /IMPL/GRAPH IN RADIOSS STARTER'
491 CALL spmd_mstop()
492 RETURN
493 ENDIF
494C
495 RETURN
496 END SUBROUTINE imp_sol_init
497!||====================================================================
498!|| imp_restarcp ../engine/source/implicit/imp_sol_init.F
499!||--- called by ------------------------------------------------------
500!|| resol ../engine/source/engine/resol.f
501!||--- calls -----------------------------------------------------
502!|| cp_dm ../engine/source/implicit/produt_v.F
503!|| cp_real ../engine/source/implicit/produt_v.F
504!|| zeror ../engine/source/system/zero.F
505!||--- uses -----------------------------------------------------
506!|| impbufdef_mod ../engine/share/modules/impbufdef_mod.F
507!||====================================================================
508 SUBROUTINE imp_restarcp(X,V,VR,GEO,IGEO,DMCP,IMPBUF_TAB)
509C-----------------------------------------------
510C M o d u l e s
511C-----------------------------------------------
512 USE impbufdef_mod
513C-----------------------------------------------
514C I m p l i c i t T y p e s
515C-----------------------------------------------
516#include "implicit_f.inc"
517C-----------------------------------------------
518C C o m m o n B l o c k s
519C-----------------------------------------------
520#include "com01_c.inc"
521#include "com04_c.inc"
522#include "param_c.inc"
523#include "impl1_c.inc"
524C-----------------------------------------------
525C D u m m y A r g u m e n t s
526C-----------------------------------------------
527 INTEGER
528 . IGEO(NPROPGI,*)
529C REAL
530 my_real
531 . GEO(NPROPG,*),X(*),DMCP(*),V(*),VR(*)
532 TYPE (IMPBUF_STRUCT_) :: IMPBUF_TAB
533C-----------------------------------------------
534C L o c a l V a r i a b l e s
535C-----------------------------------------------
536 INTEGER I,J,NNDL
537C-----------------------------------------------
538 IF (IDYNA==0) then
539 CALL zeror(v,numnod)
540 IF (iroddl/=0) CALL zeror(vr,numnod)
541 CALL cp_dm(numgeo,geo,igeo,dmcp,2)
542 ENDIF
543C ----------------for restart---
544 IF (ismdisp>0) THEN
545 nndl=3*numnod
546 CALL cp_real(nndl,impbuf_tab%X_A,x)
547 END IF
548C
549 RETURN
550 END SUBROUTINE imp_restarcp
my_real
#define my_real
Definition cppsort.cpp:32
imp_trans0
subroutine imp_trans0(r_imp, nr)
Definition wrrest.F:399
zero1
subroutine zero1(r, n)
Definition hm_read_rbe3.F:1495
imp_frii
subroutine imp_frii(ninter)
Definition imp_fri.F:456
imp_sol_init
subroutine imp_sol_init(geo, npby, lpby, itab, ipari, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, ixtg1, ixs10, ixs20, ixs16, iparg, elbuf, nint7, nbintc, x, dmcp, fr_elem, iad_elem, fr_i2m, iad_i2m, nprw, num_imp1, num_impl, monvol, igrsurf, fr_mv, ipm, igeo, iad_rby, fr_rby, sh4tree, sh3tree, irbe3, lrbe3, fr_rbe3m, iad_rbe3m, irbe2, lrbe2, ibfv, vel, elbuf_tab, iframe, intbuf_tab, nddl0, nnzk0, impbuf_tab)
Definition imp_sol_init.F:63
imp_restarcp
subroutine imp_restarcp(x, v, vr, geo, igeo, dmcp, impbuf_tab)
Definition imp_sol_init.F:509
fil_span0
subroutine fil_span0(nrbyac, irbyac, npby, iddl, ndof, nddl)
Definition ind_glob_k.F:7432
dim_span
subroutine dim_span(nn, nddl, iadk, jdik, l_nz, ndmax)
Definition ind_glob_k.F:7015
dim_int7
subroutine dim_int7(ninter, ipari, intbuf_tab, nnmax)
Definition ind_glob_k.F:4554
dim_glob_k
subroutine dim_glob_k(geo, npby, lpby, itab, nrbyac, irbyac, nint2, iint2, ipari, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, ixtg1, ixs10, ixs20, ixs16, iparg, ndof, iddl, nddl, nnzk, elbuf, inloc, lsize, fr_elem, iad_elem, fr_i2m, iad_i2m, nprw, nmonv, imonv, monvol, igrsurf, fr_mv, ipm, igeo, iad_rby, fr_rby, sh4tree, sh3tree, irbe3, lrbe3, fr_rbe3m, iad_rbe3m, irbe2, lrbe2, ibfv, vel, elbuf_tab, iframe, intbuf_tab)
Definition ind_glob_k.F:3385
ind_glob_k
subroutine ind_glob_k(npby, lpby, itab, nrbyac, irbyac, nsc, isij, nmc, imij, nss, iss, nint2, iint2, ipari, intbuf_tab, nsc2, isij2, nss2, iss2, iparg, elbuf, elbuf_tab, ixs, ixq, ixc, ixt, ixp, ixr, ixtg, ixtg1, ixs10, ixs20, ixs16, iddl, ndof, iadk, jdik, nddl, nnzk, nnmax, nkine, inloc, nkmax, nrowk, icok, icokm, nmc2, imij2, irk, npn, npp, fr_elem, iad_elem, ipm, igeo, irbe3, lrbe3, iss3, fr_i2m, iad_i2m, fr_rbe3m, iad_rbe3m, irbe2, lrbe2, isb2, nsrb2)
Definition ind_glob_k.F:3748
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
groupdef_mod
Definition groupdef_mod.F:668
impbufdef_mod
Definition impbufdef_mod.F:36
impbufdef_mod::s_lt_k
integer s_lt_k
Definition impbufdef_mod.F:77
impbufdef_mod::s_lb0
integer s_lb0
Definition impbufdef_mod.F:81
impbufdef_mod::s_ikc
integer s_ikc
Definition impbufdef_mod.F:59
impbufdef_mod::s_jdim
integer s_jdim
Definition impbufdef_mod.F:63
impbufdef_mod::s_ddr
integer s_ddr
Definition impbufdef_mod.F:89
impbufdef_mod::s_iadk
integer s_iadk
Definition impbufdef_mod.F:42
impbufdef_mod::s_elbuf_c
integer s_elbuf_c
Definition impbufdef_mod.F:85
impbufdef_mod::s_jdik
integer s_jdik
Definition impbufdef_mod.F:43
impbufdef_mod::s_ikud
integer s_ikud
Definition impbufdef_mod.F:60
impbufdef_mod::s_acr
integer s_acr
Definition impbufdef_mod.F:98
impbufdef_mod::s_cand_n
integer s_cand_n
Definition impbufdef_mod.F:70
impbufdef_mod::s_dg
integer s_dg
Definition impbufdef_mod.F:92
impbufdef_mod::s_x_a
integer s_x_a
Definition impbufdef_mod.F:90
impbufdef_mod::s_diag_k
integer s_diag_k
Definition impbufdef_mod.F:76
impbufdef_mod::s_w_ddl
integer s_w_ddl
Definition impbufdef_mod.F:61
impbufdef_mod::s_fext
integer s_fext
Definition impbufdef_mod.F:91
impbufdef_mod::s_bkud
integer s_bkud
Definition impbufdef_mod.F:82
impbufdef_mod::s_ikinw
integer s_ikinw
Definition impbufdef_mod.F:66
impbufdef_mod::s_diag_m
integer s_diag_m
Definition impbufdef_mod.F:78
impbufdef_mod::s_iddli
integer s_iddli
Definition impbufdef_mod.F:65
impbufdef_mod::s_bufmat_c
integer s_bufmat_c
Definition impbufdef_mod.F:86
impbufdef_mod::s_ac
integer s_ac
Definition impbufdef_mod.F:97
impbufdef_mod::s_cand_e
integer s_cand_e
Definition impbufdef_mod.F:71
impbufdef_mod::s_lb
integer s_lb
Definition impbufdef_mod.F:80
impbufdef_mod::s_dd
integer s_dd
Definition impbufdef_mod.F:88
impbufdef_mod::s_inbuf_c
integer s_inbuf_c
Definition impbufdef_mod.F:67
impbufdef_mod::s_bufin_c
integer s_bufin_c
Definition impbufdef_mod.F:96
impbufdef_mod::s_lt_m
integer s_lt_m
Definition impbufdef_mod.F:79
impbufdef_mod::s_x_c
integer s_x_c
Definition impbufdef_mod.F:87
impbufdef_mod::s_dgr0
integer s_dgr0
Definition impbufdef_mod.F:95
impbufdef_mod::s_iadm
integer s_iadm
Definition impbufdef_mod.F:62
impbufdef_mod::s_ndofi
integer s_ndofi
Definition impbufdef_mod.F:64
impbufdef_mod::s_dgr
integer s_dgr
Definition impbufdef_mod.F:93
impbufdef_mod::s_dg0
integer s_dg0
Definition impbufdef_mod.F:94
impbufdef_mod::s_dr_imp
integer s_dr_imp
Definition impbufdef_mod.F:84
impbufdef_mod::s_indsubt
integer s_indsubt
Definition impbufdef_mod.F:72
impbufdef_mod::s_d_imp
integer s_d_imp
Definition impbufdef_mod.F:83
message_mod
Definition message_mod.F:1257
cp_real
subroutine cp_real(n, x, xc)
Definition produt_v.F:871
buf_dim1
subroutine buf_dim1(l1, lt)
Definition produt_v.F:973
cp_dm
subroutine cp_dm(numgeo, geo, igeo, dmcp, iflag)
Definition produt_v.F:2563
resol
subroutine resol(timers, element, nodes, coupling, af, iaf, iskwn, neth, ipart, nom_opt, kxx, ixx, ixtg, ixs, ixq, ixt, ixp, ixr, ifill, mat_elem, ims, npc, ibcl, ibfv, idum, las, laccelm, nnlink, lnlink, iparg, dd_iad, igrv, iexlnk, ipari, iconx, npby, lpby, lrivet, nstrf, ljoint, nodpor, monvol, ilink, llink, linale, neflsw, nnflsw, icut, cluster, itask, inoise, thke, damp, pm, skews, geo, eani, bufmat, bufgeo, bufsf, w, veul, fill, dfill, alph, wb, dsave, asave, msnf, tf, forc, vel, fsav, fzero, xlas, accelm, agrv, fr_wave, failwave, parts0, elbuf, sensors, rby, rivet, secbuf, volmon, lambda, wa, fv, partsav, uwa, val2, phi, segvar, r, crflsw, flsw, xcut, tani, secfcum, bufnois, idata, rdata, iframe, kxsp, ixsp, nod2sp, ispsym, ispcond, xframe, spbuf, xspsym, vspsym, pv, fsavd, ibvel, lbvel, wasph, w16, isphio, lprtsph, lonfsph, vsphio, fbvel, lagbuf, ibcslag, iactiv, dampr, gjbufi, gjbufr, rbmpc, ibmpc, sphveln, nbrcvois, nbsdvois, lnrcvois, lnsdvois, nercvois, nesdvois, lercvois, lesdvois, npsegcom, lsegcom, nporgeo, ixtg1, npbyl, lpbyl, rbyl, igeo, ipm, madprt, madsh4, madsh3, madsol, madnod, madfail, iad_rby, fr_rby, iad_rby2, fr_rby2, iad_i2m, fr_i2m, addcni2, procni2, iadi2, fr_mv, iadmv2, fr_ll, fr_rl, iadcj, fr_cj, fr_sec, iad_sec, iad_cut, fr_cut, rg_cut, newfront, fr_mad, fxbipm, fxbrpm, fxbnod, fxbmod, fxbglm, fxbcpm, fxbcps, fxblm, fxbfls, fxbdls, fxbdep, fxbvit, fxbacc, fxbelm, fxbsig, fxbgrvi, fxbgrvr, eigipm, eigibuf, eigrpm, lnodpor, fr_i18, graphe, iflow, rflow, lgrav, dd_r2r, fasolfr, fr_lagf, llagf, icontact, rcontact, sh4tree, sh3tree, ipadmesh, padmesh, msc, mstg, inc, intg, ptg, iskwp, nskwp, isensp, nsensp, iaccp, naccp, ipart_state, acontact, pcontact, factiv, sh4trim, sh3trim, mscnd, incnd, ibfflux, fbfflux, rbym, irbym, lnrbym, icodrbym, ibcv, fconv, ibftemp, fbftemp, iad_rbym, fr_rbym, weight_rm, ms_ply, zi_ply, inod_pxfem, iel_pxfem, iadc_pxfem, adsky_pxfem, icode_ply, icodt_ply, iskew_ply, admsms, madclnod, nom_sect, mcpc, mcptg, dmelc, dmeltg, mssa, dmels, mstr, dmeltr, msp, dmelp, msrt, dmelrt, ibcr, fradia, res_sms, table, irbe2, lrbe2, iad_rbe2, fr_rbe2, phie, msf, procne_pxfem, iadsdp_pxfem, iadrcp_pxfem, icfield, lcfield, cfield, msz2, diag_sms, iloadp, lloadp, loadp, inod_crk, iel_crk, iadc_crk, adsky_crk, cne_crk, procne_crk, iadsdp_crk, iadrcp_crk, ibufssg_io, ibc_ply, dmint2, ibordnode, elbuf_tab, por, nodedge, iad_edge, fr_edge, fr_nbedge, crknodiad, lgauge, gauge, igaup, ngaup, nodlevxf, dd_r2r_elem, nodglobxfe, sph2sol, sol2sph, irst, dmsph, wagap, xfem_tab, elcutc, nodenr, kxfenod2elc, enrtag, rthbu f, kxig3d, ixig3d, knot, wige, wsmcomp, stack, cputime_mp_glob, cputime_mp, tab_ump, poin_ump, sol2sph_typ, irunn_bis, addcsrect, iad_frnor, fr_nor, procnor, iad_fredg, fr_edg, drape_sh4n, drape_sh3n, tab_mat, nativ0_sms, multi_fvm, segquadfr, ms_2d, h3d_data, subsets, igrnod, igrbric, igrquad, igrsh4n, igrsh3n, igrtruss, igrbeam, igrspring, igrpart, igrsurf, forneqs, nloc_dmg, iskwp_l, knotlocpc, knotlocel, pinch_data, tag_skins6, ale_connectivity, xcell, xface, ne_nercvois, ne_nesdvois, ne_lercvois, ne_lesdvois, ibcscyc, lbcscyc, t_monvol, id_global_vois, face_vois, dynain_data, fcont_max, ebcs_tab, diffusion, kloadpinter, loadpinter, dgaploadint, drapeg, user_windows, output, interfaces, dt, loads, python, dpl0cld, vel0cld, ndamp_vrel, id_damp_vrel, fr_damp_vrel, ndamp_vrel_rbyg, names_and_titles, unitab, liflow, lrflow, glob_therm, pblast, rbe3, rwall)
Definition resol.F:648
spmd_mstop
subroutine spmd_mstop()
Definition spmd_mstop.F:33
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
zeror
subroutine zeror(a, n)
Definition zero.F:39