1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
|
!
! This file is part of MUMPS 4.9.2, built on Thu Nov 5 07:05:08 UTC 2009
!
!
! This version of MUMPS is provided to you free of charge. It is public
! domain, based on public domain software developed during the Esprit IV
! European project PARASOL (1996-1999) by CERFACS, ENSEEIHT-IRIT and RAL.
! Since this first public domain version in 1999, the developments are
! supported by the following institutions: CERFACS, CNRS, INPT(ENSEEIHT)-
! IRIT, and INRIA.
!
! Current development team includes Patrick Amestoy, Alfredo Buttari,
! Abdou Guermouche, Jean-Yves L'Excellent, Bora Ucar.
!
! Up-to-date copies of the MUMPS package can be obtained
! from the Web pages:
! http://mumps.enseeiht.fr/ or http://graal.ens-lyon.fr/MUMPS
!
!
! THIS MATERIAL IS PROVIDED AS IS, WITH ABSOLUTELY NO WARRANTY
! EXPRESSED OR IMPLIED. ANY USE IS AT YOUR OWN RISK.
!
!
! User documentation of any code that uses this software can
! include this complete notice. You can acknowledge (using
! references [1] and [2]) the contribution of this package
! in any scientific publication dependent upon the use of the
! package. You shall use reasonable endeavours to notify
! the authors of the package of this publication.
!
! [1] P. R. Amestoy, I. S. Duff, J. Koster and J.-Y. L'Excellent,
! A fully asynchronous multifrontal solver using distributed dynamic
! scheduling, SIAM Journal of Matrix Analysis and Applications,
! Vol 23, No 1, pp 15-41 (2001).
!
! [2] P. R. Amestoy and A. Guermouche and J.-Y. L'Excellent and
! S. Pralet, Hybrid scheduling for the parallel solution of linear
! systems. Parallel Computing Vol 32 (2), pp 136-156 (2006).
!
TYPE DMUMPS_ROOT_STRUC
SEQUENCE
INTEGER MBLOCK, NBLOCK, NPROW, NPCOL
INTEGER MYROW, MYCOL
INTEGER ROOT_SIZE, TOT_ROOT_SIZE
INTEGER :: CNTXT_BLACS
INTEGER, DIMENSION(:), POINTER :: RG2L_ROW
INTEGER, DIMENSION(:), POINTER :: RG2L_COL
INTEGER , DIMENSION(:), POINTER :: IPIV
INTEGER, DIMENSION( 9 ) :: DESCRIPTOR, DESCB
LOGICAL yes, gridinit_done
INTEGER LPIV
! Used to access Schur easily from root structure
DOUBLE PRECISION, DIMENSION(:), POINTER :: SCHUR_POINTER
INTEGER SCHUR_MLOC, SCHUR_NLOC, SCHUR_LLD
!
! Data for nullspace/QR
!
DOUBLE PRECISION, DIMENSION(:), POINTER :: QR_TAU
DOUBLE PRECISION QR_RCOND
!
! Givens rotations
!
INTEGER MAXG, GIND
DOUBLE PRECISION, DIMENSION(:),POINTER::GROW, GCOS, GSIN
!
! RRRLU data
!
INTEGER ELG_MAX,NULL_MAX
INTEGER ELIND,EUIND,NLUPDATE,NUUPDATE
INTEGER,DIMENSION(:),POINTER::PERM_ROW,PERM_COL
INTEGER,DIMENSION(:),POINTER::ELROW, EUROW, PTREL, PTREU
DOUBLE PRECISION, DIMENSION(:), POINTER :: ELELG, EUELG, DL
!
END TYPE DMUMPS_ROOT_STRUC
|
