2! This file is part of
MUMPS 5.5.1, released
3! on Tue
Jul 12 13:17:24 UTC 2022
9! This version of
MUMPS is provided to you free of charge. It is
10! released under
the CeCILL-C license
11! (see doc/CeCILL-C_V1-en.txt, doc/CeCILL-C_V1-fr.txt, and
14! This file includes various internal datastructures
15! passed through
the main MUMPS structure between successive
16! phases of
the solver. The
main one is
root information
for
17!
the multifrontal tree.
18 TYPE SMUMPS_ROOT_STRUC
20 INTEGER :: MBLOCK, NBLOCK, NPROW, NPCOL
21 INTEGER :: MYROW, MYCOL
22 INTEGER :: SCHUR_MLOC, SCHUR_NLOC, SCHUR_LLD
24 INTEGER :: ROOT_SIZE, TOT_ROOT_SIZE
25! descriptor
for scalapack
26 INTEGER, DIMENSION( 9 ) :: DESCRIPTOR
27 INTEGER :: CNTXT_BLACS, LPIV, rootpad0
28 INTEGER, DIMENSION(:), POINTER :: RG2L_ROW
29 INTEGER, DIMENSION(:), POINTER :: RG2L_COL
30 INTEGER , DIMENSION(:), POINTER :: IPIV, rootpad1
31! Centralized master of
root
32 REAL, DIMENSION(:), POINTER :: RHS_CNTR_MASTER_ROOT
33! Used to access Schur easily from
root structure
34 REAL, DIMENSION(:), POINTER :: SCHUR_POINTER
35!
for try_null_space preprocessing constant only:
36 REAL, DIMENSION(:), POINTER :: QR_TAU, rootpad2
38! case of scalapack
root: to store
RHS in 2D block cyclic
39! format compatible with
root distribution
40 REAL, DIMENSION(:,:), POINTER :: RHS_ROOT, rootpad
41!
for try_nullspace preprocessing constant only:
42 REAL :: QR_RCOND, rootpad3
43 LOGICAL :: yes, gridinit_done
44!
for SVD on
root (#define try_null_space)
45 REAL, DIMENSION(:,:), POINTER :: SVD_U, SVD_VT
46!
for RR on
root (#define try_null_space)
47 REAL, DIMENSION(:), POINTER :: SINGULAR_VALUES
48 INTEGER :: NB_SINGULAR_VALUES,rootpad4
50 END
TYPE SMUMPS_ROOT_STRUC
52 TYPE SMUMPS_L0OMPFAC_T
54 REAL, POINTER, DIMENSION(:) :: A
56 END
TYPE SMUMPS_L0OMPFAC_T
!This file is part of MUMPS
!This file is part of released !on Tue ENS INP !Mumps Technologies
!This file is part of released !on Tue ENS INP Toulouse
!This file is part of released !on Tue ENS Lyon
!This file is part of released !on Tue Jul
!This file is part of released !on Tue ENS INP Inria
!This file is part of released !on Tue CNRS
end diagonal values have been computed in the(sparse) matrix id.SOL
we set the rigth hand side id RHS
for(i8=*sizetab-1;i8 >=0;i8--)
int main(int argc, char *argv[])