OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i25sto.F File Reference
#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "lockon.inc"
#include "lockoff.inc"

Go to the source code of this file.

Functions/Subroutines

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)
subroutine i25s1s2 (llt, nsn, eshift, prov_n, prov_e, nbinflg, mbinflg, pene)

Function/Subroutine Documentation

◆ i25s1s2()

subroutine i25s1s2 ( integer llt,
integer nsn,
integer eshift,
integer, dimension(*) prov_n,
integer, dimension(*) prov_e,
integer, dimension(*) nbinflg,
integer, dimension(*) mbinflg,
pene )

Definition at line 140 of file i25sto.F.

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
integer function bitget(i, n)
Definition bitget.F:37
#define my_real
Definition cppsort.cpp:32
integer i24iremp
Definition tri7box.F:423
integer, dimension(:,:), allocatable irem
Definition tri7box.F:339

◆ i25sto()

subroutine i25sto ( integer j_stok,
integer, dimension(4,*) irect,
x,
integer, dimension(*) nsv,
integer ii_stok,
integer, dimension(*) cand_n,
integer, dimension(*) cand_e,
integer mulnsn,
integer noint,
marge,
integer i_mem,
integer, dimension(mvsiz) prov_n,
integer, dimension(mvsiz) prov_e,
integer eshift,
v,
integer nsn,
integer nrtm,
gap_s,
gap_m,
curv_max,
integer nin,
pene_old,
integer, dimension(*) nbinflg,
integer, dimension(*) mbinflg,
integer ilev,
integer, dimension(*) msegtyp,
integer, dimension(*) itab,
integer igap,
gap_s_l,
gap_m_l,
integer, dimension(*) icodt,
integer, dimension(*) iskew,
intent(in) drad,
intent(in) dgapload )

Definition at line 34 of file i25sto.F.

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
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
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
subroutine i25s1s2(llt, nsn, eshift, prov_n, prov_e, nbinflg, mbinflg, pene)
Definition i25sto.F:142
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
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