OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i23tri.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!|| i23tri ../starter/source/interfaces/inter3d1/i23tri.F
25!||--- called by ------------------------------------------------------
26!|| i23buc1 ../starter/source/interfaces/inter3d1/i23buc3.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../starter/source/output/message/message.F
29!|| i23cor3t ../starter/source/interfaces/inter3d1/i23cor3t.F
30!|| i7cmp3 ../starter/source/interfaces/inter3d1/i7cmp3.F
31!|| i7dst3 ../starter/source/interfaces/inter3d1/i7dst3.F
32!|| i7dstk ../starter/source/interfaces/inter3d1/i7dstk.F
33!|| i7pen3 ../starter/source/interfaces/inter3d1/i7pen3.f
34!||--- uses -----------------------------------------------------
35!|| message_mod ../starter/share/message_module/message_mod.f
36!||====================================================================
37 SUBROUTINE i23tri(
38 1 BPE ,PE ,BPN ,PN ,ADD ,
39 2 IRECT ,X ,NB_NC ,NB_EC ,XYZM ,
40 3 I_ADD ,NSV ,I_AMAX,XMAX ,YMAX ,
41 4 ZMAX ,MAXSIZ,I_STOK,I_MEM ,NB_N_B,
42 5 CAND_N ,CAND_E,NSN ,NOINT ,TZINF ,
43 6 MAXBOX ,MINBOX,J_STOK,MSR ,
44 7 MULTIMP,ITAB ,GAP ,GAP_S ,IGAP ,
45 8 GAPMIN,GAPMAX,MARGE ,ID ,TITR ,
46 9 GAP_M ,PROV_N,PROV_E ,IX1,IX2,
47 1 IX3 ,IX4 ,NSVG ,X1 ,X2 ,
48 2 X3 ,X4 ,Y1 ,Y2 ,Y3 ,
49 3 Y4 ,Z1 ,Z2 ,Z3 ,Z4 ,
50 4 XI ,YI ,ZI ,X0 ,Y0 ,
51 5 Z0 ,NX1 ,NY1 ,NZ1,NX2,
52 6 NY2 ,NZ2 ,NX3 ,NY3,NZ3,
53 7 NX4 ,NY4 ,NZ4 ,P1 ,P2 ,
54 8 P3 ,P4 ,LB1 ,LB2,LB3,
55 9 LB4 ,LC1 ,LC2 ,LC3,LC4,
56 1 N11 ,N21 ,N31 ,PENE)
57 USE message_mod
58 USE names_and_titles_mod , ONLY : nchartitle
59C-----------------------------------------------
60C I m p l i c i t T y p e s
61C-----------------------------------------------
62#include "implicit_f.inc"
63C-----------------------------------------------
64C G l o b a l P a r a m e t e r s
65C-----------------------------------------------
66#include "mvsiz_p.inc"
67C-----------------------------------------------
68C C o m m o n B l o c k s
69C-----------------------------------------------
70#include "com01_c.inc"
71#include "com04_c.inc"
72#include "param_c.inc"
73#include "vect07_c.inc"
74C-----------------------------------------------
75C D u m m y A r g u m e n t s
76C-----------------------------------------------
77 INTEGER NB_NC,NB_EC,I_ADD,MAXSIZ,I_STOK,J_STOK,I_MEM
78 INTEGER I_BID, I_AMAX,NB_N_B, NOINT, NSN,MULTIMP, IGAP
79 INTEGER ADD(2,0:*),IRECT(4,*),BPE(*),PE(*),BPN(*),PN(*)
80 INTEGER NSV(*),CAND_N(*),CAND_E(*), ITAB(*), MSR(*)
81C REAL
82 my_real
83 . X(3,*),XYZM(6,*),TZINF,DBUC,
84 . MAXBOX,MINBOX, XMAX, YMAX, ZMAX,
85 . GAP, GAP_S(*), GAP_M(*),
86 . GAPMIN, GAPMAX, MARGE, GAPSMX, BGAPSMX
87 INTEGER ID
88 CHARACTER(LEN=NCHARTITLE) :: TITR
89 INTEGER, DIMENSION(MVSIZ), INTENT(INOUT) ::PROV_N,PROV_E
90 INTEGER, DIMENSION(MVSIZ), INTENT(INOUT) :: IX1,IX2,IX3,IX4,NSVG
91 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: X1,X2,X3,X4
92 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: Y1,Y2,Y3,Y4
93 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: Z1,Z2,Z3,Z4
94 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: XI,YI,ZI
95 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: X0,Y0,Z0
96 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: n11,n21,n31,pene
97 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx1,ny1,nz1
98 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx2,ny2,nz2
99 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx3,ny3,nz3
100 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx4,ny4,nz4
101 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: p1,p2,p3,p4
102 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: lb1,lb2,lb3,lb4
103 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: lc1,lc2,lc3,lc4
104C-----------------------------------------------
105C L o c a l V a r i a b l e s
106C-----------------------------------------------
107 INTEGER NB_NCN,NB_ECN,ADDNN,ADDNE,IPOS,I,IP,J
108 INTEGER INF,SUP,DIR,N1,N2,N3,N4,NN,NE
109C REAL
110 my_real
111 . bid,dx,dy,dz,dsup,seuil,xmx,xmn,gapsmax,
112 . gapv(mvsiz)
113C-----------------------------------------------
114C ROLE DE LA ROUTINE:
115C ===================
116C CLASSE LES ELETS DE BPE ET LES NOEUDS DE BPN EN TWO ZONES
117C > OU < A UNE FRONTIERE ICI DETERMINEE ET SORT LE TOUT
118C DANS bpe,hpe, et bpn,hpn
119C-----------------------------------------------
120C D u m m y A r g u m e n t s
121C
122C NOM DESCRIPTION E/S
123C
124C BPE TABLEAU DES FACETTES A TRIER E/S
125C ET DU RESULTAT COTE MAX
126C PE TABLEAU DES FACETTES S
127C RESULTAT COTE MIN
128C BPN TABLEAU DES NOEUDS A TRIER E/S
129C ET DU RESULTAT COTE MAX
130C PN TABLEAU DES NOEUDS S
131C RESULTAT COTE MIN
132C ADD(2,*) TABLEAU DES ADRESSES E/S
133C 1.......ADRESSES NOEUDS
134C 2.......ADRESSES ELEMENTS
135C ZYZM(6,*) TABLEAU DES XYZMIN E/S
136C 1.......XMIN BOITE
137C 2.......YMIN BOITE
138C 3.......ZMIN BOITE
139C 4.......XMAX BOITE
140C 5.......YMAX BOITE
141C 6.......ZMAX BOITE
142C IRECT(4,*) TABLEAU DES CONEC FACETTES E
143C X(3,*) COORDONNEES NODALES E
144C NB_NC NOMBRE DE NOEUDS CANDIDATS E/S
145C NB_EC NOMBRE D'ELTS CANDIDATS E/S
146C I_ADD POSITION DANS LE TAB DES ADRESSES E/S
147C NSV NOS SYSTEMES DES NOEUDS E
148C XMAX plus grande abcisse existante E
149C XMAX plus grande ordonn. existante E
150C XMAX plus grande cote existante E
151C MAXSIZ TAILLE MEMOIRE MAX POSSIBLE E
152C I_STOK niveau de stockage des couples
153C candidats impact E/S
154C CAND_N boites resultats noeuds
155C CAND_E adresses des boites resultat elements
156C MULTIMP*NSN TAILLE MAX ADMISE MAINTENANT POUR LES
157C COUPLES NOEUDS,ELT CANDIDATS
158C NOINT NUMERO USER DE L'INTERFACE
159C TZINF TAILLE ZONE INFLUENCE
160C MAXBOX TAILLE MAX BUCKET
161C MINBOX TAILLE MIN BUCKET
162C=======================================================================
163C
164C
165C 1- TEST ARRET = BOITE VIDE
166C BOITE TROP PETITE
167C BOITE NE CONTENANT QU'ONE NOEUD
168C PLUS DE MEMOIRE DISPONIBLE
169C
170C-----------------------------------------------------------
171C
172C IF(MEMX>ADD(2,1)+NB_EC)THEN
173C WRITE(ISTDO,*)' *******MEM MAX=',MEMX
174C MEMX=-1
175C ELSEIF(MEMX/=-1)THEN
176C MEMX=ADD(2,1)+NB_EC
177C ENDIF
178C--------------------TEST SUR BOITE VIDES--------------
179C
180 IF(nb_ec==0.OR.nb_nc==0) THEN
181C write(6,*)" BOITE VIDE"
182C IL FAUT COPIER LES BAS DES PILES DANS BAS_DE_PILE CORRESPONDANTS
183C AVANT DE REDESCENDRE DANS LA BRANCHE MITOYENNE
184 CALL i7dstk(i_add,nb_nc,nb_ec,add,bpn,pn,bpe,pe)
185 RETURN
186 ENDIF
187C
188C-------------------TEST SUR FIN DE BRANCHE / MEMOIRE DEPASSEE------------
189C
190 dx = xyzm(4,i_add) - xyzm(1,i_add)
191 dy = xyzm(5,i_add) - xyzm(2,i_add)
192 dz = xyzm(6,i_add) - xyzm(3,i_add)
193 dsup= max(dx,dy,dz)
194C
195 IF(add(2,1)+nb_ec>=maxsiz) THEN
196C PLUS DE PLACE DANS LA PILE DES ELEMENTS BOITES TROP PETITES
197 IF ( nb_n_b == numnod) THEN
198C WRITE(IOUT,*)'***ERROR INFINITE LOOP DETECTED '
199C WRITE(ISTDO,*)'***ERROR INFINITE LOOP DETECTED '
200C CALL ARRET(2)
201 IF (istamping == 1) THEN
202 CALL ancmsg(msgid=775,
203 . msgtype=msgerror,
204 . anmode=aninfo,
205 . i1=id,
206 . c1=titr)
207 ELSE
208 CALL ancmsg(msgid=685,
209 . msgtype=msgerror,
210 . anmode=aninfo,
211 . i1=id,
212 . c1=titr)
213 ENDIF
214 ENDIF
215 i_mem = 1
216 RETURN
217 ENDIF
218 IF(dsup<minbox.OR.
219 . nb_nc<=nb_n_b.AND.dsup<maxbox.OR.
220 . nb_nc<=nb_n_b.AND.nb_ec==1) THEN
221C
222C write(6,*)" NOUVELLE BOITE "
223C 1- STOCKAGE DU OU DES NOEUD CANDIDAT ET DES ELTS CORRESP.
224C VIRER LES INUTILES
225 DO 20 i=1,nb_ec
226 ne = bpe(i)
227 n1=irect(1,ne)
228 n2=irect(2,ne)
229 n3=irect(3,ne)
230 n4=irect(4,ne)
231 DO 20 j=1,nb_nc
232 nn=nsv(bpn(j))
233 IF(nn/=n1.AND.nn/=n2.AND.nn/=n3.AND.nn/=n4) THEN
234 j_stok = j_stok + 1
235 prov_n(j_stok) = bpn(j)
236 prov_e(j_stok) = ne
237 IF(j_stok==nvsiz) THEN
238 lft = 1
239 llt = nvsiz
240 nft = 0
241 j_stok = 0
242 CALL i23cor3t(x ,irect,nsv ,prov_e ,prov_n,
243 2 gapv ,igap ,gap ,gap_s,gapmin ,
244 3 gapmax,msr ,gap_m ,ix1 ,ix2 ,
245 4 ix3 ,ix4 ,nsvg ,x1 ,x2 ,
246 5 x3 ,x4 ,y1 ,y2 ,y3 ,
247 6 y4 ,z1 ,z2 ,z3 ,z4 ,
248 7 xi ,yi ,zi )
249 CALL i7dst3(ix3,ix4,x1 ,x2 ,x3 ,
250 1 x4 ,y1 ,y2 ,y3 ,y4 ,
251 2 z1 ,z2 ,z3 ,z4 ,xi ,
252 3 yi ,zi ,x0 ,y0 ,z0 ,
253 4 nx1,ny1,nz1,nx2,ny2,
254 5 nz2,nx3,ny3,nz3,nx4,
255 6 ny4,nz4,p1 ,p2 ,p3 ,
256 7 p4 ,lb1,lb2,lb3,lb4,
257 8 lc1,lc2,lc3,lc4,llt)
258 CALL i7pen3(marge,gapv,n11,n21,n31,
259 1 pene ,nx1 ,ny1,nz1,nx2,
260 2 ny2 ,nz2 ,nx3,ny3,nz3,
261 3 nx4 ,ny4 ,nz4,p1 ,p2 ,
262 4 p3 ,p4,llt)
263 IF(i_stok+nvsiz<multimp*nsn) THEN
264 CALL i7cmp3(i_stok,cand_e ,cand_n,1,pene,
265 1 prov_n,prov_e)
266 ELSE
267 i_bid = 0
268 CALL i7cmp3(i_bid,cand_e,cand_n,0,pene,
269 1 prov_n,prov_e)
270 IF(i_stok+i_bid<multimp*nsn) THEN
271 CALL i7cmp3(i_stok,cand_e,cand_n,1,pene,
272 1 prov_n,prov_e)
273 ELSE
274 i_mem = 2
275C WRITE(IOUT,*)' **INFO: INTERFACE NB:',NOINT
276C WRITE(IOUT,*)' TOO MANY POSSIBLE IMPACTS'
277C WRITE(IOUT,*)' SIZE OF INFLUENCE ZONE IS'
278C WRITE(IOUT,*)' MULTIPLIED BY 0.75'
279 CALL ancmsg(msgid=103,
280 . msgtype=msgwarning,
281 . anmode=aninfo_blind_2,
282 . i1=id,
283 . c1=titr)
284 RETURN
285 ENDIF
286 ENDIF
287 ENDIF
288C write(6,*)"Noeud candidat",BPN(J)
289C write(6,*)"Element candidat",NE
290 ENDIF
291 20 CONTINUE
292C IL FAUT COPIER LES BAS DES PILES DANS BAS_DE_PILE CORRESPONDANTS
293C AVANT DE REDESCENDRE DANS LA BRANCHE MITOYENNE
294 CALL i7dstk(i_add,nb_nc,nb_ec,add,bpn,pn,bpe,pe)
295 RETURN
296 ENDIF
297C
298C-----------------------------------------------------------
299C
300C
301C 2- PHASE DE TRI SUR LA MEDIANE SELON LA + GDE DIRECTION
302C
303C
304C-----------------------------------------------------------
305C
306C
307C 1- DETERMINER LA DIRECTION A DIVISER X,Y OU Z
308C
309 dir = 1
310 IF(dy==dsup) THEN
311 dir = 2
312 ELSE IF(dz==dsup) THEN
313 dir = 3
314 ENDIF
315 seuil =(xyzm(dir+3,i_add)+xyzm(dir,i_add))/2
316C
317C 2- DIVISER LES NOEUDS EN TWO ZONES
318C
319 nb_ncn= 0
320 addnn= add(1,1)
321 inf = 0
322 sup = 0
323 IF(igap==0)THEN
324 DO i=1,nb_nc
325 IF(x(dir,nsv(bpn(i)))<seuil) THEN
326C ON STOCKE DANS LE BAS DE LA PILE BP
327 addnn = addnn + 1
328 pn(addnn) = bpn(i)
329 inf = 1
330 ELSE
331 nb_ncn = nb_ncn + 1
332 bpn(nb_ncn) = bpn(i)
333C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPN
334 sup = 1
335 ENDIF
336 END DO
337 ELSE
338 gapsmx = zero
339 bgapsmx = zero
340 DO i=1,nb_nc
341 IF(x(dir,nsv(bpn(i)))<seuil) THEN
342C ON STOCKE DANS LE BAS DE LA PILE BP
343 addnn = addnn + 1
344 pn(addnn) = bpn(i)
345 gapsmx = max(gapsmx,gap_s(bpn(i)))
346 inf = 1
347 ELSE
348C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPN
349 nb_ncn = nb_ncn + 1
350 bpn(nb_ncn) = bpn(i)
351 bgapsmx = max(bgapsmx,gap_s(bpn(i)))
352 sup = 1
353 ENDIF
354 END DO
355 END IF
356CC
357CC 3- DIVISER LES ELEMENTS
358CC
359C NB_ECN= 0
360C ADDNE= ADD(2,1)
361C SEUILI = SEUIL-TZINF
362C SEUILS = SEUIL+TZINF
363C DO 85 I=1,NB_EC
364C INF = 0
365C SUP = 0
366C DO 80 J=1,4
367C IP = IRECT(J,BPE(I))
368C IF(X(DIR,IP)<SEUILS) THEN
369C INF = 1
370C IF(SUP==1) GOTO 81
371C ENDIF
372C IF(X(DIR,IP)>=SEUILI) THEN
373C SUP = 1
374C IF(INF==1) GOTO 81
375C ENDIF
376C 80 CONTINUE
377C
378C 81 CONTINUE
379C IF(INF==1) THEN
380C ON STOCKE DANS LE BAS DE LA PILE BP
381C ADDNE = ADDNE + 1
382C PE(ADDNE) = BPE(I)
383C ENDIF
384C IF(SUP==1) THEN
385C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPE
386C NB_ECN = NB_ECN + 1
387C BPE(NB_ECN) = BPE(I)
388C ENDIF
389C 85 CONTINUE
390C
391C 3- DIVISER LES ELEMENTS
392C
393C 2 LIGNES PROV POUR TEST
394C INF = 1
395C SUP = 1
396C
397 nb_ecn= 0
398 addne= add(2,1)
399 IF(igap==0)THEN
400 DO i=1,nb_ec
401 xmx=max(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
402 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
403 . + tzinf
404 xmn=min(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
405 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
406 . - tzinf
407 IF(xmn<seuil.AND.inf==1) THEN
408C ON STOCKE DANS LE BAS DE LA PILE BP
409 addne = addne + 1
410 pe(addne) = bpe(i)
411 ENDIF
412 IF(xmx>=seuil.AND.sup==1) THEN
413C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPE
414 nb_ecn = nb_ecn + 1
415 bpe(nb_ecn) = bpe(i)
416 ENDIF
417 ENDDO
418 ELSE
419 DO i=1,nb_ec
420 xmn=min(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
421 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
422 - -sqrt(three)*min(max(gapsmx+gap_m(bpe(i)),gapmin),gapmax)
423 - -marge
424 IF(xmn<seuil.AND.inf==1) THEN
425C ON STOCKE DANS LE BAS DE LA PILE BP
426 addne = addne + 1
427 pe(addne) = bpe(i)
428 ENDIF
429 xmx=max(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
430 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
431 + +sqrt(three)*min(max(bgapsmx+gap_m(bpe(i)),gapmin),gapmax)
432 + +marge
433 IF(xmx>=seuil.AND.sup==1) THEN
434C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPE
435 nb_ecn = nb_ecn + 1
436 bpe(nb_ecn) = bpe(i)
437 ENDIF
438 ENDDO
439 END IF
440C
441C 4- REMPLIR LES TABLEAUX D'ADRESSES
442C
443 add(1,2) = addnn
444 add(2,2) = addne
445C-----on remplit les min de la boite suivante et les max de la courante
446C (i.e. seuil est un max pour la courante)
447C on va redescendre et donc on definit une nouvelle boite
448C on remplit les max de la nouvelle boite
449C initialises dans i7buc1 a 1.E30 comme ca on recupere
450C soit XMAX soit le max de la boite
451 xyzm(1,i_add+1) = xyzm(1,i_add)
452 xyzm(2,i_add+1) = xyzm(2,i_add)
453 xyzm(3,i_add+1) = xyzm(3,i_add)
454 xyzm(4,i_add+1) = xyzm(4,i_add)
455 xyzm(5,i_add+1) = xyzm(5,i_add)
456 xyzm(6,i_add+1) = xyzm(6,i_add)
457 xyzm(dir,i_add+1) = seuil
458 xyzm(dir+3,i_add) = seuil
459C
460 nb_nc = nb_ncn
461 nb_ec = nb_ecn
462C on incremente le niveau de descente avant de sortir
463 i_add = i_add + 1
464 IF(i_add>=1000) THEN
465C TROP NIVEAUX PILE ON VAS LES PRENDRE PLUS GRANDES...
466 IF ( nb_n_b == numnod) THEN
467C WRITE(IOUT,*)'***COMPUTATION STOPPED WHILE INFINITELY LOOPING'
468C WRITE(ISTDO,*)'***COMPUTATION STOPPED WHILE INFINITELY LOOPING'
469C CALL ARRET(2)
470 IF (istamping == 1) THEN
471 CALL ancmsg(msgid=775,
472 . msgtype=msgerror,
473 . anmode=aninfo,
474 . i1=id,
475 . c1=titr)
476 ELSE
477 CALL ancmsg(msgid=685,
478 . msgtype=msgerror,
479 . anmode=aninfo,
480 . i1=id,
481 . c1=titr)
482 ENDIF
483 ENDIF
484 i_mem = 1
485 RETURN
486 ENDIF
487C
488C ce return n'est atteint que dans le cas ok = 0
489 RETURN
490 END
i23tri
subroutine i23tri(bpe, pe, bpn, pn, add, irect, x, nb_nc, nb_ec, xyzm, i_add, nsv, i_amax, xmax, ymax, zmax, maxsiz, i_stok, i_mem, nb_n_b, cand_n, cand_e, nsn, noint, tzinf, maxbox, minbox, j_stok, msr, multimp, itab, gap, gap_s, igap, gapmin, gapmax, marge, id, titr, gap_m, prov_n, prov_e, ix1, ix2, ix3, ix4, nsvg, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, x0, y0, z0, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, n11, n21, n31, pene)
Definition i23tri.F:57
i7cmp3
subroutine i7cmp3(i_stok, cand_e, cand_n, iflag, pene, prov_n, prov_e)
Definition i7cmp3.F:82
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
message_mod
Definition message_mod.F:1249
names_and_titles_mod
Definition names_and_titles_mod.F:997
names_and_titles_mod::nchartitle
integer, parameter nchartitle
Definition names_and_titles_mod.F:1003
i23cor3t
subroutine i23cor3t(x, irect, nsv, cand_e, cand_n, gapv, igap, gap, gap_s, gapmin, gapmax, msr, gap_m, ix1, ix2, ix3, ix4, nsvg, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi)
Definition i23cor3t.F:36
i7dst3
subroutine i7dst3(ix3, ix4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, x0, y0, z0, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, lb1, lb2, lb3, lb4, lc1, lc2, lc3, lc4, last)
Definition i7dst3.F:46
i7dstk
subroutine i7dstk(i_add, nb_nc, nb_ec, add, bpn, pn, bpe, pe)
Definition i7dstk.F:34
i7pen3
subroutine i7pen3(marge, gapv, n1, n2, n3, pene, nx1, ny1, nz1, nx2, ny2, nz2, nx3, ny3, nz3, nx4, ny4, nz4, p1, p2, p3, p4, last)
Definition i7pen3.F:43
ancmsg
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
Definition message.F:889
starter
program starter
Definition starter.F:39