OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i25sto.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!|| i25sto ../engine/source/interfaces/intsort/i25sto.F
25!||--- called by ------------------------------------------------------
26!|| i25trivox ../engine/source/interfaces/intsort/i25trivox.f
27!||--- calls -----------------------------------------------------
28!|| i25cor3t ../engine/source/interfaces/intsort/i25cor3t.F
29!|| i25pen3 ../engine/source/interfaces/intsort/i25pen3.F
30!|| i25s1s2 ../engine/source/interfaces/intsort/i25sto.F
31!||--- uses -----------------------------------------------------
32!|| tri7box ../engine/share/modules/tri7box.F
33!||====================================================================
34 SUBROUTINE i25sto(
35 1 J_STOK,IRECT ,X ,NSV ,II_STOK,
36 2 CAND_N,CAND_E ,MULNSN,NOINT ,MARGE ,
37 3 I_MEM ,PROV_N ,PROV_E,ESHIFT,V ,
38 4 NSN ,NRTM ,GAP_S ,GAP_M ,CURV_MAX,NIN ,
39 5 PENE_OLD,NBINFLG,MBINFLG,ILEV ,MSEGTYP,
40 6 ITAB ,IGAP,GAP_S_L,GAP_M_L,ICODT,ISKEW,
41 7 DRAD ,DGAPLOAD)
42C============================================================================
43 USE tri7box
44C-----------------------------------------------
45C I m p l i c i t T y p e s
46C-----------------------------------------------
47#include "implicit_f.inc"
48#include "comlock.inc"
49C-----------------------------------------------
50C G l o b a l P a r a m e t e r s
51C-----------------------------------------------
52#include "mvsiz_p.inc"
53C-----------------------------------------------
54C D u m m y A r g u m e n t s
55C-----------------------------------------------
56 INTEGER I_MEM, NSN, NRTM, NIN,NBINFLG(*),MBINFLG(*),ILEV, IGAP
57 INTEGER J_STOK,MULNSN,NOINT,ESHIFT
58 INTEGER IRECT(4,*),NSV(*),CAND_N(*),CAND_E(*)
59 INTEGER PROV_N(MVSIZ),PROV_E(MVSIZ),II_STOK,MSEGTYP(*),
60 . ITAB(*),ICODT(*),ISKEW(*)
61C REAL
62 my_real , INTENT(IN) :: DGAPLOAD ,DRAD
63 my_real
64 . x(3,*), v(3,*), gap_s(*), gap_m(*),
65 . marge, curv_max(*),pene_old(5,nsn),
66 . gap_s_l(*), gap_m_l(*)
67C-----------------------------------------------
68C L o c a l V a r i a b l e s
69C-----------------------------------------------
70 INTEGER I,K_STOK,I_STOK,N,NE,J,ITYPE,ISH
71 INTEGER IX1(MVSIZ), IX2(MVSIZ), IX3(MVSIZ), IX4(MVSIZ),
72 . ETYP(MVSIZ), IBC(MVSIZ)
73C REAL
74 my_real
75 . X1(MVSIZ), X2(MVSIZ), X3(MVSIZ), X4(MVSIZ),
76 . Y1(MVSIZ), Y2(MVSIZ), Y3(MVSIZ), Y4(MVSIZ),
77 . z1(mvsiz), z2(mvsiz), z3(mvsiz), z4(mvsiz),
78 . xi(mvsiz), yi(mvsiz), zi(mvsiz), stif(mvsiz),
79 . pene(mvsiz), gapv(mvsiz)
80 DATA itype/25/
81C-----------------------------------------------
82 CALL i25cor3t( j_stok ,x ,irect ,nsv ,prov_e ,
83 1 prov_n ,x1 ,x2 ,
84 2 x3 ,x4 ,y1 ,y2 ,y3 ,
85 3 y4 ,z1 ,z2 ,z3 ,z4 ,
86 4 xi ,yi ,zi ,stif ,ix1 ,
87 5 ix2 ,ix3 ,ix4 ,nsn ,gap_s ,
88 6 gap_m ,gapv ,curv_max,itype ,nin ,
89 7 v ,igap ,gap_s_l ,gap_m_l,msegtyp,
90 8 etyp ,icodt ,iskew,ibc ,drad ,
91 9 dgapload)
92C-----------------------------------------------
93 CALL i25pen3( j_stok ,marge ,x1 ,x2 ,x3 ,
94 . x4 ,y1 ,y2 ,y3 ,y4 ,
95 . z1 ,z2 ,z3 ,z4 ,xi ,
96 . yi ,zi ,pene ,ix1 ,ix2 ,
97 . ix3 ,ix4 ,gapv ,nrtm ,etyp ,
98 . ibc )
99C-----------------------------------------------
100 IF (ilev==2)
101 . CALL i25s1s2(j_stok,nsn,eshift,prov_n,prov_e,
102 . nbinflg,mbinflg,pene)
103C-----------------------------------------------
104 k_stok = 0
105 DO i=1,j_stok
106 IF(pene(i)/=zero) THEN
107 k_stok = k_stok + 1
108 END IF
109 ENDDO
110 IF(k_stok==0)RETURN
111C
112#include "lockon.inc"
113 i_stok = ii_stok
114 IF(i_stok+k_stok>mulnsn) THEN
115 i_mem = 2
116#include "lockoff.inc"
117 RETURN
118 ENDIF
119 ii_stok = i_stok + k_stok
120#include "lockoff.inc"
121 DO i=1,j_stok
122 IF(pene(i)/=zero)THEN
123 i_stok = i_stok + 1
124 cand_n(i_stok) = prov_n(i)
125 cand_e(i_stok) = prov_e(i)+eshift
126 ENDIF
127 ENDDO
128C-----------------------------------------------
129 RETURN
130 END
131!||====================================================================
132!|| i25s1s2 ../engine/source/interfaces/intsort/i25sto.F
133!||--- called by ------------------------------------------------------
134!|| i25sto ../engine/source/interfaces/intsort/i25sto.f
135!||--- calls -----------------------------------------------------
136!|| bitget ../engine/source/interfaces/intsort/i20sto.F
137!||--- uses -----------------------------------------------------
138!|| tri7box ../engine/share/modules/tri7box.F
139!||====================================================================
140 SUBROUTINE i25s1s2(LLT ,NSN ,ESHIFT,PROV_N,PROV_E,
141 . NBINFLG,MBINFLG,PENE)
142C-----------------------------------------------
143C M o d u l e s
144C-----------------------------------------------
145 USE tri7box
146C-----------------------------------------------
147C I m p l i c i t T y p e s
148C-----------------------------------------------
149#include "implicit_f.inc"
150C-----------------------------------------------
151C D u m m y A r g u m e n t s
152C-----------------------------------------------
153 INTEGER LLT,NSN,ESHIFT,PROV_N(*),PROV_E(*),NBINFLG(*),MBINFLG(*)
154C REAL
155C-----------------------------------------------
156C L o c a l V a r i a b l e s
157C-----------------------------------------------
158 INTEGER I,N,NE,IMS1,IMS2,ISS1,ISS2
159C REAL
160C-----------------------------------------------
161 my_real
162 . pene(*)
163 INTEGER BITGET
164 EXTERNAL BITGET
165C=======================================================================
166 DO i=1,llt
167 n = prov_n(i)
168 ne = prov_e(i)+eshift
169 ims1 = bitget(mbinflg(ne),0)
170 ims2 = bitget(mbinflg(ne),1)
171 IF(n <= nsn) THEN
172 iss1 = bitget(nbinflg(n),0)
173 iss2 = bitget(nbinflg(n),1)
174 ELSE
175 iss1 = bitget(irem(i24iremp+6,n-nsn),0)
176 iss2 = bitget(irem(i24iremp+6,n-nsn),1)
177 ENDIF
178 IF(.NOT.((ims1 == 1 .and. iss2==1).or.
179 . (ims2 == 1 .and. iss1==1)))THEN
180 pene(i)=zero
181 ENDIF
182 ENDDO
183C
184 RETURN
185 END
my_real
#define my_real
Definition cppsort.cpp:32
i25s1s2
subroutine i25s1s2(llt, nsn, eshift, prov_n, prov_e, nbinflg, mbinflg, pene)
Definition i25sto.F:142
i25sto
subroutine i25sto(j_stok, irect, x, nsv, ii_stok, cand_n, cand_e, mulnsn, noint, marge, i_mem, prov_n, prov_e, eshift, v, nsn, nrtm, gap_s, gap_m, curv_max, nin, pene_old, nbinflg, mbinflg, ilev, msegtyp, itab, igap, gap_s_l, gap_m_l, icodt, iskew, drad, dgapload)
Definition i25sto.F:42
i25trivox
subroutine i25trivox(nsn, nsnr, isznsnr, i_mem, vmaxdt, irect, x, stf, stfn, xyzm, nsv, ii_stok, cand_n, eshift, cand_e, mulnsn, noint, v, bgapsmx, voxel, nbx, nby, nbz, pmax_gap, nrtm, gap_s, gap_m, marge, curv_max, nin, itask, pene_old, itab, nbinflg, mbinflg, ilev, msegtyp, flagremnode, kremnod, remnod, igap, gap_s_l, gap_m_l, icodt, iskew, drad, dgapload)
Definition i25trivox.F:46
tri7box::i24iremp
integer i24iremp
Definition tri7box.F:423
tri7box::irem
integer, dimension(:,:), allocatable irem
Definition tri7box.F:339
i25cor3t
subroutine i25cor3t(jlt, x, irect, nsv, cand_e, cand_n, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, ix1, ix2, ix3, ix4, nsn, nrtm, marge, gap_s, gap_m, gapv, ityp, igap, gap_s_l, gap_m_l, edge_l2, msegtyp, etyp, icode, iskew, ibc, drad, dgapload, nrtmt)
Definition i25cor3t.F:40
i25pen3
subroutine i25pen3(jlt, cand_n, cand_e, penmin, penmax, x1, x2, x3, x4, x0, y1, y2, y3, y4, y0, z1, z2, z3, z4, z0, xi, yi, zi, nsn, ix1, ix2, ix3, ix4, nsvg, nrtm, mseglo, gaps, irect, irtlm, time_s, pene_old, itab, msegtyp, isharp, nnx, nny, nnz, gap_nm, mvoisn, gapmxl, ivis2, ibound, vtx_bisector, ilev, inacti)
Definition i25pen3.F:43