OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
integrator.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!|| integratorl ../engine/source/implicit/integrator.F
25!||====================================================================
26 SUBROUTINE integratorl(NODFT,NODLT,D ,DR ,
27 1 X ,V ,VR ,A ,AR )
28C-----------------------------------------------
29C I m p l i c i t T y p e s
30C-----------------------------------------------
31#include "implicit_f.inc"
32C-----------------------------------------------
33C C o m m o n B l o c k s
34C-----------------------------------------------
35#include "com01_c.inc"
36#include "com08_c.inc"
37C-----------------------------------------------
38C D u m m y A r g u m e n t s
39C-----------------------------------------------
40 INTEGER NODFT,NODLT
41C REAL
42 my_real
43 . d(3,*),dr(3,*),v(3,*),vr(3,*),a(3,*),ar(3,*),x(3,*)
44C-----------------------------------------------
45C L o c a l V a r i a b l e s
46C-----------------------------------------------
47 INTEGER I,K
48 my_real
49 . dti
50C------- integrateur ----------
51 dti = one/dt2
52 DO i=nodft,nodlt
53 v(1,i)=d(1,i)*dti
54 v(2,i)=d(2,i)*dti
55 v(3,i)=d(3,i)*dti
56 a(1,i)=zero
57 a(2,i)=zero
58 a(3,i)=zero
59 ENDDO
60 IF (iroddl/=0) THEN
61 DO i=nodft,nodlt
62 vr(1,i)=dr(1,i)*dti
63 vr(2,i)=dr(2,i)*dti
64 vr(3,i)=dr(3,i)*dti
65 ar(1,i)=zero
66 ar(2,i)=zero
67 ar(3,i)=zero
68 ENDDO
69 ENDIF
70C--------------------------------------------
71 RETURN
72 END
73!||====================================================================
74!|| integrator ../engine/source/implicit/integrator.F
75!||--- calls -----------------------------------------------------
76!|| inte_dyna ../engine/source/implicit/imp_dyna.F
77!|| integrator0 ../engine/source/implicit/integrator.f
78!||====================================================================
79 SUBROUTINE integrator(NODFT,NODLT ,NDT ,D ,DR ,
80 1 X ,V ,VR ,A ,AR )
81C-----------------------------------------------
82C I m p l i c i t T y p e s
83C-----------------------------------------------
84#include "implicit_f.inc"
85C-----------------------------------------------
86C C o m m o n B l o c k s
87C-----------------------------------------------
88#include "com01_c.inc"
89#include "com08_c.inc"
90#include "impl1_c.inc"
91C-----------------------------------------------
92C D u m m y A r g u m e n t s
93C-----------------------------------------------
94 INTEGER NODFT,NODLT,NDT
95C REAL
96 my_real
97 . d(3,*),dr(3,*),v(3,*),vr(3,*),a(3,*),ar(3,*),x(3,*)
98C-----------------------------------------------
99C L o c a l V a r i a b l e s
100C-----------------------------------------------
101 INTEGER I,K
102 my_real
103 . dti
104C------- integrateur ----------
105 IF (idyna>0) CALL inte_dyna(nodft,nodlt ,
106 . d ,dr ,v ,vr )
107 dti = one/dt2/ndt
108 DO i=nodft,nodlt
109 v(1,i)=dti*d(1,i)
110 v(2,i)=dti*d(2,i)
111 v(3,i)=dti*d(3,i)
112 ENDDO
113 IF (iroddl/=0) THEN
114 DO i=nodft,nodlt
115 vr(1,i)=dti*dr(1,i)
116 vr(2,i)=dti*dr(2,i)
117 vr(3,i)=dti*dr(3,i)
118 ENDDO
119 ENDIF
120C
121 CALL integrator0(nodft,nodlt,v ,x )
122 DO i=nodft,nodlt
123 a(1,i)=zero
124 a(2,i)=zero
125 a(3,i)=zero
126 ENDDO
127 IF (iroddl/=0) THEN
128 DO i=nodft,nodlt
129 ar(1,i)=zero
130 ar(2,i)=zero
131 ar(3,i)=zero
132 ENDDO
133 ENDIF
134C--------------------------------------------
135 RETURN
136 END
137!||====================================================================
138!|| integrator0 ../engine/source/implicit/integrator.F
139!||--- called by ------------------------------------------------------
140!|| integrator ../engine/source/implicit/integrator.F
141!|| integrator_hp ../engine/source/implicit/integrator.f
142!||====================================================================
143 SUBROUTINE integrator0(NODFT,NODLT,V ,X )
144C-----------------------------------------------
145C I m p l i c i t T y p e s
146C-----------------------------------------------
147#include "implicit_f.inc"
148C-----------------------------------------------
149C C o m m o n B l o c k s
150C-----------------------------------------------
151#include "com08_c.inc"
152C-----------------------------------------------
153C D u m m y A r g u m e n t s
154C-----------------------------------------------
155 INTEGER NODFT,NODLT
156C REAL
157 my_real
158 . v(3,*),x(3,*)
159C-----------------------------------------------
160C L o c a l V a r i a b l e s
161C-----------------------------------------------
162 INTEGER I
163C------- integrateur de base---------
164 DO i=nodft,nodlt
165 x(1,i)=x(1,i)+v(1,i)*dt2
166 x(2,i)=x(2,i)+v(2,i)*dt2
167 x(3,i)=x(3,i)+v(3,i)*dt2
168 ENDDO
169C--------------------------------------------
170 RETURN
171 END
172!||====================================================================
173!|| integrator1 ../engine/source/implicit/integrator.F
174!||====================================================================
175 SUBROUTINE integrator1(NODFT,NODLT,D ,X )
176C-----------------------------------------------
177C I m p l i c i t T y p e s
178C-----------------------------------------------
179#include "implicit_f.inc"
180C-----------------------------------------------
181C D u m m y A r g u m e n t s
182C-----------------------------------------------
183 INTEGER NODFT,NODLT
184C REAL
185 my_real
186 . d(3,*),x(3,*)
187C-----------------------------------------------
188C L o c a l V a r i a b l e s
189C-----------------------------------------------
190 INTEGER I,K
191C------- integrateur de base---------
192 DO i=nodft,nodlt
193 x(1,i)=x(1,i)+d(1,i)
194 x(2,i)=x(2,i)+d(2,i)
195 x(3,i)=x(3,i)+d(3,i)
196 ENDDO
197C--------------------------------------------
198 RETURN
199 END
200!||====================================================================
201!|| integrator2 ../engine/source/implicit/integrator.F
202!||--- called by ------------------------------------------------------
203!|| al_constraint1 ../engine/source/implicit/nl_solv.F
204!|| al_constrainth1 ../engine/source/implicit/nl_solv.F
205!||====================================================================
206 SUBROUTINE integrator2(NODFT,NODLT,IDDL,NDOF,IKC,D,DR,DD,DDR)
207C-----------------------------------------------
208C I m p l i c i t T y p e s
209C-----------------------------------------------
210#include "implicit_f.inc"
211C-----------------------------------------------
212C C o m m o n B l o c k s
213C-----------------------------------------------
214#include "com01_c.inc"
215C-----------------------------------------------
216C D u m m y A r g u m e n t s
217C-----------------------------------------------
218 INTEGER NODFT,NODLT ,IDDL(*),IKC(*) ,NDOF(*)
219C REAL
220 my_real
221 . d(3,*),dr(3,*),dd(3,*),ddr(3,*)
222C-----------------------------------------------
223C L o c a l V a r i a b l e s
224C-----------------------------------------------
225 INTEGER I,J,ID
226 DO I=nodft,nodlt
227 DO j=1,3
228 id = iddl(i)+j
229 IF (ikc(id)==0) d(j,i)=d(j,i)+dd(j,i)
230 ENDDO
231 ENDDO
232 IF (iroddl/=0) THEN
233 DO i=nodft,nodlt
234 IF (ndof(i)>3) THEN
235 DO j=1,3
236 id = iddl(i)+j+3
237 IF (ikc(id)==0) dr(j,i)=dr(j,i)+ddr(j,i)
238 ENDDO
239 ENDIF
240 ENDDO
241 ENDIF
242C--------------------------------------------
243 RETURN
244 END
245!||====================================================================
246!|| frac_d ../engine/source/implicit/integrator.F
247!||====================================================================
248 SUBROUTINE frac_d(NODFT,NODLT,IDDL,NDOF,IKC,D,DR,FAC)
249C-----------------------------------------------
250C I m p l i c i t T y p e s
251C-----------------------------------------------
252#include "implicit_f.inc"
253C-----------------------------------------------
254C C o m m o n B l o c k s
255C-----------------------------------------------
256#include "com01_c.inc"
257C-----------------------------------------------
258C D u m m y A r g u m e n t s
259C-----------------------------------------------
260 INTEGER NODFT,NODLT ,IDDL(*),IKC(*) ,NDOF(*)
261C REAL
262 my_real
263 . fac,d(3,*),dr(3,*)
264C-----------------------------------------------
265C L o c a l V a r i a b l e s
266C-----------------------------------------------
267 INTEGER I,J,ID
268C-----------------------------------------------
269 DO i=nodft,nodlt
270 DO j=1,3
271 id = iddl(i)+j
272 IF (ikc(id)==0) d(j,i)=d(j,i)*fac
273 ENDDO
274 ENDDO
275 IF (iroddl/=0) THEN
276 DO i=nodft,nodlt
277 IF (ndof(i)>3) THEN
278 DO j=1,3
279 id = iddl(i)+j+3
280 IF (ikc(id)==0) dr(j,i)=dr(j,i)*fac
281 ENDDO
282 ENDIF
283 ENDDO
284 ENDIF
285C--------------------------------------------
286 RETURN
287 END
288!||====================================================================
289!|| frac_dd ../engine/source/implicit/integrator.F
290!||--- called by ------------------------------------------------------
291!|| al_constraint1 ../engine/source/implicit/nl_solv.F
292!|| al_constrainth1 ../engine/source/implicit/nl_solv.F
293!||====================================================================
294 SUBROUTINE frac_dd(NODFT,NODLT,IDDL,NDOF,IKC,D,DR,DD,DDR,FAC)
295C-----------------------------------------------
296C I m p l i c i t T y p e s
297C-----------------------------------------------
298#include "implicit_f.inc"
299C-----------------------------------------------
300C C o m m o n B l o c k s
301C-----------------------------------------------
302#include "com01_c.inc"
303C-----------------------------------------------
304C D u m m y A r g u m e n t s
305C-----------------------------------------------
306 INTEGER NODFT,NODLT ,IDDL(*),IKC(*) ,NDOF(*)
307C REAL
308 my_real
309 . fac,d(3,*),dr(3,*),dd(3,*),ddr(3,*)
310C-----------------------------------------------
311C L o c a l V a r i a b l e s
312C-----------------------------------------------
313 INTEGER I,J,ID
314C-----------------------------------------------
315 DO i=nodft,nodlt
316 DO j=1,3
317 id = iddl(i)+j
318 IF (ikc(id)==0) d(j,i)=d(j,i)+fac*dd(j,i)
319 ENDDO
320 ENDDO
321 IF (iroddl/=0) THEN
322 DO i=nodft,nodlt
323 IF (ndof(i)>3) THEN
324 DO j=1,3
325 id = iddl(i)+j+3
326 IF (ikc(id)==0) dr(j,i)=dr(j,i)+fac*ddr(j,i)
327 ENDDO
328 ENDIF
329 ENDDO
330 ENDIF
331C--------------------------------------------
332 RETURN
333 END
334C-----_HP: Hybrid SMP // inside
335!||====================================================================
336!|| integratorl_hp ../engine/source/implicit/integrator.F
337!||--- called by ------------------------------------------------------
338!|| imp_solv ../engine/source/implicit/imp_solv.F
339!||--- calls -----------------------------------------------------
340!|| imp_smpini ../engine/source/implicit/imp_solv.F
341!||====================================================================
342 SUBROUTINE integratorl_hp(D ,DR ,
343 1 X ,V ,VR ,A ,AR )
344C-----------------------------------------------
345C I m p l i c i t T y p e s
346C-----------------------------------------------
347#include "implicit_f.inc"
348C-----------------------------------------------
349C C o m m o n B l o c k s
350C-----------------------------------------------
351#include "com01_c.inc"
352#include "com04_c.inc"
353#include "com08_c.inc"
354C-----------------------------------------------
355C D u m m y A r g u m e n t s
356C-----------------------------------------------
357C REAL
358 my_real
359 . d(3,*),dr(3,*),v(3,*),vr(3,*),a(3,*),ar(3,*),x(3,*)
360C-----------------------------------------------
361C L o c a l V a r i a b l e s
362C-----------------------------------------------
363 INTEGER NODFT,NODLT,ITSK
364 INTEGER I,K
365 my_real
366 . dti
367C------- integrateur ----------
368 dti = one/dt2
369!$OMP PARALLEL PRIVATE(ITSK,NODFT ,NODLT,I)
370 CALL imp_smpini(itsk ,nodft ,nodlt ,numnod )
371 DO i=nodft,nodlt
372 v(1,i)=d(1,i)*dti
373 v(2,i)=d(2,i)*dti
374 v(3,i)=d(3,i)*dti
375 a(1,i)=zero
376 a(2,i)=zero
377 a(3,i)=zero
378 ENDDO
379 IF (iroddl/=0) THEN
380 DO i=nodft,nodlt
381 vr(1,i)=dr(1,i)*dti
382 vr(2,i)=dr(2,i)*dti
383 vr(3,i)=dr(3,i)*dti
384 ar(1,i)=zero
385 ar(2,i)=zero
386 ar(3,i)=zero
387 ENDDO
388 ENDIF
389!$OMP END PARALLEL
390C--------------------------------------------
391 RETURN
392 END
393!||====================================================================
394!|| integrator_hp ../engine/source/implicit/integrator.F
395!||--- called by ------------------------------------------------------
396!|| imp_solv ../engine/source/implicit/imp_solv.F
397!||--- calls -----------------------------------------------------
398!|| imp_smpini ../engine/source/implicit/imp_solv.F
399!|| inte_dyna ../engine/source/implicit/imp_dyna.F
400!|| integrator0 ../engine/source/implicit/integrator.F
401!||====================================================================
402 SUBROUTINE integrator_hp(NDT ,D ,DR ,
403 1 X ,V ,VR ,A ,AR )
404C-----------------------------------------------
405C I m p l i c i t T y p e s
406C-----------------------------------------------
407#include "implicit_f.inc"
408C-----------------------------------------------
409C C o m m o n B l o c k s
410C-----------------------------------------------
411#include "com01_c.inc"
412#include "com08_c.inc"
413#include "impl1_c.inc"
414#include "com04_c.inc"
415C-----------------------------------------------
416C D u m m y A r g u m e n t s
417C-----------------------------------------------
418 INTEGER NDT
419C REAL
420 my_real
421 . D(3,*),DR(3,*),V(3,*),VR(3,*),A(3,*),AR(3,*),X(3,*)
422C-----------------------------------------------
423C L o c a l V a r i a b l e s
424C-----------------------------------------------
425 INTEGER NODFT,NODLT,ITSK
426 INTEGER I,K
427 my_real
428 . DTI
429C------- integrateur ----------
430 dti = one/dt2/ndt
431!$OMP PARALLEL PRIVATE(ITSK,NODFT ,NODLT,I)
432 CALL imp_smpini(itsk ,nodft ,nodlt ,numnod )
433 IF (idyna>0) CALL inte_dyna(nodft,nodlt ,
434 . d ,dr ,v ,vr )
435 DO i=nodft,nodlt
436 v(1,i)=dti*d(1,i)
437 v(2,i)=dti*d(2,i)
438 v(3,i)=dti*d(3,i)
439 ENDDO
440 IF (iroddl/=0) THEN
441 DO i=nodft,nodlt
442 vr(1,i)=dti*dr(1,i)
443 vr(2,i)=dti*dr(2,i)
444 vr(3,i)=dti*dr(3,i)
445 ENDDO
446 ENDIF
447C
448 CALL integrator0(nodft,nodlt,v ,x )
449 DO i=nodft,nodlt
450 a(1,i)=zero
451 a(2,i)=zero
452 a(3,i)=zero
453 ENDDO
454 IF (iroddl/=0) THEN
455 DO i=nodft,nodlt
456 ar(1,i)=zero
457 ar(2,i)=zero
458 ar(3,i)=zero
459 ENDDO
460 ENDIF
461!$OMP END PARALLEL
462C--------------------------------------------
463 RETURN
464 END
465!||====================================================================
466!|| integrator1_hp ../engine/source/implicit/integrator.F
467!||--- called by ------------------------------------------------------
468!|| imp_solv ../engine/source/implicit/imp_solv.F
469!||--- calls -----------------------------------------------------
470!|| imp_smpini ../engine/source/implicit/imp_solv.F
471!||====================================================================
472 SUBROUTINE integrator1_hp(D ,X )
473C-----------------------------------------------
474C I m p l i c i t T y p e s
475C-----------------------------------------------
476#include "implicit_f.inc"
477C-----------------------------------------------
478C C o m m o n B l o c k s
479C-----------------------------------------------
480#include "com04_c.inc"
481C-----------------------------------------------
482C D u m m y A r g u m e n t s
483C-----------------------------------------------
484C REAL
485 my_real
486 . d(3,*),x(3,*)
487C-----------------------------------------------
488C L o c a l V a r i a b l e s
489C-----------------------------------------------
490 INTEGER NODFT,NODLT ,ITSK
491 INTEGER I,K
492C------- integrateur de base---------
493!$OMP PARALLEL PRIVATE(ITSK,NODFT ,NODLT,I)
494 CALL imp_smpini(itsk ,nodft ,nodlt ,numnod )
495 DO i=nodft,nodlt
496 x(1,i)=x(1,i)+d(1,i)
497 x(2,i)=x(2,i)+d(2,i)
498 x(3,i)=x(3,i)+d(3,i)
499 ENDDO
500!$OMP END PARALLEL
501C--------------------------------------------
502 RETURN
503 END
504!||====================================================================
505!|| integrator2_hp ../engine/source/implicit/integrator.F
506!||--- called by ------------------------------------------------------
507!|| al_constraint1_hp ../engine/source/implicit/nl_solv.F
508!|| nl_solv ../engine/source/implicit/nl_solv.F
509!||--- calls -----------------------------------------------------
510!|| imp_smpini ../engine/source/implicit/imp_solv.F
511!||====================================================================
512 SUBROUTINE integrator2_hp(IDDL,NDOF,IKC,D,DR,DD,DDR)
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"
522C-----------------------------------------------
523C D u m m y A r g u m e n t s
524C-----------------------------------------------
525 INTEGER IDDL(*),IKC(*) ,NDOF(*)
526C REAL
527 my_real
528 . D(3,*),DR(3,*),DD(3,*),DDR(3,*)
529C-----------------------------------------------
530C L o c a l V a r i a b l e s
531C-----------------------------------------------
532 INTEGER NODFT,NODLT ,ITSK
533 INTEGER I,J,ID
534C-----------------------------------------------
535!$OMP PARALLEL PRIVATE(ITSK,NODFT ,NODLT,ID,I,J)
536 CALL imp_smpini(itsk ,nodft ,nodlt ,numnod )
537 DO i=nodft,nodlt
538 DO j=1,3
539 id = iddl(i)+j
540 IF (ikc(id)==0) d(j,i)=d(j,i)+dd(j,i)
541 ENDDO
542 ENDDO
543 IF (iroddl/=0) THEN
544 DO i=nodft,nodlt
545 IF (ndof(i)>3) THEN
546 DO j=1,3
547 id = iddl(i)+j+3
548 IF (ikc(id)==0) dr(j,i)=dr(j,i)+ddr(j,i)
549 ENDDO
550 ENDIF
551 ENDDO
552 ENDIF
553!$OMP END PARALLEL
554C--------------------------------------------
555 RETURN
556 END
557!||====================================================================
558!|| frac_d_hp ../engine/source/implicit/integrator.f
559!||--- called by ------------------------------------------------------
560!|| nl_solv ../engine/source/implicit/nl_solv.f
561!||--- calls -----------------------------------------------------
562!|| imp_smpini ../engine/source/implicit/imp_solv.F
563!||====================================================================
564 SUBROUTINE frac_d_hp(IDDL,NDOF,IKC,D,DR,FAC)
565C-----------------------------------------------
566C I m p l i c i t T y p e s
567C-----------------------------------------------
568#include "implicit_f.inc"
569C-----------------------------------------------
570C C o m m o n B l o c k s
571C-----------------------------------------------
572#include "com01_c.inc"
573#include "com04_c.inc"
574C-----------------------------------------------
575C D u m m y A r g u m e n t s
576C-----------------------------------------------
577 INTEGER IDDL(*),IKC(*) ,NDOF(*)
578C REAL
579 my_real
580 . FAC,D(3,*),DR(3,*)
581C-----------------------------------------------
582C L o c a l V a r i a b l e s
583C-----------------------------------------------
584 INTEGER NODFT,NODLT ,ITSK
585 INTEGER I,J,ID
586C-----------------------------------------------
587!$OMP PARALLEL PRIVATE(ITSK,NODFT ,NODLT,ID,I,J)
588 CALL imp_smpini(itsk ,nodft ,nodlt ,numnod )
589 DO i=nodft,nodlt
590 DO j=1,3
591 id = iddl(i)+j
592 IF (ikc(id)==0) d(j,i)=d(j,i)*fac
593 ENDDO
594 ENDDO
595 IF (iroddl/=0) THEN
596 DO i=nodft,nodlt
597 IF (ndof(i)>3) THEN
598 DO j=1,3
599 id = iddl(i)+j+3
600 IF (ikc(id)==0) dr(j,i)=dr(j,i)*fac
601 ENDDO
602 ENDIF
603 ENDDO
604 ENDIF
605!$OMP END PARALLEL
606C--------------------------------------------
607 RETURN
608 END
609!||====================================================================
610!|| frac_dd_hp ../engine/source/implicit/integrator.f
611!||--- called by ------------------------------------------------------
612!|| al_constraint1_hp ../engine/source/implicit/nl_solv.F
613!|| nl_solv ../engine/source/implicit/nl_solv.F
614!||--- calls -----------------------------------------------------
615!|| imp_smpini ../engine/source/implicit/imp_solv.F
616!||====================================================================
617 SUBROUTINE frac_dd_hp(IDDL,NDOF,IKC,D,DR,DD,DDR,FAC)
618C-----------------------------------------------
619C I m p l i c i t T y p e s
620C-----------------------------------------------
621#include "implicit_f.inc"
622C-----------------------------------------------
623C C o m m o n B l o c k s
624C-----------------------------------------------
625#include "com01_c.inc"
626#include "com04_c.inc"
627C-----------------------------------------------
628C D u m m y A r g u m e n t s
629C-----------------------------------------------
630 INTEGER IDDL(*),IKC(*) ,NDOF(*)
631C REAL
632 my_real
633 . FAC,D(3,*),DR(3,*),DD(3,*),DDR(3,*)
634C-----------------------------------------------
635C L o c a l V a r i a b l e s
636C-----------------------------------------------
637 INTEGER NODFT,NODLT ,ITSK
638 INTEGER I,J,ID
639C-----------------------------------------------
640!$OMP PARALLEL PRIVATE(ITSK,NODFT ,NODLT,I,J,ID)
641 CALL imp_smpini(itsk ,nodft ,nodlt ,numnod )
642 DO i=nodft,nodlt
643 DO j=1,3
644 id = iddl(i)+j
645 IF (ikc(id)==0) d(j,i)=d(j,i)+fac*dd(j,i)
646 ENDDO
647 ENDDO
648 IF (iroddl/=0) THEN
649 DO i=nodft,nodlt
650 IF (ndof(i)>3) THEN
651 DO j=1,3
652 id = iddl(i)+j+3
653 IF (ikc(id)==0) dr(j,i)=dr(j,i)+fac*ddr(j,i)
654 ENDDO
655 ENDIF
656 ENDDO
657 ENDIF
658!$OMP END PARALLEL
659C--------------------------------------------
660 RETURN
661 END
my_real
#define my_real
Definition cppsort.cpp:32
inte_dyna
subroutine inte_dyna(nodft, nodlt, d, dr, v, vr)
Definition imp_dyna.F:186
imp_smpini
subroutine imp_smpini(itsk, n1ftsk, n1ltsk, n1)
Definition imp_solv.F:6895
integratorl
subroutine integratorl(nodft, nodlt, d, dr, x, v, vr, a, ar)
Definition integrator.F:28
integrator1
subroutine integrator1(nodft, nodlt, d, x)
Definition integrator.F:176
integratorl_hp
subroutine integratorl_hp(d, dr, x, v, vr, a, ar)
Definition integrator.F:344
integrator
subroutine integrator(nodft, nodlt, ndt, d, dr, x, v, vr, a, ar)
Definition integrator.F:81
integrator2
subroutine integrator2(nodft, nodlt, iddl, ndof, ikc, d, dr, dd, ddr)
Definition integrator.F:207
integrator1_hp
subroutine integrator1_hp(d, x)
Definition integrator.F:473
frac_d
subroutine frac_d(nodft, nodlt, iddl, ndof, ikc, d, dr, fac)
Definition integrator.F:249
frac_dd
subroutine frac_dd(nodft, nodlt, iddl, ndof, ikc, d, dr, dd, ddr, fac)
Definition integrator.F:295
frac_dd_hp
subroutine frac_dd_hp(iddl, ndof, ikc, d, dr, dd, ddr, fac)
Definition integrator.F:618
frac_d_hp
subroutine frac_d_hp(iddl, ndof, ikc, d, dr, fac)
Definition integrator.F:565
integrator0
subroutine integrator0(nodft, nodlt, v, x)
Definition integrator.F:144
integrator_hp
subroutine integrator_hp(ndt, d, dr, x, v, vr, a, ar)
Definition integrator.F:404
integrator2_hp
subroutine integrator2_hp(iddl, ndof, ikc, d, dr, dd, ddr)
Definition integrator.F:513
nl_solv
subroutine nl_solv(nddl, iddl, ndof, ikc, d, dr, nnz, iadk, jdik, diag_k, lt_k, f, nddli, iadi, jdii, diag_i, lt_i, itok, iadm, jdim, diag_m, lt_m, r02, dd, ddr, itask0, it, itc, ru0, rold, idiv, inprint, icprec, istop, e02, de0, eimp, inloc, nddl0, ls, u02, gap, itab, fr_elem, iad_elem, w_ddl, a, ar, v, ms, x, ipari, intbuf_tab, num_imp, ns_imp, ne_imp, nsrem, nsl, icont, graphe, fac_k, ipiv_k, nk, nmonv, imonv, monvol, igrsurf, fr_mv, volmon, ibfv, skew, xframe, mumps_par, cddlp, ind_imp, nbintc, intlist, newfront, isendto, irecvfrom, irbe3, lrbe3, ndiv, icont0, isign, fext, dg, dgr, dg0, dgr0, rfext, ls1, nodft, nodlt, irbe2, lrbe2, idiv0, relres, anew_stif)
Definition nl_solv.F:74