00001
00002
00003
00004
00005
00006 #if defined(_WIN32) || defined(_WIN64)
00007 #define pardiso_ PARDISO
00008 #else
00009 #define PARDISO pardiso_
00010 #endif
00011
00012 extern "C" { int PARDISO (void *, int *, int *, int *, int *, int *,
00013 double *, int *, int *, int *, int *, int *,
00014 int *, double *, double *, int *);
00015 int pardisoinit_(void *, int *, int *);
00016
00017 }
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031 class spasol
00032 {
00033 public:
00034
00035 spasol ();
00036 ~spasol ();
00037
00038 void symbfact (double *a,long *ci,long *adr,long ndof);
00039 void numfact (double *a,long *ci,long *adr,long ndof);
00040 void backsubst (double *a,long *ci,long *adr,long ndof,double *x,double *y);
00041
00042
00043
00044
00045
00046
00047 int nRows;
00048
00049 int nCols;
00050
00051 int nNonZeros;
00052
00053 int nrhs;
00054
00055 int num_procs;
00056
00057
00058
00059
00060
00061
00062
00063
00064 int mtype;
00065
00066
00067
00068
00069 void *pt[64];
00070
00071
00072 int iparm[64];
00073 int maxfct, mnum, phase, error, msglvl;
00074
00075 double ddum;
00076 int idum;
00077
00078
00079 };
00080