OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sfac_asm_master_m.F
Go to the documentation of this file.
1C
2C This file is part of MUMPS 5.5.1, released
3C on Tue Jul 12 13:17:24 UTC 2022
4C
5C
6C Copyright 1991-2022 CERFACS, CNRS, ENS Lyon, INP Toulouse, Inria,
7C Mumps Technologies, University of Bordeaux.
8C
9C This version of MUMPS is provided to you free of charge. It is
10C released under the CeCILL-C license
11C (see doc/CeCILL-C_V1-en.txt, doc/CeCILL-C_V1-fr.txt, and
12C https://cecill.info/licences/Licence_CeCILL-C_V1-en.html)
13C
14 MODULE smumps_fac_asm_master_m
15 CONTAINS
16 SUBROUTINE smumps_fac_asm_niv1( COMM_LOAD, ASS_IRECV,
17 & N, INODE, IW, LIW, A, LA, INFO, ND,
18 & FILS, FRERE, DAD, MAXFRW, root,
19 & OPASSW, OPELIW, PTRIST, PTLUST, PTRFAC, PTRAST,
20 & STEP, PIMASTER, PAMASTER,PTRARW,
21 & PTRAIW, ITLOC, RHS_MUMPS, NSTEPS, SON_LEVEL2,
22 & COMP, LRLU, IPTRLU, IWPOS, IWPOSCB, POSFAC, LRLUS, LRLUSM,
23 & ICNTL, KEEP,KEEP8,DKEEP,INTARR,LINTARR,DBLARR,LDBLARR,
24 &
25 & NSTK_S,PROCNODE_STEPS, SLAVEF, COMM,MYID,
26 & BUFR, LBUFR, LBUFR_BYTES, NBFIN, IPOOL, LPOOL, LEAF,
27 & PERM,
28 & ISTEP_TO_INIV2, TAB_POS_IN_PERE, JOBASS, ETATASS
29 & , LRGROUPS
30 & , MUMPS_TPS_ARR, SMUMPS_TPS_ARR, L0_OMP_MAPPING
31 & )
32!$ USE OMP_LIB
33 USE mumps_tps_m
34 USE smumps_tps_m
35 USE smumps_dynamic_memory_m, ONLY : smumps_dm_set_ptr,
36 & smumps_dm_free_block
37 USE mumps_build_sort_index_m
38 USE smumps_buf
39 USE smumps_load
40 USE smumps_lr_core, ONLY : is_front_blr_candidate,
41 & smumps_blr_asm_niv1
42 USE smumps_lr_data_m, ONLY : smumps_blr_init_front,
43 & smumps_blr_save_nfs4father
44 USE smumps_struc_def, ONLY : smumps_root_struc
45 USE smumps_ana_lr, ONLY : get_cut
46 USE smumps_lr_core, ONLY : max_cluster
47 USE mumps_lr_common, ONLY : compute_blr_vcs
48 IMPLICIT NONE
49 TYPE (SMUMPS_ROOT_STRUC) :: root
50 INTEGER COMM_LOAD, ASS_IRECV
51 INTEGER IZERO
52 parameter(izero=0)
53 INTEGER N, NSTEPS
54 INTEGER(8) LA, LRLU, LRLUS, LRLUSM, IPTRLU, POSFAC
55 INTEGER KEEP(500), ICNTL(60)
56 INTEGER(8) KEEP8(150)
57 REAL DKEEP(230)
58 INTEGER, INTENT(INOUT) :: INFO(2)
59 INTEGER INODE,MAXFRW,
60 & iwposcb, comp
61 INTEGER, TARGET :: IWPOS, LIW
62 TYPE (MUMPS_TPS_T), TARGET, OPTIONAL :: MUMPS_TPS_ARR(:)
63 TYPE (SMUMPS_TPS_T), TARGET, OPTIONAL :: SMUMPS_TPS_ARR(:)
64 INTEGER, INTENT(IN), OPTIONAL :: L0_OMP_MAPPING(:)
65 INTEGER IDUMMY(1)
66 INTEGER, PARAMETER :: LIDUMMY = 1
67 INTEGER, TARGET :: IW(LIW)
68 INTEGER(8), INTENT(IN) :: PTRARW(N), PTRAIW(N)
69 INTEGER ITLOC(N+KEEP(253)),
70 & nd(keep(28)), perm(n),
71 & fils(n), frere(keep(28)), dad(keep(28)),
72 & ptrist(keep(28)), ptlust(keep(28)),
73 & step(n), pimaster(keep(28))
74 REAL :: RHS_MUMPS(KEEP(255))
75 INTEGER(8) :: PTRFAC(KEEP(28)), PTRAST(KEEP(28)),
76 & pamaster(keep(28))
77 INTEGER COMM, NBFIN, SLAVEF, MYID
78 INTEGER ISTEP_TO_INIV2(KEEP(71)),
79 & tab_pos_in_pere(slavef+2,max(1,keep(56)))
80 INTEGER JOBASS,ETATASS
81 LOGICAL SON_LEVEL2
82 REAL, TARGET :: A(LA)
83 INTEGER, INTENT(IN) :: LRGROUPS(N)
84 DOUBLE PRECISION OPASSW, OPELIW
85 INTEGER(8), INTENT(IN) :: LINTARR, LDBLARR
86 REAL DBLARR(LDBLARR)
87 INTEGER INTARR(LINTARR)
88 INTEGER LPOOL, LEAF
89 INTEGER LBUFR, LBUFR_BYTES
90 INTEGER IPOOL( LPOOL )
91 INTEGER NSTK_S(KEEP(28))
92 INTEGER PROCNODE_STEPS(KEEP(28))
93 INTEGER BUFR( LBUFR )
94 LOGICAL PACKED_CB, IS_CB_LR
95 INTEGER, EXTERNAL :: MUMPS_TYPENODE
96 INTEGER, EXTERNAL :: MUMPS_PROCNODE
97 include 'mpif.h'
98 INTEGER :: IERR
99 INTEGER :: STATUS(MPI_STATUS_SIZE)
100!$ INTEGER :: NOMP
101 include 'mumps_headers.h'
102 INTEGER LP, HS, HF
103 LOGICAL LPOK
104 INTEGER NBPANELS_L, NBPANELS_U
105 INTEGER IN,NUMSTK,NASS,ISON,IFSON,NASS1,IELL
106 INTEGER NFS4FATHER, ESTIM_NFS4FATHER_ATSON
107 INTEGER IFATH
108 INTEGER(8) NFRONT8, LAELL8, LAELL_REQ8, ITMP8, KEEP8TMPCOPY
109 INTEGER NFRONT,NFRONT_EFF,ISTCHK,ISTCHK_CB_RIGHT,LSTK,LREQ
110 INTEGER LREQ_OOC
111 INTEGER :: SON_XXS, SON_XXLR, SON_XXG
112 INTEGER(8) LSTK8, SIZFR8
113 LOGICAL :: IS_DYNAMIC_CB
114 INTEGER(8) :: DYN_SIZE
115 INTEGER SIZFI, NCB
116 INTEGER NCOLS, NROWS, LDA_SON
117 INTEGER NELIM, IORG, IBROT
118#if ! defined(ZERO_TRIANGLE)
119 INTEGER(8) :: NUMROWS, JJ3
120#endif
121 INTEGER :: TOPDIAG
122!$ INTEGER :: CHUNK
123!$ INTEGER(8) :: CHUNK8
124 INTEGER JPOS,ICT11
125 INTEGER IJROW,NBCOL,NUMORG,IOLDPS
126 INTEGER(8) IACHK, POSELT, LAPOS2, IACHK_ini
127 INTEGER(8) APOS, APOS2, APOS3, POSEL1, ICT12
128 INTEGER(8) :: JJ2, ICT13
129 INTEGER(8) :: JK8, J18, J28, J38, J48, JJ8
130 INTEGER(8) :: AINPUT8
131 INTEGER :: K1, K2, K3, KK, KK1
132 INTEGER :: J253
133 INTEGER NSLAVES, NSLSON, NPIVS,NPIV_ANA,NPIV
134 INTEGER PTRCOL, ISLAVE, PDEST,LEVEL
135 INTEGER ISON_IN_PLACE
136 LOGICAL SKIP_TOP_STACK
137 INTEGER ISON_TOP
138 INTEGER(8) SIZE_ISON_TOP8, DYN_SIZE_ISON_TOP8
139 LOGICAL RESET_TO_ZERO, RISK_OF_SAME_POS,
140 & risk_of_same_pos_this_line
141!$ LOGICAL OMP_PARALLEL_FLAG
142 LOGICAL LEVEL1, NIV1
143 INTEGER TROW_SIZE
144 INTEGER INDX, FIRST_INDEX, SHIFT_INDEX
145 INTEGER PARPIV_T1
146 LOGICAL BLOCKING, SET_IRECV, MESSAGE_RECEIVED
147 INTEGER LRSTATUS
148 LOGICAL COMPRESS_PANEL, LR_ACTIVATED, COMPRESS_CB,
149 & oocwrite_compatible_with_blr
150 INTEGER :: ITHREAD
151 INTEGER, POINTER :: SON_IWPOS, SON_LIW
152 INTEGER, POINTER, DIMENSION(:) :: SON_IW
153 REAL, POINTER, DIMENSION(:) :: SON_A
154 INTEGER NCBSON
155 LOGICAL SAME_PROC
156 INTEGER, POINTER, DIMENSION(:) :: BEGS_BLR
157 INTEGER :: NB_BLR, NPARTSCB, NPARTSASS, MAXI_CLUSTER,
158 & ibcksz2, minsize
159 INTRINSIC real
160 REAL ZERO
161 parameter( zero = 0.0e0 )
162 INTEGER NELT, LPTRAR
163 EXTERNAL mumps_inssarbr
164 LOGICAL MUMPS_INSSARBR
165 LOGICAL SSARBR
166 DOUBLE PRECISION FLOP1,FLOP1_EFF
167 EXTERNAL mumps_in_or_root_ssarbr
168 LOGICAL MUMPS_IN_OR_ROOT_SSARBR
169!$ NOMP = OMP_GET_MAX_THREADS()
170 lp = icntl(1)
171 lpok = ((lp.GT.0).AND.(icntl(4).GE.1))
172 nelt = 1
173 lptrar = n
174 nfs4father = -1
175 packed_cb = .false.
176 is_cb_lr = .false.
177 in = inode
178 level = mumps_typenode(procnode_steps(step(inode)),keep(199))
179 IF (level.NE.1) THEN
180 WRITE(*,*) 'INTERNAL ERROR 1 in SMUMPS_FAC_ASM_NIV1 '
181 CALL mumps_abort()
182 END IF
183 nslaves = 0
184 hf = 6 + nslaves + keep(ixsz)
185 IF (jobass.EQ.0) THEN
186 etatass= 0
187 ELSE
188 etatass= 2
189 ioldps = ptlust(step(inode))
190 nfront = iw(ioldps + keep(ixsz))
191 nass1 = iabs(iw(ioldps + 2 + keep(ixsz)))
192 ict11 = ioldps + hf - 1 + nfront
193 ssarbr=mumps_inssarbr(procnode_steps(step(inode)),
194 & keep(199))
195 numorg = 0
196 DO WHILE (in.GT.0)
197 numorg = numorg + 1
198 in = fils(in)
199 ENDDO
200 numstk = 0
201 ifson = -in
202 ison = ifson
203 IF (ison .NE. 0) THEN
204 DO WHILE (ison .GT. 0)
205 numstk = numstk + 1
206 ison = frere(step(ison))
207 ENDDO
208 ENDIF
209 GOTO 123
210 ENDIF
211 numorg = 0
212 DO WHILE (in.GT.0)
213 numorg = numorg + 1
214 in = fils(in)
215 END DO
216 npiv_ana=numorg
217 nsteps = nsteps + 1
218 numstk = 0
219 nass = 0
220 ifson = -in
221 ison = ifson
222 IF (ison .NE. 0) THEN
223 DO WHILE (ison .GT. 0)
224 numstk = numstk + 1
225 son_iw => iw
226 IF (keep(400).GT.0) THEN
227 IF (present(l0_omp_mapping)) THEN
228 ithread=l0_omp_mapping(step(ison))
229 IF (ithread .NE.0) THEN
230 son_iw=>mumps_tps_arr(ithread)%IW
231 ENDIF
232 ENDIF
233 ENDIF
234 nass = nass + son_iw(pimaster(step(ison))+1+keep(ixsz))
235 ison = frere(step(ison))
236 END DO
237 ENDIF
238 nfront = nd(step(inode)) + nass + keep(253)
239 nass1 = nass + numorg
240 CALL is_front_blr_candidate(inode, 1, nd(step(inode)),
241 & numorg, keep(486),
242 & keep(489), keep(490), keep(491), keep(492),
243 & keep(20), keep(60), dad(step(inode)), keep(38),
244 & lrstatus, n, lrgroups)
245 IF (dad(step(inode)).NE.0) THEN
246 IF ( mumps_procnode(procnode_steps(step(dad(step(inode)))),
247 & keep(199) )
248 & .NE. myid
249 & .AND.
250 & mumps_typenode(procnode_steps(step(dad(step(inode)))),
251 & keep(199))
252 & .EQ.1
253 & ) THEN
254 IF (lrstatus.EQ.1 .OR. lrstatus.EQ.3) THEN
255 lrstatus = lrstatus-1
256 ENDIF
257 ENDIF
258 ENDIF
259 compress_panel = (lrstatus.GE.2)
260 compress_cb = ((lrstatus.EQ.1).OR.
261 & (lrstatus.EQ.3))
262 lr_activated = (lrstatus.GT.0)
263 IF (compress_cb.AND.(.NOT.compress_panel)) THEN
264 compress_panel = .true.
265 lrstatus = 3
266 ENDIF
267 oocwrite_compatible_with_blr =
268 & ( .NOT.lr_activated.OR.(.NOT.compress_panel).OR.
269 & (keep(486).NE.2)
270 & )
271 lreq_ooc = 0
272 IF (keep(201).EQ.1.AND.oocwrite_compatible_with_blr) THEN
273 CALL smumps_ooc_get_pp_sizes(keep(50), nfront, nfront, nass1,
274 & nbpanels_l, nbpanels_u, lreq_ooc)
275 ENDIF
276 lreq = hf + 2 * nfront + lreq_ooc
277 IF ((iwpos + lreq -1) .GT. iwposcb) THEN
278 CALL smumps_compre_new(n, keep,
279 & iw, liw, a, la,
280 & lrlu, iptrlu,
281 & iwpos, iwposcb, ptrist, ptrast,
282 & step, pimaster, pamaster, lrlus,
283 & keep(ixsz), comp, dkeep(97), myid, slavef,
284 & procnode_steps, dad)
285 IF (lrlu .NE. lrlus) THEN
286 IF (lpok) THEN
287 WRITE(lp, * ) 'INTERNAL ERROR 2 after compress '
288 WRITE(lp, * ) 'IN SMUMPS_FAC_ASM_NIV1 '
289 WRITE(lp, * ) 'LRLU,LRLUS=', lrlu,lrlus
290 ENDIF
291 GOTO 270
292 END IF
293 IF ((iwpos + lreq -1) .GT. iwposcb) GOTO 270
294 END IF
295 ioldps = iwpos
296 iwpos = iwpos + lreq
297 ison_top = -9999
298 ison_in_place = -9999
299 size_ison_top8 = 0_8
300 IF (keep(234).NE.0) THEN
301 IF ( iwposcb .NE. liw ) THEN
302 IF ( iwposcb+iw(iwposcb+1+xxi).NE.liw) THEN
303 ison = iw( iwposcb + 1 + xxn )
304 IF ( dad( step( ison ) ) .EQ. inode .AND.
305 & mumps_typenode(procnode_steps(step(ison)),keep(199))
306 & .EQ. 1 )
307 & THEN
308 ison_top = ison
309 CALL mumps_geti8(size_ison_top8,iw(iwposcb + 1 + xxr))
310 CALL mumps_geti8(dyn_size_ison_top8, iw(iwposcb + 1 + xxd))
311 IF (dyn_size_ison_top8 .EQ. 0_8) THEN
312 IF (lrlu .LT. int(nfront,8) * int(nfront,8)) THEN
313 ison_in_place = ison
314 ENDIF
315 ENDIF
316 END IF
317 END IF
318 END IF
319 END IF
320 niv1 = .true.
321 IF (.NOT. present(mumps_tps_arr).AND.
322 & .NOT. present(l0_omp_mapping) ) THEN
323 CALL mumps_build_sort_index(
324 & myid, inode, n, ioldps, hf, lp, lpok,
325 & nfront, nfront_eff, perm, dad,
326 & nass1, nass, numstk, numorg, iwposcb, iwpos,
327 & ifson, step, pimaster, ptrist, ptraiw, iw, liw,
328 & intarr, lintarr, itloc, fils, frere,
329 & son_level2, niv1, keep, keep8, info(1),
330 & ison_in_place,
331 & procnode_steps, slavef, idummy, lidummy )
332 ELSE
333 CALL mumps_build_sort_index(
334 & myid, inode, n, ioldps, hf, lp, lpok,
335 & nfront, nfront_eff, perm, dad,
336 & nass1, nass, numstk, numorg, iwposcb, iwpos,
337 & ifson, step, pimaster, ptrist, ptraiw, iw, liw,
338 & intarr, lintarr, itloc, fils, frere,
339 & son_level2, niv1, keep, keep8, info(1),
340 & ison_in_place,
341 & procnode_steps, slavef, idummy, lidummy
342 & , mumps_tps_arr, l0_omp_mapping )
343 ENDIF
344 IF (info(1).LT.0) GOTO 300
345 IF (nfront_eff.NE.nfront) THEN
346 IF (nfront.GT.nfront_eff) THEN
347 IF(mumps_in_or_root_ssarbr(procnode_steps(step(inode)),
348 & keep(199)))THEN
349 npiv=nass1-(nfront_eff-nd(step(inode)))
350 CALL mumps_get_flops_cost(nd(step(inode))+keep(253),
351 & npiv,npiv,
352 & keep(50),1,flop1)
353 npiv=npiv_ana
354 CALL mumps_get_flops_cost(nd(step(inode))+keep(253),
355 & npiv,npiv,
356 & keep(50),1,flop1_eff)
357 CALL smumps_load_update(0,.false.,flop1-flop1_eff,
358 & keep,keep8)
359 ENDIF
360 iwpos = iwpos - ((2*nfront)-(2*nfront_eff))
361 nfront = nfront_eff
362 lreq = hf + 2 * nfront + lreq_ooc
363 ELSE
364 IF (lpok) THEN
365 WRITE(lp,*) ' INTERNAL ERROR 3 ',
366 & ' IN SMUMPS_FAC_ASM_NIV1 ',
367 & ' NFRONT, NFRONT_EFF = ',
368 & nfront, nfront_eff
369 ENDIF
370 GOTO 270
371 ENDIF
372 ENDIF
373 IF (keep(201).EQ.1.AND.keep(50).NE.1.AND.
374 & oocwrite_compatible_with_blr) THEN
375 CALL smumps_ooc_pp_set_ptr(keep(50),
376 & nbpanels_l, nbpanels_u, nass1,
377 & ioldps + hf + 2 * nfront, iw, liw)
378 ENDIF
379 ncb = nfront - nass1
380 maxfrw = max0(maxfrw, nfront)
381 ict11 = ioldps + hf - 1 + nfront
382 CALL smumps_set_parpivt1 ( inode, nfront, nass1, keep,
383 & lr_activated, parpiv_t1)
384 nfront8=int(nfront,8)
385 laell8 = nfront8 * nfront8
386 IF(parpiv_t1.NE.0) THEN
387 laell8 = laell8+int(nass1,8)
388 ENDIF
389 laell_req8 = laell8
390 IF ( ison_in_place > 0 ) THEN
391 laell_req8 = laell8 - size_ison_top8
392 ENDIF
393 skip_top_stack = (ison_in_place.GT.0)
394 CALL smumps_get_size_needed
395 & (0, laell_req8, skip_top_stack,
396 & keep(1), keep8(1),
397 & n,iw,liw,a,la,
398 & lrlu,iptrlu,iwpos,iwposcb,
399 & ptrist,ptrast,
400 & step, pimaster,pamaster,lrlus,
401 & keep(ixsz), comp, dkeep(97), myid,
402 & slavef, procnode_steps, dad,
403 & info(1), info(2))
404 IF (info(1).LT.0) GOTO 490
405 lrlu = lrlu - laell8
406 lrlus = lrlus - laell8 + size_ison_top8
407 lrlusm = min( lrlus, lrlusm )
408 itmp8 = laell8 - size_ison_top8
409 IF (keep(405).EQ.0) THEN
410 keep8(69) = keep8(69) + itmp8
411 keep8(68) = max(keep8(69), keep8(68))
412 ELSE
413!$OMP ATOMIC CAPTURE
414 keep8(69) = keep8(69) + itmp8
415 keep8tmpcopy = keep8(69)
416!$OMP END ATOMIC
417!$OMP ATOMIC UPDATE
418 keep8(68) = max(keep8(68), keep8tmpcopy)
419!$OMP END ATOMIC
420 ENDIF
421 poselt = posfac
422 posfac = posfac + laell8
423 ssarbr=mumps_inssarbr(procnode_steps(step(inode)),keep(199))
424 CALL smumps_load_mem_update(ssarbr,.false.,
425 & la-lrlus,
426 & 0_8,
427 & laell8-size_ison_top8,
428 & keep,keep8,
429 & lrlus)
430 IF (keep(405).EQ.0) keep(429)= keep(429)+1
431#if defined(ZERO_TRIANGLE)
432 lapos2 = min(poselt + laell8 - 1_8, iptrlu)
433 a(poselt:lapos2) = zero
434#else
435 IF ( keep(50) .eq. 0 .OR. nfront .LT. keep(63) ) THEN
436 lapos2 = min(poselt + laell8 - 1_8, iptrlu)
437!$ CHUNK8=int(KEEP(361),8)
438!$OMP PARALLEL DO PRIVATE(JJ8) SCHEDULE(STATIC, CHUNK8)
439!$OMP& IF (LAPOS2 - POSELT + 1_8 > int(KEEP(361),8) .AND. NOMP .GT. 1)
440#if defined(__ve__)
441!NEC$ IVDEP
442#endif
443 DO jj8 = poselt, lapos2
444 a( jj8 ) = zero
445 ENDDO
446!$OMP END PARALLEL DO
447 ELSE
448 topdiag = max(keep(7), keep(8), keep(218))-1
449 IF (lr_activated) THEN
450 NULLIFY(begs_blr)
451 CALL get_cut(iw(ioldps+hf:ioldps+hf+nfront-1), nass1,
452 & ncb, lrgroups, npartscb,
453 & npartsass, begs_blr)
454 nb_blr = npartsass + npartscb
455 CALL max_cluster(begs_blr,nb_blr,maxi_cluster)
456 DEALLOCATE(begs_blr)
457 CALL compute_blr_vcs(keep(472), ibcksz2, keep(488), nass1)
458 minsize = int(ibcksz2 / 2)
459 topdiag = max(2*minsize + maxi_cluster-1,topdiag)
460 ENDIF
461 IF (etatass.EQ.1) THEN
462 IF (keep(234).NE.0) THEN
463 WRITE(*,*)
464 & .EQ."Internal error: ETATASS1 and IN-PLACE ACTIVATED"
465 CALL mumps_abort()
466 ENDIF
467#if defined(__ve__)
468!NEC$ IVDEP
469#endif
470!$ CHUNK = max( KEEP(360)/2, (NFRONT+NOMP-1) / NOMP )
471!$OMP PARALLEL DO PRIVATE(APOS, JJ3) SCHEDULE( STATIC, CHUNK )
472!$OMP& IF (NFRONT8 - 1_8 > KEEP(360))
473 DO jj8 = 0_8, nfront8 - 1_8
474 jj3 = min(jj8+topdiag,int(nass1-1,8))
475 apos = poselt + jj8 * nfront8
476 a(apos:apos+jj3) = zero
477 END DO
478!$OMP END PARALLEL DO
479 ELSE
480 numrows = min(nfront8, (iptrlu-poselt) / nfront8 )
481!$ CHUNK = max(KEEP(360)/2,
482!$ & ( ((int(numrows)+nomp-1) / nomp + 2) / 3) )
483#if defined(__ve__)
484!NEC$ IVDEP
485#endif
486!$OMP PARALLEL DO PRIVATE(APOS,JJ3) SCHEDULE(STATIC, CHUNK )
487!$OMP& IF (int(NUMROWS - 1) .GT. KEEP(360) .AND. NOMP .GT. 1)
488 DO jj8 = 0_8, numrows - 1_8
489 apos = poselt + jj8 * nfront8
490 jj3 = min( nfront8 - 1_8, jj8 + topdiag )
491 a(apos:apos + jj3) = zero
492 ENDDO
493!$OMP END PARALLEL DO
494 IF( numrows .LT. nfront8 ) THEN
495 apos = poselt + nfront8*numrows
496 a(apos : min(iptrlu,apos+numrows)) = zero
497 ENDIF
498 ENDIF
499 END IF
500#endif
501 ptrast(step(inode)) = poselt
502 ptrfac(step(inode)) = poselt
503 ptlust(step(inode)) = ioldps
504 iw(ioldps+xxi) = lreq
505 CALL mumps_storei8(laell8,iw(ioldps+xxr))
506 CALL mumps_storei8(0_8,iw(ioldps+xxd))
507 iw(ioldps+xxs) = -9999
508 iw(ioldps+xxn) = -99999
509 iw(ioldps+xxp) = -99999
510 iw(ioldps+xxa) = -99999
511 iw(ioldps+xxf) = -99999
512 iw(ioldps+xxlr) = lrstatus
513 iw(ioldps + keep(ixsz)) = nfront
514 iw(ioldps + keep(ixsz)+ 1) = 0
515 iw(ioldps + keep(ixsz) + 2) = -nass1
516 iw(ioldps + keep(ixsz) + 3) = -nass1
517 iw(ioldps + keep(ixsz) + 4) = step(inode)
518 iw(ioldps + keep(ixsz) + 5) = nslaves
519 IF (lr_activated.AND.
520 & (keep(480).NE.0
521 & .OR.
522 & (
523 & (keep(486).EQ.2)
524 & )
525 & .OR.compress_cb
526 & )) THEN
527 CALL smumps_blr_init_front (iw(ioldps+xxf), info,
528 & mtk405=keep(405))
529 IF (info(1).LT.0) GOTO 500
530 ENDIF
531 estim_nfs4father_atson = -9999
532 IF (keep(219).NE.0.AND.keep(50).EQ.2) THEN
533 ifath = dad( step( inode) )
534 IF (ifath.NE.0) THEN
535 IF (compress_cb.AND.
536 & mumps_typenode(procnode_steps(step(ifath)),keep(199))
537 & .EQ. 2 ) THEN
538 ioldps = ptlust(step(inode))
539 CALL smumps_compute_estim_nfs4father (
540 & n, inode, ifath, fils, perm, keep,
541 & ioldps, hf, iw, liw, nfront, nass1,
542 & estim_nfs4father_atson
543 & )
544 CALL smumps_blr_save_nfs4father ( iw(ioldps+xxf),
545 & estim_nfs4father_atson )
546 IF (info(1).LT.0) GOTO 500
547 ENDIF
548 ENDIF
549 ENDIF
550 123 CONTINUE
551 IF (numstk.NE.0) THEN
552 IF (ison_top > 0) THEN
553 ison = ison_top
554 ELSE
555 ison = ifson
556 ENDIF
557 DO 220 iell = 1, numstk
558 istchk = pimaster(step(ison))
559 son_iw => iw
560 son_liw => liw
561 son_iwpos => iwpos
562 son_a => a
563 ithread = 0
564 IF (keep(400).GT.0) THEN
565 IF (present(l0_omp_mapping)) THEN
566 ithread=l0_omp_mapping(step(ison))
567 IF (ithread .NE.0) THEN
568 son_liw => mumps_tps_arr(ithread)%LIW
569 son_iw => mumps_tps_arr(ithread)%IW
570 son_iwpos => mumps_tps_arr(ithread)%IWPOS
571 son_a => smumps_tps_arr(ithread)%A
572 ENDIF
573 ENDIF
574 ENDIF
575 lstk = son_iw(istchk + keep(ixsz))
576 lstk8 = int(lstk,8)
577 nelim = son_iw(istchk + keep(ixsz) + 1)
578 npivs = son_iw(istchk + keep(ixsz) + 3)
579 IF ( npivs .LT. 0 ) npivs = 0
580 nslson = son_iw(istchk + keep(ixsz) + 5)
581 hs = 6 + keep(ixsz) + nslson
582 ncols = npivs + lstk
583 same_proc = (istchk.LT.son_iwpos)
584 IF ( same_proc ) THEN
585 istchk_cb_right = ptrist(step(ison))
586 ELSE
587 istchk_cb_right = istchk
588 ENDIF
589 son_xxs = son_iw(istchk_cb_right+xxs)
590 son_xxlr = son_iw(istchk_cb_right+xxlr)
591 son_xxg = son_iw(istchk_cb_right+xxg)
592 packed_cb = ( son_xxs .EQ. s_cb1comp )
593 is_cb_lr = ( son_xxlr.EQ.1 .OR. son_xxlr.EQ.3 )
594 & .AND. (keep(489).EQ.1.OR.keep(489).EQ.3)
595 level1 = nslson.EQ.0
596 IF (.NOT.same_proc) THEN
597 nrows = son_iw( istchk + keep(ixsz) + 2)
598 ELSE
599 nrows = ncols
600 ENDIF
601 sizfi = hs + nrows + ncols
602 k1 = istchk + hs + nrows + npivs
603 IF ( .NOT. level1 .AND. nelim.EQ.0 ) GOTO 205
604 IF (level1 .AND. .NOT. is_cb_lr) THEN
605 k2 = k1 + lstk - 1
606 IF (packed_cb) THEN
607 sizfr8 = (lstk8*(lstk8+1_8)/2_8)
608 ELSE
609 sizfr8 = lstk8*lstk8
610 ENDIF
611 ELSE
612 IF ( keep(50).eq.0 ) THEN
613 sizfr8 = int(nelim,8) * lstk8
614 ELSE
615 IF (packed_cb) THEN
616 sizfr8 = int(nelim,8) * int(nelim+1,8)/2_8
617 ELSE
618 sizfr8 = int(nelim,8) * int(nelim,8)
619 ENDIF
620 END IF
621 k2 = k1 + nelim - 1
622 ENDIF
623 IF (jobass.EQ.0) THEN
624 IF (level1 .AND. .NOT. is_cb_lr) THEN
625 IF (keep(50).EQ.0) THEN
626 opassw = opassw + lstk8*lstk8
627 ELSE
628 opassw = opassw + lstk8*(lstk8+1)/2_8
629 ENDIF
630 ELSE
631 IF (keep(50).EQ.0) THEN
632 opassw = opassw + int(nelim,8)*lstk8
633 ELSE
634 opassw = opassw + int(nelim,8)*int(nelim,8)/2_8
635 ENDIF
636 ENDIF
637 ENDIF
638 CALL mumps_geti8(dyn_size, son_iw(istchk_cb_right+xxd))
639 is_dynamic_cb = dyn_size .GT. 0_8
640 IF ( is_dynamic_cb ) THEN
641 CALL smumps_dm_set_ptr( pamaster(step(ison)), dyn_size,
642 & son_a )
643 iachk = 1_8
644 ELSE
645 iachk = pamaster(step(ison))
646 ENDIF
647 IF (is_cb_lr .AND. level1) THEN
648 posel1 = ptrast(step(inode))
649 CALL smumps_blr_asm_niv1 (a, la,
650 & posel1, nfront, nass1, son_iw(istchk+xxf),
651 & son_iw, son_liw,
652 & lstk, nelim, k1, k1+lstk-1, keep(50),
653 & keep, keep8, opassw)
654 ENDIF
655 IF ( keep(50) .eq. 0 ) THEN
656 posel1 = ptrast(step(inode)) - nfront8
657 IF (nfront .EQ. lstk.AND. ison.EQ.ison_in_place
658 & .AND.iachk + sizfr8 - 1_8 .EQ. posfac - 1_8 ) THEN
659 GOTO 205
660 ENDIF
661 IF (k2.GE.k1) THEN
662 reset_to_zero = (iachk .LT. posfac .AND.
663 & ison.EQ.ison_in_place)
664 risk_of_same_pos = iachk + sizfr8 - 1_8 .EQ. posfac - 1_8
665 & .AND. ison.EQ.ison_in_place
666 risk_of_same_pos_this_line = .false.
667 iachk_ini = iachk
668!$ OMP_PARALLEL_FLAG = (RESET_TO_ZERO.EQV..FALSE.).AND.
669!$ & ((K2-K1).GT.KEEP(360))
670!$omp parallel IF(omp_parallel_flag) private(apos, kk1, jj2,iachk)
671!$OMP& FIRSTPRIVATE(RISK_OF_SAME_POS_THIS_LINE,RESET_TO_ZERO)
672!$OMP DO
673 DO 170 kk = k1, k2
674 apos = posel1 + int(son_iw(kk),8) * int(nfront,8)
675 iachk = iachk_ini + int(kk-k1,8)*int(lstk,8)
676 IF (reset_to_zero) THEN
677 IF (risk_of_same_pos) THEN
678 IF (kk.EQ.k2) THEN
679 risk_of_same_pos_this_line =
680 & (ison .EQ. ison_in_place)
681 & .AND. ( apos + int(son_iw(k1+lstk-1)-1,8).EQ.
682 & iachk+int(lstk-1,8) )
683 ENDIF
684 ENDIF
685 IF ((iachk .GE. posfac).AND.(kk>k1))THEN
686 reset_to_zero =.false.
687 ENDIF
688 IF (risk_of_same_pos_this_line) THEN
689 DO kk1 = 1, lstk
690 jj2 = apos + int(son_iw(k1 + kk1 - 1) - 1,8)
691 IF ( iachk+int(kk1-1,8) .NE. jj2 ) THEN
692 a(jj2) = a(iachk + int(kk1 - 1,8))
693 a(iachk + int(kk1 -1,8)) = zero
694 ENDIF
695 ENDDO
696 ELSE
697#if defined(__ve__)
698!NEC$ IVDEP
699#endif
700 DO kk1 = 1, lstk
701 jj2 = apos + int(son_iw(k1+kk1-1),8) - 1_8
702 a(jj2) = a(iachk + int(kk1 - 1,8))
703 a(iachk + int(kk1 -1,8)) = zero
704 ENDDO
705 ENDIF
706 ELSE
707#if defined(__ve__)
708!NEC$ IVDEP
709#endif
710 DO kk1 = 1, lstk
711 jj2 = apos + int(son_iw(k1+kk1-1),8) - 1_8
712 a(jj2) = a(jj2) + son_a(iachk + int(kk1 - 1,8))
713 ENDDO
714 ENDIF
715 170 CONTINUE
716!$OMP END DO
717!$OMP END PARALLEL
718 END IF
719 ELSE
720 IF (level1 .AND. .NOT. is_cb_lr) THEN
721 lda_son = lstk
722 ELSE
723 lda_son = nelim
724 ENDIF
725 IF (ison .EQ. ison_in_place) THEN
726 CALL smumps_ldlt_asm_niv12_ip(a, la,
727 & ptrast(step( inode )), nfront, nass1,
728 & iachk, lda_son, sizfr8,
729 & son_iw( k1 ), k2 - k1 + 1, nelim, etatass,
730 & packed_cb)
731 ELSE
732 IF (sizfr8 .GT. 0) THEN
733 CALL smumps_ldlt_asm_niv12(a, la, son_a(iachk),
734 & ptrast(step( inode )), nfront, nass1,
735 & lda_son, sizfr8,
736 & son_iw( k1 ), k2 - k1 + 1, nelim, etatass,
737 & packed_cb
738!$ & , KEEP(360)
739 & )
740 ENDIF
741 ENDIF
742 ENDIF
743 205 IF (level1) THEN
744 IF (same_proc) istchk = ptrist(step(ison))
745 IF ((same_proc).AND.etatass.NE.1) THEN
746 IF (keep(50).NE.0) THEN
747 k2 = k1 + lstk - 1
748#if defined(__ve__)
749!NEC$ IVDEP
750#endif
751 DO kk = k1, k2
752 son_iw(kk) = son_iw(kk - nrows)
753 ENDDO
754 ELSE
755 k2 = k1 + lstk - 1
756 k3 = k1 + nelim
757#if defined(__ve__)
758!NEC$ IVDEP
759#endif
760 DO kk = k3, k2
761 son_iw(kk) = son_iw(kk - nrows)
762 ENDDO
763 IF (nelim .NE. 0) THEN
764 k3 = k3 - 1
765#if defined(__ve__)
766!NEC$ IVDEP
767#endif
768 DO kk = k1, k3
769 jpos = son_iw(kk) + ict11
770 son_iw(kk) = iw(jpos)
771 ENDDO
772 ENDIF
773 ENDIF
774 ENDIF
775 IF (etatass.NE.1) THEN
776 IF ( same_proc ) THEN
777 ptrist(step(ison)) = -99999999
778 ELSE
779 pimaster(step( ison )) = -99999999
780 ENDIF
781 IF (ithread .EQ. 0) THEN
782 CALL smumps_free_block_cb_static(
783 & ssarbr, myid, n, istchk,
784 & iw, liw, lrlu, lrlus, iptrlu,
785 & iwposcb, la, keep,keep8,
786 & (ison .EQ. ison_top)
787 & )
788 ELSE
789 CALL mumps_load_disable()
790 CALL smumps_free_block_cb_static(
791 & ssarbr, myid, n, istchk,
792 & mumps_tps_arr(ithread)%IW(1),
793 & mumps_tps_arr(ithread)%LIW,
794 & mumps_tps_arr(ithread)%LRLU,
795 & mumps_tps_arr(ithread)%LRLUS,
796 & mumps_tps_arr(ithread)%IPTRLU,
797 & mumps_tps_arr(ithread)%IWPOSCB,
798 & mumps_tps_arr(ithread)%LA, keep,keep8, .false.
799 & )
800 CALL mumps_load_enable()
801 ENDIF
802 IF (is_dynamic_cb) THEN
803 CALL smumps_dm_free_block(son_xxg,
804 & son_a, sizfr8,
805 & keep(405).EQ.1, keep8 )
806 ENDIF
807 ENDIF
808 ELSE
809 pdest = istchk + 6 + keep(ixsz)
810 ncbson = lstk - nelim
811 ptrcol = istchk + hs + nrows + npivs + nelim
812 DO islave = 0, nslson-1
813 IF (iw(pdest+islave).EQ.myid) THEN
814 CALL mumps_bloc2_get_slave_info(
815 & keep, keep8, ison, step, n, slavef,
816 & istep_to_iniv2, tab_pos_in_pere,
817 & islave+1, ncbson,
818 & nslson,
819 & trow_size, first_index )
820 shift_index = first_index - 1
821 indx = ptrcol + shift_index
822 CALL smumps_maplig( comm_load, ass_irecv,
823 & bufr, lbufr, lbufr_bytes,
824 & inode, ison, nslaves, idummy,
825 & nfront, nass1, nfs4father,
826 & trow_size, iw( indx ),
827 & procnode_steps,
828 & slavef, posfac, iwpos, iwposcb, iptrlu, lrlu,
829 & lrlus, n, iw, liw, a, la,
830 & ptrist, ptlust, ptrfac, ptrast, step,
831 & pimaster, pamaster, nstk_s, comp,
832 & info(1), info(2), myid, comm, perm, ipool, lpool,
833 & leaf, nbfin, icntl, keep, keep8, dkeep, root,
834 & opassw, opeliw, itloc, rhs_mumps, fils, dad,
835 & ptrarw, ptraiw, intarr, dblarr, nd, frere,
836 & lptrar, nelt, iw, iw,
837 &
838 & istep_to_iniv2, tab_pos_in_pere, lrgroups
839 & )
840 IF ( info(1) .LT. 0 ) GOTO 500
841 EXIT
842 ENDIF
843 ENDDO
844 IF (pimaster(step(ison)).GT.0) THEN
845 ierr = -1
846 DO WHILE (ierr.EQ.-1)
847 ptrcol = pimaster(step(ison)) + hs + nrows + npivs + nelim
848 pdest = pimaster(step(ison)) + 6 + keep(ixsz)
849 CALL smumps_buf_send_maplig(
850 & inode, nfront, nass1, nfs4father,
851 & ison, myid,
852 & izero, idummy, iw(ptrcol), ncbson,
853 & comm, ierr, iw(pdest), nslson, slavef,
854 & keep, keep8, step, n,
855 & istep_to_iniv2, tab_pos_in_pere
856 & )
857 IF (ierr.EQ.-1) THEN
858 blocking = .false.
859 set_irecv = .true.
860 message_received = .false.
861 CALL smumps_try_recvtreat( comm_load, ass_irecv,
862 & blocking, set_irecv, message_received,
863 & mpi_any_source, mpi_any_tag,
864 & status,
865 & bufr, lbufr, lbufr_bytes, procnode_steps, posfac,
866 & iwpos, iwposcb, iptrlu,
867 & lrlu, lrlus, n, iw, liw, a, la,
868 & ptrist, ptlust, ptrfac,
869 & ptrast, step, pimaster, pamaster, nstk_s, comp,
870 & info(1), info(2), comm,
871 & perm,
872 & ipool, lpool, leaf,
873 & nbfin, myid, slavef,
874 & root, opassw, opeliw, itloc, rhs_mumps,
875 & fils, dad, ptrarw, ptraiw,
876 & intarr, dblarr, icntl, keep, keep8,dkeep, nd, frere,
877 & lptrar, nelt, iw, iw,
878 & istep_to_iniv2, tab_pos_in_pere, .true., lrgroups )
879 IF ( info(1) .LT. 0 ) GOTO 500
880 ENDIF
881 ENDDO
882 IF (ierr .EQ. -2) GOTO 290
883 IF (ierr .EQ. -3) GOTO 295
884 ENDIF
885 ENDIF
886 ison = frere(step(ison))
887 IF (ison .LE. 0) THEN
888 ison = ifson
889 ENDIF
890 220 CONTINUE
891 END IF
892 IF (etatass.EQ.2) GOTO 500
893 poselt = ptrast(step(inode))
894 ibrot = inode
895 DO 260 iorg = 1, numorg
896 jk8 = ptraiw(ibrot)
897 ainput8 = ptrarw(ibrot)
898 jj8 = jk8 + 1_8
899 j18 = jj8 + 1_8
900 j28 = j18 + intarr(jk8)
901 j38 = j28 + 1
902 j48 = j28 - intarr(jj8)
903 ijrow = intarr(j18)
904 ict12 = poselt + int(ijrow - nfront - 1,8)
905#if defined(__ve__)
906 IF ( keep(265).NE. 0 ) THEN
907!NEC$ IVDEP
908#endif
909 DO jj8 = j18, j28
910 apos2 = ict12 + int(intarr(jj8),8) * nfront8
911 a(apos2) = a(apos2) + dblarr(ainput8)
912 ainput8 = ainput8 + 1_8
913 ENDDO
914#if defined(__ve__)
915 ELSE
916 DO jj8 = j18, j28
917 apos2 = ict12 + int(intarr(jj8),8) * nfront8
918 a(apos2) = a(apos2) + dblarr(ainput8)
919 ainput8 = ainput8 + 1_8
920 ENDDO
921 ENDIF
922#endif
923 IF (j38 .LE. j48) THEN
924 ict13 = poselt + int(ijrow - 1,8) * nfront8
925 nbcol = int(j48 - j38 + 1_8)
926#if defined(__ve__)
927 IF ( keep(265) .NE. 0 ) THEN
928!NEC$ IVDEP
929#endif
930 DO jj8 = 1_8, int(nbcol,8)
931 apos3 = ict13 + int(intarr(j38 + jj8 - 1_8) - 1_8,8)
932 a(apos3) = a(apos3) + dblarr(ainput8 + jj8 - 1_8)
933 ENDDO
934#if defined(__ve__)
935 ELSE
936 DO jj8 = 1_8, int(nbcol,8)
937 apos3 = ict13 + int(intarr(j38 + jj8 - 1_8) - 1_8,8)
938 a(apos3) = a(apos3) + dblarr(ainput8 + jj8 - 1_8)
939 ENDDO
940 ENDIF
941#endif
942 ENDIF
943 IF (keep(50).EQ.0) THEN
944 DO j253=1, keep(253)
945 apos = poselt+
946 & int(ijrow-1,8) * nfront8 +
947 & int(nfront-keep(253)+j253-1,8)
948 a(apos) = a(apos) + rhs_mumps( (j253-1)*keep(254)+ibrot )
949 ENDDO
950 ELSE
951 DO j253=1, keep(253)
952 apos = poselt+
953 & int(nfront-keep(253)+j253-1,8) * nfront8 +
954 & int(ijrow-1,8)
955 a(apos) = a(apos) + rhs_mumps( (j253-1)*keep(254)+ibrot )
956 ENDDO
957 ENDIF
958 ibrot = fils(ibrot)
959 260 CONTINUE
960 IF (parpiv_t1.NE.0.AND.(.NOT.son_level2)) THEN
961 ioldps = ptlust(step(inode))
962 CALL smumps_parpivt1_set_nvschur_max (
963 & n, inode, iw, liw, a, la, keep, perm,
964 & ioldps, poselt,
965 & nfront, nass1, lr_activated, parpiv_t1, nass)
966 ENDIF
967 GOTO 500
968 270 CONTINUE
969 info(1) = -8
970 info(2) = lreq
971 IF (lpok) THEN
972 WRITE( lp, * )
973 &' failure in INTEGER ALLOCATION DURING SMUMPS_FAC_ASM'
974 ENDIF
975 GOTO 490
976 290 CONTINUE
977 IF (LPOK) THEN
978 WRITE( LP, * )
979 & ' FAILURE, SEND BUFFER TOO SMALL DURING SMUMPS_FAC_ASM'
980 ENDIF
981 INFO(1) = -17
982 LREQ = NCBSON + 6+NSLSON+KEEP(IXSZ)
983 INFO(2) = LREQ * KEEP( 34 )
984 GOTO 490
985 295 CONTINUE
986 IF (LPOK) THEN
987 WRITE( LP, * )
988 & ' FAILURE, RECV BUFFER TOO SMALL DURING SMUMPS_FAC_ASM'
989 ENDIF
990 INFO(1) = -20
991 LREQ = NCBSON + 6+NSLSON+KEEP(IXSZ)
992 INFO(2) = LREQ * KEEP( 34 )
993 GOTO 490
994 300 CONTINUE
995.EQ. IF( INFO(1)-13 ) THEN
996 IF (LPOK) THEN
997 WRITE( LP, * )
998 & ' FAILURE IN INTEGER DYNAMIC ALLOCATION DURING SMUMPS_FAC_ASM'
999 ENDIF
1000 INFO(2) = NUMSTK + 1
1001 ENDIF
1002 490 CONTINUE
1003.EQ. IF ( KEEP(405) 0 ) THEN
1004 CALL SMUMPS_BDC_ERROR( MYID, SLAVEF, COMM, KEEP )
1005 ENDIF
1006 500 CONTINUE
1007 RETURN
1008 END SUBROUTINE SMUMPS_FAC_ASM_NIV1
1009 SUBROUTINE SMUMPS_FAC_ASM_NIV2(COMM_LOAD, ASS_IRECV,
1010 & N, INODE, IW, LIW, A, LA, INFO,
1011 & ND, FILS, FRERE, DAD,
1012 & CAND,
1013 & ISTEP_TO_INIV2, TAB_POS_IN_PERE,
1014 & MAXFRW, root,
1015 & OPASSW, OPELIW, PTRIST, PTLUST, PTRFAC,
1016 & PTRAST, STEP, PIMASTER, PAMASTER, PTRARW, NSTK_S,
1017 & PTRAIW, ITLOC, RHS_MUMPS, NSTEPS,
1018 & COMP, LRLU, IPTRLU, IWPOS, IWPOSCB, POSFAC, LRLUS,
1019 & ICNTL, KEEP, KEEP8,DKEEP,INTARR,LINTARR,DBLARR,LDBLARR,
1020 & PROCNODE_STEPS, SLAVEF, COMM,MYID,
1021 & BUFR, LBUFR, LBUFR_BYTES, NBFIN, LEAF, IPOOL, LPOOL,
1022 & PERM, MEM_DISTRIB
1023 & , LRGROUPS
1024 & )
1025!$ USE OMP_LIB
1026 USE MUMPS_BUILD_SORT_INDEX_M
1027 USE SMUMPS_BUF
1028 USE SMUMPS_LOAD
1029 USE SMUMPS_LR_CORE, ONLY : IS_FRONT_BLR_CANDIDATE
1030 USE SMUMPS_STRUC_DEF, ONLY : SMUMPS_ROOT_STRUC
1031 USE SMUMPS_DYNAMIC_MEMORY_M, ONLY : SMUMPS_DM_SET_PTR,
1032 & SMUMPS_DM_IS_DYNAMIC
1033 USE SMUMPS_ANA_LR, ONLY : GET_CUT
1034 USE SMUMPS_LR_CORE, ONLY : MAX_CLUSTER
1035 USE MUMPS_LR_COMMON, ONLY : COMPUTE_BLR_VCS
1036 IMPLICIT NONE
1037 TYPE (SMUMPS_ROOT_STRUC) :: root
1038 INTEGER COMM_LOAD, ASS_IRECV
1039 INTEGER N,LIW,NSTEPS, NBFIN
1040 INTEGER KEEP(500), ICNTL(60)
1041 INTEGER(8) KEEP8(150)
1042 REAL DKEEP(230)
1043 INTEGER(8) :: LRLUS, LRLU, IPTRLU, POSFAC, LA
1044 INTEGER, INTENT(INOUT) :: INFO(2)
1045 INTEGER INODE, MAXFRW, LPOOL, LEAF,
1046 & IWPOS, IWPOSCB, COMP, SLAVEF
1047 REAL, TARGET :: A(LA)
1048 INTEGER, intent(in) :: LRGROUPS(N)
1049 DOUBLE PRECISION OPASSW, OPELIW
1050 INTEGER COMM, MYID, LBUFR, LBUFR_BYTES
1051 INTEGER, DIMENSION(0:SLAVEF - 1) :: MEM_DISTRIB
1052 INTEGER IPOOL(LPOOL)
1053 INTEGER(8) :: PTRAST(KEEP(28))
1054 INTEGER(8) :: PTRFAC(KEEP(28))
1055 INTEGER(8) :: PAMASTER(KEEP(28))
1056 INTEGER(8), INTENT(IN) :: PTRARW(N), PTRAIW(N)
1057 INTEGER IW(LIW), ITLOC(N+KEEP(253)),
1058 & ND(KEEP(28)),
1059 & FILS(N), FRERE(KEEP(28)), DAD (KEEP(28)),
1060 & PTRIST(KEEP(28)), PTLUST(KEEP(28)),
1061 & STEP(N),
1062 & PIMASTER(KEEP(28)),
1063 & NSTK_S(KEEP(28)), PERM(N)
1064 REAL :: RHS_MUMPS(KEEP(255))
1065 INTEGER CAND(SLAVEF+1, max(1,KEEP(56)))
1066 INTEGER ISTEP_TO_INIV2(KEEP(71)),
1067 & TAB_POS_IN_PERE(SLAVEF+2,max(1,KEEP(56)))
1068 INTEGER PROCNODE_STEPS(KEEP(28)), BUFR(LBUFR)
1069 INTEGER(8), INTENT(IN) :: LINTARR,LDBLARR
1070 REAL DBLARR(LDBLARR)
1071 INTEGER INTARR(LINTARR)
1072 INCLUDE 'mpif.h'
1073 INTEGER :: IERR
1074 INTEGER :: STATUS(MPI_STATUS_SIZE)
1075!$ INTEGER :: NOMP
1076 INTEGER LP, HS, HF, HF_OLD, NCBSON, NSLAVES_OLD
1077 LOGICAL LPOK
1078 INTEGER NCBSON_MAX
1079 INTEGER IN,NUMSTK,NASS,ISON,IFSON,NASS1,IELL
1080 INTEGER :: IBC_SOURCE
1081 REAL, DIMENSION(:), POINTER :: SON_A
1082 INTEGER :: MAXWASTEDPROCS
1083 PARAMETER (MAXWASTEDPROCS=1)
1084 INTEGER NFS4FATHER, ESTIM_NFS4FATHER_ATSON
1085 INTEGER IFATH
1086 INTEGER I
1087 INTEGER NFRONT,NFRONT_EFF,ISTCHK,ISTCHK_CB_RIGHT,LSTK,LREQ
1088 INTEGER :: SON_XXS
1089 INTEGER(8) :: LAELL8
1090 INTEGER LREQ_OOC
1091 INTEGER NBPANELS_L, NBPANELS_U
1092 LOGICAL PACKED_CB, IS_CB_LR
1093 INTEGER(8) :: LCB
1094 LOGICAL :: IS_DYNAMIC_CB
1095 INTEGER(8) :: DYN_SIZE
1096 INTEGER NCB
1097 INTEGER MP
1098 INTEGER :: K1, K2, KK, KK1
1099 INTEGER :: J253
1100 INTEGER(8) :: JK8, AINPUT8, J18, J28, J38, J48, JJ8
1101 INTEGER(8) :: LAPOS2, JJ2, JJ3
1102 INTEGER(8) :: ICT13
1103 INTEGER(8) :: IACHK, APOS, APOS2, POSELT, ICT12, POSEL1
1104#if ! defined(ZERO_TRIANGLE)
1105 INTEGER :: TOPDIAG
1106#endif
1107!$ INTEGER :: CHUNK
1108!$ INTEGER(8) :: CHUNK8
1109 INTEGER NELIM,NPIVS,NCOLS,NROWS,
1110 & IORG
1111 INTEGER LDAFS, LDA_SON, IJROW, IBROT
1112 INTEGER NBCOL,NUMORG,IOLDPS, NUMORG_SPLIT
1113 INTEGER NSLAVES, NSLSON
1114 INTEGER NBLIG, PTRCOL, PTRROW, PDEST
1115 INTEGER PDEST1(1)
1116 INTEGER :: ISLAVE
1117 INTEGER TYPESPLIT
1118 INTEGER ISON_IN_PLACE
1119 LOGICAL IS_ofType5or6, SPLIT_MAP_RESTART
1120 INTEGER NMB_OF_CAND, NMB_OF_CAND_ORIG
1121 LOGICAL SAME_PROC, NIV1, SON_LEVEL2
1122 LOGICAL BLOCKING, SET_IRECV, MESSAGE_RECEIVED
1123 INTEGER TROW_SIZE, INDX, FIRST_INDEX, SHIFT_INDEX
1124 INTEGER LRSTATUS
1125 LOGICAL COMPRESS_PANEL, LR_ACTIVATED, COMPRESS_CB,
1126 & OOCWRITE_COMPATIBLE_WITH_BLR
1127 INTEGER IZERO
1128 INTEGER IDUMMY(1)
1129 PARAMETER( IZERO = 0 )
1130 INTEGER MUMPS_PROCNODE, MUMPS_TYPENODE, MUMPS_TYPESPLIT
1131 EXTERNAL MUMPS_PROCNODE, MUMPS_TYPENODE, MUMPS_TYPESPLIT
1132 REAL ZERO
1133 REAL RZERO
1134 PARAMETER( RZERO = 0.0E0 )
1135 PARAMETER( ZERO = 0.0E0 )
1136 INTEGER NELT, LPTRAR
1137 logical :: force_cand
1138 INTEGER ETATASS
1139 INCLUDE 'mumps_headers.h'
1140 INTEGER(8) :: APOSMAX
1141 REAL MAXARR
1142 INTEGER INIV2, SIZE_TMP_SLAVES_LIST, allocok
1143 INTEGER NCB_SPLIT, SIZE_LIST_SPLIT, NBSPLIT
1144 INTEGER, ALLOCATABLE, DIMENSION(:) :: TMP_SLAVES_LIST, COPY_CAND
1145 INTEGER, ALLOCATABLE, DIMENSION(:) :: SONROWS_PER_ROW
1146 INTEGER, POINTER, DIMENSION(:) :: BEGS_BLR
1147 INTEGER :: NB_BLR, NPARTSCB, NPARTSASS, MAXI_CLUSTER,
1148 & IBCKSZ2, MINSIZE
1149!$ NOMP = OMP_GET_MAX_THREADS()
1150 MP = ICNTL(2)
1151 LP = ICNTL(1)
1152.GT..AND..GE. LPOK = ((LP0)(ICNTL(4)1))
1153 IS_ofType5or6 = .FALSE.
1154 PACKED_CB = .FALSE.
1155 ETATASS = 0
1156 IN = INODE
1157 NSTEPS = NSTEPS + 1
1158 KEEP(429) = KEEP(429)+1
1159 NUMORG = 0
1160.GT. DO WHILE (IN0)
1161 NUMORG = NUMORG + 1
1162 IN = FILS(IN)
1163 ENDDO
1164 NUMSTK = 0
1165 NASS = 0
1166 IFSON = -IN
1167 ISON = IFSON
1168 NCBSON_MAX = 0
1169 NELT = 1
1170 LPTRAR = 1
1171.GT. DO WHILE (ISON 0)
1172 NUMSTK = NUMSTK + 1
1173.AND. IF ( KEEP(48)==5
1174 & MUMPS_TYPENODE(PROCNODE_STEPS(STEP(ISON)),
1175.EQ. & KEEP(199)) 1) THEN
1176 NCBSON_MAX =
1177 & max(NCBSON_MAX,IW(PIMASTER(STEP(ISON))+KEEP(IXSZ)))
1178 ENDIF
1179 NASS = NASS + IW(PIMASTER(STEP(ISON)) + 1 + KEEP(IXSZ))
1180 ISON = FRERE(STEP(ISON))
1181 ENDDO
1182 NFRONT = ND(STEP(INODE)) + NASS + KEEP(253)
1183 NASS1 = NASS + NUMORG
1184 NCB = NFRONT - NASS1
1185 CALL IS_FRONT_BLR_CANDIDATE(INODE, 2, NFRONT, NASS1, KEEP(486),
1186 & KEEP(489), KEEP(490), KEEP(491), KEEP(492),
1187 & KEEP(20), KEEP(60), DAD(STEP(INODE)), KEEP(38),
1188 & LRSTATUS, N, LRGROUPS)
1189.GE. COMPRESS_PANEL = (LRSTATUS2)
1190.EQ..OR. COMPRESS_CB = ((LRSTATUS1)
1191.EQ. & (LRSTATUS3))
1192.GT. LR_ACTIVATED = (LRSTATUS0)
1193.AND..NOT. IF (COMPRESS_CB(COMPRESS_PANEL)) THEN
1194 COMPRESS_PANEL = .TRUE.
1195 LRSTATUS = 3
1196 ENDIF
1197 OOCWRITE_COMPATIBLE_WITH_BLR =
1198.NOT..OR..NOT..OR. & ( LR_ACTIVATED(COMPRESS_PANEL)
1199.NE. & (KEEP(486)2)
1200 & )
1201.eq..or..eq. IF((KEEP(24)0)(KEEP(24)1)) then
1202 force_cand=.FALSE.
1203 ELSE
1204.eq. force_cand=(mod(KEEP(24),2)0)
1205 end if
1206 TYPESPLIT = MUMPS_TYPESPLIT (PROCNODE_STEPS(STEP(INODE)),
1207 & KEEP(199))
1208.EQ..OR..EQ. IS_ofType5or6 = (TYPESPLIT5 TYPESPLIT6)
1209 ISTCHK = PIMASTER(STEP(IFSON))
1210 PDEST = ISTCHK + 6 + KEEP(IXSZ)
1211 NSLSON = IW(ISTCHK + KEEP(IXSZ) + 5)
1212 SPLIT_MAP_RESTART = .FALSE.
1213 IF (force_cand) THEN
1214 INIV2 = ISTEP_TO_INIV2( STEP( INODE ))
1215 NMB_OF_CAND = CAND( SLAVEF+1, INIV2 )
1216 NMB_OF_CAND_ORIG = NMB_OF_CAND
1217 SIZE_TMP_SLAVES_LIST = NMB_OF_CAND
1218 IF (IS_ofType5or6) THEN
1219 DO I=NMB_OF_CAND+1,SLAVEF
1220.LT. IF ( CAND( I, INIV2 )0) EXIT
1221 NMB_OF_CAND = NMB_OF_CAND +1
1222 ENDDO
1223 SIZE_TMP_SLAVES_LIST = NSLSON-1
1224.EQ. IF (INODE-999999) THEN
1225 SPLIT_MAP_RESTART = .TRUE.
1226 ENDIF
1227 ENDIF
1228.AND. IF (IS_ofType5or6SPLIT_MAP_RESTART) THEN
1229 TYPESPLIT = 4
1230 IS_ofType5or6 = .FALSE.
1231 SIZE_TMP_SLAVES_LIST = NMB_OF_CAND
1232 CAND (SLAVEF+1, INIV2) = SIZE_TMP_SLAVES_LIST
1233 ENDIF
1234 ELSE
1235 INIV2 = 1
1236 SIZE_TMP_SLAVES_LIST = SLAVEF - 1
1237 NMB_OF_CAND = SLAVEF - 1
1238 NMB_OF_CAND_ORIG = SLAVEF - 1
1239 ENDIF
1240 ALLOCATE(TMP_SLAVES_LIST(SIZE_TMP_SLAVES_LIST),stat=allocok)
1241 IF (allocok > 0 ) THEN
1242 GOTO 265
1243 ENDIF
1244.EQ. IF ( (TYPESPLIT4)
1245.OR..EQ..OR..EQ. & (TYPESPLIT5)(TYPESPLIT6)
1246 & ) THEN
1247.EQ. IF (TYPESPLIT4) THEN
1248 ALLOCATE(COPY_CAND(SLAVEF+1),stat=allocok)
1249 IF (allocok > 0 ) THEN
1250 GOTO 245
1251 ENDIF
1252 CALL SMUMPS_SPLIT_PREP_PARTITION (
1253 & INODE, STEP, N, SLAVEF,
1254 & PROCNODE_STEPS, KEEP, DAD, FILS,
1255 & CAND(1,INIV2), ICNTL, COPY_CAND,
1256 & NBSPLIT, NUMORG_SPLIT, TMP_SLAVES_LIST(1),
1257 & SIZE_TMP_SLAVES_LIST
1258 & )
1259 NCB_SPLIT = NCB-NUMORG_SPLIT
1260 SIZE_LIST_SPLIT = SIZE_TMP_SLAVES_LIST - NBSPLIT
1261 CALL SMUMPS_LOAD_SET_PARTITION( NCBSON_MAX, SLAVEF, KEEP,KEEP8,
1262 & ICNTL, COPY_CAND,
1263 & MEM_DISTRIB(0), NCB_SPLIT, NFRONT, NSLAVES,
1264 & TAB_POS_IN_PERE(1,INIV2),
1265 & TMP_SLAVES_LIST(NBSPLIT+1),
1266 & SIZE_LIST_SPLIT,INODE
1267 & )
1268 DEALLOCATE (COPY_CAND)
1269 CALL SMUMPS_SPLIT_POST_PARTITION (
1270 & INODE, STEP, N, SLAVEF, NBSPLIT, NCB,
1271 & PROCNODE_STEPS, KEEP, DAD, FILS,
1272 & ICNTL,
1273 & TAB_POS_IN_PERE(1,INIV2),
1274 & NSLAVES
1275 & )
1276 IF (SPLIT_MAP_RESTART) THEN
1277 IS_ofType5or6 = .TRUE.
1278 TYPESPLIT = MUMPS_TYPESPLIT (PROCNODE_STEPS(STEP(INODE)),
1279 & KEEP(199))
1280 CAND( SLAVEF+1, INIV2 ) = NMB_OF_CAND_ORIG
1281 ENDIF
1282 ELSE
1283 ISTCHK = PIMASTER(STEP(IFSON))
1284 PDEST = ISTCHK + 6 + KEEP(IXSZ)
1285 NSLSON = IW(ISTCHK + KEEP(IXSZ) + 5)
1286.EQ. IF (KEEP(376) 1) THEN
1287 NFRONT = IW( PIMASTER(STEP(IFSON)) + KEEP(IXSZ))
1288 ENDIF
1289 CALL SMUMPS_SPLIT_PROPAGATE_PARTI (
1290 & INODE, TYPESPLIT, IFSON,
1291 & CAND(1,INIV2), NMB_OF_CAND_ORIG,
1292 & IW(PDEST), NSLSON,
1293 & STEP, N, SLAVEF,
1294 & PROCNODE_STEPS, KEEP, DAD, FILS,
1295 & ICNTL, ISTEP_TO_INIV2, INIV2,
1296 & TAB_POS_IN_PERE, NSLAVES,
1297 & TMP_SLAVES_LIST,
1298 & SIZE_TMP_SLAVES_LIST
1299 & )
1300 ENDIF
1301 ELSE
1302 CALL SMUMPS_LOAD_SET_PARTITION( NCBSON_MAX, SLAVEF, KEEP,KEEP8,
1303 & ICNTL, CAND(1,INIV2),
1304 & MEM_DISTRIB(0), NCB, NFRONT, NSLAVES,
1305 & TAB_POS_IN_PERE(1,INIV2),
1306 & TMP_SLAVES_LIST,
1307 & SIZE_TMP_SLAVES_LIST,INODE
1308 & )
1309 ENDIF
1310 HF = NSLAVES + 6 + KEEP(IXSZ)
1311 LREQ_OOC = 0
1312.EQ..AND. IF (KEEP(201)1OOCWRITE_COMPATIBLE_WITH_BLR) THEN
1313 CALL SMUMPS_OOC_GET_PP_SIZES(KEEP(50), NASS1, NFRONT, NASS1,
1314 & NBPANELS_L, NBPANELS_U, LREQ_OOC)
1315 ENDIF
1316 LREQ = HF + 2 * NFRONT + LREQ_OOC
1317.GT. IF ((IWPOS + LREQ -1) IWPOSCB) THEN
1318 CALL SMUMPS_COMPRE_NEW(N, KEEP,
1319 & IW, LIW, A, LA,
1320 & LRLU, IPTRLU,
1321 & IWPOS, IWPOSCB, PTRIST, PTRAST,
1322 & STEP, PIMASTER, PAMASTER,
1323 & LRLUS,KEEP(IXSZ),
1324 & COMP, DKEEP(97), MYID, SLAVEF,
1325 & PROCNODE_STEPS, DAD)
1326.NE. IF (LRLU LRLUS) THEN
1327 IF (LPOK) THEN
1328 WRITE(LP, * ) 'PB compress SMUMPS_FAC_ASM_NIV2 ',
1329 & 'LRLU,LRLUS=',LRLU,LRLUS
1330 ENDIF
1331 GOTO 270
1332 ENDIF
1333.GT. IF ((IWPOS + LREQ -1) IWPOSCB) GOTO 270
1334 ENDIF
1335 IOLDPS = IWPOS
1336 IWPOS = IWPOS + LREQ
1337 NIV1 = .FALSE.
1338 ALLOCATE(SONROWS_PER_ROW(NFRONT-NASS1), stat=allocok)
1339 IF (allocok > 0) THEN
1340 GOTO 275
1341 ENDIF
1342 ISON_IN_PLACE = -9999
1343 CALL MUMPS_BUILD_SORT_INDEX(
1344 & MYID, INODE, N, IOLDPS, HF, LP, LPOK,
1345 & NFRONT, NFRONT_EFF, PERM, DAD,
1346 & NASS1, NASS, NUMSTK, NUMORG, IWPOSCB, IWPOS,
1347 & IFSON, STEP, PIMASTER, PTRIST, PTRAIW, IW, LIW,
1348 & INTARR, LINTARR, ITLOC, FILS, FRERE,
1349 & SON_LEVEL2, NIV1, KEEP,KEEP8, INFO(1),
1350 & ISON_IN_PLACE,
1351 & PROCNODE_STEPS, SLAVEF, SONROWS_PER_ROW,
1352 & NFRONT-NASS1)
1353.LT. IF (INFO(1)0) GOTO 250
1354.NE. IF ( NFRONT NFRONT_EFF ) THEN
1355 IF (
1356.EQ..OR..EQ. & (TYPESPLIT5) (TYPESPLIT6)) THEN
1357 WRITE(*,*) ' internal error 1 in fac_ass due to splitting ',
1358 & ' inode, nfront, nfront_eff =', INODE, NFRONT, NFRONT_EFF
1359 WRITE(*,*) ' splitting not yet ready for that'
1360 CALL MUMPS_ABORT()
1361 ENDIF
1362.GT. IF (NFRONTNFRONT_EFF) THEN
1363 NCB = NFRONT_EFF - NASS1
1364 NSLAVES_OLD = NSLAVES
1365 HF_OLD = HF
1366.EQ. IF (TYPESPLIT4) THEN
1367 ALLOCATE(COPY_CAND(SLAVEF+1),stat=allocok)
1368 IF (allocok > 0 ) THEN
1369 GOTO 245
1370 ENDIF
1371 CALL SMUMPS_SPLIT_PREP_PARTITION (
1372 & INODE, STEP, N, SLAVEF,
1373 & PROCNODE_STEPS, KEEP, DAD, FILS,
1374 & CAND(1,INIV2), ICNTL, COPY_CAND,
1375 & NBSPLIT, NUMORG_SPLIT, TMP_SLAVES_LIST(1),
1376 & SIZE_TMP_SLAVES_LIST
1377 & )
1378 NCB_SPLIT = NCB-NUMORG_SPLIT
1379 SIZE_LIST_SPLIT = SIZE_TMP_SLAVES_LIST - NBSPLIT
1380 CALL SMUMPS_LOAD_SET_PARTITION( NCBSON_MAX,
1381 & SLAVEF, KEEP,KEEP8,
1382 & ICNTL, COPY_CAND,
1383 & MEM_DISTRIB(0), NCB_SPLIT, NFRONT_EFF, NSLAVES,
1384 & TAB_POS_IN_PERE(1,INIV2),
1385 & TMP_SLAVES_LIST(NBSPLIT+1),
1386 & SIZE_LIST_SPLIT,INODE
1387 & )
1388 DEALLOCATE (COPY_CAND)
1389 CALL SMUMPS_SPLIT_POST_PARTITION (
1390 & INODE, STEP, N, SLAVEF, NBSPLIT, NCB,
1391 & PROCNODE_STEPS, KEEP, DAD, FILS,
1392 & ICNTL,
1393 & TAB_POS_IN_PERE(1,INIV2),
1394 & NSLAVES
1395 & )
1396 ELSE
1397 CALL SMUMPS_LOAD_SET_PARTITION( NCBSON_MAX,
1398 & SLAVEF, KEEP, KEEP8, ICNTL,
1399 & CAND(1,INIV2),
1400 & MEM_DISTRIB(0), NCB, NFRONT_EFF, NSLAVES,
1401 & TAB_POS_IN_PERE(1,INIV2),
1402 & TMP_SLAVES_LIST, SIZE_TMP_SLAVES_LIST,INODE
1403 & )
1404 ENDIF
1405 HF = NSLAVES + 6 + KEEP(IXSZ)
1406 IWPOS = IWPOS - ((2*NFRONT)-(2*NFRONT_EFF)) -
1407 & (NSLAVES_OLD - NSLAVES)
1408.NE. IF (NSLAVES_OLD NSLAVES) THEN
1409 IF (NSLAVES_OLD > NSLAVES) THEN
1410 DO KK=0,2*NFRONT_EFF-1
1411 IW(IOLDPS+HF+KK)=IW(IOLDPS+HF_OLD+KK)
1412 ENDDO
1413 ELSE
1414 IF (IWPOS - 1 > IWPOSCB ) GOTO 270
1415 DO KK=2*NFRONT_EFF-1, 0, -1
1416 IW(IOLDPS+HF+KK) = IW(IOLDPS+HF_OLD+KK)
1417 ENDDO
1418 END IF
1419 END IF
1420 NFRONT = NFRONT_EFF
1421 LREQ = HF + 2 * NFRONT + LREQ_OOC
1422 ELSE
1423 IF (LPOK) THEN
1424 WRITE(LP,*) MYID,': internal error 2 ',
1425 & ' in smumps_fac_asm_niv2 , inode=',
1426 & INODE, ' nfront, nfront_eff=', NFRONT, NFRONT_EFF
1427 ENDIF
1428 GOTO 270
1429 ENDIF
1430 ENDIF
1431.EQ..AND..NE..AND. IF (KEEP(201)1KEEP(50)1
1432 & OOCWRITE_COMPATIBLE_WITH_BLR) THEN
1433 CALL SMUMPS_OOC_PP_SET_PTR(KEEP(50),
1434 & NBPANELS_L, NBPANELS_U, NASS1,
1435 & IOLDPS + HF + 2 * NFRONT, IW, LIW)
1436 ENDIF
1437 MAXFRW = max0(MAXFRW, NFRONT)
1438 PTLUST(STEP(INODE)) = IOLDPS
1439 IW(IOLDPS+KEEP(IXSZ)) = NFRONT
1440 IW(IOLDPS + 1+KEEP(IXSZ)) = 0
1441 IW(IOLDPS + 2+KEEP(IXSZ)) = -NASS1
1442 IW(IOLDPS + 3+KEEP(IXSZ)) = -NASS1
1443 IW(IOLDPS + 4+KEEP(IXSZ)) = STEP(INODE)
1444 IW(IOLDPS+5+KEEP(IXSZ)) = NSLAVES
1445 IW(IOLDPS+6+KEEP(IXSZ):IOLDPS+5+KEEP(IXSZ)+NSLAVES)=
1446 & TMP_SLAVES_LIST(1:NSLAVES)
1447 ESTIM_NFS4FATHER_ATSON = -9999
1448.NE..AND..EQ. IF (KEEP(219)0KEEP(50)2) THEN
1449 IFATH = DAD( STEP( INODE) )
1450.NE. IF (IFATH0) THEN
1451.AND. IF (COMPRESS_CB
1452 & MUMPS_TYPENODE(PROCNODE_STEPS(STEP(IFATH)),KEEP(199))
1453.EQ. & 2 ) THEN
1454 IOLDPS = PTLUST(STEP(INODE))
1455 CALL SMUMPS_COMPUTE_ESTIM_NFS4FATHER (
1456 & N, INODE, IFATH, FILS, PERM, KEEP,
1457 & IOLDPS, HF, IW, LIW, NFRONT, NASS1,
1458 & ESTIM_NFS4FATHER_ATSON
1459 & )
1460 ENDIF
1461 ENDIF
1462 ENDIF
1463 CALL SMUMPS_LOAD_MASTER_2_ALL(MYID, SLAVEF, COMM_LOAD,
1464 & TAB_POS_IN_PERE(1,ISTEP_TO_INIV2(STEP(INODE))),
1465 & NASS1, KEEP, KEEP8, IW(IOLDPS+6+KEEP(IXSZ)), NSLAVES,INODE)
1466.EQ. IF(KEEP(86)1)THEN
1467.eq. IF(mod(KEEP(24),2)0)THEN
1468 CALL SMUMPS_LOAD_SEND_MD_INFO(SLAVEF,
1469 & CAND(SLAVEF+1,INIV2),
1470 & CAND(1,INIV2),
1471 & TAB_POS_IN_PERE(1,ISTEP_TO_INIV2(STEP(INODE))),
1472 & NASS1, KEEP,KEEP8, TMP_SLAVES_LIST,
1473 & NSLAVES,INODE)
1474.EQ..OR..EQ. ELSEIF((KEEP(24)0)(KEEP(24)1))THEN
1475 CALL SMUMPS_LOAD_SEND_MD_INFO(SLAVEF,
1476 & SLAVEF-1,
1477 & TMP_SLAVES_LIST,
1478 & TAB_POS_IN_PERE(1,ISTEP_TO_INIV2(STEP(INODE))),
1479 & NASS1, KEEP,KEEP8, TMP_SLAVES_LIST,
1480 & NSLAVES,INODE)
1481 ENDIF
1482 ENDIF
1483 DEALLOCATE(TMP_SLAVES_LIST)
1484.EQ. IF (KEEP(50)0) THEN
1485 LAELL8 = int(NASS1,8) * int(NFRONT,8)
1486 LDAFS = NFRONT
1487 ELSE
1488 LAELL8 = int(NASS1,8)*int(NASS1,8)
1489.NE..AND..EQ. IF(KEEP(219)0KEEP(50) 2)
1490 & LAELL8 = LAELL8+int(NASS1,8)
1491 LDAFS = NASS1
1492 ENDIF
1493 CALL SMUMPS_GET_SIZE_NEEDED
1494 & (0, LAELL8, .FALSE.,
1495 & KEEP(1), KEEP8(1),
1496 & N,IW,LIW,A,LA,
1497 & LRLU,IPTRLU,IWPOS,IWPOSCB,
1498 & PTRIST,PTRAST,
1499 & STEP, PIMASTER,PAMASTER,LRLUS,
1500 & KEEP(IXSZ), COMP, DKEEP(97), MYID,
1501 & SLAVEF, PROCNODE_STEPS, DAD,
1502 & INFO(1), INFO(2))
1503.LT. IF (INFO(1)0) GOTO 490
1504 LRLU = LRLU - LAELL8
1505 LRLUS = LRLUS - LAELL8
1506 KEEP8(67) = min(LRLUS, KEEP8(67))
1507 KEEP8(69) = KEEP8(69) + LAELL8
1508 KEEP8(68) = max(KEEP8(69), KEEP8(68))
1509 POSELT = POSFAC
1510 PTRAST(STEP(INODE)) = POSELT
1511 PTRFAC(STEP(INODE)) = POSELT
1512 POSFAC = POSFAC + LAELL8
1513 IW(IOLDPS+XXI) = LREQ
1514 CALL MUMPS_STOREI8(LAELL8,IW(IOLDPS+XXR))
1515 CALL MUMPS_STOREI8(0_8,IW(IOLDPS+XXD))
1516 IW(IOLDPS+XXS) = -9999
1517 IW(IOLDPS+XXN) = -99999
1518 IW(IOLDPS+XXP) = -99999
1519 IW(IOLDPS+XXA) = -99999
1520 IW(IOLDPS+XXF) = -99999
1521 IW(IOLDPS+XXLR)= LRSTATUS
1522 IW(IOLDPS+XXG) = MemNotPinned
1523 CALL SMUMPS_LOAD_MEM_UPDATE(.FALSE.,.FALSE.,LA-LRLUS,0_8,LAELL8,
1524 & KEEP,KEEP8,LRLUS)
1525 POSEL1 = POSELT - int(LDAFS,8)
1526#if defined(ZERO_TRIANGLE)
1527 LAPOS2 = POSELT + LAELL8 - 1_8
1528 A(POSELT:LAPOS2) = ZERO
1529#else
1530.eq..OR..lt. IF ( KEEP(50) 0 LDAFS KEEP(63) ) THEN
1531 LAPOS2 = POSELT + LAELL8 - 1_8
1532!$ CHUNK8 = int(KEEP(361),8)
1533!$OMP PARALLEL DO PRIVATE(JJ8) SCHEDULE(STATIC, CHUNK8)
1534.AND..GT.!$OMP& IF (LAPOS2 - POSELT > int(KEEP(361),8) NOMP 1)
1535#if defined(__ve__)
1536!NEC$ IVDEP
1537#endif
1538 DO JJ8 = POSELT, LAPOS2
1539 A(JJ8) = ZERO
1540 ENDDO
1541!$OMP END PARALLEL DO
1542 ELSE
1543 TOPDIAG = max(KEEP(7), KEEP(8))-1
1544 IF (LR_ACTIVATED) THEN
1545 NULLIFY(BEGS_BLR)
1546 CALL GET_CUT(IW(IOLDPS+HF:IOLDPS+HF+NFRONT-1), NASS1,
1547 & 0, LRGROUPS, NPARTSCB,
1548 & NPARTSASS, BEGS_BLR)
1549 NB_BLR = NPARTSASS + NPARTSCB
1550 CALL MAX_CLUSTER(BEGS_BLR,NB_BLR,MAXI_CLUSTER)
1551 DEALLOCATE(BEGS_BLR)
1552 CALL COMPUTE_BLR_VCS(KEEP(472), IBCKSZ2, KEEP(488), NASS1)
1553 MINSIZE = int(IBCKSZ2 / 2)
1554 TOPDIAG = max(2*MINSIZE + MAXI_CLUSTER-1, TOPDIAG)
1555 ENDIF
1556!$ CHUNK = max(KEEP(360)/2,
1557!$ & ( (( LDAFS + NOMP -1 ) / NOMP + 2) / 3) )
1558 APOS = POSELT
1559#if defined(__ve__)
1560!NEC$ IVDEP
1561#endif
1562!$OMP PARALLEL DO PRIVATE(APOS,JJ3) SCHEDULE(STATIC, CHUNK)
1563.GT..AND..GT.!$OMP& IF (LDAFS - 1 KEEP(360) NOMP 1)
1564 DO JJ8 = 0_8, int(LDAFS-1,8)
1565 APOS = POSELT + JJ8 * int(LDAFS,8)
1566 JJ3 = min( int(LDAFS,8) - 1_8, JJ8 + TOPDIAG )
1567 A(APOS:APOS+JJ3) = ZERO
1568 END DO
1569!$OMP END PARALLEL DO
1570.NE..AND..EQ. IF (KEEP(219)0KEEP(50)2) THEN
1571 APOSMAX = POSELT + int(NASS1,8)*int(NASS1,8)
1572 A(APOSMAX:APOSMAX+int(LDAFS-1,8))=ZERO
1573 ENDIF
1574 END IF
1575#endif
1576.NE..AND..NE. IF ((NUMSTK0)(NASS0)) THEN
1577 ISON = IFSON
1578 DO 220 IELL = 1, NUMSTK
1579 ISTCHK = PIMASTER(STEP(ISON))
1580 NELIM = IW(ISTCHK + KEEP(IXSZ) + 1)
1581.EQ. IF (NELIM0) GOTO 210
1582 LSTK = IW(ISTCHK + KEEP(IXSZ))
1583 NPIVS = IW(ISTCHK + 3+KEEP(IXSZ))
1584.LT. IF (NPIVS0) NPIVS=0
1585 NSLSON = IW(ISTCHK + 5+KEEP(IXSZ))
1586 HS = 6 + NSLSON + KEEP(IXSZ)
1587 NCOLS = NPIVS + LSTK
1588.LT. SAME_PROC = (ISTCHKIWPOS)
1589 IF ( SAME_PROC ) THEN
1590 ISTCHK_CB_RIGHT=PTRIST(STEP(ISON))
1591 ELSE
1592 ISTCHK_CB_RIGHT=ISTCHK
1593 ENDIF
1594 SON_XXS = IW(ISTCHK_CB_RIGHT + XXS)
1595.EQ. PACKED_CB = ( SON_XXS S_CB1COMP )
1596.NOT. IF (SAME_PROC) THEN
1597 NROWS = IW(ISTCHK + KEEP(IXSZ) + 2)
1598 ELSE
1599 NROWS = NCOLS
1600 ENDIF
1601.EQ. IF (KEEP(50)0) THEN
1602 LDA_SON = LSTK
1603 LCB = int(NELIM,8)*int(LSTK,8)
1604 ELSE
1605.EQ. IF (NSLSON0) THEN
1606 IF (SAME_PROC) THEN
1607.EQ..OR. IS_CB_LR = IW(ISTCHK_CB_RIGHT+XXLR) 1
1608.EQ. & IW(ISTCHK_CB_RIGHT+XXLR) 3
1609 IF (IS_CB_LR) THEN
1610 LDA_SON = NELIM
1611 ELSE
1612 LDA_SON = LSTK
1613 ENDIF
1614 ELSE
1615 LDA_SON = LSTK
1616 ENDIF
1617 ELSE
1618 LDA_SON = NELIM
1619 ENDIF
1620 IF (PACKED_CB) THEN
1621 LCB = (int(NELIM,8)*int(NELIM+1,8))/2_8
1622 ELSE
1623 LCB = int(LDA_SON,8)*int(NELIM,8)
1624 ENDIF
1625 ENDIF
1626.EQ. IF (KEEP(50) 0) THEN
1627 OPASSW = OPASSW + dble(LCB)
1628 ELSE
1629 OPASSW = OPASSW + int(NELIM,8)*int(NELIM+1,8)/2_8
1630 ENDIF
1631 IS_DYNAMIC_CB =
1632 & SMUMPS_DM_IS_DYNAMIC(IW(ISTCHK_CB_RIGHT+XXD:
1633 & ISTCHK_CB_RIGHT+XXD+1))
1634 IF ( IS_DYNAMIC_CB ) THEN
1635 CALL MUMPS_GETI8(DYN_SIZE, IW(ISTCHK_CB_RIGHT+XXD))
1636 CALL SMUMPS_DM_SET_PTR( PAMASTER(STEP(ISON)), DYN_SIZE,
1637 & SON_A )
1638 IACHK = 1_8
1639 ELSE
1640 IACHK = PAMASTER(STEP(ISON))
1641 SON_A=>A
1642 ENDIF
1643 K1 = ISTCHK + HS + NROWS + NPIVS
1644 K2 = K1 + NELIM - 1
1645.eq. IF (KEEP(50)0) THEN
1646 IF (IS_ofType5or6) THEN
1647 APOS = POSELT
1648 DO JJ8 = 1_8, int(NELIM,8)*int(LSTK,8)
1649 A(APOS+JJ8-1_8) = A(APOS+JJ8-1_8) + SON_A(IACHK+JJ8-1_8)
1650 ENDDO
1651 ELSE
1652 DO 170 KK = K1, K2
1653 APOS = POSEL1 + int(IW(KK),8) * int(LDAFS,8)
1654 DO 160 KK1 = 1, LSTK
1655 JJ2 = APOS + int(IW(K1 + KK1 - 1),8) - 1_8
1656 A(JJ2) = A(JJ2) + SON_A(IACHK + int(KK1 - 1,8))
1657 160 CONTINUE
1658 IACHK = IACHK + int(LSTK,8)
1659 170 CONTINUE
1660 ENDIF
1661 ELSE
1662.GT. IF (LCB 0) THEN
1663 CALL SMUMPS_LDLT_ASM_NIV12(A, LA, SON_A(IACHK),
1664 & POSELT, LDAFS, NASS1,
1665 & LDA_SON, LCB,
1666 & IW( K1 ), NELIM, NELIM, ETATASS,
1667 & PACKED_CB
1668!$ & , KEEP(360)
1669 & )
1670 ENDIF
1671 ENDIF
1672 210 ISON = FRERE(STEP(ISON))
1673 220 CONTINUE
1674 ENDIF
1675 IBROT = INODE
1676 APOSMAX = POSELT + int(NASS1,8)*int(NASS1,8)
1677 DO 260 IORG = 1, NUMORG
1678 JK8 = PTRAIW(IBROT)
1679 AINPUT8 = PTRARW(IBROT)
1680 JJ8 = JK8 + 1_8
1681 J18 = JJ8 + 1_8
1682 J28 = J18 + INTARR(JK8)
1683 J38 = J28 + 1_8
1684 J48 = J28 - INTARR(JJ8)
1685 IJROW = INTARR(J18)
1686 ICT12 = POSELT + int(IJROW - 1 - LDAFS, 8)
1687 MAXARR = RZERO
1688 DO JJ8 = J18, J28
1689.NE. IF (KEEP(219)0) THEN
1690.LE. IF (INTARR(JJ8)NASS1) THEN
1691 APOS2 = ICT12 + int(INTARR(JJ8),8) * int(LDAFS,8)
1692 A(APOS2) = A(APOS2) + DBLARR(AINPUT8)
1693.EQ. ELSEIF (KEEP(50)2) THEN
1694 MAXARR = max(MAXARR,abs(DBLARR(AINPUT8)))
1695 ENDIF
1696 ELSE
1697.LE. IF (INTARR(JJ8)NASS1) THEN
1698 APOS2 = ICT12 + int(INTARR(JJ8),8) * int(LDAFS,8)
1699 A(APOS2) = A(APOS2) + DBLARR(AINPUT8)
1700 ENDIF
1701 ENDIF
1702 AINPUT8 = AINPUT8 + 1_8
1703 ENDDO
1704.NE..AND..EQ. IF(KEEP(219)0KEEP(50) 2) THEN
1705 A(APOSMAX+int(IJROW-1,8)) = MAXARR
1706 ENDIF
1707.GT. IF (J38 J48) GOTO 255
1708 ICT13 = POSELT + int(IJROW - 1,8) * int(LDAFS,8)
1709 NBCOL = int(J48 - J38 + 1_8)
1710 DO JJ8 = 1_8, int(NBCOL,8)
1711 JJ3 = ICT13 + int(INTARR(J38 + JJ8 - 1_8),8) - 1_8
1712 A(JJ3) = A(JJ3) + DBLARR(AINPUT8 + JJ8 - 1_8)
1713 ENDDO
1714 255 CONTINUE
1715.EQ. IF (KEEP(50)0) THEN
1716 DO J253 = 1, KEEP(253)
1717 APOS = POSELT +
1718 & int(IJROW-1,8) * int(LDAFS,8) +
1719 & int(LDAFS-KEEP(253)+J253-1,8)
1720 A(APOS) = A(APOS) + RHS_MUMPS( (J253-1)*KEEP(254)+IBROT )
1721 ENDDO
1722 ENDIF
1723 IBROT = FILS(IBROT)
1724 260 CONTINUE
1725 PTRCOL = IOLDPS + HF + NFRONT
1726 PTRROW = IOLDPS + HF + NASS1
1727 PDEST = IOLDPS + 6 + KEEP(IXSZ)
1728 IBC_SOURCE = MYID
1729 DO ISLAVE = 1, NSLAVES
1730 CALL MUMPS_BLOC2_GET_SLAVE_INFO(
1731 & KEEP,KEEP8, INODE, STEP, N, SLAVEF,
1732 & ISTEP_TO_INIV2, TAB_POS_IN_PERE,
1733 & ISLAVE, NCB,
1734 & NSLAVES,
1735 & NBLIG, FIRST_INDEX )
1736 SHIFT_INDEX = FIRST_INDEX - 1
1737 IERR = -1
1738.EQ. DO WHILE (IERR -1)
1739.eq. IF ( KEEP(50) 0 ) THEN
1740 NBCOL = NFRONT
1741 CALL SMUMPS_BUF_SEND_DESC_BANDE( INODE,
1742 & sum(SONROWS_PER_ROW(FIRST_INDEX:FIRST_INDEX+NBLIG-1)),
1743 & NBLIG, IW(PTRROW), NBCOL, IW(PTRCOL), NASS1,
1744 & IZERO, IDUMMY,
1745 & NSLAVES,
1746 & ESTIM_NFS4FATHER_ATSON,
1747 & IW(PDEST), IBC_SOURCE, NFRONT, COMM, KEEP, IERR
1748 & , IW(IOLDPS+XXLR)
1749 & )
1750 ELSE
1751 NBCOL = NASS1+SHIFT_INDEX+NBLIG
1752 CALL SMUMPS_BUF_SEND_DESC_BANDE( INODE,
1753 & sum(SONROWS_PER_ROW(FIRST_INDEX:FIRST_INDEX+NBLIG-1)),
1754 & NBLIG, IW(PTRROW), NBCOL, IW(PTRCOL), NASS1,
1755 & NSLAVES-ISLAVE,
1756 & IW( PTLUST(STEP(INODE))+6+KEEP(IXSZ)+ISLAVE),
1757 & NSLAVES,
1758 & ESTIM_NFS4FATHER_ATSON,
1759 & IW(PDEST), IBC_SOURCE, NFRONT, COMM, KEEP, IERR
1760 & , IW(IOLDPS+XXLR)
1761 & )
1762 ENDIF
1763.EQ. IF (IERR-1) THEN
1764 BLOCKING = .FALSE.
1765 SET_IRECV = .TRUE.
1766 MESSAGE_RECEIVED = .FALSE.
1767 CALL SMUMPS_TRY_RECVTREAT( COMM_LOAD, ASS_IRECV,
1768 & BLOCKING, SET_IRECV, MESSAGE_RECEIVED,
1769 & MPI_ANY_SOURCE, MPI_ANY_TAG,
1770 & STATUS, BUFR, LBUFR,
1771 & LBUFR_BYTES,
1772 & PROCNODE_STEPS, POSFAC, IWPOS, IWPOSCB, IPTRLU,
1773 & LRLU, LRLUS, N, IW, LIW, A, LA, PTRIST,
1774 & PTLUST, PTRFAC,
1775 & PTRAST, STEP, PIMASTER, PAMASTER, NSTK_S, COMP, INFO(1),
1776 & INFO(2), COMM,
1777 & PERM,
1778 & IPOOL, LPOOL, LEAF, NBFIN, MYID, SLAVEF,
1779 & root, OPASSW, OPELIW, ITLOC, RHS_MUMPS,
1780 & FILS, DAD, PTRARW, PTRAIW,
1781 & INTARR, DBLARR, ICNTL, KEEP,KEEP8,DKEEP, ND, FRERE,
1782 & LPTRAR, NELT, IW, IW,
1783 & ISTEP_TO_INIV2, TAB_POS_IN_PERE, .TRUE.
1784 & , LRGROUPS
1785 & )
1786.LT. IF ( INFO(1) 0 ) GOTO 500
1787 IF (MESSAGE_RECEIVED) THEN
1788 IOLDPS = PTLUST(STEP(INODE))
1789 PTRCOL = IOLDPS + HF + NFRONT
1790 PTRROW = IOLDPS + HF + NASS1 + SHIFT_INDEX
1791 ENDIF
1792 ENDIF
1793 ENDDO
1794.EQ. IF (IERR -2) GOTO 300
1795.EQ. IF (IERR -3) GOTO 305
1796 PTRROW = PTRROW + NBLIG
1797 PDEST = PDEST + 1
1798 ENDDO
1799 DEALLOCATE(SONROWS_PER_ROW)
1800.EQ. IF (NUMSTK0) GOTO 500
1801 ISON = IFSON
1802 DO IELL = 1, NUMSTK
1803 ISTCHK = PIMASTER(STEP(ISON))
1804 NELIM = IW(ISTCHK + 1 + KEEP(IXSZ))
1805 LSTK = IW(ISTCHK + KEEP(IXSZ))
1806 NPIVS = IW(ISTCHK + 3 + KEEP(IXSZ))
1807.LT. IF ( NPIVS 0 ) NPIVS = 0
1808 NSLSON = IW(ISTCHK + 5 + KEEP(IXSZ))
1809 HS = 6 + NSLSON + KEEP(IXSZ)
1810 NCOLS = NPIVS + LSTK
1811.LT. SAME_PROC = (ISTCHKIWPOS)
1812.NOT. IF (SAME_PROC) THEN
1813 NROWS = IW(ISTCHK + 2 + KEEP(IXSZ) )
1814 ELSE
1815 NROWS = NCOLS
1816 ENDIF
1817 PDEST = ISTCHK + 6 + KEEP(IXSZ)
1818 NCBSON = LSTK - NELIM
1819 PTRCOL = ISTCHK + HS + NROWS + NPIVS + NELIM
1820.NE..AND..EQ. IF (KEEP(219)0KEEP(50)2) THEN
1821 NFS4FATHER = NCBSON
1822 DO I=0,NCBSON-1
1823.GT. IF(IW(PTRCOL+I) NASS1) THEN
1824 NFS4FATHER = I
1825 EXIT
1826 ENDIF
1827 ENDDO
1828 NFS4FATHER = NFS4FATHER + NELIM
1829 ELSE
1830 NFS4FATHER = 0
1831 ENDIF
1832.EQ. IF (NSLSON0) THEN
1833 NSLSON = 1
1834 PDEST1(1) = MUMPS_PROCNODE(PROCNODE_STEPS(STEP(ISON)),
1835 & KEEP(199))
1836.EQ. IF (PDEST1(1)MYID) THEN
1837 CALL SMUMPS_MAPLIG_FILS_NIV1( COMM_LOAD, ASS_IRECV,
1838 & BUFR, LBUFR, LBUFR_BYTES,
1839 & INODE, ISON, NSLAVES,
1840 & IW( PTLUST(STEP(INODE)) + 6 +KEEP(IXSZ)),
1841 & NFRONT, NASS1, NFS4FATHER, NCBSON, IW( PTRCOL ),
1842 & PROCNODE_STEPS,
1843 & SLAVEF, POSFAC, IWPOS, IWPOSCB, IPTRLU, LRLU,
1844 & LRLUS, N, IW, LIW, A, LA,
1845 & PTRIST, PTLUST, PTRFAC, PTRAST, STEP,
1846 & PIMASTER, PAMASTER, NSTK_S, COMP,
1847 & INFO(1), INFO(2), MYID, COMM, PERM,
1848 & IPOOL, LPOOL, LEAF,
1849 & NBFIN, ICNTL, KEEP, KEEP8, DKEEP, root,
1850 & OPASSW, OPELIW,
1851 & ITLOC, RHS_MUMPS, FILS, DAD,
1852 & PTRARW, PTRAIW, INTARR, DBLARR,
1853 & ND, FRERE, LPTRAR, NELT, IW, IW,
1854 & ISTEP_TO_INIV2, TAB_POS_IN_PERE,
1855 & LRGROUPS )
1856.LT. IF ( INFO(1) 0 ) GOTO 500
1857 ELSE
1858 IERR = -1
1859.EQ. DO WHILE (IERR-1)
1860 PTRCOL = PIMASTER(STEP(ISON)) + HS + NROWS + NPIVS + NELIM
1861 CALL SMUMPS_BUF_SEND_MAPLIG(
1862 & INODE, NFRONT,NASS1,NFS4FATHER,
1863 & ISON, MYID,
1864 & NSLAVES, IW( PTLUST(STEP(INODE)) + 6 +KEEP(IXSZ)),
1865 & IW(PTRCOL), NCBSON,
1866 & COMM, IERR, PDEST1, NSLSON, SLAVEF,
1867 & KEEP,KEEP8, STEP, N,
1868 & ISTEP_TO_INIV2, TAB_POS_IN_PERE
1869 & )
1870.EQ. IF (IERR-1) THEN
1871 BLOCKING = .FALSE.
1872 SET_IRECV = .TRUE.
1873 MESSAGE_RECEIVED = .FALSE.
1874 CALL SMUMPS_TRY_RECVTREAT( COMM_LOAD, ASS_IRECV,
1875 & BLOCKING, SET_IRECV, MESSAGE_RECEIVED,
1876 & MPI_ANY_SOURCE, MPI_ANY_TAG,
1877 & STATUS, BUFR, LBUFR, LBUFR_BYTES,
1878 & PROCNODE_STEPS, POSFAC, IWPOS, IWPOSCB, IPTRLU,
1879 & LRLU, LRLUS, N, IW, LIW, A, LA, PTRIST,
1880 & PTLUST, PTRFAC,
1881 & PTRAST, STEP, PIMASTER, PAMASTER, NSTK_S, COMP, INFO(1),
1882 & INFO(2), COMM,
1883 & PERM,
1884 & IPOOL, LPOOL, LEAF, NBFIN, MYID, SLAVEF,
1885 & root,OPASSW, OPELIW, ITLOC, RHS_MUMPS, FILS, DAD,
1886 & PTRARW, PTRAIW,
1887 & INTARR, DBLARR, ICNTL, KEEP,KEEP8,DKEEP, ND, FRERE,
1888 & LPTRAR,
1889 & NELT, IW, IW,
1890 & ISTEP_TO_INIV2, TAB_POS_IN_PERE, .TRUE.
1891 & , LRGROUPS
1892 & )
1893.LT. IF ( INFO(1) 0 ) GOTO 500
1894 ENDIF
1895 ENDDO
1896.EQ. IF (IERR -2) GOTO 290
1897.EQ. IF (IERR -3) GOTO 295
1898 ENDIF
1899 ELSE
1900.GT. IF (PIMASTER(STEP(ISON))0) THEN
1901 IERR = -1
1902.EQ. DO WHILE (IERR-1)
1903 PTRCOL = PIMASTER(STEP(ISON)) + HS + NROWS + NPIVS + NELIM
1904 PDEST = PIMASTER(STEP(ISON)) + 6 + KEEP(IXSZ)
1905 CALL SMUMPS_BUF_SEND_MAPLIG(
1906 & INODE, NFRONT, NASS1, NFS4FATHER,
1907 & ISON, MYID,
1908 & NSLAVES, IW(PTLUST(STEP(INODE))+6+KEEP(IXSZ)),
1909 & IW(PTRCOL), NCBSON,
1910 & COMM, IERR, IW(PDEST), NSLSON, SLAVEF,
1911 & KEEP,KEEP8, STEP, N,
1912 & ISTEP_TO_INIV2, TAB_POS_IN_PERE
1913 & )
1914.EQ. IF (IERR-1) THEN
1915 BLOCKING = .FALSE.
1916 SET_IRECV = .TRUE.
1917 MESSAGE_RECEIVED = .FALSE.
1918 CALL SMUMPS_TRY_RECVTREAT( COMM_LOAD, ASS_IRECV,
1919 & BLOCKING, SET_IRECV, MESSAGE_RECEIVED,
1920 & MPI_ANY_SOURCE, MPI_ANY_TAG,
1921 & STATUS, BUFR, LBUFR,
1922 & LBUFR_BYTES,
1923 & PROCNODE_STEPS, POSFAC, IWPOS, IWPOSCB, IPTRLU,
1924 & LRLU, LRLUS, N, IW, LIW, A, LA, PTRIST,
1925 & PTLUST, PTRFAC,
1926 & PTRAST, STEP, PIMASTER, PAMASTER, NSTK_S, COMP, INFO(1),
1927 & INFO(2), COMM,
1928 & PERM,
1929 & IPOOL, LPOOL, LEAF, NBFIN, MYID, SLAVEF,
1930 & root,OPASSW, OPELIW, ITLOC, RHS_MUMPS,
1931 & FILS, DAD, PTRARW, PTRAIW,
1932 & INTARR, DBLARR, ICNTL, KEEP,KEEP8,DKEEP, ND, FRERE,
1933 & LPTRAR, NELT, IW, IW,
1934 & ISTEP_TO_INIV2, TAB_POS_IN_PERE, .TRUE.
1935 & , LRGROUPS
1936 & )
1937.LT. IF ( INFO(1) 0 ) GOTO 500
1938 ENDIF
1939 ENDDO
1940.EQ. IF (IERR -2) GOTO 290
1941.EQ. IF (IERR -3) GOTO 295
1942 ENDIF
1943 DO ISLAVE = 0, NSLSON-1
1944.EQ. IF (IW(PDEST+ISLAVE)MYID) THEN
1945 CALL MUMPS_BLOC2_GET_SLAVE_INFO(
1946 & KEEP,KEEP8, ISON, STEP, N, SLAVEF,
1947 & ISTEP_TO_INIV2, TAB_POS_IN_PERE,
1948 & ISLAVE+1, NCBSON,
1949 & NSLSON,
1950 & TROW_SIZE, FIRST_INDEX )
1951 SHIFT_INDEX = FIRST_INDEX - 1
1952 INDX = PTRCOL + SHIFT_INDEX
1953 CALL SMUMPS_MAPLIG( COMM_LOAD, ASS_IRECV,
1954 & BUFR, LBUFR, LBUFR_BYTES,
1955 & INODE, ISON, NSLAVES,
1956 & IW( PTLUST(STEP(INODE))+6+KEEP(IXSZ)),
1957 & NFRONT, NASS1,NFS4FATHER,
1958 & TROW_SIZE, IW( INDX ),
1959 & PROCNODE_STEPS,
1960 & SLAVEF, POSFAC, IWPOS, IWPOSCB, IPTRLU, LRLU,
1961 & LRLUS, N, IW, LIW, A, LA,
1962 & PTRIST, PTLUST, PTRFAC, PTRAST, STEP,
1963 & PIMASTER, PAMASTER, NSTK_S, COMP, INFO(1), INFO(2),
1964 & MYID, COMM, PERM, IPOOL, LPOOL, LEAF,
1965 & NBFIN, ICNTL, KEEP,KEEP8,DKEEP, root,
1966 & OPASSW, OPELIW, ITLOC, RHS_MUMPS, FILS, DAD,
1967 & PTRARW, PTRAIW,
1968 & INTARR, DBLARR, ND, FRERE, LPTRAR, NELT, IW,
1969 & IW,
1970 &
1971 & ISTEP_TO_INIV2, TAB_POS_IN_PERE, LRGROUPS)
1972.LT. IF ( INFO(1) 0 ) GOTO 500
1973 EXIT
1974 ENDIF
1975 ENDDO
1976 ENDIF
1977 ISON = FRERE(STEP(ISON))
1978 ENDDO
1979 GOTO 500
1980 250 CONTINUE
1981.EQ. IF (INFO(1)-13) THEN
1982 IF (LPOK) THEN
1983 WRITE( LP, * )
1984 & ' failure in INTEGER DYNAMIC ALLOCATION DURING
1985 & SMUMPS_FAC_ASM_NIV2'
1986 ENDIF
1987 INFO(2) = NUMSTK + 1
1988 ENDIF
1989 GOTO 490
1990 245 CONTINUE
1991 IF (LPOK) THEN
1992 WRITE( LP, * ) ' FAILURE ALLOCATING COPY_CAND',
1993 & ' DURING SMUMPS_FAC_ASM_NIV2'
1994 ENDIF
1995 INFO(1) = -13
1996 INFO(2) = SLAVEF+1
1997 GOTO 490
1998 265 CONTINUE
1999 IF (LPOK) THEN
2000 WRITE( LP, * ) ' FAILURE ALLOCATING TMP_SLAVES_LIST',
2001 & ' DURING SMUMPS_FAC_ASM_NIV2'
2002 ENDIF
2003 INFO(1) = -13
2004 INFO(2) = SIZE_TMP_SLAVES_LIST
2005 GOTO 490
2006 270 CONTINUE
2007 INFO(1) = -8
2008 INFO(2) = LREQ
2009 IF (LPOK) THEN
2010 WRITE( LP, * )
2011 & ' FAILURE IN INTEGER ALLOCATION DURING SMUMPS_FAC_ASM_NIV2'
2012 ENDIF
2013 GOTO 490
2014 275 CONTINUE
2015 IF (LPOK) THEN
2016 WRITE( LP, * ) ' FAILURE ALLOCATING SONROWS_PER_ROW',
2017 & ' DURING SMUMPS_FAC_ASM_NIV2'
2018 ENDIF
2019 INFO(1) = -13
2020 INFO(2) = NFRONT-NASS1
2021 GOTO 490
2022 290 CONTINUE
2023.GT..AND..GE. IF ((ICNTL(1) 0) (ICNTL(4) 1)) THEN
2024 LP = ICNTL(1)
2025 WRITE( LP, * )
2026 &' FAILURE, SEND BUFFER TOO SMALL (1) DURING SMUMPS_FAC_ASM_NIV2'
2027 ENDIF
2028 INFO(1) = -17
2029 LREQ = NCBSON + 6 + NSLSON+KEEP(IXSZ)
2030 INFO(2) = LREQ * KEEP( 34 )
2031 GOTO 490
2032 295 CONTINUE
2033.GT..AND..GE. IF ((ICNTL(1) 0) (ICNTL(4) 1)) THEN
2034 LP = ICNTL(1)
2035 WRITE( LP, * )
2036 &' FAILURE, RECV BUFFER TOO SMALL (1) DURING SMUMPS_FAC_ASM_NIV2'
2037 ENDIF
2038 INFO(1) = -20
2039 LREQ = NCBSON + 6 + NSLSON+KEEP(IXSZ)
2040 INFO(2) = LREQ * KEEP( 34 )
2041 GOTO 490
2042 300 CONTINUE
2043.GT..AND..GE. IF ((ICNTL(1) 0) (ICNTL(4) 1)) THEN
2044 LP = ICNTL(1)
2045 WRITE( LP, * )
2046 &' FAILURE, SEND BUFFER TOO SMALL (2)',
2047 &' DURING SMUMPS_FAC_ASM_NIV2'
2048 ENDIF
2049 INFO(1) = -17
2050 LREQ = NBLIG + NBCOL + 4 + KEEP(IXSZ)
2051 INFO(2) = LREQ * KEEP( 34 )
2052 GOTO 490
2053 305 CONTINUE
2054.GT..AND..GE. IF ((ICNTL(1) 0) (ICNTL(4) 1)) THEN
2055 LP = ICNTL(1)
2056 WRITE( LP, * )
2057 &' FAILURE, RECV BUFFER TOO SMALL (2)',
2058 &' DURING SMUMPS_FAC_ASM_NIV2'
2059 ENDIF
2060 INFO(1) = -20
2061 LREQ = NBLIG + NBCOL + 4 + KEEP(IXSZ)
2062 INFO(2) = LREQ * KEEP( 34 )
2063 GOTO 490
2064 490 CALL SMUMPS_BDC_ERROR( MYID, SLAVEF, COMM, KEEP )
2065 500 CONTINUE
2066 RETURN
2067 END SUBROUTINE SMUMPS_FAC_ASM_NIV2
2068 END MODULE SMUMPS_FAC_ASM_MASTER_M
mumps_abort
#define mumps_abort
Definition VE_Metis.h:25
mumps_get_flops_cost
subroutine mumps_get_flops_cost(nfront, npiv, nass, keep50, level, cost)
Definition estim_flops.F:74
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
for
for(i8=*sizetab-1;i8 >=0;i8--)
Definition mumps_common.c:91
mumps_bloc2_get_slave_info
subroutine mumps_bloc2_get_slave_info(keep, keep8, inode, step, n, slavef, istep_to_iniv2, tab_pos_in_pere islave, ncb, nslaves, size, first_index)
Definition mumps_type2_blocking.F:367
mumps_build_sort_index_m
Definition fac_asm_build_sort_index_m.F:14
mumps_build_sort_index_m::mumps_build_sort_index
subroutine mumps_build_sort_index(myid, inode, n, ioldps, hf, lp, lpok, nfront, nfront_eff, perm, dad, nass1, nass, numstk, numorg, iwposcb, iwpos, ifson, step, pimaster, ptrist, ptraiw, iw, liw, intarr, lintarr, itloc, fils, frere_steps, son_level2, niv1, keep, keep8, iflag, ison_in_place, procnode_steps, slavef, sonrows_per_row, lsonrows_per_row, mumps_tps_arr, l0_omp_mapping)
Definition fac_asm_build_sort_index_m.F:27
mumps_lr_common
Definition lr_common.F:14
mumps_lr_common::compute_blr_vcs
subroutine compute_blr_vcs(k472, ibcksz, maxsize, nass)
Definition lr_common.F:18
mumps_tps_m
Definition omp_tps_common_m.F:14
smumps_ana_lr
Definition sana_lr.F:14
smumps_ana_lr::get_cut
subroutine get_cut(iwr, nass, ncb, lrgroups, npartscb, npartsass, cut)
Definition sana_lr.F:25
smumps_buf
Definition smumps_comm_buffer.F:14
smumps_buf::smumps_buf_send_maplig
subroutine, public smumps_buf_send_maplig(inode, nfront, nass1, nfs4father, ison, myid, nslaves, slaves_pere, trow, ncbson, comm, ierr, dest, ndest, slavef, keep, keep8, step, n, istep_to_iniv2, tab_pos_in_per)
Definition smumps_comm_buffer.F:1649
smumps_dynamic_memory_m
Definition sfac_mem_dynamic.F:14
smumps_dynamic_memory_m::smumps_dm_free_block
subroutine smumps_dm_free_block(xxg_status, dynptr, sizfr8, atomic_updates, keep8)
Definition sfac_mem_dynamic.F:433
smumps_dynamic_memory_m::smumps_dm_set_ptr
subroutine smumps_dm_set_ptr(address, sizfr8, cbptr)
Definition sfac_mem_dynamic.F:417
smumps_fac_asm_master_m
Definition sfac_asm_master_m.F:14
smumps_fac_asm_master_m::smumps_fac_asm_niv1
subroutine smumps_fac_asm_niv1(comm_load, ass_irecv, n, inode, iw, liw, a, la, info, nd, fils, frere, dad, maxfrw, root, opassw, opeliw, ptrist, ptlust, ptrfac, ptrast, step, pimaster, pamaster, ptrarw, ptraiw, itloc, rhs_mumps, nsteps, son_level2, comp, lrlu, iptrlu, iwpos, iwposcb, posfac, lrlus, lrlusm, icntl, keep, keep8, dkeep, intarr, lintarr, dblarr, ldblarr nstk_s, procnode_steps, slavef, comm, myid, bufr, lbufr, lbufr_bytes, nbfin, ipool, lpool, leaf, perm, istep_to_iniv2, tab_pos_in_pere, jobass, etatass, lrgroups, mumps_tps_arr, smumps_tps_arr, l0_omp_mapping)
Definition sfac_asm_master_m.F:32
smumps_load
Definition smumps_load.F:14
smumps_load::smumps_load_mem_update
subroutine, public smumps_load_mem_update(ssarbr, process_bande_arg, mem_value, new_lu, inc_mem_arg, keep, keep8, lrlus)
Definition smumps_load.F:920
smumps_load::mumps_load_disable
subroutine, public mumps_load_disable()
Definition smumps_load.F:135
smumps_load::smumps_load_update
subroutine, public smumps_load_update(check_flops, process_bande, inc_load, keep, keep8)
Definition smumps_load.F:824
smumps_load::mumps_load_enable
subroutine, public mumps_load_enable()
Definition smumps_load.F:130
smumps_lr_core
Definition slr_core.F:18
smumps_lr_core::max_cluster
subroutine max_cluster(cut, cut_size, maxi_cluster)
Definition slr_core.F:1304
smumps_lr_core::smumps_blr_asm_niv1
subroutine smumps_blr_asm_niv1(a, la, posel1, nfront, nass1, iwhandler, son_iw, liw, lstk, nelim, k1, k2, sym, keep, keep8, opassw)
Definition slr_core.F:1400
smumps_lr_core::is_front_blr_candidate
subroutine is_front_blr_candidate(inode, niv, nfront, nass, blron, k489, k490, k491, k492, k20, k60, idad, k38, lrstatus, n, lrgroups)
Definition slr_core.F:45
smumps_lr_data_m
Definition smumps_lr_data_m.F:14
smumps_lr_data_m::smumps_blr_init_front
subroutine, public smumps_blr_init_front(iwhandler, info, mtk405)
Definition smumps_lr_data_m.F:175
smumps_lr_data_m::smumps_blr_save_nfs4father
subroutine, public smumps_blr_save_nfs4father(iwhandler, nfs4father)
Definition smumps_lr_data_m.F:856
smumps_struc_def
Definition smumps_struc_def.F:14
smumps_tps_m
Definition somp_tps_m.F:14
comp
int comp(int a, int b)
Definition set_surface_lines.cpp:102
smumps_ldlt_asm_niv12
subroutine smumps_ldlt_asm_niv12(a, la, son_a, iafath, nfront, nass1, ncols, lcb, iw, nrows, nelim, etatass, cb_is_compressed)
Definition sfac_asm.F:406
smumps_ldlt_asm_niv12_ip
subroutine smumps_ldlt_asm_niv12_ip(a, la, iafath, nfront, nass1, iacb, ncols, lcb, iw, nrows, nelim, etatass, cb_is_compressed)
Definition sfac_asm.F:318
smumps_parpivt1_set_nvschur_max
subroutine smumps_parpivt1_set_nvschur_max(n, inode, iw, liw, a, la, keep, perm, ioldps, poselt, nfront, nass1, lr_activated, parpiv_t1, nb_postponed)
Definition sfac_asm.F:950
smumps_set_parpivt1
subroutine smumps_set_parpivt1(inode, nfront, nass1, keep, lr_activated, parpiv_t1)
Definition sfac_asm.F:788
smumps_get_size_needed
subroutine smumps_get_size_needed(sizei_needed, sizer_needed, skip_top_stack, keep, keep8, n, iw, liw, a, la, lrlu, iptrlu, iwpos, iwposcb, ptrist, ptrast, step, pimaster, pamaster, lrlus, xsize, comp, acc_time, myid, slavef, procnode_steps, dad, iflag, ierror)
Definition sfac_mem_compress_cb.F:437
smumps_compre_new
subroutine smumps_compre_new(n, keep, iw, liw, a, la, lrlu, iptrlu, iwpos, iwposcb, ptrist, ptrast, step, pimaster, pamaster, lrlus, xsize, comp, acc_time, myid, slavef, procnode_steps, dad)
Definition sfac_mem_compress_cb.F:120
smumps_free_block_cb_static
subroutine smumps_free_block_cb_static(ssarbr, myid, n, iposblock, iw, liw, lrlu, lrlus, iptrlu, iwposcb, la, keep, keep8, in_place_stats)
Definition sfac_mem_free_block_cb.F:20
smumps_maplig
recursive subroutine smumps_maplig(comm_load, ass_irecv, bufr, lbufr, lbufr_bytes inode_pere, ison, nslaves_pere, list_slaves_pere, nfront_pere, nass_pere, nfs4father, lmap, trow, procnode_steps, slavef, posfac, iwpos, iwposcb, iptrlu, lrlu, lrlus, n, iw, liw, a, la, ptrist, ptlust, ptrfac, ptrast, step, pimaster, pamaster, nstk, comp, iflag, ierror, myid, comm, perm, ipool, lpool, leaf, nbfin, icntl, keep, keep8, dkeep, root, opassw, opeliw, itloc, rhs_mumps, fils, dad, ptrarw, ptraiw, intarr, dblarr, nd, frere, lptrar, nelt, frtptr, frtelt, istep_to_iniv2, tab_pos_in_pere, lrgroups)
Definition sfac_process_maprow.F:34
smumps_try_recvtreat
recursive subroutine smumps_try_recvtreat(comm_load, ass_irecv, blocking, set_irecv, message_received, msgsou, msgtag, status, bufr, lbufr, lbufr_bytes, procnode_steps, posfac, iwpos, iwposcb, iptrlu, lrlu, lrlus, n, iw, liw, a, la, ptrist, ptlust, ptrfac, ptrast, step, pimaster, pamaster, nstk_s, comp, iflag, ierror, comm, perm, ipool, lpool, leaf, nbfin, myid, slavef root, opassw, opeliw, itloc, rhs_mumps, fils, dad, ptrarw, ptraiw, intarr, dblarr, icntl, keep, keep8, dkeep, nd, frere, lptrar, nelt, frtptr, frtelt istep_to_iniv2, tab_pos_in_pere, stack_right_authorized, lrgroups)
Definition sfac_process_message.F:560
smumps_ooc_get_pp_sizes
subroutine smumps_ooc_get_pp_sizes(k50, nbrow_l, nbcol_u, nass, nbpanels_l, nbpanels_u, lreq)
Definition sooc_panel_piv.F:303
smumps_ooc_pp_set_ptr
subroutine smumps_ooc_pp_set_ptr(k50, nbpanels_l, nbpanels_u, nass, ipos, iw, liw)
Definition sooc_panel_piv.F:195
smumps_compute_estim_nfs4father
subroutine smumps_compute_estim_nfs4father(n, inode, ifath, fils, perm, keep, ioldps, hf, iw, liw, nfront, nass1, estim_nfs4father_atson)
Definition stools.F:1612
mumps_in_or_root_ssarbr
logical function mumps_in_or_root_ssarbr(procinfo_inode, k199)
Definition tools_common.F:214
mumps_storei8
subroutine mumps_storei8(i8, int_array)
Definition tools_common.F:871
mumps_inssarbr
logical function mumps_inssarbr(procinfo_inode, k199)
Definition tools_common.F:201
mumps_typenode
integer function mumps_typenode(procinfo_inode, k199)
Definition tools_common.F:90
mumps_procnode
integer function mumps_procnode(procinfo_inode, k199)
Definition tools_common.F:151
mumps_geti8
subroutine mumps_geti8(i8, int_array)
Definition tools_common.F:888