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
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
|
#ifndef BOOST_NUMERIC_BINDINGS_MUMPS_464_DMUMPS_C_HPP
#define BOOST_NUMERIC_BINDINGS_MUMPS_464_DMUMPS_C_HPP
/*
THIS FILE IS PART OF MUMPS VERSION 4.6.4
This Version was built on Thu Jan 11 13:32:35 2007
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, ENSEEIHT-IRIT, and
INRIA.
Main contributors are Patrick Amestoy, Iain Duff, Abdou Guermouche,
Jacko Koster, Jean-Yves L'Excellent, and Stephane Pralet.
Up-to-date copies of the MUMPS package can be obtained
from the Web pages http://www.enseeiht.fr/apo/MUMPS/
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], [2], and [3] 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 and J.-Y. L'Excellent,
Multifrontal parallel distributed symmetric and unsymmetric solvers,
in Comput. Methods in Appl. Mech. Eng., 184, 501-520 (2000).
[2] 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).
[3] 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).
*/
/* $Id: dmumps_c.hpp 39206 2007-09-12 06:58:12Z karlmeerbergen $ */
/* Mostly written in march 2002 (JYL) */
// This file is modified by Karl Meerbergen for C++ users
/* Complex datatypes */
/* Next line defines DMUMPS_INT, DMUMPS_DOUBLE and DMUMPS_DOUBLE2 */
#include "dmumps_prec.h"
/*
* Definition of the (simplified)
* MUMPS C structure
*/
typedef struct
{
DMUMPS_INT sym, par, job;
DMUMPS_INT comm_fortran; /* Fortran communicator */
DMUMPS_INT icntl[40];
DMUMPS_DOUBLE2 cntl[5];
DMUMPS_INT n;
DMUMPS_INT nz_alloc; /* used in matlab interface to decide if
we free + malloc when we have large variation */
/* Assembled entry */
DMUMPS_INT nz; DMUMPS_INT *irn; DMUMPS_INT *jcn; DMUMPS_DOUBLE *a;
/* Distributed entry */
DMUMPS_INT nz_loc; DMUMPS_INT *irn_loc; DMUMPS_INT *jcn_loc; DMUMPS_DOUBLE *a_loc;
/* Element entry */
DMUMPS_INT nelt; DMUMPS_INT *eltptr; DMUMPS_INT *eltvar; DMUMPS_DOUBLE *a_elt;
/* Ordering, if given by user */
DMUMPS_INT *perm_in;
/* Orderings returned to user */
/* symmetric permutation */
DMUMPS_INT *sym_perm;
/* column permutation */
DMUMPS_INT *uns_perm;
/* Scaling (input only in this version) */
DMUMPS_DOUBLE *colsca; DMUMPS_DOUBLE *rowsca;
/* RHS, solution, ouptput data and statistics */
DMUMPS_DOUBLE *rhs, *rhs_sparse, *sol_loc;
DMUMPS_INT *irhs_sparse, *irhs_ptr, *isol_loc;
DMUMPS_INT nrhs, lrhs, nz_rhs, lsol_loc;
DMUMPS_INT schur_mloc, schur_nloc, schur_lld;
DMUMPS_INT mblock, nblock, nprow, npcol;
DMUMPS_INT info[40],infog[40];
DMUMPS_DOUBLE2 rinfo[20], rinfog[20];
/* Null space */
DMUMPS_INT deficiency; DMUMPS_DOUBLE * nullspace; DMUMPS_INT * mapping;
/* Schur */
DMUMPS_INT size_schur; DMUMPS_INT *listvar_schur; DMUMPS_DOUBLE *schur;
/* Internal parameters */
DMUMPS_INT instance_number;
/* For out-of-core */
char ooc_tmpdir[151];
char ooc_prefix[151];
} DMUMPS_STRUC_C;
#define MUMPS_CALL
#if defined(_WIN32)
/*
* Next line May be needed depending on your Windows environment:
* #define MUMPS_CALL __stdcall
*/
#endif
extern "C" {
void dmumps_c(DMUMPS_STRUC_C * dmumps_par);
}
#endif
|
