00001
00063 #ifndef __SUPERLU_sSP_DEFS
00064 #define __SUPERLU_sSP_DEFS
00065
00066
00067
00068
00069
00070
00071
00072 #ifdef _CRAY
00073 #include <fortran.h>
00074 #include <string.h>
00075 #endif
00076
00077
00078 typedef int int_t;
00079
00080 #include <math.h>
00081 #include <limits.h>
00082 #include "slu_Cnames.h"
00083 #include "supermatrix.h"
00084 #include "slu_util.h"
00085
00086
00087
00088 typedef struct {
00089 int *xsup;
00090 int *supno;
00091 int *lsub;
00092 int *xlsub;
00093 float *lusup;
00094 int *xlusup;
00095 float *ucol;
00096 int *usub;
00097 int *xusub;
00098 int nzlmax;
00099 int nzumax;
00100 int nzlumax;
00101 int n;
00102 LU_space_t MemModel;
00103 int num_expansions;
00104 ExpHeader *expanders;
00105 LU_stack_t stack;
00106 } GlobalLU_t;
00107
00108
00109
00110
00111 #ifdef __cplusplus
00112 extern "C" {
00113 #endif
00114
00116 extern void
00117 sgssv(superlu_options_t *, SuperMatrix *, int *, int *, SuperMatrix *,
00118 SuperMatrix *, SuperMatrix *, SuperLUStat_t *, int *);
00119 extern void
00120 sgssvx(superlu_options_t *, SuperMatrix *, int *, int *, int *,
00121 char *, float *, float *, SuperMatrix *, SuperMatrix *,
00122 void *, int, SuperMatrix *, SuperMatrix *,
00123 float *, float *, float *, float *,
00124 mem_usage_t *, SuperLUStat_t *, int *);
00125
00126 extern void
00127 sgsisv(superlu_options_t *, SuperMatrix *, int *, int *, SuperMatrix *,
00128 SuperMatrix *, SuperMatrix *, SuperLUStat_t *, int *);
00129 extern void
00130 sgsisx(superlu_options_t *, SuperMatrix *, int *, int *, int *,
00131 char *, float *, float *, SuperMatrix *, SuperMatrix *,
00132 void *, int, SuperMatrix *, SuperMatrix *, float *, float *,
00133 mem_usage_t *, SuperLUStat_t *, int *);
00134
00135
00137 extern void
00138 sCreate_CompCol_Matrix(SuperMatrix *, int, int, int, float *,
00139 int *, int *, Stype_t, Dtype_t, Mtype_t);
00140 extern void
00141 sCreate_CompRow_Matrix(SuperMatrix *, int, int, int, float *,
00142 int *, int *, Stype_t, Dtype_t, Mtype_t);
00143 extern void
00144 sCopy_CompCol_Matrix(SuperMatrix *, SuperMatrix *);
00145 extern void
00146 sCreate_Dense_Matrix(SuperMatrix *, int, int, float *, int,
00147 Stype_t, Dtype_t, Mtype_t);
00148 extern void
00149 sCreate_SuperNode_Matrix(SuperMatrix *, int, int, int, float *,
00150 int *, int *, int *, int *, int *,
00151 Stype_t, Dtype_t, Mtype_t);
00152 extern void
00153 sCopy_Dense_Matrix(int, int, float *, int, float *, int);
00154
00155 extern void countnz (const int, int *, int *, int *, GlobalLU_t *);
00156 extern void ilu_countnz (const int, int *, int *, GlobalLU_t *);
00157 extern void fixupL (const int, const int *, GlobalLU_t *);
00158
00159 extern void sallocateA (int, int, float **, int **, int **);
00160 extern void sgstrf (superlu_options_t*, SuperMatrix*,
00161 int, int, int*, void *, int, int *, int *,
00162 SuperMatrix *, SuperMatrix *, SuperLUStat_t*, int *);
00163 extern int ssnode_dfs (const int, const int, const int *, const int *,
00164 const int *, int *, int *, GlobalLU_t *);
00165 extern int ssnode_bmod (const int, const int, const int, float *,
00166 float *, GlobalLU_t *, SuperLUStat_t*);
00167 extern void spanel_dfs (const int, const int, const int, SuperMatrix *,
00168 int *, int *, float *, int *, int *, int *,
00169 int *, int *, int *, int *, GlobalLU_t *);
00170 extern void spanel_bmod (const int, const int, const int, const int,
00171 float *, float *, int *, int *,
00172 GlobalLU_t *, SuperLUStat_t*);
00173 extern int scolumn_dfs (const int, const int, int *, int *, int *, int *,
00174 int *, int *, int *, int *, int *, GlobalLU_t *);
00175 extern int scolumn_bmod (const int, const int, float *,
00176 float *, int *, int *, int,
00177 GlobalLU_t *, SuperLUStat_t*);
00178 extern int scopy_to_ucol (int, int, int *, int *, int *,
00179 float *, GlobalLU_t *);
00180 extern int spivotL (const int, const double, int *, int *,
00181 int *, int *, int *, GlobalLU_t *, SuperLUStat_t*);
00182 extern void spruneL (const int, const int *, const int, const int,
00183 const int *, const int *, int *, GlobalLU_t *);
00184 extern void sreadmt (int *, int *, int *, float **, int **, int **);
00185 extern void sGenXtrue (int, int, float *, int);
00186 extern void sFillRHS (trans_t, int, float *, int, SuperMatrix *,
00187 SuperMatrix *);
00188 extern void sgstrs (trans_t, SuperMatrix *, SuperMatrix *, int *, int *,
00189 SuperMatrix *, SuperLUStat_t*, int *);
00190
00191 extern void sgsitrf (superlu_options_t*, SuperMatrix*, int, int, int*,
00192 void *, int, int *, int *, SuperMatrix *, SuperMatrix *,
00193 SuperLUStat_t*, int *);
00194 extern int sldperm(int, int, int, int [], int [], float [],
00195 int [], float [], float []);
00196 extern int ilu_ssnode_dfs (const int, const int, const int *, const int *,
00197 const int *, int *, GlobalLU_t *);
00198 extern void ilu_spanel_dfs (const int, const int, const int, SuperMatrix *,
00199 int *, int *, float *, float *, int *, int *,
00200 int *, int *, int *, int *, GlobalLU_t *);
00201 extern int ilu_scolumn_dfs (const int, const int, int *, int *, int *,
00202 int *, int *, int *, int *, int *,
00203 GlobalLU_t *);
00204 extern int ilu_scopy_to_ucol (int, int, int *, int *, int *,
00205 float *, int, milu_t, double, int,
00206 float *, int *, GlobalLU_t *, float *);
00207 extern int ilu_spivotL (const int, const double, int *, int *, int, int *,
00208 int *, int *, int *, double, milu_t,
00209 float, GlobalLU_t *, SuperLUStat_t*);
00210 extern int ilu_sdrop_row (superlu_options_t *, int, int, double,
00211 int, int *, double *, GlobalLU_t *,
00212 float *, float *, int);
00213
00214
00217 extern void sgsequ (SuperMatrix *, float *, float *, float *,
00218 float *, float *, int *);
00219 extern void slaqgs (SuperMatrix *, float *, float *, float,
00220 float, float, char *);
00221 extern void sgscon (char *, SuperMatrix *, SuperMatrix *,
00222 float, float *, SuperLUStat_t*, int *);
00223 extern float sPivotGrowth(int, SuperMatrix *, int *,
00224 SuperMatrix *, SuperMatrix *);
00225 extern void sgsrfs (trans_t, SuperMatrix *, SuperMatrix *,
00226 SuperMatrix *, int *, int *, char *, float *,
00227 float *, SuperMatrix *, SuperMatrix *,
00228 float *, float *, SuperLUStat_t*, int *);
00229
00230 extern int sp_strsv (char *, char *, char *, SuperMatrix *,
00231 SuperMatrix *, float *, SuperLUStat_t*, int *);
00232 extern int sp_sgemv (char *, float, SuperMatrix *, float *,
00233 int, float, float *, int);
00234
00235 extern int sp_sgemm (char *, char *, int, int, int, float,
00236 SuperMatrix *, float *, int, float,
00237 float *, int);
00238 extern float slamch_(char *);
00239
00240
00242 extern int sLUMemInit (fact_t, void *, int, int, int, int, int,
00243 float, SuperMatrix *, SuperMatrix *,
00244 GlobalLU_t *, int **, float **);
00245 extern void sSetRWork (int, int, float *, float **, float **);
00246 extern void sLUWorkFree (int *, float *, GlobalLU_t *);
00247 extern int sLUMemXpand (int, int, MemType, int *, GlobalLU_t *);
00248
00249 extern float *floatMalloc(int);
00250 extern float *floatCalloc(int);
00251 extern int smemory_usage(const int, const int, const int, const int);
00252 extern int sQuerySpace (SuperMatrix *, SuperMatrix *, mem_usage_t *);
00253 extern int ilu_sQuerySpace (SuperMatrix *, SuperMatrix *, mem_usage_t *);
00254
00256 extern void sreadhb(int *, int *, int *, float **, int **, int **);
00257 extern void sreadrb(int *, int *, int *, float **, int **, int **);
00258 extern void sreadtriple(int *, int *, int *, float **, int **, int **);
00259 extern void sCompRow_to_CompCol(int, int, int, float*, int*, int*,
00260 float **, int **, int **);
00261 extern void sfill (float *, int, float);
00262 extern void sinf_norm_error (int, SuperMatrix *, float *);
00263 extern void PrintPerf (SuperMatrix *, SuperMatrix *, mem_usage_t *,
00264 float, float, float *, float *, char *);
00265 extern float sqselect(int, float *, int);
00266
00267
00269 extern void sPrint_CompCol_Matrix(char *, SuperMatrix *);
00270 extern void sPrint_SuperNode_Matrix(char *, SuperMatrix *);
00271 extern void sPrint_Dense_Matrix(char *, SuperMatrix *);
00272 extern void sprint_lu_col(char *, int, int, int *, GlobalLU_t *);
00273 extern int print_double_vec(char *, int, double *);
00274 extern void check_tempv(int, float *);
00275
00276 #ifdef __cplusplus
00277 }
00278 #endif
00279
00280 #endif
00281