OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i17xsave.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/.
23C | I17MAIN_TRI /int17/i17main_lag.F
24!||====================================================================
25!|| i17xsave ../engine/source/interfaces/int17/i17xsave.F
26!||--- called by ------------------------------------------------------
27!|| i17main_tri ../engine/source/interfaces/int17/i17main_pena.F
28!||--- uses -----------------------------------------------------
29!|| element_mod ../common_source/modules/elements/element_mod.F90
30!||====================================================================
31 SUBROUTINE i17xsave(
32 1 X ,V ,A ,IXS ,IXS16 ,
33 2 NMES ,NME_T ,ESH_T ,NELEM ,NELES ,
34 3 EMINXM,ITASK ,XSAV ,XMIN ,YMIN ,
35 4 ZMIN ,XMAX ,YMAX ,ZMAX )
36 use element_mod , only : nixs
37C saving of xsav and calculation of domain bounds
38C-----------------------------------------------
39C I m p l i c i t T y p e s
40C-----------------------------------------------
41#include "implicit_f.inc"
42#include "comlock.inc"
43C-----------------------------------------------
44C C o m m o n B l o c k s
45C-----------------------------------------------
46#include "com04_c.inc"
47#include "com08_c.inc"
48#include "task_c.inc"
49C-----------------------------------------------
50C D u m m y A r g u m e n t s
51C-----------------------------------------------
52 INTEGER NMES,NME_T,ESH_T,ITASK,
53 . NELEM(*),NELES(*),IXS(NIXS,*),IXS16(8,*)
54 my_real
55 . XMAX, YMAX, ZMAX, XMIN, YMIN, ZMIN,
56 . x(3,*), v(3,*), a(3,*), eminxm(6,*), xsav(3,*)
57C-----------------------------------------------
58C L o c a l V a r i a b l e s
59C-----------------------------------------------
60 INTEGER NMES_F, NMES_L, I, J, K
61C-----------------------------------------------
62C S o u r c e L i n e s
63C-----------------------------------------------
64C
65 nmes_f = 1 + itask*nmes / nthread
66 nmes_l = (itask+1)*nmes / nthread
67C
68 xmin=ep30
69 xmax=-ep30
70 ymin=ep30
71 ymax=-ep30
72 zmin=ep30
73 zmax=-ep30
74C
75 DO k=1,8
76 DO i=1+esh_t,nme_t+esh_t
77 j=ixs(k+1,nelem(i))
78 xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
79 xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
80 xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
81 j=ixs16(k,nelem(i)-numels8-numels10-numels20)
82 xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
83 xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
84 xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
85C
86 xmin = min( xmin , eminxm(1,i) )
87 ymin = min( ymin , eminxm(2,i) )
88 zmin = min( zmin , eminxm(3,i) )
89 xmax = max( xmax , eminxm(4,i) )
90 ymax = max( ymax , eminxm(5,i) )
91 zmax = max( zmax , eminxm(6,i) )
92 ENDDO
93 DO i=nmes_f,nmes_l
94 j=ixs(k+1,neles(i))
95 xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
96 xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
97 xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
98 j=ixs16(k,neles(i)-numels8-numels10-numels20)
99 xsav(1,j) = x(1,j)+dt2*(v(1,j)+dt12*a(1,j))
100 xsav(2,j) = x(2,j)+dt2*(v(2,j)+dt12*a(2,j))
101 xsav(3,j) = x(3,j)+dt2*(v(3,j)+dt12*a(3,j))
102 ENDDO
103 ENDDO
104C
105 RETURN
106 END
i17xsave
subroutine i17xsave(x, v, a, ixs, ixs16, nmes, nme_t, esh_t, nelem, neles, eminxm, itask, xsav, xmin, ymin, zmin, xmax, ymax, zmax)
Definition i17xsave.F:36
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21