OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i7xsave.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 | I10MAIN_TRI /intsort/i10main_tri.F
24C | I11MAIN_TRI /intsort/i11main_tri.F
25C | I24MAIN_TRI /intsort/i24main_tri.F
26!||====================================================================
27!|| i7xsave ../engine/source/interfaces/intsort/i7xsave.F
28!||--- called by ------------------------------------------------------
29!|| i23main_tri ../engine/source/interfaces/intsort/i23main_tri.F
30!|| i24main_tri ../engine/source/interfaces/intsort/i24main_tri.f
31!|| i25main_tri ../engine/source/interfaces/intsort/i25main_tri.F
32!|| i7main_tri ../engine/source/interfaces/intsort/i7main_tri.F
33!||====================================================================
34 SUBROUTINE i7xsave(
35 1 X ,NSV ,MSR ,NSN ,NMN ,
36 2 ITASK ,XSAV ,XMIN ,YMIN ,ZMIN ,
37 3 XMAX ,YMAX ,ZMAX ,C_MAX ,CURV_MAX,
38 4 ICURV ,IRECT,NRTM_T,SX ,SY ,
39 5 SZ ,SX2 ,SY2 ,SZ2 ,NMN_L)
40C sauvegarde des XSAV et calcul borne domaine
41C-----------------------------------------------
42C I m p l i c i t T y p e s
43C-----------------------------------------------
44#include "implicit_f.inc"
45#include "comlock.inc"
46C-----------------------------------------------
47C C o m m o n B l o c k s
48C-----------------------------------------------
49#include "com04_c.inc"
50#include "task_c.inc"
51#include "impl1_c.inc"
52#include "ige3d_c.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56 INTEGER NSN,NMN,ITASK, ICURV, NRTM_T,
57 . NSV(*),MSR(*), IRECT(4,*),NMN_L
58 my_real
59 . XMAX, YMAX, ZMAX, XMIN, YMIN, ZMIN, C_MAX,
60 . SX, SY, SZ, SX2, SY2, SZ2,
61 . x(3,*), xsav(3,*), curv_max(*)
62C-----------------------------------------------
63C L o c a l V a r i a b l e s
64C-----------------------------------------------
65 INTEGER NSNF,NMNF,NSNL,NMNL,I, J, II
66 my_real
67 . XXX, YYY, ZZZ
68 INTEGER :: INCONV_L
69C-----------------------------------------------
70C S o u r c e L i n e s
71C-----------------------------------------------
72C
73 nsnf=1+itask*nsn/nthread
74 nsnl=(itask+1)*nsn/nthread
75 nmnf=1+itask*nmn/nthread
76 nmnl=(itask+1)*nmn/nthread
77C
78 xmin=ep30
79 xmax=-ep30
80 ymin=ep30
81 ymax=-ep30
82 zmin=ep30
83 zmax=-ep30
84C
85 sx=zero
86 sy=zero
87 sz=zero
88 sx2=zero
89 sy2=zero
90 sz2=zero
91 nmn_l = 0
92C
93 inconv_l = inconv
94 IF(nsn+nmn < numnod+numfakenodigeo)THEN
95 IF(inconv_l == 1) THEN
96 DO i=nsnf,nsnl
97 j=nsv(i)
98 IF(j>0.AND.j<=(numnod+numfakenodigeo))THEN
99 xsav(1,i)=x(1,j)
100 xsav(2,i)=x(2,j)
101 xsav(3,i)=x(3,j)
102 END IF
103 END DO
104 ENDIF
105 DO i=nmnf,nmnl
106 ii = i+nsn
107 j=msr(i)
108C msr < 0 <=> shooting nodes
109 IF(j>0) THEN
110 xmin= min(xmin,x(1,j))
111 ymin= min(ymin,x(2,j))
112 zmin= min(zmin,x(3,j))
113 xmax= max(xmax,x(1,j))
114 ymax= max(ymax,x(2,j))
115 zmax= max(zmax,x(3,j))
116 IF (inconv_l==1) THEN
117 xsav(1,ii)=x(1,j)
118 xsav(2,ii)=x(2,j)
119 xsav(3,ii)=x(3,j)
120 END IF
121 sx=sx+x(1,j)
122 sy=sy+x(2,j)
123 sz=sz+x(3,j)
124 sx2=sx2+x(1,j)**2
125 sy2=sy2+x(2,j)**2
126 sz2=sz2+x(3,j)**2
127 nmn_l= nmn_l + 1
128 ENDIF
129 END DO
130 ELSE
131C
132 IF(inconv_l == 1) THEN
133 DO i=nsnf,nsnl
134 j=nsv(i)
135 IF(j>0.AND.j<=(numnod+numfakenodigeo))THEN
136 xsav(1,j)=x(1,j)
137 xsav(2,j)=x(2,j)
138 xsav(3,j)=x(3,j)
139 END IF
140 END DO
141 ENDIF
142 DO i=nmnf,nmnl
143 j=msr(i)
144C msr < 0 <=> shooting nodes
145 IF(j>0) THEN
146 xmin= min(xmin,x(1,j))
147 ymin= min(ymin,x(2,j))
148 zmin= min(zmin,x(3,j))
149 xmax= max(xmax,x(1,j))
150 ymax= max(ymax,x(2,j))
151 zmax= max(zmax,x(3,j))
152 IF (inconv_l==1) THEN
153 xsav(1,j)=x(1,j)
154 xsav(2,j)=x(2,j)
155 xsav(3,j)=x(3,j)
156 END IF
157 sx=sx+x(1,j)
158 sy=sy+x(2,j)
159 sz=sz+x(3,j)
160 sx2=sx2+x(1,j)**2
161 sy2=sy2+x(2,j)**2
162 sz2=sz2+x(3,j)**2
163 nmn_l= nmn_l + 1
164 END IF
165 END DO
166 ENDIF
167C
168 c_max = zero
169 IF(icurv/=0)THEN
170 DO i=1,nrtm_t
171 xxx=max(x(1,irect(1,i)),x(1,irect(2,i)),
172 . x(1,irect(3,i)),x(1,irect(4,i)))
173 . -min(x(1,irect(1,i)),x(1,irect(2,i)),
174 . x(1,irect(3,i)),x(1,irect(4,i)))
175 yyy=max(x(2,irect(1,i)),x(2,irect(2,i)),
176 . x(2,irect(3,i)),x(2,irect(4,i)))
177 . -min(x(2,irect(1,i)),x(2,irect(2,i)),
178 . x(2,irect(3,i)),x(2,irect(4,i)))
179 zzz=max(x(3,irect(1,i)),x(3,irect(2,i)),
180 . x(3,irect(3,i)),x(3,irect(4,i)))
181 . -min(x(3,irect(1,i)),x(3,irect(2,i)),
182 . x(3,irect(3,i)),x(3,irect(4,i)))
183 curv_max(i) = half * max(xxx,yyy,zzz)
184 c_max = max(c_max,curv_max(i))
185 ENDDO
186 ELSE
187 DO i=1,nrtm_t
188 curv_max(i)=zero
189 ENDDO
190 ENDIF
191C
192 RETURN
193 END
i24main_tri
subroutine i24main_tri(timers, ipari, x, v, intbuf_tab, ms, nin, itask, mwag, weight, isendto, ircvfrom, retri, iad_elem, fr_elem, itab, kinet, temp, nrtm_t, renum, nsnfiold, eshift, num_imp, ind_imp, nodnx_sms, h3d_data, t2main_sms, forneqs, t2fac_sms, parameters, intheat, idt_therm, nodadt_therm)
Definition i24main_tri.F:60
i7xsave
subroutine i7xsave(x, nsv, msr, nsn, nmn, itask, xsav, xmin, ymin, zmin, xmax, ymax, zmax, c_max, curv_max, icurv, irect, nrtm_t, sx, sy, sz, sx2, sy2, sz2, nmn_l)
Definition i7xsave.F:40
ymax
subroutine ymax(idn, fac, npc, pld, stiffmin, stiffmax, stiffini, stiffavg)
Definition law100_upd.F:272
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21