OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i7tri.F File Reference
#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"
#include "com04_c.inc"
#include "vect07_c.inc"

Go to the source code of this file.

Functions/Subroutines

subroutine i7tri (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, stf, stfn, j_stok, multimp, istf, itab, gap, gap_s, gap_m, igap, gapmin, gapmax, marge, gap_s_l, gap_m_l, id, titr, ix1, ix2, ix3, ix4, nsvg, prov_n, prov_e, n11, n12, n13, pene, 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, stif)

Function/Subroutine Documentation

◆ i7tri()

subroutine i7tri ( integer, dimension(*) bpe,
integer, dimension(*) pe,
integer, dimension(*) bpn,
integer, dimension(*) pn,
integer, dimension(2,0:*) add,
integer, dimension(4,*) irect,
x,
integer nb_nc,
integer nb_ec,
xyzm,
integer i_add,
integer, dimension(*) nsv,
integer i_amax,
xmax,
ymax,
zmax,
integer maxsiz,
integer i_stok,
integer i_mem,
integer nb_n_b,
integer, dimension(*) cand_n,
integer, dimension(*) cand_e,
integer nsn,
integer noint,
tzinf,
maxbox,
minbox,
stf,
stfn,
integer j_stok,
integer multimp,
integer istf,
integer, dimension(*) itab,
gap,
gap_s,
gap_m,
integer igap,
gapmin,
gapmax,
marge,
gap_s_l,
gap_m_l,
integer id,
character(len=nchartitle) titr,
integer, dimension(mvsiz), intent(inout) ix1,
integer, dimension(mvsiz), intent(inout) ix2,
integer, dimension(mvsiz), intent(inout) ix3,
integer, dimension(mvsiz), intent(inout) ix4,
integer, dimension(mvsiz), intent(inout) nsvg,
integer, dimension(mvsiz), intent(inout) prov_n,
integer, dimension(mvsiz), intent(inout) prov_e,
intent(inout) n11,
intent(inout) n12,
intent(inout) n13,
intent(inout) pene,
intent(inout) x1,
intent(inout) x2,
intent(inout) x3,
intent(inout) x4,
intent(inout) y1,
intent(inout) y2,
intent(inout) y3,
intent(inout) y4,
intent(inout) z1,
intent(inout) z2,
intent(inout) z3,
intent(inout) z4,
intent(inout) xi,
intent(inout) yi,
intent(inout) zi,
intent(inout) x0,
intent(inout) y0,
intent(inout) z0,
intent(inout) nx1,
intent(inout) ny1,
intent(inout) nz1,
intent(inout) nx2,
intent(inout) ny2,
intent(inout) nz2,
intent(inout) nx3,
intent(inout) ny3,
intent(inout) nz3,
intent(inout) nx4,
intent(inout) ny4,
intent(inout) nz4,
intent(inout) p1,
intent(inout) p2,
intent(inout) p3,
intent(inout) p4,
intent(inout) lb1,
intent(inout) lb2,
intent(inout) lb3,
intent(inout) lb4,
intent(inout) lc1,
intent(inout) lc2,
intent(inout) lc3,
intent(inout) lc4,
intent(inout) stif )

Definition at line 37 of file i7tri.F.

58
59C-----------------------------------------------
60C M o d u l e s
61C-----------------------------------------------
62 USE message_mod
63 USE names_and_titles_mod , ONLY : nchartitle
64C-----------------------------------------------
65C I m p l i c i t T y p e s
66C-----------------------------------------------
67#include "implicit_f.inc"
68C-----------------------------------------------
69C G l o b a l P a r a m e t e r s
70C-----------------------------------------------
71#include "mvsiz_p.inc"
72#include "param_c.inc"
73C-----------------------------------------------
74C C o m m o n B l o c k s
75C-----------------------------------------------
76#include "com04_c.inc"
77#include "vect07_c.inc"
78C-----------------------------------------------
79C D u m m y A r g u m e n t s
80C-----------------------------------------------
81 INTEGER NB_NC,NB_EC,I_ADD,MAXSIZ,I_STOK,J_STOK,I_MEM,ISTF
82 INTEGER I_BID, I_AMAX,NB_N_B, NOINT, NSN,MULTIMP, IGAP
83 INTEGER ADD(2,0:*),IRECT(4,*),BPE(*),PE(*),BPN(*),PN(*)
84 INTEGER NSV(*),CAND_N(*),CAND_E(*), ITAB(*)
85C REAL
86 my_real
87 . x(3,*),xyzm(6,*),tzinf,dbuc,stf(*),stfn(*),
88 . maxbox,minbox, xmax, ymax, zmax,
89 . gap, gap_s(*), gap_m(*),
90 . gapmin, gapmax, marge, gapsmx, bgapsmx,
91 . gap_s_l(*),gap_m_l(*)
92 INTEGER ID
93 CHARACTER(LEN=NCHARTITLE) :: TITR
94 INTEGER, DIMENSION(MVSIZ), INTENT(INOUT) :: PROV_N,PROV_E,NSVG
95 INTEGER, DIMENSION(MVSIZ), INTENT(INOUT) :: IX1,IX2,IX3,IX4
96 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: n11,n12,n13,pene
97 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: x1,x2,x3,x4
98 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: y1,y2,y3,y4
99 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: z1,z2,z3,z4
100 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: xi,yi,zi
101 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: x0,y0,z0
102 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx1,ny1,nz1
103 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx2,ny2,nz2
104 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx3,ny3,nz3
105 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: nx4,ny4,nz4
106 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: p1,p2,p3,p4
107 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: lb1,lb2,lb3,lb4
108 my_real, DIMENSION(MVSIZ), INTENT(INOUT) :: lc1,lc2,lc3,lc4,stif
109C-----------------------------------------------
110C L o c a l V a r i a b l e s
111C-----------------------------------------------
112 INTEGER NB_NCN,NB_ECN,ADDNN,ADDNE,IPOS,I,IP,J
113 INTEGER INF,SUP,DIR,N1,N2,N3,N4,NN,NE
114 my_real
115 . bid,dx,dy,dz,dsup,seuil,xmx,xmn,gapsmax,
116 . gapv(mvsiz)
117C-----------------------------------------------
118C ROLE DE LA ROUTINE:
119C ===================
120C CLASSE LES ELETS DE BPE ET LES NOEUDS DE BPN EN TWO ZONES
121C > OU < A UNE FRONTIERE ICI DETERMINEE ET SORT LE TOUT
122C DANS bpe,hpe, et bpn,hpn
123C-----------------------------------------------
124C D u m m y A r g u m e n t s
125C
126C NOM DESCRIPTION E/S
127C
128C BPE TABLEAU DES FACETTES A TRIER E/S
129C ET DU RESULTAT COTE MAX
130C PE TABLEAU DES FACETTES S
131C RESULTAT COTE MIN
132C BPN TABLEAU DES NOEUDS A TRIER E/S
133C ET DU RESULTAT COTE MAX
134C PN TABLEAU DES NOEUDS S
135C RESULTAT COTE MIN
136C ADD(2,*) TABLEAU DES ADRESSES E/S
137C 1.......ADRESSES NOEUDS
138C 2.......ADRESSES ELEMENTS
139C ZYZM(6,*) TABLEAU DES XYZMIN E/S
140C 1.......XMIN BOITE
141C 2.......YMIN BOITE
142C 3.......ZMIN BOITE
143C 4.......XMAX BOITE
144C 5.......YMAX BOITE
145C 6.......ZMAX BOITE
146C IRECT(4,*) TABLEAU DES CONEC FACETTES E
147C X(3,*) COORDONNEES NODALES E
148C NB_NC NOMBRE DE NOEUDS CANDIDATS E/S
149C NB_EC NOMBRE D'ELTS CANDIDATS E/S
150C I_ADD POSITION DANS LE TAB DES ADRESSES E/S
151C NSV NOS SYSTEMES DES NOEUDS E
152C XMAX plus grande abcisse existante E
153C XMAX plus grande ordonn. existante E
154C XMAX plus grande cote existante E
155C MAXSIZ TAILLE MEMOIRE MAX POSSIBLE E
156C I_STOK niveau de stockage des couples
157C candidats impact E/S
158C CAND_N boites resultats noeuds
159C CAND_E adresses des boites resultat elements
160C MULTIMP*NSN TAILLE MAX ADMISE MAINTENANT POUR LES
161C COUPLES NOEUDS,ELT CANDIDATS
162C NOINT NUMERO USER DE L'INTERFACE
163C TZINF TAILLE ZONE INFLUENCE
164C MAXBOX TAILLE MAX BUCKET
165C MINBOX TAILLE MIN BUCKET
166C=======================================================================
167C
168C
169C 1- TEST ARRET = BOITE VIDE
170C BOITE TROP PETITE
171C BOITE NE CONTENANT QU'ONE NOEUD
172C PLUS DE MEMOIRE DISPONIBLE
173C
174C-----------------------------------------------------------
175C
176C IF(MEMX>ADD(2,1)+NB_EC)THEN
177C WRITE(ISTDO,*)' *******MEM MAX=',MEMX
178C MEMX=-1
179C ELSEIF(MEMX/=-1)THEN
180C MEMX=ADD(2,1)+NB_EC
181C ENDIF
182C--------------------TEST SUR BOITE VIDES--------------
183C
184 IF(nb_ec==0.OR.nb_nc==0) THEN
185C write(6,*)" BOITE VIDE"
186C IL FAUT COPIER LES BAS DES PILES DANS BAS_DE_PILE CORRESPONDANTS
187C AVANT DE REDESCENDRE DANS LA BRANCHE MITOYENNE
188 CALL i7dstk(i_add,nb_nc,nb_ec,add,bpn,pn,bpe,pe)
189 RETURN
190 ENDIF
191C
192C-------------------TEST SUR FIN DE BRANCHE / MEMOIRE DEPASSEE------------
193C
194 dx = xyzm(4,i_add) - xyzm(1,i_add)
195 dy = xyzm(5,i_add) - xyzm(2,i_add)
196 dz = xyzm(6,i_add) - xyzm(3,i_add)
197 dsup= max(dx,dy,dz)
198C
199 IF(add(2,1)+nb_ec>=maxsiz) THEN
200C PLUS DE PLACE DANS LA PILE DES ELEMENTS BOITES TROP PETITES
201 IF ( nb_n_b == numnod) THEN
202 CALL ancmsg(msgid=83,
203 . msgtype=msgerror,
204 . anmode=aninfo,
205 . i1=id,
206 . c1=titr)
207 ENDIF
208 i_mem = 1
209 RETURN
210 ENDIF
211 IF(dsup<minbox.OR.
212 . nb_nc<=nb_n_b.AND.dsup<maxbox.OR.
213 . nb_nc<=nb_n_b.AND.nb_ec==1) THEN
214C
215C write(6,*)" NOUVELLE BOITE "
216C 1- STOCKAGE DU OU DES NOEUD CANDIDAT ET DES ELTS CORRESP.
217C VIRER LES INUTILES
218 DO 20 i=1,nb_ec
219 ne = bpe(i)
220 n1=irect(1,ne)
221 n2=irect(2,ne)
222 n3=irect(3,ne)
223 n4=irect(4,ne)
224 DO 20 j=1,nb_nc
225 nn=nsv(bpn(j))
226 IF(nn/=n1.AND.nn/=n2.AND.nn/=n3.AND.nn/=n4) THEN
227 j_stok = j_stok + 1
228 prov_n(j_stok) = bpn(j)
229 prov_e(j_stok) = ne
230 IF(j_stok==nvsiz) THEN
231 lft = 1
232 llt = nvsiz
233 nft = 0
234 j_stok = 0
235 CALL i7cor3(x ,irect,nsv ,prov_e ,prov_n,
236 2 stf ,stfn ,gapv ,igap ,gap ,
237 3 gap_s,gap_m,istf ,gapmin ,gapmax,
238 4 gap_s_l,gap_m_l ,zero ,ix1 ,ix2 ,
239 5 ix3 ,ix4 ,nsvg,x1 ,x2 ,
240 6 x3 ,x4 ,y1 ,y2 ,y3 ,
241 7 y4 ,z1 ,z2 ,z3 ,z4 ,
242 8 xi ,yi ,zi ,stif ,zero ,
243 9 llt)
244 CALL i7dst3(ix3,ix4,x1 ,x2 ,x3 ,
245 1 x4 ,y1 ,y2 ,y3 ,y4 ,
246 2 z1 ,z2 ,z3 ,z4 ,xi ,
247 3 yi ,zi ,x0 ,y0 ,z0 ,
248 4 nx1,ny1,nz1,nx2,ny2,
249 5 nz2,nx3,ny3,nz3,nx4,
250 6 ny4,nz4,p1 ,p2 ,p3 ,
251 7 p4 ,lb1,lb2,lb3,lb4,
252 8 lc1,lc2,lc3,lc4,llt)
253 CALL i7pen3(marge,gapv,n11,n12,n13 ,
254 1 pene ,nx1 ,ny1,nz1,nx2,
255 2 ny2 ,nz2 ,nx3,ny3,nz3,
256 3 nx4 ,ny4 ,nz4,p1 ,p2 ,
257 4 p3 ,p4,llt)
258 IF(i_stok+nvsiz<multimp*nsn) THEN
259 CALL i7cmp3(i_stok,cand_e ,cand_n,1,pene,
260 1 prov_n,prov_e)
261 ELSE
262 i_bid = 0
263 CALL i7cmp3(i_bid,cand_e,cand_n,0,pene,
264 1 prov_n,prov_e)
265 IF(i_stok+i_bid<multimp*nsn) THEN
266 CALL i7cmp3(i_stok,cand_e,cand_n,1,pene,
267 1 prov_n,prov_e)
268 ELSE
269 i_mem = 2
270c CALL ANSTCKI(NOINT)
271c CALL ANCWARN(103,ANINFO_BLIND_2)
272 RETURN
273 ENDIF
274 ENDIF
275 ENDIF
276C write(6,*)"Noeud candidat",BPN(J)
277C write(6,*)"Element candidat",NE
278 ENDIF
279 20 CONTINUE
280C IL FAUT COPIER LES BAS DES PILES DANS BAS_DE_PILE CORRESPONDANTS
281C AVANT DE REDESCENDRE DANS LA BRANCHE MITOYENNE
282 CALL i7dstk(i_add,nb_nc,nb_ec,add,bpn,pn,bpe,pe)
283 RETURN
284 ENDIF
285C
286C-----------------------------------------------------------
287C
288C
289C 2- PHASE DE TRI SUR LA MEDIANE SELON LA + GDE DIRECTION
290C
291C
292C-----------------------------------------------------------
293C
294C
295C 1- DETERMINER LA DIRECTION A DIVISER X,Y OU Z
296C
297 dir = 1
298 IF(dy==dsup) THEN
299 dir = 2
300 ELSE IF(dz==dsup) THEN
301 dir = 3
302 ENDIF
303 seuil =(xyzm(dir+3,i_add)+xyzm(dir,i_add))/2
304C
305C 2- DIVISER LES NOEUDS EN TWO ZONES
306C
307 nb_ncn= 0
308 addnn= add(1,1)
309 inf = 0
310 sup = 0
311 IF(igap==0)THEN
312 DO i=1,nb_nc
313 IF(x(dir,nsv(bpn(i)))<seuil) THEN
314C ON STOCKE DANS LE BAS DE LA PILE BP
315 addnn = addnn + 1
316 pn(addnn) = bpn(i)
317 inf = 1
318 ELSE
319 nb_ncn = nb_ncn + 1
320 bpn(nb_ncn) = bpn(i)
321C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPN
322 sup = 1
323 ENDIF
324 END DO
325 ELSE
326 gapsmx = zero
327 bgapsmx = zero
328 DO i=1,nb_nc
329 IF(x(dir,nsv(bpn(i)))<seuil) THEN
330C ON STOCKE DANS LE BAS DE LA PILE BP
331 addnn = addnn + 1
332 pn(addnn) = bpn(i)
333 gapsmx = max(gapsmx,gap_s(bpn(i)))
334 inf = 1
335 ELSE
336C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPN
337 nb_ncn = nb_ncn + 1
338 bpn(nb_ncn) = bpn(i)
339 bgapsmx = max(bgapsmx,gap_s(bpn(i)))
340 sup = 1
341 ENDIF
342 END DO
343 END IF
344
345CC
346CC 3- DIVISER LES ELEMENTS
347CC
348C NB_ECN= 0
349C ADDNE= ADD(2,1)
350C SEUILI = SEUIL-TZINF
351C SEUILS = SEUIL+TZINF
352C DO 85 I=1,NB_EC
353C INF = 0
354C SUP = 0
355C DO 80 J=1,4
356C IP = IRECT(J,BPE(I))
357C IF(X(DIR,IP)<SEUILS) THEN
358C INF = 1
359C IF(SUP==1) GOTO 81
360C ENDIF
361C IF(X(DIR,IP)>=SEUILI) THEN
362C SUP = 1
363C IF(INF==1) GOTO 81
364C ENDIF
365C 80 CONTINUE
366C
367C 81 CONTINUE
368C IF(INF==1) THEN
369C ON STOCKE DANS LE BAS DE LA PILE BP
370C ADDNE = ADDNE + 1
371C PE(ADDNE) = BPE(I)
372C ENDIF
373C IF(SUP==1) THEN
374C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPE
375C NB_ECN = NB_ECN + 1
376C BPE(NB_ECN) = BPE(I)
377C ENDIF
378C 85 CONTINUE
379C
380C 3- DIVISER LES ELEMENTS
381C
382C 2 LIGNES PROV POUR TEST
383C INF = 1
384C SUP = 1
385C
386 nb_ecn= 0
387 addne= add(2,1)
388 IF(igap==0)THEN
389 DO i=1,nb_ec
390 xmx = max(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
391 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
392 . + tzinf
393 xmn = min(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
394 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
395 . - tzinf
396 IF(xmn<seuil.AND.inf==1) THEN
397C ON STOCKE DANS LE BAS DE LA PILE BP
398 addne = addne + 1
399 pe(addne) = bpe(i)
400 ENDIF
401 IF(xmx>=seuil.AND.sup==1) THEN
402C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPE
403 nb_ecn = nb_ecn + 1
404 bpe(nb_ecn) = bpe(i)
405 ENDIF
406 ENDDO
407 ELSE
408 DO i=1,nb_ec
409 xmn = min(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
410 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
411 . - max(min(gapsmx+gap_m(bpe(i)),gapmax),gapmin)-marge
412 IF(xmn<seuil.AND.inf==1) THEN
413C ON STOCKE DANS LE BAS DE LA PILE BP
414 addne = addne + 1
415 pe(addne) = bpe(i)
416 ENDIF
417 xmx = max(x(dir,irect(1,bpe(i))),x(dir,irect(2,bpe(i))),
418 . x(dir,irect(3,bpe(i))),x(dir,irect(4,bpe(i))))
419 . + max(min(bgapsmx+gap_m(bpe(i)),gapmax),gapmin)+marge
420 IF(xmx>=seuil.AND.sup==1) THEN
421C ON STOCKE EN ECRASANT PROGRESSIVEMENT BPE
422 nb_ecn = nb_ecn + 1
423 bpe(nb_ecn) = bpe(i)
424 ENDIF
425 ENDDO
426 END IF
427C
428C 4- REMPLIR LES TABLEAUX D'ADRESSES
429C
430 add(1,2) = addnn
431 add(2,2) = addne
432C-----on remplit les min de la boite suivante et les max de la courante
433C (i.e. seuil est un max pour la courante)
434C on va redescendre et donc on definit une nouvelle boite
435C on remplit les max de la nouvelle boite
436C initialises dans i7buc1 a 1.E30 comme ca on recupere
437C soit XMAX soit le max de la boite
438 xyzm(1,i_add+1) = xyzm(1,i_add)
439 xyzm(2,i_add+1) = xyzm(2,i_add)
440 xyzm(3,i_add+1) = xyzm(3,i_add)
441 xyzm(4,i_add+1) = xyzm(4,i_add)
442 xyzm(5,i_add+1) = xyzm(5,i_add)
443 xyzm(6,i_add+1) = xyzm(6,i_add)
444 xyzm(dir,i_add+1) = seuil
445 xyzm(dir+3,i_add) = seuil
446C
447 nb_nc = nb_ncn
448 nb_ec = nb_ecn
449C on incremente le niveau de descente avant de sortir
450 i_add = i_add + 1
451 IF(i_add>=1000) THEN
452C TROP NIVEAUX PILE ON VAS LES PRENDRE PLUS GRANDES...
453 IF ( nb_n_b == numnod) THEN
454 CALL ancmsg(msgid=83,
455 . msgtype=msgerror,
456 . anmode=aninfo,
457 . i1=id,
458 . c1=titr)
459 ENDIF
460 i_mem = 1
461 RETURN
462 ENDIF
463C
464C ce return n'est atteint que dans le cas ok = 0
465 RETURN
my_real
#define my_real
Definition cppsort.cpp:32
i7cmp3
subroutine i7cmp3(i_stok, cand_e, cand_n, iflag, pene, prov_n, prov_e)
Definition i7cmp3.F:82
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
id
initmumps id
Definition multiplerhs_example.m:4
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
i7cor3
subroutine i7cor3(x, irect, nsv, cand_e, cand_n, stf, stfn, gapv, igap, gap, gap_s, gap_m, istf, gapmin, gapmax, gap_s_l, gap_m_l, drad, ix1, ix2, ix3, ix4, nsvg, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, xi, yi, zi, stif, dgapload, last)
Definition i7cor3.F:43
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