OpenRadioss 2025.1.11
OpenRadioss project
Loading...
Searching...
No Matches
sfac_front_LDLT_type1.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_fac1_ldlt_m
15 CONTAINS
16 SUBROUTINE smumps_fac1_ldlt( N, INODE, IW, LIW, A, LA,
17 & IOLDPS, POSELT, IFLAG, IERROR,
18 & UU, NNEGW, NPVW, NB22T1W, NBTINYW,
19 & DET_EXPW, DET_MANTW, DET_SIGNW,
20 & KEEP,KEEP8,
21 & MYID, SEUIL, AVOID_DELAYED, ETATASS,
22 & DKEEP,PIVNUL_LIST,LPN_LIST, IWPOS
23 & , LRGROUPS
24 & , PERM
25 & )
26 USE smumps_fac_front_aux_m
27 USE smumps_ooc
28 USE smumps_fac_lr
29 USE smumps_lr_type
30 USE smumps_lr_stats
31 USE smumps_ana_lr, ONLY : get_cut
32 USE smumps_lr_data_m
33#if defined(BLR_MT)
34#endif
35!$ USE OMP_LIB
36 IMPLICIT NONE
37 INTEGER(8) :: LA, POSELT
38 INTEGER N, INODE, LIW, IFLAG, IERROR
39 INTEGER, intent(inout) :: NNEGW, NPVW, NB22T1W, NBTINYW
40 INTEGER, intent(inout) :: DET_EXPW, DET_SIGNW
41 REAL, intent(inout) :: DET_MANTW
42 INTEGER MYID, IOLDPS
43 INTEGER KEEP( 500 )
44 INTEGER(8) KEEP8(150)
45 REAL UU, SEUIL
46 REAL A( LA )
47 INTEGER, TARGET :: IW( LIW )
48 INTEGER, intent(in) :: PERM(N)
49 LOGICAL AVOID_DELAYED
50 INTEGER ETATASS, IWPOS
51 INTEGER LPN_LIST
52 INTEGER PIVNUL_LIST(LPN_LIST)
53 REAL DKEEP(230)
54 INTEGER :: LRGROUPS(N)
55 INTEGER INOPV, IFINB, NFRONT, NPIV, IBEG_BLOCK, IEND_BLOCK
56 INTEGER NASS, NBKJIB_ORIG, XSIZE
57 INTEGER :: LDA
58 REAL UUTEMP
59 LOGICAL STATICMODE
60 REAL SEUIL_LOC
61 LOGICAL IS_MAXFROMM_AVAIL
62 INTEGER PIVOT_OPTION
63 INTEGER LRTRSM_OPTION
64 INTEGER LAST_ROW, FIRST_ROW
65 REAL MAXFROMM
66 INTEGER(8) :: LAFAC
67 INTEGER LIWFAC, STRAT, NextPiv2beWritten, IFLAG_OOC,
68 & idummy, pp_first2swap_l, pp_lastpivrptrfilled
69 TYPE(io_block) :: MonBloc
70 LOGICAL LAST_CALL
71 INTEGER PARPIV_T1, OFFSET
72 INTEGER NFS4FATHER
73 REAL, ALLOCATABLE, DIMENSION(:) :: M_ARRAY
74 LOGICAL LASTBL
75 INTEGER CURRENT_BLR
76 LOGICAL LR_ACTIVATED
77 LOGICAL COMPRESS_CB, COMPRESS_PANEL
78 LOGICAL OOCWRITE_COMPATIBLE_WITH_BLR, OOC_EFFECTIVE_ON_FRONT,
79 & ooc_eff_and_write_bypanel
80 INTEGER K473_LOC
81 INTEGER INFO_TMP(2), MAXI_RANK
82 INTEGER FIRST_BLOCK, LAST_BLOCK
83 INTEGER HF, NPARTSASS, NPARTSCB, NB_BLR
84 INTEGER MAXI_CLUSTER, LWORK, NELIM, NELIM_LOC
85 TYPE(lrb_type), POINTER, DIMENSION(:,:) :: CB_LRB
86 INTEGER, POINTER, DIMENSION(:) :: PTDummy
87 TYPE(lrb_type), POINTER, DIMENSION(:) :: ACC_LUA
88 INTEGER, POINTER, DIMENSION(:) :: BEGS_BLR
89 TYPE(lrb_type), POINTER, DIMENSION(:) :: BLR_L
90 REAL, POINTER, DIMENSION(:) :: DIAG
91 INTEGER, POINTER, DIMENSION(:) :: BEGS_BLR_TMP
92 TYPE(lrb_type), POINTER, DIMENSION(:) :: BLR_PANEL
93 INTEGER :: DIAGSIZ_STA, DIAGSIZ_DYN, DIAGPOS, I, IP, MEM, MEM_TOT
94 INTEGER(8) :: POSELT_DIAG
95 REAL, ALLOCATABLE :: WORK(:), TAU(:)
96 INTEGER, ALLOCATABLE :: JPVT(:)
97 REAL,ALLOCATABLE :: RWORK(:)
98 REAL, ALLOCATABLE :: BLOCK(:,:)
99 INTEGER :: allocok,J
100 INTEGER :: OMP_NUM
101 INTEGER :: II,JJ
102 INTEGER(8) :: UPOS, LPOS, DPOS
103 REAL :: ONE, MONE, ZERO
104 parameter(one = 1.0e0, mone=-1.0e0)
105 parameter(zero=0.0e0)
106 INTEGER :: MY_NUM
107 TYPE(lrb_type), POINTER, DIMENSION(:) :: NEXT_BLR_L
108 INTEGER, POINTER, DIMENSION(:) :: BEGS_BLR_STATIC
109 INTEGER :: NVSCHUR, NVSCHUR_K253, IROW_L
110 include 'mumps_headers.h'
111 INTEGER NBLR_ORIG, IBEG_BLR, IEND_BLR
112 INTEGER Inextpiv
113 INTEGER PIVSIZ,IWPOSP2
114 is_maxfromm_avail = .false.
115 IF (keep(206).GE.1) THEN
116 inextpiv = 1
117 ELSE
118 inextpiv = 0
119 ENDIF
120 inopv = 0
121 IF(keep(97) .EQ. 0) THEN
122 staticmode = .false.
123 ELSE
124 staticmode = .true.
125 ENDIF
126 uutemp=uu
127 IF (avoid_delayed) THEN
128 staticmode = .true.
129 seuil_loc = max(seuil,epsilon(seuil))
130 ELSE
131 seuil_loc = seuil
132 ENDIF
133 lafac = -9999_8
134 xsize = keep(ixsz)
135 nfront = iw(ioldps+xsize)
136 lda = nfront
137 nass = iabs(iw(ioldps+2+xsize))
138 iw(ioldps+3+xsize) = -99999
139 lr_activated= .false.
140 compress_panel = .false.
141 compress_cb = .false.
142 NULLIFY(ptdummy)
143 NULLIFY(begs_blr)
144 NULLIFY(cb_lrb)
145 NULLIFY(acc_lua)
146 NULLIFY(blr_l)
147 NULLIFY(begs_blr_tmp)
148 NULLIFY(blr_panel)
149 NULLIFY(diag)
150 compress_panel = (iw(ioldps+xxlr).GE.2)
151 compress_cb = ((iw(ioldps+xxlr).EQ.1).OR.
152 & (iw(ioldps+xxlr).EQ.3))
153 lr_activated = (iw(ioldps+xxlr).GT.0)
154 IF (compress_cb.AND.(.NOT.compress_panel)) THEN
155 compress_panel = .true.
156 k473_loc = 1
157 ELSE
158 k473_loc = keep(473)
159 ENDIF
160 oocwrite_compatible_with_blr =
161 & ( .NOT.lr_activated.OR.(.NOT.compress_panel).OR.
162 & (keep(486).NE.2)
163 & )
164 ooc_effective_on_front= ((keep(201).EQ.1).AND.
165 & oocwrite_compatible_with_blr)
166 CALL smumps_set_parpivt1 ( inode, nfront, nass, keep,
167 & lr_activated, parpiv_t1)
168 lrtrsm_option = keep(475)
169 pivot_option = keep(468)
170 IF (uutemp.EQ.zero) THEN
171 pivot_option = 0
172 ELSE IF (parpiv_t1.NE.0) THEN
173 pivot_option = min(pivot_option,2)
174 ENDIF
175 IF (lr_activated) THEN
176 IF (lrtrsm_option.EQ.3) THEN
177 pivot_option = min(pivot_option,1)
178 ELSEIF (lrtrsm_option.EQ.2) THEN
179 pivot_option = min(pivot_option, 2)
180 ENDIF
181 ENDIF
182 IF (pivot_option.LE.1) THEN
183 parpiv_t1 = 0
184 ENDIF
185 IF (nass.LT.keep(4)) THEN
186 nbkjib_orig = nass
187 ELSE IF (nass .GT. keep(3)) THEN
188 nbkjib_orig = min( keep(6), nass )
189 ELSE
190 nbkjib_orig = min( keep(5), nass )
191 ENDIF
192 IF (.not.lr_activated) THEN
193 nblr_orig = keep(420)
194 ELSE
195 nblr_orig = -9999
196 ENDIF
197 IF ((keep(114).EQ.1) .AND.
198 & (keep(116).GT.0) .AND. ((nfront-nass-keep(253)).GT.0)
199 & ) THEN
200 irow_l = ioldps+6+xsize+nass
201 CALL smumps_get_size_schur_in_front (
202 & n,
203 & nfront-nass-keep(253),
204 & keep(116),
205 & iw(irow_l), perm,
206 & nvschur )
207 ELSE
208 nvschur = 0
209 ENDIF
210 iend_block = 0
211 iend_blr = 0
212 current_blr = 0
213 lastbl = .false.
214 CALL mumps_geti8(lafac,iw(ioldps+xxr))
215 liwfac = iw(ioldps+xxi)
216 IF (ooc_effective_on_front) THEN
217 idummy = -8765
218 nextpiv2bewritten = 1
219 pp_first2swap_l = nextpiv2bewritten
220 monbloc%LastPanelWritten_L = 0
221 pp_lastpivrptrfilled = 0
222 monbloc%INODE = inode
223 monbloc%MASTER = .true.
224 monbloc%Typenode = 1
225 monbloc%NROW = nfront
226 monbloc%NCOL = nfront
227 monbloc%NFS = nass
228 monbloc%Last = .false.
229 monbloc%LastPiv = -77777
230 monbloc%INDICES =>
231 & iw(ioldps+6+nfront+xsize:
232 & ioldps+5+nfront+xsize+nfront)
233 ENDIF
234 IF (lr_activated) THEN
235 IF (keep(405) .EQ. 1) THEN
236!$OMP ATOMIC UPDATE
237 cnt_nodes = cnt_nodes + 1
238!$OMP END ATOMIC
239 ELSE
240 cnt_nodes = cnt_nodes + 1
241 ENDIF
242 ELSE IF (keep(486).NE.0) THEN
243 ENDIF
244 ooc_eff_and_write_bypanel = ( (pivot_option.GE.3) .AND.
245 & ooc_effective_on_front )
246 hf = 6 + iw(ioldps+5+xsize)+xsize
247 IF (lr_activated) THEN
248 CALL get_cut(iw(ioldps+hf:ioldps+hf+nfront-1), nass,
249 & nfront-nass, lrgroups, npartscb,
250 & npartsass, begs_blr)
251 CALL regrouping2(begs_blr, npartsass, nass, npartscb,
252 & nfront-nass, keep(488), .false., keep(472))
253 nb_blr = npartsass + npartscb
254 call max_cluster(begs_blr,nb_blr,maxi_cluster)
255 maxi_rank = keep(479)*maxi_cluster
256 lwork = maxi_cluster*maxi_cluster
257 omp_num = 1
258#if defined(BLR_MT)
259!$ OMP_NUM = OMP_GET_MAX_THREADS()
260#endif
261 ALLOCATE(block(maxi_cluster, omp_num*maxi_cluster),
262 & rwork(2*maxi_cluster*omp_num),
263 & tau(maxi_cluster*omp_num),
264 & jpvt(maxi_cluster*omp_num),
265 & work(lwork*omp_num),stat=allocok)
266 IF (allocok > 0) THEN
267 iflag = -13
268 ierror = omp_num*(lwork + maxi_cluster*(maxi_cluster+4))
269 GOTO 490
270 ENDIF
271 ALLOCATE(acc_lua(omp_num),stat=allocok)
272 IF (allocok > 0) THEN
273 iflag = -13
274 ierror = omp_num
275 GOTO 490
276 ENDIF
277 IF (keep(480).GE.3) THEN
278 DO my_num=1,omp_num
279 CALL alloc_lrb(acc_lua(my_num), maxi_rank,
280 & maxi_cluster, maxi_cluster, .true.,
281 & iflag, ierror, keep8)
282 IF (iflag.LT.0) GOTO 500
283 acc_lua(my_num)%K = 0
284 ENDDO
285 ENDIF
286 ENDIF
287 IF (lr_activated.AND.(keep(480).NE.0
288 & .OR.
289 & (
290 & (keep(486).EQ.2)
291 & )
292 & .OR.compress_cb
293 & )) THEN
294 info_tmp(1) = iflag
295 info_tmp(2) = ierror
296 IF (iflag.LT.0) GOTO 500
297 CALL smumps_blr_save_init(iw(ioldps+xxf),
298 & .true.,
299 & .false.,
300 & .false.,
301 & npartsass,
302 & begs_blr, ptdummy,
303 & huge(npartsass),
304 & info_tmp)
305 iflag = info_tmp(1)
306 ierror = info_tmp(2)
307 IF (iflag.LT.0) GOTO 500
308 ENDIF
309 IF (compress_cb.AND.npartscb.GT.0) THEN
310 allocate(cb_lrb(npartscb,npartscb),stat=allocok)
311 IF (allocok > 0) THEN
312 iflag = -13
313 ierror = npartscb*npartscb
314 GOTO 490
315 ENDIF
316 DO ii=1,npartscb
317 DO jj=1,npartscb
318 cb_lrb(ii,jj)%M=0
319 cb_lrb(ii,jj)%N=0
320 NULLIFY(cb_lrb(ii,jj)%Q)
321 NULLIFY(cb_lrb(ii,jj)%R)
322 cb_lrb(ii,jj)%ISLR = .false.
323 ENDDO
324 ENDDO
325 CALL smumps_blr_save_cb_lrb(iw(ioldps+xxf),cb_lrb)
326 ENDIF
327 DO WHILE (iend_blr < nass )
328 current_blr = current_blr + 1
329 ibeg_blr = iw(ioldps+1+keep(ixsz)) + 1
330 IF (.NOT. lr_activated) THEN
331 iend_blr = min(iend_blr + nblr_orig, nass)
332 ELSE
333 iend_blr = begs_blr(current_blr+1)-1
334 begs_blr( current_blr ) = ibeg_blr
335 IF ( iend_blr - ibeg_blr + 1 .GT. maxi_cluster ) THEN
336 maxi_cluster = iend_blr - ibeg_blr + 1
337 lwork = maxi_cluster*maxi_cluster
338 DEALLOCATE(block, work, rwork, tau, jpvt)
339 ALLOCATE(block(maxi_cluster, omp_num*maxi_cluster),
340 & rwork(2*maxi_cluster*omp_num),
341 & tau(maxi_cluster*omp_num),
342 & jpvt(maxi_cluster*omp_num),
343 & work(lwork*omp_num),stat=allocok)
344 IF (allocok > 0) THEN
345 iflag = -13
346 ierror = omp_num*(lwork + maxi_cluster*(maxi_cluster+4))
347 GOTO 490
348 ENDIF
349 IF (keep(480).GE.3) THEN
350 DO my_num=1,omp_num
351 CALL dealloc_lrb(acc_lua(my_num), keep8, keep(34))
352 CALL alloc_lrb(acc_lua(my_num), maxi_rank,
353 & maxi_cluster, maxi_cluster, .true.,
354 & iflag, ierror, keep8)
355 IF (iflag.LT.0) GOTO 500
356 acc_lua(my_num)%K = 0
357 ENDDO
358 ENDIF
359 ENDIF
360 IF (keep(480).GE.5) THEN
361 IF (current_blr.EQ.1) THEN
362 ALLOCATE(blr_l(nb_blr-current_blr),stat=allocok)
363 IF (allocok > 0) THEN
364 iflag = -13
365 ierror = nb_blr-current_blr
366 GOTO 490
367 ENDIF
368 IF (nb_blr.GT.current_blr) THEN
369 blr_l(1:nb_blr-current_blr)%ISLR=.false.
370 CALL smumps_blr_save_panel_loru (
371 & iw(ioldps+xxf),
372 & 0,
373 & current_blr, blr_l)
374 ENDIF
375 ELSE
376 IF (nb_blr.GT.current_blr) THEN
377 CALL smumps_blr_retrieve_panel_loru(
378 & iw(ioldps+xxf),
379 & 0,
380 & current_blr, blr_l)
381 ENDIF
382 ENDIF
383 IF (current_blr.LT.npartsass) THEN
384 ALLOCATE(next_blr_l(nb_blr-current_blr-1),stat=allocok)
385 IF (allocok > 0) THEN
386 iflag = -13
387 ierror = nb_blr-current_blr-1
388 GOTO 490
389 ENDIF
390 IF (nb_blr.GT.current_blr+1) THEN
391 CALL smumps_blr_save_panel_loru (
392 & iw(ioldps+xxf),
393 & 0,
394 & current_blr+1, next_blr_l)
395 ENDIF
396 ENDIF
397 ELSE
398 ALLOCATE(blr_l(nb_blr-current_blr),stat=allocok)
399 IF (allocok > 0) THEN
400 iflag = -13
401 ierror = nb_blr-current_blr
402 GOTO 490
403 ENDIF
404 ENDIF
405 ENDIF
406 IF (lr_activated) THEN
407 ENDIF
408 DO WHILE (iend_block < iend_blr )
409 ibeg_block = iw(ioldps+1+keep(ixsz)) + 1
410 IF (keep(405).EQ.0) THEN
411 keep(425)=max(keep(425),iend_block-ibeg_block)
412 ELSE
413!$OMP ATOMIC UPDATE
414 keep(425)=max(keep(425),iend_block-ibeg_block)
415!$OMP END ATOMIC
416 ENDIF
417 iend_block = min(iend_block + nbkjib_orig, iend_blr)
418 50 CONTINUE
419 CALL smumps_fac_i_ldlt(nfront,nass,inode,
420 & ibeg_block, iend_block,
421 & iw,liw,a,la,
422 & inopv, nnegw, nb22t1w, nbtinyw,
423 & det_expw, det_mantw, det_signw,
424 & iflag,ioldps,poselt,uutemp,
425 & seuil_loc,keep,keep8,pivsiz,
426 & dkeep(1),pivnul_list(1),lpn_list, xsize,
427 & pp_first2swap_l, monbloc%LastPanelWritten_L,
428 & pp_lastpivrptrfilled, maxfromm, is_maxfromm_avail,
429 & pivot_option, iend_blr, inextpiv,
430 & ooc_effective_on_front,
431 & nvschur, parpiv_t1, lr_activated
432 & )
433 IF (iflag.LT.0) GOTO 500
434 IF (inopv.EQ.1) THEN
435 IF(staticmode) THEN
436 inopv = -1
437 GOTO 50
438 ENDIF
439 lastbl = .true.
440 ELSE IF ( inopv.LE.0 ) THEN
441 inopv = 0
442 npvw = npvw + pivsiz
443 nvschur_k253 = 0
444 IF (pivot_option.GE.3) THEN
445 last_row = nfront
446 nvschur_k253 = nvschur + keep(253)
447 ELSEIF (pivot_option.EQ.2) THEN
448 last_row = nass
449 ELSE
450 last_row = iend_blr
451 ENDIF
452 CALL smumps_fac_mq_ldlt(iend_block,
453 & nfront, nass, iw(ioldps+1+xsize),
454 & inode,a,la,
455 & lda,
456 & poselt,ifinb,
457 & pivsiz, maxfromm,
458 & is_maxfromm_avail, (uutemp.NE.0.0e0),
459 & parpiv_t1,
460 & last_row, iend_blr, nvschur_k253,
461 & lr_activated
462 & )
463 IF(pivsiz .EQ. 2) THEN
464 iwposp2 = ioldps+iw(ioldps+1+xsize)+6
465 iw(iwposp2+nfront+xsize) =
466 & -iw(iwposp2+nfront+xsize)
467 ENDIF
468 iw(ioldps+1+xsize) = iw(ioldps+1+xsize) + pivsiz
469 IF (ifinb.EQ.0) THEN
470 GOTO 50
471 ELSE IF (ifinb.EQ.-1) THEN
472 lastbl = .true.
473 ENDIF
474 ENDIF
475 IF ( ooc_eff_and_write_bypanel ) THEN
476 monbloc%Last = lastbl
477 monbloc%LastPiv= iw(ioldps+1+xsize)
478 last_call=.false.
479 CALL smumps_ooc_io_lu_panel(
480 & strat_try_write,
481 & typef_l, a(poselt),
482 & lafac, monbloc, nextpiv2bewritten, idummy,
483 & iw(ioldps), liwfac,
484 & myid, keep8(31), iflag_ooc,last_call )
485 IF (iflag_ooc < 0 ) THEN
486 iflag=iflag_ooc
487 GOTO 500
488 ENDIF
489 ENDIF
490 npiv = iw(ioldps+1+xsize)
491 IF ( iend_blr .GT. iend_block ) THEN
492 IF (pivot_option.GE.3) THEN
493 last_row = nfront
494 ELSEIF (pivot_option.EQ.2) THEN
495 last_row = nass
496 ELSE
497 last_row = iend_blr
498 ENDIF
499 CALL smumps_fac_sq_ldlt(ibeg_block,iend_block,
500 & npiv, nfront,nass,inode,a,la,
501 & lda, poselt,
502 & keep, keep8,
503 & -6666, -6666,
504 & iend_blr, last_row,
505 & .false., .true., lr_activated,
506 & iw, liw, -6666
507 & )
508 ENDIF
509 END DO
510 npiv = iw(ioldps+1+xsize)
511 IF (.NOT. lr_activated
512 & .OR. (.NOT. compress_panel)
513 & ) THEN
514 IF (pivot_option.GE.3) THEN
515 last_row = nfront
516 ELSEIF (pivot_option.EQ.2) THEN
517 last_row = nass
518 ELSE
519 last_row = iend_blr
520 ENDIF
521 CALL smumps_fac_sq_ldlt(ibeg_blr,iend_blr,npiv,
522 & nfront,nass,inode,a,la,
523 & lda, poselt,
524 & keep, keep8,
525 & iend_blr, nass,
526 & nass, last_row,
527 & (pivot_option.LE.1), .true., lr_activated,
528 & iw, liw, ioldps+6+xsize+nfront+ibeg_blr-1)
529 ELSE
530 nelim = iend_block - npiv
531 IF (nelim .EQ. iend_blr - ibeg_blr + 1) THEN
532 IF (keep(480).GE.2
533 & .OR.
534 & (
535 & (keep(486).EQ.2)
536 & )
537 & ) THEN
538 DO j=1,nb_blr-current_blr
539 blr_l(j)%M=0
540 blr_l(j)%N=0
541 blr_l(j)%K=0
542 blr_l(j)%ISLR=.false.
543 NULLIFY(blr_l(j)%Q)
544 NULLIFY(blr_l(j)%R)
545 ENDDO
546 CALL smumps_blr_save_panel_loru (
547 & iw(ioldps+xxf),
548 & 0,
549 & current_blr, blr_l)
550 NULLIFY(blr_l)
551 IF (keep(480).GE.2 .AND. iend_blr.LT.nass) THEN
552 IF (lrtrsm_option.EQ.2) THEN
553 first_block = npartsass-current_blr
554 ELSE
555 first_block = 1
556 ENDIF
557#if defined(BLR_MT)
558!$OMP PARALLEL
559#endif
560 CALL smumps_blr_upd_panel_left_ldlt(a, la, poselt,
561 & nfront, iw(ioldps+xxf),
562 & begs_blr, current_blr, nb_blr, npartsass,
563 & nelim,
564 & iw(hf+ioldps+nfront), block,
565 & acc_lua, maxi_cluster, maxi_rank,
566 & 1, iflag, ierror,
567 & keep(481), dkeep(11), keep(466), keep(477),
568 & keep(480), keep(479), keep(478), keep(476),
569 & keep(483), keep8, first_block=first_block)
570#if defined(BLR_MT)
571!$OMP END PARALLEL
572#endif
573 IF (iflag.LT.0) GOTO 500
574 ENDIF
575 ENDIF
576 IF (keep(486).EQ.3) THEN
577 IF (keep(480).EQ.0) THEN
578 DEALLOCATE(blr_l)
579 NULLIFY(blr_l)
580 ENDIF
581 ENDIF
582 GOTO 100
583 ENDIF
584 IF (pivot_option.GE.3) THEN
585 first_row = nfront
586 ELSEIF (pivot_option.EQ.2) THEN
587 first_row = nass
588 ELSE
589 first_row = iend_blr
590 ENDIF
591 IF (lrtrsm_option.EQ.3) THEN
592 last_row = iend_blr
593 ELSEIF (lrtrsm_option.EQ.2) THEN
594 last_row = nass
595 ELSE
596 last_row = nfront
597 ENDIF
598 IF ((iend_blr.LT.nfront) .AND. (last_row-first_row.GT.0)) THEN
599 CALL smumps_fac_sq_ldlt(ibeg_blr, iend_blr,
600 & npiv, nfront, nass,
601 & inode, a, la, lda, poselt,
602 & keep, keep8,
603 & first_row, last_row,
604 & -6666, -6666,
605 & .true., .false., lr_activated,
606 & iw, liw, ioldps+6+xsize+nfront+ibeg_blr-1)
607 ENDIF
608#if defined(BLR_MT)
609#endif
610#if defined(BLR_MT)
611!$omp parallel private(upos,lpos,dpos,offset)
612!$OMP& FIRSTPRIVATE(FIRST_BLOCK,LAST_BLOCK)
613#endif
614 CALL smumps_compress_panel(a, la, poselt, iflag, ierror,
615 & nfront,
616 & begs_blr, nb_blr, dkeep(8), keep(466), k473_loc, blr_l,
617 & current_blr,
618 & 'V', work, tau, jpvt, lwork, rwork,
619 & block, maxi_cluster, nelim,
620 & .false., 0, 0,
621 & 1, keep(483), keep8,
622 & k480=keep(480)
623 & )
624#if defined(BLR_MT)
625!$OMP BARRIER
626#endif
627 IF (iflag.LT.0) GOTO 400
628 IF (pivot_option.LT.3) THEN
629 IF (lrtrsm_option.GE.2) THEN
630 IF (pivot_option.LE.1.AND.lrtrsm_option.EQ.3) THEN
631 first_block = current_blr+1
632 ELSE
633 first_block = npartsass+1
634 ENDIF
635 CALL smumps_blr_panel_lrtrsm(a, la, poselt, nfront,
636 & ibeg_blr, nb_blr, blr_l,
637 & current_blr, first_block, nb_blr,
638 & 1, 1, 0,
639 & .false.,
640 & iw, offset_iw=ioldps+6+xsize+nfront+ibeg_blr-1)
641#if defined(BLR_MT)
642!$OMP BARRIER
643#endif
644 ENDIF
645 IF (nelim.GT.0) THEN
646 IF (pivot_option.LE.1) THEN
647 first_block = current_blr+1
648 ELSE
649 first_block = npartsass+1
650 ENDIF
651 lpos = poselt
652 & +int(begs_blr(current_blr+1)-1-nelim,8)*int(nfront,8)
653 & +int(begs_blr(current_blr)-1,8)
654 dpos = poselt
655 & +int(begs_blr(current_blr)-1,8)*int(nfront,8)
656 & +int(begs_blr(current_blr)-1,8)
657 offset=ioldps+6+xsize+nfront+ibeg_blr-1
658 upos = poselt+int(begs_blr(current_blr)-1,8)*int(nfront,8)
659 & +int(begs_blr(current_blr+1)-1-nelim,8)
660#if defined(BLR_MT)
661!$OMP SINGLE
662#endif
663 CALL smumps_fac_ldlt_copyscale_u( nelim, 1,
664 & keep(424), nfront, npiv-ibeg_blr+1,
665 & liw, iw, offset, la, a, poselt, lpos, upos, dpos)
666#if defined(BLR_MT)
667!$OMP END SINGLE
668#endif
669 lpos = poselt
670 & +int(begs_blr(current_blr+1)-1,8)*int(nfront,8)
671 & +int(begs_blr(current_blr+1)-1-nelim,8)
672 CALL smumps_blr_upd_nelim_var_l(
673 & a, la, upos, a, la, lpos,
674 & iflag, ierror, nfront, nfront,
675 & begs_blr, current_blr, blr_l, nb_blr,
676 & first_block, nelim, 'N')
677 ENDIF
678 ENDIF
679 IF (iflag.LT.0) GOTO 400
680#if defined(BLR_MT)
681!$OMP MASTER
682#endif
683 IF (keep(480).NE.0
684 & .OR.
685 & (
686 & (keep(486).EQ.2)
687 & )
688 & ) THEN
689 IF (keep(480).LT.5) THEN
690 CALL smumps_blr_save_panel_loru (
691 & iw(ioldps+xxf),
692 & 0,
693 & current_blr, blr_l)
694 ENDIF
695 ENDIF
696#if defined(BLR_MT)
697!$OMP END MASTER
698!$OMP BARRIER
699#endif
700 IF (keep(480).GE.2) THEN
701 IF (iend_blr.LT.nass) THEN
702 IF (lrtrsm_option.EQ.2) THEN
703 first_block = npartsass-current_blr
704 ELSE
705 first_block = 1
706 ENDIF
707 CALL smumps_blr_upd_panel_left_ldlt(a, la, poselt,
708 & nfront, iw(ioldps+xxf),
709 & begs_blr, current_blr, nb_blr, npartsass,
710 & nelim,
711 & iw(hf+ioldps+nfront), block,
712 & acc_lua, maxi_cluster, maxi_rank,
713 & 1, iflag, ierror,
714 & keep(481), dkeep(11), keep(466), keep(477),
715 & keep(480), keep(479), keep(478), keep(476),
716 & keep(483), keep8, first_block=first_block)
717 ENDIF
718 ELSE
719 CALL smumps_blr_update_trailing_ldlt(a, la, poselt,
720 & iflag, ierror, nfront,
721 & begs_blr, nb_blr, current_blr, blr_l, nelim,
722 & iw(hf+ioldps+nfront+ibeg_blr-1), block,
723 & maxi_cluster, npiv,
724 & 1,
725 & keep(481), dkeep(11), keep(466), keep(477)
726 & )
727 ENDIF
728#if defined(BLR_MT)
729!$OMP BARRIER
730#endif
731 IF (iflag.LT.0) GOTO 400
732 IF (lrtrsm_option.GE.2) THEN
733 IF (lrtrsm_option.EQ.2) THEN
734 first_block = npartsass+1
735 ELSE
736 first_block = current_blr+1
737 ENDIF
738 IF (keep(486).NE.2) THEN
739 last_block = nb_blr
740 ELSEIF(uu.GT.0) THEN
741 last_block = npartsass
742 ELSE
743 last_block = current_blr
744 ENDIF
745 CALL smumps_decompress_panel(a, la, poselt, nfront, nfront,
746 & .true.,
747 & begs_blr(current_blr),
748 & begs_blr(current_blr+1), nb_blr, blr_l, current_blr, 'V',
749 & 1,
750 & beg_i_in=first_block, end_i_in=last_block)
751 ENDIF
752 400 CONTINUE
753#if defined(blr_mt)
754!$OMP END PARALLEL
755#endif
756 IF (iflag.LT.0) GOTO 500
757 IF (keep(486).EQ.3) THEN
758 IF (keep(480).EQ.0) THEN
759 CALL dealloc_blr_panel(blr_l, nb_blr-current_blr, keep8,
760 & keep(34))
761 DEALLOCATE(blr_l)
762 ELSE
763 NULLIFY(next_blr_l)
764 ENDIF
765 ENDIF
766 NULLIFY(blr_l)
767 ENDIF
768 IF ( ooc_eff_and_write_bypanel ) THEN
769 monbloc%Last = lastbl
770 monbloc%LastPiv= npiv
771 last_call=.false.
772 CALL smumps_ooc_io_lu_panel(
773 & strat_try_write,
774 & typef_l, a(poselt),
775 & lafac, monbloc, nextpiv2bewritten, idummy, iw(ioldps),
776 & liwfac, myid, keep8(31), iflag_ooc,last_call )
777 IF (iflag_ooc < 0 ) THEN
778 iflag=iflag_ooc
779 GOTO 500
780 ENDIF
781 ENDIF
782 100 CONTINUE
783 END DO
784 IF (lr_activated) THEN
785 ibeg_blr = iw(ioldps+1+keep(ixsz)) + 1
786 begs_blr( current_blr + 1 ) = ibeg_blr
787 IF (
788 & (keep(486).EQ.2)
789 & ) THEN
790 CALL smumps_blr_retrieve_begsblr_sta(iw(ioldps+xxf),
791 & begs_blr_static)
792 IF (uu.GT.0) THEN
793 allocate(begs_blr_tmp(nb_blr+1),stat=allocok)
794 IF (allocok > 0) THEN
795 iflag = -13
796 ierror = nb_blr+1
797 GOTO 500
798 ENDIF
799 DO j=1,nb_blr+1
800 begs_blr_tmp(j) = begs_blr_static(j)
801 ENDDO
802 ENDIF
803 ENDIF
804 mem_tot = 0
805#if defined(BLR_MT)
806!$OMP PARALLEL
807!$OMP& PRIVATE(IP, NELIM_LOC, BLR_PANEL)
808#endif
809 IF ( (keep(486).EQ.2)
810 & ) THEN
811#if defined(BLR_MT)
812!$OMP DO PRIVATE(DIAG, DIAGSIZ_STA, DIAGSIZ_DYN, DIAGPOS, POSELT_DIAG,
813!$OMP& MEM, allocok)
814!$OMP& REDUCTION(+:MEM_TOT)
815#endif
816 DO ip=1,npartsass
817 IF (iflag.LT.0) cycle
818 diagsiz_dyn = begs_blr(ip+1)-begs_blr(ip)
819 diagsiz_sta = begs_blr_static(ip+1)-begs_blr(ip)
820 mem = diagsiz_dyn*diagsiz_sta
821 mem_tot = mem_tot + mem
822 ALLOCATE(diag(mem),stat=allocok)
823 IF (allocok > 0) THEN
824 iflag = -13
825 ierror = mem
826 cycle
827 ENDIF
828 diagpos = 1
829 poselt_diag = poselt + int(begs_blr(ip)-1,8)*int(nfront,8)
830 & + int(begs_blr(ip)-1,8)
831 DO i=1,diagsiz_sta
832 diag(diagpos:diagpos+diagsiz_dyn-1) =
833 & a(poselt_diag:poselt_diag+int(diagsiz_dyn-1,8))
834 diagpos = diagpos + diagsiz_dyn
835 poselt_diag = poselt_diag + int(nfront,8)
836 ENDDO
837 CALL smumps_blr_save_diag_block(
838 & iw(ioldps+xxf),
839 & ip, diag)
840 ENDDO
841#if defined(BLR_MT)
842!$OMP ENDDO
843!$OMP SINGLE
844#endif
845 CALL mumps_dm_fac_upd_dyn_memcnts(int(mem_tot,8),
846 & (keep(405).NE.0), keep8, iflag, ierror, .true., .true.)
847#if defined(BLR_MT)
848!$OMP END SINGLE
849#endif
850 IF (iflag.LT.0) GOTO 447
851 IF (uu.GT.0) THEN
852 DO ip=1,npartsass
853 nelim_loc = begs_blr_tmp(ip+1)-begs_blr(ip+1)
854 CALL smumps_blr_retrieve_panel_loru(
855 & iw(ioldps+xxf), 0, ip, blr_panel)
856#if defined(BLR_MT)
857!$OMP SINGLE
858#endif
859 CALL dealloc_blr_panel(blr_panel, npartsass-ip, keep8,
860 & keep(34))
861#if defined(BLR_MT)
862!$OMP END SINGLE
863#endif
864 CALL smumps_compress_panel(a, la, poselt, iflag,
865 & ierror, nfront, begs_blr_tmp,
866 & nb_blr, dkeep(8), keep(466), k473_loc,
867 & blr_panel, ip,
868 & 'V', work, tau, jpvt, lwork, rwork,
869 & block, maxi_cluster, nelim_loc,
870 & .false., 0, 0,
871 & 1, keep(483), keep8,
872 & end_i_in=npartsass, frswap=.true.
873 & )
874#if defined(BLR_MT)
875!$OMP BARRIER
876#endif
877 IF (iflag.LT.0) GOTO 445
878#if defined(BLR_MT)
879!$OMP SINGLE
880#endif
881 begs_blr_tmp(ip+1) = begs_blr(ip+1)
882#if defined(BLR_MT)
883!$OMP END SINGLE
884#endif
885 ENDDO
886#if defined(BLR_MT)
887!$OMP BARRIER
888#endif
889 445 CONTINUE
890 ENDIF
891 447 CONTINUE
892 ENDIF
893 IF (keep(480) .GE. 2) THEN
894#if defined(BLR_MT)
895!$OMP SINGLE
896#endif
897 CALL smumps_blr_retrieve_begsblr_sta(iw(ioldps+xxf),
898 & begs_blr_static)
899#if defined(BLR_MT)
900!$OMP END SINGLE
901#endif
902 CALL smumps_blr_upd_cb_left_ldlt(a, la, poselt, nfront,
903 & begs_blr_static, begs_blr, npartscb, npartsass, nass,
904 & iw(ioldps+xxf),
905 & iw(hf+ioldps+nfront), block,
906 & acc_lua, maxi_cluster, maxi_rank,
907 & 1, iflag, ierror,
908 & keep(481), dkeep(11), keep(466), keep(477),
909 & keep(480), keep(479), keep(478), keep(476),
910 & keep(484), keep8)
911#if defined(BLR_MT)
912!$OMP BARRIER
913#endif
914 END IF
915 IF (iflag.LT.0) GOTO 450
916#if defined(BLR_MT)
917!$OMP MASTER
918#endif
919 IF (compress_cb
920 & .OR.
921 & (
922 & (keep(486).EQ.2)
923 & )
924 & ) THEN
925 CALL smumps_blr_save_begs_blr_dyn(iw(ioldps+xxf),
926 & begs_blr)
927 ENDIF
928 IF (compress_cb) THEN
929 iend_blr = begs_blr(current_blr+2)
930 IF ( iend_blr - ibeg_blr + 1 .GT. maxi_cluster ) THEN
931 maxi_cluster = iend_blr - ibeg_blr + 1
932 lwork = maxi_cluster*maxi_cluster
933 DEALLOCATE(block, work, rwork, tau, jpvt)
934 ALLOCATE(block(maxi_cluster, omp_num*maxi_cluster),
935 & rwork(2*maxi_cluster*omp_num),
936 & tau(maxi_cluster*omp_num),
937 & jpvt(maxi_cluster*omp_num),
938 & work(lwork*omp_num),stat=allocok)
939 IF (allocok > 0) THEN
940 iflag = -13
941 ierror = omp_num*(lwork + maxi_cluster*(maxi_cluster+4))
942 ENDIF
943 ENDIF
944 ENDIF
945#if defined(BLR_MT)
946!$OMP END MASTER
947!$OMP BARRIER
948#endif
949 IF (iflag.LT.0) GOTO 450
950 IF (compress_cb) THEN
951#if defined(BLR_MT)
952!$OMP MASTER
953#endif
954 nfs4father = -9999
955 IF ( (keep(219).NE.0).AND.(keep(50).EQ.2) ) THEN
956 CALL smumps_blr_retrieve_nfs4father ( iw(ioldps+xxf),
957 & nfs4father )
958 IF (nfs4father.GE.0) nfs4father = nfs4father + nelim
959 ENDIF
960 ALLOCATE(m_array(max(nfs4father,1)), stat=allocok)
961 IF ( allocok.GT.0 ) THEN
962 iflag = -13
963 ierror = max(nfs4father,1)
964 ENDIF
965#if defined(BLR_MT)
966!$omp END master
967!$OMP BARRIER
968#endif
969 IF (iflag.LT.0) GOTO 448
970 CALL smumps_compress_cb(a, la, poselt, nfront,
971 & begs_blr, begs_blr, npartscb, npartscb, npartsass,
972 & nfront-nass, nfront-nass, inode,
973 & iw(ioldps+xxf), 2, 1, iflag, ierror,
974 & dkeep(12), keep(466), keep(484), keep(489), cb_lrb,
975 & work, tau, jpvt, lwork, rwork, block,
976 & maxi_cluster, keep8,
977 & nfs4father, npiv, nvschur+keep(253), keep(1),
978 & m_array=m_array,
979 & nelim=nelim )
980#if defined(BLR_MT)
981!$OMP BARRIER
982#endif
983 IF (iflag.LT.0) GOTO 448
984#if defined(BLR_MT)
985!$OMP MASTER
986#endif
987 IF ( (keep(219).NE.0).AND.(keep(50).EQ.2).AND.
988 & nfs4father.GT.0 ) THEN
989 info_tmp(1) = iflag
990 info_tmp(2) = ierror
991 CALL smumps_blr_save_m_array( iw(ioldps+xxf),
992 & m_array, info_tmp)
993 iflag = info_tmp(1)
994 ierror = info_tmp(2)
995 ENDIF
996 DEALLOCATE(m_array)
997#if defined(BLR_MT)
998!$OMP END MASTER
999!$OMP BARRIER
1000#endif
1001 448 CONTINUE
1002 ENDIF
1003 450 CONTINUE
1004#if defined(BLR_MT)
1005!$OMP END PARALLEL
1006#endif
1007 IF ( (
1008 & (keep(486).EQ.2)
1009 & )
1010 & .AND.uu.GT.0) THEN
1011 deallocate(begs_blr_tmp)
1012 ENDIF
1013 IF (iflag.LT.0) GOTO 500
1014 CALL upd_mry_lu_fr(nass, nfront-nass, 1, nass-npiv)
1015 DO ip=1,npartsass
1016 CALL smumps_blr_retrieve_panel_loru(
1017 & iw(ioldps+xxf), 0, ip, blr_panel)
1018 CALL upd_mry_lu_lrgain(blr_panel, nb_blr-ip
1019 & )
1020 ENDDO
1021 CALL upd_flop_facto_fr(nfront, nass, npiv, 2, 1)
1022 ENDIF
1023 IF (.NOT. compress_panel) THEN
1024 CALL smumps_fac_t_ldlt(nfront,nass,iw,liw,a,la,
1025 & lda, ioldps,poselt, keep,keep8,
1026 & (pivot_option.NE.3), etatass,
1027 & typef_l, lafac, monbloc, nextpiv2bewritten,
1028 & liwfac, myid, iflag, ioldps+6+xsize+nfront, inode )
1029 ENDIF
1030 IF (keep(486).NE.0) THEN
1031 IF (.NOT.lr_activated) THEN
1032 CALL upd_flop_frfronts(nfront, npiv, nass, 1, 1)
1033 ENDIF
1034 ENDIF
1035 IF (ooc_effective_on_front) THEN
1036 strat = strat_write_max
1037 monbloc%Last = .true.
1038 monbloc%LastPiv = iw(ioldps+1+xsize)
1039 last_call = .true.
1040 CALL smumps_ooc_io_lu_panel
1041 & ( strat, typef_l,
1042 & a(poselt), lafac, monbloc,
1043 & nextpiv2bewritten, idummy,
1044 & iw(ioldps), liwfac,
1045 & myid, keep8(31), iflag_ooc,last_call )
1046 IF (iflag_ooc < 0 ) THEN
1047 iflag=iflag_ooc
1048 GOTO 500
1049 ENDIF
1050 CALL smumps_ooc_pp_tryrelease_space (iwpos,
1051 & ioldps, iw, liw, monbloc , nfront, keep)
1052 ENDIF
1053 GOTO 600
1054 490 CONTINUE
1055 500 CONTINUE
1056 600 CONTINUE
1057 IF (lr_activated) THEN
1058 IF (allocated(rwork)) DEALLOCATE(rwork)
1059 IF (allocated(work)) DEALLOCATE(work)
1060 IF (allocated(tau)) deallocate(tau)
1061 IF (allocated(jpvt)) deallocate(jpvt)
1062 IF (allocated(block)) deallocate(block)
1063 IF (associated(acc_lua)) THEN
1064 IF (keep(480).GE.3) THEN
1065 DO my_num=1,omp_num
1066 CALL dealloc_lrb(acc_lua(my_num), keep8, keep(34))
1067 ENDDO
1068 ENDIF
1069 DEALLOCATE(acc_lua)
1070 ENDIF
1071 IF (associated(begs_blr)) THEN
1072 DEALLOCATE(begs_blr)
1073 NULLIFY(begs_blr)
1074 ENDIF
1075 ENDIF
1076 IF (lr_activated.AND.keep(480).NE.0) THEN
1077 IF (.NOT.
1078 & (
1079 & (keep(486).EQ.2)
1080 & )
1081 & ) THEN
1082 CALL smumps_blr_free_all_panels(iw(ioldps+xxf), 0,
1083 & keep8, keep(34))
1084 ENDIF
1085 ENDIF
1086 IF (lr_activated) THEN
1087 IF (.NOT.
1088 & (
1089 & (keep(486).EQ.2)
1090 & )
1091 & .AND. .NOT.compress_cb) THEN
1092 CALL smumps_blr_end_front(iw(ioldps+xxf), iflag, keep8,
1093 & keep(34), mtk405=keep(405))
1094 ENDIF
1095 ENDIF
1096 RETURN
1097 END SUBROUTINE smumps_fac1_ldlt
1098 END MODULE smumps_fac1_ldlt_m
min
#define min(a, b)
Definition macros.h:20
max
#define max(a, b)
Definition macros.h:21
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_fac1_ldlt_m
Definition sfac_front_LDLT_type1.F:14
smumps_fac1_ldlt_m::smumps_fac1_ldlt
subroutine smumps_fac1_ldlt(n, inode, iw, liw, a, la, ioldps, poselt, iflag, ierror, uu, nnegw, npvw, nb22t1w, nbtinyw, det_expw, det_mantw, det_signw, keep, keep8, myid, seuil, avoid_delayed, etatass, dkeep, pivnul_list, lpn_list, iwpos, lrgroups, perm)
Definition sfac_front_LDLT_type1.F:26
smumps_fac_front_aux_m
Definition sfac_front_aux.F:14
smumps_fac_front_aux_m::smumps_fac_i_ldlt
subroutine smumps_fac_i_ldlt(nfront, nass, inode, ibeg_block, iend_block, iw, liw, a, la, inopv, nnegw, nb22t1w, nbtinyw, det_expw, det_mantw, det_signw, iflag, ioldps, poselt, uu, seuil, keep, keep8, pivsiz, dkeep, pivnul_list, lpn_list, xsize, pp_first2swap_l, pp_lastpanelondisk, pp_lastpivrptrindexfilled, maxfromm, is_maxfromm_avail, pivot_option, iend_blr, inextpiv, ooc_effective_on_front, nvschur, parpiv_t1, lr_activated)
Definition sfac_front_aux.F:1102
smumps_fac_front_aux_m::smumps_fac_mq_ldlt
subroutine smumps_fac_mq_ldlt(iend_block, nfront, nass, npiv, inode, a, la, lda, poselt, ifinb, pivsiz, maxfromm, is_maxfromm_avail, is_max_useful, parpiv_t1, last_row, iend_blr, nvschur_k253, lr_activated)
Definition sfac_front_aux.F:1609
smumps_fac_front_aux_m::smumps_fac_sq_ldlt
subroutine smumps_fac_sq_ldlt(ibeg_block, iend_block, npiv, nfront, nass, inode, a, la, lda, poselt, keep, keep8, first_row_trsm, last_row_trsm, last_col_gemm, last_row_gemm, call_trsm, call_gemm, lr_activated, iw, liw, offset_iw)
Definition sfac_front_aux.F:1928
smumps_fac_front_aux_m::smumps_fac_ldlt_copyscale_u
subroutine smumps_fac_ldlt_copyscale_u(irowmax, irowmin, sizecopy, lda, ncols, liw, iw, offset_iw, la, a, poselt, a_lpos, a_upos, a_dpos)
Definition sfac_front_aux.F:2212
smumps_fac_front_aux_m::smumps_get_size_schur_in_front
subroutine smumps_get_size_schur_in_front(n, ncb, size_schur, row_indices, perm, nvschur)
Definition sfac_front_aux.F:2505
smumps_fac_lr
Definition sfac_lr.F:14
smumps_fac_lr::smumps_compress_panel
subroutine smumps_compress_panel(a, la, poselt, iflag, ierror, nfront, begs_blr, nb_blr, toleps, tol_opt, k473, blr_panel, current_blr, dir, work, tau, jpvt, lwork, rwork, block, maxi_cluster, nelim, lbandslave, npiv, ishift, niv, kpercent, keep8, k480, beg_i_in, end_i_in, frswap)
Definition sfac_lr.F:2199
smumps_fac_lr::smumps_blr_upd_cb_left_ldlt
subroutine smumps_blr_upd_cb_left_ldlt(a, la, poselt, nfront, begs_blr, begs_blr_dyn, nb_incb, nb_inasm, nass, iwhandler, iw2, block, acc_lua, maxi_cluster, maxi_rank, niv, iflag, ierror, midblk_compress, toleps, tol_opt, kpercent_rmb, k480, k479, k478, kpercent_lua, kpercent, keep8)
Definition sfac_lr.F:1130
smumps_fac_lr::smumps_decompress_panel
subroutine smumps_decompress_panel(a, la, poselt, lda11, lda21, copy_dense_blocks, begs_blr_diag, begs_blr_first_offdiag, nb_blr, blr_panel, current_blr, dir, decomp_timer, beg_i_in, end_i_in, only_nelim_in, cbasm_tofix_in)
Definition sfac_lr.F:1754
smumps_fac_lr::smumps_blr_panel_lrtrsm
subroutine smumps_blr_panel_lrtrsm(a, la, poselt, nfront, ibeg_block, nb_blr, blr_loru, current_blr, first_block, last_block, niv, sym, loru, lbandslave, iw, offset_iw, nass)
Definition sfac_lr.F:2437
smumps_fac_lr::smumps_blr_upd_panel_left_ldlt
subroutine smumps_blr_upd_panel_left_ldlt(a, la, poselt, nfront, iwhandler, begs_blr, current_blr, nb_blr, npartsass, nelim, iw2, block, acc_lua, maxi_cluster, maxi_rank, niv, iflag, ierror, midblk_compress, toleps, tol_opt, kpercent_rmb, k480, k479, k478, kpercent_lua, kpercent, keep8, first_block)
Definition sfac_lr.F:447
smumps_fac_lr::smumps_blr_update_trailing_ldlt
subroutine smumps_blr_update_trailing_ldlt(a, la, poselt, iflag, ierror, nfront, begs_blr, nb_blr, current_blr, blr_l, nelim, iw2, block, maxi_cluster, npiv, niv, midblk_compress, toleps, tol_opt, kpercent)
Definition sfac_lr.F:24
smumps_fac_lr::smumps_blr_upd_nelim_var_l
subroutine smumps_blr_upd_nelim_var_l(a_u, la_u, upos, a_l, la_l, lpos, iflag, ierror, ldu, ldl, begs_blr_l, current_blr, blr_l, nb_blr_l, first_block, nelim, utrans)
Definition sfac_lr.F:259
smumps_lr_data_m
Definition smumps_lr_data_m.F:14
smumps_lr_data_m::smumps_blr_save_begs_blr_dyn
subroutine, public smumps_blr_save_begs_blr_dyn(iwhandler, begs_blr_dynamic)
Definition smumps_lr_data_m.F:588
smumps_lr_data_m::smumps_blr_save_panel_loru
subroutine, public smumps_blr_save_panel_loru(iwhandler, loru, ipanel, lrb_panel)
Definition smumps_lr_data_m.F:512
smumps_lr_data_m::smumps_blr_retrieve_panel_loru
subroutine, public smumps_blr_retrieve_panel_loru(iwhandler, loru, ipanel, thelrbpanel)
Definition smumps_lr_data_m.F:750
smumps_lr_data_m::smumps_blr_save_init
subroutine, public smumps_blr_save_init(iwhandler, issym, ist2, isslave, nb_panels, begs_blr_l, begs_blr_col, nb_accesses_init, info)
Definition smumps_lr_data_m.F:235
smumps_lr_data_m::smumps_blr_save_cb_lrb
subroutine, public smumps_blr_save_cb_lrb(iwhandler, cb_lrb)
Definition smumps_lr_data_m.F:532
smumps_lr_data_m::smumps_blr_save_diag_block
subroutine, public smumps_blr_save_diag_block(iwhandler, ipanel, d)
Definition smumps_lr_data_m.F:547
smumps_lr_data_m::smumps_blr_retrieve_nfs4father
subroutine, public smumps_blr_retrieve_nfs4father(iwhandler, nfs4father)
Definition smumps_lr_data_m.F:868
smumps_lr_data_m::smumps_blr_retrieve_begsblr_sta
subroutine, public smumps_blr_retrieve_begsblr_sta(iwhandler, begs_blr_static)
Definition smumps_lr_data_m.F:622
smumps_lr_stats
Definition slr_stats.F:14
smumps_lr_stats::cnt_nodes
integer cnt_nodes
Definition slr_stats.F:23
smumps_lr_type
Definition slr_type.F:14
smumps_lr_type::dealloc_blr_panel
subroutine dealloc_blr_panel(blr_panel, iend, keep8, k34, ibeg_in)
Definition slr_type.F:56
smumps_lr_type::dealloc_lrb
subroutine dealloc_lrb(lrb_out, keep8, k34)
Definition slr_type.F:25
smumps_ooc
Definition smumps_ooc.F:14
smumps_ooc::smumps_ooc_io_lu_panel
subroutine, public smumps_ooc_io_lu_panel(strat, typefile, afac, lafac, monbloc, lnextpiv2bewritten, unextpiv2bewritten, iw, liwfac, myid, filesize, ierr, last_call)
Definition smumps_ooc.F:3032
smumps_set_parpivt1
subroutine smumps_set_parpivt1(inode, nfront, nass1, keep, lr_activated, parpiv_t1)
Definition sfac_asm.F:788
smumps_ooc_pp_tryrelease_space
subroutine smumps_ooc_pp_tryrelease_space(iwpos, ioldps, iw, liw, monbloc, nfront, keep)
Definition sooc_panel_piv.F:235
smumps_lr_type::lrb_type
Definition slr_type.F:16
mumps_geti8
subroutine mumps_geti8(i8, int_array)
Definition tools_common.F:888
mumps_dm_fac_upd_dyn_memcnts
subroutine mumps_dm_fac_upd_dyn_memcnts(mem_count_allocated, atomic_updates, keep8, iflag, ierror, k69upd, k71upd)
Definition tools_common.F:1476