OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
i11buce.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!|| i11buce_vox ../engine/source/interfaces/intsort/i11buce.F
25!||--- called by ------------------------------------------------------
26!|| i11main_tri ../engine/source/interfaces/intsort/i11main_tri.F
27!||--- calls -----------------------------------------------------
28!|| ancmsg ../engine/source/output/message/message.F
29!|| arret ../engine/source/system/arret.F
30!|| i11trivox ../engine/source/interfaces/intsort/i11trivox.F
31!|| my_barrier ../engine/source/system/machine.F
32!||--- uses -----------------------------------------------------
33!|| message_mod ../engine/share/message_module/message_mod.F
34!|| tri11 ../engine/share/modules/tri11_mod.F
35!|| tri7box ../engine/share/modules/tri7box.F
36!||====================================================================
37 SUBROUTINE i11buce_vox(
38 1 X ,IRECTS ,IRECTM ,NRTS ,NMN ,
39 2 NRTM ,NSN ,CAND_M ,CAND_S ,MAXGAP ,
40 3 NOINT ,II_STOK ,TZINF ,MAXBOX ,MINBOX ,
41 4 NB_N_B ,ESHIFT ,ILD ,BMINMA ,NCONTACT,
42 6 ADDCM ,CHAINE ,NIN ,ITAB ,NRTSR ,
43 7 NCONT ,GAP_S , STIFS ,PENIS ,IGAP ,
44 8 STIFM ,IAUTO , I_MEM ,ITASK ,IFORM ,
45 9 IFPEN ,DRAD , GAP_M , GAP_S_L,
46 1 GAP_M_L ,GAPMIN, BGAPSMX, GAP,
47 2 FLAGREMNODE,KREMNODE,REMNODE,DGAPLOAD)
48C-----------------------------------------------
49C M o d u l e s
50C-----------------------------------------------
51 USE message_mod
52 USE tri7box
53 USE tri11
54
55C============================================================================
56C I m p l i c i t T y p e s
57C-----------------------------------------------
58#include "implicit_f.inc"
59#include "comlock.inc"
60C-----------------------------------------------
61C C o m m o n B l o c k s
62C-----------------------------------------------
63#include "units_c.inc"
64#include "warn_c.inc"
65C-----------------------------------------------
66C D u m m y A r g u m e n t s
67C-----------------------------------------------
68 INTEGER NMN, NRTM, NSN, NOINT,IDT,NRTS, NIN, NRTSR,
69 . IGAP,IAUTO, I_MEM, ITASK
70 INTEGER IRECTS(2,*),IRECTM(2,*),ADDCM(*),CHAINE(2,*)
71 INTEGER CAND_M(*),CAND_S(*),IFPEN(*),FLAGREMNODE,KREMNODE(*),REMNODE(*)
72 INTEGER ESHIFT,ILD,NB_N_B, NCONTACT, NCONT, ITAB(*),
73 . IFORM,II_STOK
74C REAL
75 my_real
76 . TZINF,MAXBOX,MINBOX,BMINMA(6),BGAPSMX
77 my_real
78 . maxgap,gapmin,gap
79 my_real , INTENT(IN) :: dgapload,drad
80 my_real
81 . x(3,*),stifs(*),penis(2,*),stifm(*),
82 . gap_s(*),gap_m(*),gap_s_l(*),gap_m_l(*)
83C-----------------------------------------------
84C L o c a l V a r i a b l e s
85C-----------------------------------------------
86 INTEGER I_ADD_MAX,ISZNSNR
87 parameter(i_add_max = 1001)
88C
89 INTEGER I, J, N1, N2, I_ADD, MAXSIZ,JJ,
90 . add(2,i_add_max), n
91 my_real
92 . xyzm(6,i_add_max-1), marge, aaa
93 my_real
94 . dd,dd1,marge_st,dx1,dy1,dz1
95 INTEGER :: L
96 INTEGER NB_OLD(2,I_ADD_MAX+1)
97 INTEGER NBX,NBY,NBZ
98 INTEGER (KIND=8) :: NBX8,NBY8,NBZ8,RES8,LVOXEL8
99C-----------------------------------------------
100C definition from TRI7BOX module
101C-----------------------------------------------
102C
103C-----------------------------------------------
104C S o u r c e L i n e s
105C-----------------------------------------------
106C
107C----- TRI PAR BOITES
108C
109C-----------------------------------------------
110C SI ON A PAS ASSEZ DE MEMOIRE POUR LES PILES ON RECOMMENCE LE TRI
111C EN INCREMENTANT LE NB_N_B (NOMBRE DE NOEUDS PAR BOITE FINIE)
112C
113C POINTEUR NOM TAILLE
114C P1........Elt Bas Pile NRTM
115C P2........Elt PILE 3*NRTM
116C P21.......Elt Bas Pile NRTS
117C P22.......Elt PILE 3*NRTS
118 maxsiz = 3*(max(nrtm,nrts+nrtsr)+100)
119C
120C-----INITIALISATION DES ADRESSES ET X,Y,Z
121C
122C ADDE ADDN X Y Z
123C 1 1 XMIN YMIN ZMIN
124C 1 1 XMAX YMAX ZMAX
125C
126 add(1,1) = 0
127 add(2,1) = 0
128 add(1,2) = 0
129 add(2,2) = 0
130 i_add = 1
131 xyzm(1,i_add) = bminma(4)
132 xyzm(2,i_add) = bminma(5)
133 xyzm(3,i_add) = bminma(6)
134 xyzm(4,i_add) = bminma(1)
135 xyzm(5,i_add) = bminma(2)
136 xyzm(6,i_add) = bminma(3)
137 i_mem = 0
138C
139 IF (iform /= 2) THEN
140 isznsnr = 0
141 DO i=1,nrtm
142 addcm(i)=0
143 ENDDO
144 ELSE
145 isznsnr = nrtsr
146 ENDIF
147C
148C
149C-----DEBUT DE LA PHASE DE TRI
150
151 marge = tzinf - max(maxgap+dgapload,drad)
152
153 IF( nmn /= 0 ) THEN
154 aaa = sqrt(nmn /
155 . ((bminma(1)-bminma(4))*(bminma(2)-bminma(5))
156 . +(bminma(2)-bminma(5))*(bminma(3)-bminma(6))
157 . +(bminma(3)-bminma(6))*(bminma(1)-bminma(4))))
158 ELSE
159 aaa = 0
160 ENDIF
161
162 aaa = 0.75*aaa
163
164 nbx = nint(aaa*(bminma(1)-bminma(4)))
165 nby = nint(aaa*(bminma(2)-bminma(5)))
166 nbz = nint(aaa*(bminma(3)-bminma(6)))
167 nbx = max(nbx,1)
168 nby = max(nby,1)
169 nbz = max(nbz,1)
170
171 nbx8=nbx
172 nby8=nby
173 nbz8=nbz
174 res8=(nbx8+2)*(nby8+2)*(nbz8+2)
175 lvoxel8 = lvoxel
176
177 IF(res8 > lvoxel8)THEN
178 aaa = lvoxel
179 aaa = aaa/((nbx8+2)*(nby8+2)*(nbz8+2))
180 aaa = aaa**(third)
181 nbx = int((nbx+2)*aaa)-2
182 nby = int((nby+2)*aaa)-2
183 nbz = int((nbz+2)*aaa)-2
184 nbx = max(nbx,1)
185 nby = max(nby,1)
186 nbz = max(nbz,1)
187 nbx8 = nbx
188 nby8 = nby
189 nbz8 = nbz
190 res8=(nbx8+2)*(nby8+2)*(nbz8+2)
191 END IF
192C If still too many voxels:
193C Reduction of the number of voxels
194 IF(res8 > lvoxel8) THEN
195 nbx = min(100,max(nbx8,1))
196 nby = min(100,max(nby8,1))
197 nbz = min(100,max(nbz8,1))
198 END IF
199
200 DO i=inivoxel,(nbx+2)*(nby+2)*(nbz+2)
201 voxel1(i)=0
202 ENDDO
203 inivoxel = max(inivoxel,(nbx+2)*(nby+2)*(nbz+2)+1)
204
205 !print *, "voxel search"
206
207 CALL i11trivox(
208 1 irects ,irectm ,x ,nrtm ,nrtsr ,
209 2 xyzm ,ii_stok ,cand_s ,cand_m ,ncontact,
210 3 noint ,tzinf ,i_mem ,eshift ,addcm ,
211 4 chaine ,nrts ,itab ,stifs ,stifm ,
212 5 iauto ,voxel1 ,nbx ,nby ,nbz ,
213 6 itask ,ifpen ,iform ,gapmin ,drad ,
214 7 marge ,gap_s ,gap_m ,gap_s_l, gap_m_l,
215 8 bgapsmx, igap ,gap ,flagremnode,kremnode,
216 9 remnode,dgapload )
217
218 CALL my_barrier
219
220 100 CONTINUE
221
222C
223C I_MEM = 1 ==> N/A
224C I_MEM = 2 ==> PAS ASSEZ DE MEMOIRE CANDIDATS
225C I_MEM = 3 ==> N/A
226
227
228C === WRITE CANDIDATES IN fort.[900 + ISPMD] files
229c IF ( ) THEN
230c DO I = 1, II_STOK
231c IF(CAND_S(I) <= NRTS) THEN
232c WRITE(900+ISPMD,*) NOINT,ITAB(IRECTS(1,CAND_S(I))),
233c . ITAB(IRECTS(2,CAND_S(I))),ITAB(IRECTM(1,CAND_M(I))),ITAB(IRECTM(2,CAND_M(I)))
234c ELSE
235c WRITE(900+ISPMD,*) -NOINT,ITAB(IRECTS(1,CAND_S(I)))
236c . ,ITAB(IRECTS(2,CAND_S(I))),ITAB(IRECTM(1,CAND_M(I))),ITAB(IRECTM(2,CAND_M(I)))
237c ENDIF
238c ENDDO
239c CALL FLUSH(900+ISPMD)
240c STOP
241c ENDIF
242
243
244
245 IF (i_mem == 2) RETURN
246
247 IF(i_mem==1)THEN
248 nb_n_b = nb_n_b + 1
249 IF ( nb_n_b > max(nrtm,nrts)) THEN
250 CALL ancmsg(msgid=85,anmode=aninfo,
251 . i1=noint)
252 CALL arret(2)
253 ENDIF
254 ild = 1
255 ELSEIF(i_mem==2) THEN
256 IF(debug(1)>=1) THEN
257 iwarn = iwarn+1
258#include "lockon.inc"
259 WRITE(istdo,*)' **WARNING INTERFACE/MEMORY'
260 WRITE(iout,*)' **WARNING INTERFACE NB:',noint
261 WRITE(iout,*)' TOO MANY POSSIBLE IMPACTS'
262 WRITE(iout,*)' SIZE OF INFLUENCE ZONE IS'
263 WRITE(iout,*)' MULTIPLIED BY 0.75'
264#include "lockoff.inc"
265 ENDIF
266 tzinf = three_over_4*tzinf
267C ne pas dimunuer la taille des boite
268C MINBOX= THREE_OVER_4*MINBOX
269C MAXBOX= THREE_OVER_4*MAXBOX
270 IF( tzinf<=max(maxgap+dgapload,drad) ) THEN
271 CALL ancmsg(msgid=98,anmode=aninfo,
272 . i1=noint,c1='(I11BUCE)')
273 CALL arret(2)
274 ENDIF
275 ild = 1
276 ELSEIF(i_mem==3)THEN
277 nb_n_b = nb_n_b + 1
278 IF ( nb_n_b > max(nrtm,nrts)) THEN
279 CALL ancmsg(msgid=99,anmode=aninfo,
280 . i1=noint,c1='(I11BUCE)')
281 CALL arret(2)
282 ENDIF
283 ild = 1
284 ENDIF
285C
286 RETURN
287 END
288
289C --- --- --- Chd|====================================================================
290C --- --- --- Chd| I11BUCE src/intsort/i11buce.F
291C --- --- --- Chd|-- called by -----------
292C --- --- --- Chd| I11MAIN_TRI src/intsort/i11main_tri.F
293C --- --- --- Chd|-- calls ---------------
294C --- --- --- Chd| ANCMSG src/message/message.F
295C --- --- --- Chd| ARRET src/sortie/arret.F
296C --- --- --- Chd| I11TRI src/intsort/i11tri.F
297C --- --- --- Chd| I11TRIVOX src/intsort/i11trivox.F
298C --- --- --- Chd| MY_BARRIER src/spe/machine.F
299C --- --- --- Chd| MESSAGE_MOD share/message_module/message_mod.F
300C --- --- --- Chd| TRI11 share/modules/tri11_mod.F
301C --- --- --- Chd| TRI7BOX share/modules/tri7box.F
302C --- --- --- Chd|====================================================================
303C --- --- --- SUBROUTINE I11BUCE(
304C --- --- --- 1 X ,IRECTS ,IRECTM ,NRTS ,NMN ,
305C --- --- --- 2 NRTM ,NSN ,CAND_M ,CAND_S ,GAP ,
306C --- --- --- 3 NOINT ,II_STOK ,TZINF ,MAXBOX ,MINBOX ,
307C --- --- --- 4 NB_N_B ,ESHIFT ,ILD ,BMINMA ,NCONTACT,
308C --- --- --- 6 ADDCM ,CHAINE ,NIN ,ITAB ,NRTSR ,
309C --- --- --- 7 NCONT ,GAP_S , STIFS ,PENIS ,IGAP ,
310C --- --- --- 8 STIFM ,IAUTO , I_MEM ,ITASK ,IFORM ,
311C --- --- --- 9 IFPEN ,DRAD )
312C --- --- --- C-----------------------------------------------
313C --- --- --- C M o d u l e s
314C --- --- --- C-----------------------------------------------
315C --- --- --- USE MESSAGE_MOD
316C --- --- --- USE TRI7BOX
317C --- --- --- USE TRI11
318C --- --- ---
319C --- --- --- C============================================================================
320C --- --- --- C I m p l i c i t T y p e s
321C --- --- --- C-----------------------------------------------
322C --- --- --- #include "implicit_f.inc"
323C --- --- --- #include "comlock.inc"
324C --- --- --- C-----------------------------------------------
325C --- --- --- C C o m m o n B l o c k s
326C --- --- --- C-----------------------------------------------
327C --- --- --- #include "units_c.inc"
328C --- --- --- #include "warn_c.inc"
329C --- --- --- #include "com01_c.inc"
330C --- --- --- #include "com04_c.inc"
331C --- --- --- #include "com08_c.inc"
332C --- --- --- #include "scr05_c.inc"
333C --- --- --- #include "task_c.inc"
334C --- --- --- #include "param_c.inc"
335C --- --- --- C-----------------------------------------------
336C --- --- --- C D u m m y A r g u m e n t s
337C --- --- --- C-----------------------------------------------
338C --- --- --- INTEGER NMN, NRTM, NSN, NOINT,IDT,NRTS, NIN, NRTSR,
339C --- --- --- . IGAP,IAUTO, I_MEM, ITASK
340C --- --- --- INTEGER IRECTS(2,*),IRECTM(2,*),ADDCM(*),CHAINE(2,*)
341C --- --- --- INTEGER CAND_M(*),CAND_S(*),IFPEN(*)
342C --- --- --- INTEGER ESHIFT,ILD,NB_N_B, NCONTACT, NCONT, ITAB(*),
343C --- --- --- . IFORM,II_STOK
344C --- --- --- C REAL
345C --- --- --- my_real
346C --- --- --- . GAP,TZINF,MAXBOX,MINBOX,BMINMA(6),DRAD
347C --- --- --- my_real
348C --- --- --- . X(3,*),GAP_S(*),STIFS(*),PENIS(2,*),STIFM(*)
349C --- --- --- C-----------------------------------------------
350C --- --- --- C G l o b a l P a r a m e t e r s
351C --- --- --- C-----------------------------------------------
352C --- --- --- #include "mvsiz_p.inc"
353C --- --- --- C-----------------------------------------------
354C --- --- --- C L o c a l V a r i a b l e s
355C --- --- --- C-----------------------------------------------
356C --- --- --- INTEGER I_ADD_MAX,ISZNSNR
357C --- --- --- PARAMETER (I_ADD_MAX = 1001)
358C --- --- --- C
359C --- --- --- INTEGER I, J, N1, N2, I_ADD, MAXSIZ,JJ,
360C --- --- --- . ADD(2,I_ADD_MAX), N
361C --- --- --- my_real
362C --- --- --- . XYZM(6,I_ADD_MAX-1), MARGE, AAA
363C --- --- --- INTEGER NB_OLD(2,I_ADD_MAX+1)
364C --- --- --- INTEGER NBX,NBY,NBZ
365C --- --- --- INTEGER (KIND=8) :: NBX8,NBY8,NBZ8,RES8,LVOXEL8
366C --- --- --- C-----------------------------------------------
367C --- --- --- C definition from TRI7BOX module
368C --- --- --- C-----------------------------------------------
369C --- --- --- C
370C --- --- --- C-----------------------------------------------
371C --- --- --- C S o u r c e L i n e s
372C --- --- --- C-----------------------------------------------
373C --- --- --- C-----------------------------------------------
374C --- --- --- C
375C --- --- --- C----- TRI PAR BOITES
376C --- --- --- C
377C --- --- --- C-----------------------------------------------
378C --- --- --- C SI ON A PAS ASSEZ DE MEMOIRE POUR LES PILES ON RECOMMENCE LE TRI
379C --- --- --- C EN INCREMENTANT LE NB_N_B (NOMBRE DE NOEUDS PAR BOITE FINIE)
380C --- --- --- C
381C --- --- --- C POINTEUR NOM TAILLE
382C --- --- --- C P1........Elt Bas Pile NRTM
383C --- --- --- C P2........Elt PILE 3*NRTM
384C --- --- --- C P21.......Elt Bas Pile NRTS
385C --- --- --- C P22.......Elt PILE 3*NRTS
386C --- --- --- MAXSIZ = 3*(MAX(NRTM,NRTS+NRTSR)+100)
387C --- --- --- C
388C --- --- --- C-----INITIALISATION DES ADRESSES ET X,Y,Z
389C --- --- --- C
390C --- --- --- C ADDE ADDN X Y Z
391C --- --- --- C 1 1 XMIN YMIN ZMIN
392C --- --- --- C 1 1 XMAX YMAX ZMAX
393C --- --- --- C
394C --- --- --- ADD(1,1) = 0
395C --- --- --- ADD(2,1) = 0
396C --- --- --- ADD(1,2) = 0
397C --- --- --- ADD(2,2) = 0
398C --- --- --- I_ADD = 1
399C --- --- --- XYZM(1,I_ADD) = BMINMA(4)
400C --- --- --- XYZM(2,I_ADD) = BMINMA(5)
401C --- --- --- XYZM(3,I_ADD) = BMINMA(6)
402C --- --- --- XYZM(4,I_ADD) = BMINMA(1)
403C --- --- --- XYZM(5,I_ADD) = BMINMA(2)
404C --- --- --- XYZM(6,I_ADD) = BMINMA(3)
405C --- --- --- I_MEM = 0
406C --- --- --- C
407C --- --- --- IF (IFORM /= 2) THEN
408C --- --- --- ISZNSNR = 0
409C --- --- --- DO I=1,NRTM
410C --- --- --- ADDCM(I)=0
411C --- --- --- ENDDO
412C --- --- --- ELSE
413C --- --- --- ISZNSNR = NRTSR
414C --- --- --- ENDIF
415C --- --- --- C
416C --- --- --- C
417C --- --- --- C-----DEBUT DE LA PHASE DE TRI
418C --- --- --- C
419C --- --- --- GOTO 50 !commenter pour debtrancher le tri voxel
420C --- --- --- print *, "bucket search"
421C --- --- --- CALL I11TRI(
422C --- --- --- 1 ADD ,
423C --- --- --- 2 IRECTS ,IRECTM ,X ,NRTM ,NRTSR ,
424C --- --- --- 3 XYZM ,I_ADD ,MAXSIZ ,II_STOK ,CAND_S ,
425C --- --- --- 4 CAND_M ,NCONTACT,NOINT ,TZINF ,MAXBOX ,
426C --- --- --- 5 MINBOX ,I_MEM ,NB_N_B ,I_ADD_MAX,ESHIFT ,
427C --- --- --- 6 ADDCM ,CHAINE ,NRTS ,ITAB ,NB_OLD ,
428C --- --- --- 7 STIFS ,STIFM ,IAUTO ,NIN ,IFPEN ,
429C --- --- --- 8 IFORM)
430C --- --- ---
431C --- --- --- !cts on debranche la routine de tri par voxel pour le moment
432C --- --- --- GOTO 100
433C --- --- ---
434C --- --- --- 50 CONTINUE
435C --- --- ---
436C --- --- --- MARGE = TZINF - MAX(GAP,DRAD)
437C --- --- ---
438C --- --- --- AAA = SQRT(NMN /
439C --- --- --- . ((BMINMA(1)-BMINMA(4))*(BMINMA(2)-BMINMA(5))
440C --- --- --- . +(BMINMA(2)-BMINMA(5))*(BMINMA(3)-BMINMA(6))
441C --- --- --- . +(BMINMA(3)-BMINMA(6))*(BMINMA(1)-BMINMA(4))))
442C --- --- ---
443C --- --- --- AAA = 0.75*AAA
444C --- --- ---
445C --- --- --- NBX = NINT(AAA*(BMINMA(1)-BMINMA(4)))
446C --- --- --- NBY = NINT(AAA*(BMINMA(2)-BMINMA(5)))
447C --- --- --- NBZ = NINT(AAA*(BMINMA(3)-BMINMA(6)))
448C --- --- --- NBX = MAX(NBX,1)
449C --- --- --- NBY = MAX(NBY,1)
450C --- --- --- NBZ = MAX(NBZ,1)
451C --- --- ---
452C --- --- --- NBX8=NBX
453C --- --- --- NBY8=NBY
454C --- --- --- NBZ8=NBZ
455C --- --- --- RES8=(NBX8+2)*(NBY8+2)*(NBZ8+2)
456C --- --- --- LVOXEL8 = LVOXEL
457C --- --- ---
458C --- --- --- IF(RES8 > LVOXEL8)THEN
459C --- --- --- AAA = LVOXEL
460C --- --- --- AAA = AAA/((NBX8+2)*(NBY8+2)*(NBZ8+2))
461C --- --- --- AAA = AAA**(THIRD)
462C --- --- --- NBX = INT((NBX+2)*AAA)-2
463C --- --- --- NBY = INT((NBY+2)*AAA)-2
464C --- --- --- NBZ = INT((NBZ+2)*AAA)-2
465C --- --- --- NBX = MAX(NBX,1)
466C --- --- --- NBY = MAX(NBY,1)
467C --- --- --- NBZ = MAX(NBZ,1)
468C --- --- --- NBX8 = NBX
469C --- --- --- NBY8 = NBY
470C --- --- --- NBZ8 = NBZ
471C --- --- --- RES8=(NBX8+2)*(NBY8+2)*(NBZ8+2)
472C --- --- --- END IF
473C --- --- ---
474C --- --- --- IF(RES8 > LVOXEL8) THEN
475C --- --- --- stop 678
476C --- --- --- END IF
477C --- --- ---
478C --- --- --- DO I=INIVOXEL,(NBX+2)*(NBY+2)*(NBZ+2)
479C --- --- --- VOXEL1(I)=0
480C --- --- --- ENDDO
481C --- --- --- INIVOXEL = MAX(INIVOXEL,(NBX+2)*(NBY+2)*(NBZ+2)+1)
482C --- --- ---
483C --- --- --- !print *, "voxel search"
484C --- --- ---
485C --- --- --- CALL I11TRIVOX(
486C --- --- --- 1 IRECTS ,IRECTM ,X ,NRTM ,NRTSR ,
487C --- --- --- 2 XYZM ,II_STOK ,CAND_S ,CAND_M ,NCONTACT,
488C --- --- --- 3 NOINT ,TZINF ,I_MEM ,ESHIFT ,ADDCM ,
489C --- --- --- 4 CHAINE ,NRTS ,ITAB ,STIFS ,STIFM ,
490C --- --- --- 5 IAUTO ,VOXEL1 ,NBX ,NBY ,NBZ ,
491C --- --- --- 6 ITASK ,IFPEN ,IFORM )
492C --- --- ---
493C --- --- --- CALL MY_BARRIER
494C --- --- ---
495C --- --- --- 100 CONTINUE
496C --- --- ---
497C --- --- --- C
498C --- --- --- C I_MEM = 1 ==> N/A
499C --- --- --- C I_MEM = 2 ==> PAS ASSEZ DE MEMOIRE CANDIDATS
500C --- --- --- C I_MEM = 3 ==> N/A
501C --- --- ---
502C --- --- --- IF (I_MEM == 2) RETURN
503C --- --- --- IF(I_MEM==1)THEN
504C --- --- --- NB_N_B = NB_N_B + 1
505C --- --- --- IF ( NB_N_B > MAX(NRTM,NRTS)) THEN
506C --- --- --- CALL ANCMSG(MSGID=85,ANMODE=ANINFO,
507C --- --- --- . I1=NOINT)
508C --- --- --- CALL ARRET(2)
509C --- --- --- ENDIF
510C --- --- --- ILD = 1
511C --- --- --- ELSEIF(I_MEM==2) THEN
512C --- --- --- IF(DEBUG(1)>=1) THEN
513C --- --- --- IWARN = IWARN+1
514C --- --- --- #include "lockon.inc"
515C --- --- --- WRITE(ISTDO,*)' **WARNING INTERFACE/MEMORY'
516C --- --- --- WRITE(IOUT,*)' **WARNING INTERFACE NB:',NOINT
517C --- --- --- WRITE(IOUT,*)' TOO MANY POSSIBLE IMPACTS'
518C --- --- --- WRITE(IOUT,*)' SIZE OF INFLUENCE ZONE IS'
519C --- --- --- WRITE(IOUT,*)' MULTIPLIED BY 0.75'
520C --- --- --- #include "lockoff.inc"
521C --- --- --- ENDIF
522C --- --- --- TZINF = THREE_OVER_4*TZINF
523C --- --- --- C ne pas dimunuer la taille des boite
524C --- --- --- C MINBOX= THREE_OVER_4*MINBOX
525C --- --- --- C MAXBOX= THREE_OVER_4*MAXBOX
526C --- --- --- IF( TZINF<=MAX(GAP,DRAD) ) THEN
527C --- --- --- CALL ANCMSG(MSGID=98,ANMODE=ANINFO,
528C --- --- --- . I1=NOINT,C1='(I11BUCE)')
529C --- --- --- CALL ARRET(2)
530C --- --- --- ENDIF
531C --- --- --- ILD = 1
532C --- --- --- ELSEIF(I_MEM==3)THEN
533C --- --- --- NB_N_B = NB_N_B + 1
534C --- --- --- IF ( NB_N_B > MAX(NRTM,NRTS)) THEN
535C --- --- --- CALL ANCMSG(MSGID=99,ANMODE=ANINFO,
536C --- --- --- . I1=NOINT,C1='(I11BUCE)')
537C --- --- --- CALL ARRET(2)
538C --- --- --- ENDIF
539C --- --- --- ILD = 1
540C --- --- --- ENDIF
541C --- --- --- C
542C --- --- --- RETURN
543C --- --- --- END
i11buce_vox
subroutine i11buce_vox(x, irects, irectm, nrts, nmn, nrtm, nsn, cand_m, cand_s, maxgap, noint, ii_stok, tzinf, maxbox, minbox, nb_n_b, eshift, ild, bminma, ncontact, addcm, chaine, nin, itab, nrtsr, ncont, gap_s, stifs, penis, igap, stifm, iauto, i_mem, itask, iform, ifpen, drad, gap_m, gap_s_l, gap_m_l, gapmin, bgapsmx, gap, flagremnode, kremnode, remnode, dgapload)
Definition i11buce.F:48
i11trivox
subroutine i11trivox(irects, irectm, x, nrtm, nrtsr, xyzm, ii_stok, cand_s, cand_m, nsn4, noint, tzinf, i_mem, eshift, addcm, chaine, nrts, itab, stfs, stfm, iauto, voxel, nbx, nby, nbz, itask, ifpen, iform, gapmin, drad, marge, gap_s, gap_m, gap_s_l, gap_m_l, bgapsmx, igap, gap, flagremnode, kremnode, remnode, dgapload)
Definition i11trivox.F:46
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
tri7box::voxel1
integer, dimension(lvoxel) voxel1
Definition tri7box.F:53
tri7box::inivoxel
integer inivoxel
Definition tri7box.F:53
tri7box::lvoxel
integer lvoxel
Definition tri7box.F:51
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
arret
subroutine arret(nn)
Definition arret.F:87
my_barrier
subroutine my_barrier
Definition machine.F:31