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
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
|
/*
*
* This file is part of MUMPS 5.8.1, released
* on Wed Jul 30 16:49:18 UTC 2025
*
*
* Copyright 1991-2025 CERFACS, CNRS, ENS Lyon, INP Toulouse, Inria,
* Mumps Technologies, University of Bordeaux.
*
* This version of MUMPS is provided to you free of charge. It is
* released under the CeCILL-C license
* (see doc/CeCILL-C_V1-en.txt, doc/CeCILL-C_V1-fr.txt, and
* https://cecill.info/licences/Licence_CeCILL-C_V1-en.html)
*
*/
/* Mostly written in march 2002 (JYL) */
#ifndef SMUMPS_C_H
#define SMUMPS_C_H
#ifdef __cplusplus
extern "C" {
#endif
#include "mumps_compat.h"
/* Next line defines MUMPS_INT, SMUMPS_COMPLEX and SMUMPS_REAL */
#include "mumps_c_types.h"
#ifndef MUMPS_VERSION
/* Protected in case headers of other arithmetics are included */
#define MUMPS_VERSION "5.8.1"
#endif
#ifndef MUMPS_VERSION_MAX_LEN
#define MUMPS_VERSION_MAX_LEN 30
#endif
/*
* Definition of the (simplified) MUMPS C structure.
* NB: SMUMPS_COMPLEX are REAL types in s and d arithmetics.
*/
typedef struct {
MUMPS_INT sym, par, job;
MUMPS_INT comm_fortran; /* Fortran communicator */
MUMPS_INT icntl[60];
MUMPS_INT keep[500];
SMUMPS_REAL cntl[15];
SMUMPS_REAL dkeep[230];
MUMPS_INT8 keep8[150];
MUMPS_INT n;
MUMPS_INT nblk;
MUMPS_INT nz_alloc; /* used in matlab interface to decide if we
free + malloc when we have large variation */
/* Assembled entry */
MUMPS_INT nz;
MUMPS_INT8 nnz;
MUMPS_INT *irn;
MUMPS_INT *jcn;
SMUMPS_COMPLEX *a;
/* Distributed entry */
MUMPS_INT nz_loc;
MUMPS_INT8 nnz_loc;
MUMPS_INT *irn_loc;
MUMPS_INT *jcn_loc;
SMUMPS_COMPLEX *a_loc;
/* Element entry */
MUMPS_INT nelt;
MUMPS_INT *eltptr;
MUMPS_INT *eltvar;
SMUMPS_COMPLEX *a_elt;
/* Matrix by blocks */
MUMPS_INT *blkptr;
MUMPS_INT *blkvar;
/* Ordering, if given by user */
MUMPS_INT *perm_in;
/* Orderings returned to user */
MUMPS_INT *sym_perm; /* symmetric permutation */
MUMPS_INT *uns_perm; /* column permutation */
/* Scaling (inout but complicated) */
SMUMPS_REAL *colsca;
SMUMPS_REAL *rowsca;
MUMPS_INT colsca_from_mumps;
MUMPS_INT rowsca_from_mumps;
/* Distributed scaling(out) */
SMUMPS_REAL *colsca_loc;
SMUMPS_REAL *rowsca_loc;
/* Info after facto */
MUMPS_INT *rowind;
MUMPS_INT *colind;
SMUMPS_COMPLEX *pivots;
/* RHS, solution, ouptput data and statistics */
SMUMPS_COMPLEX *rhs, *redrhs, *rhs_sparse, *sol_loc, *rhs_loc, *rhsintr;
MUMPS_INT *irhs_sparse, *irhs_ptr, *isol_loc, *irhs_loc, *glob2loc_rhs, *glob2loc_sol;
MUMPS_INT nrhs, lrhs, lredrhs, nz_rhs, lsol_loc, nloc_rhs, lrhs_loc, nsol_loc;
MUMPS_INT schur_mloc, schur_nloc, schur_lld;
MUMPS_INT mblock, nblock, nprow, npcol;
MUMPS_INT ld_rhsintr;
MUMPS_INT info[80],infog[80];
SMUMPS_REAL rinfo[40], rinfog[40];
/* Null space */
MUMPS_INT deficiency;
MUMPS_INT *pivnul_list;
MUMPS_INT *mapping;
SMUMPS_REAL *singular_values;
/* Schur */
MUMPS_INT size_schur;
MUMPS_INT *listvar_schur;
SMUMPS_COMPLEX *schur;
/* user workspace */
SMUMPS_COMPLEX *wk_user;
/* Version number: length=30 in FORTRAN + 1 for final \0 + 1 for alignment */
char version_number[MUMPS_VERSION_MAX_LEN + 1 + 1];
/* For out-of-core */
char ooc_tmpdir[1024];
char ooc_prefix[256];
/* To save the matrix in matrix market format */
char write_problem[1024];
MUMPS_INT lwk_user;
/* For save/restore feature */
char save_dir[1024];
char save_prefix[256];
/* Metis options */
MUMPS_INT metis_options[40];
/* Internal parameters */
MUMPS_INT instance_number;
} SMUMPS_STRUC_C;
void MUMPS_CALL
smumps_c( SMUMPS_STRUC_C * smumps_par );
#ifdef __cplusplus
}
#endif
#endif /* SMUMPS_C_H */
|
