00001
00002
00003
00004
00005
00006
00007
00008
00009
00010
00011
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039 #ifndef UmfpackGenLinSOE_h
00040 #define UmfpackGenLinSOE_h
00041
00042 #include <LinearSOE.h>
00043 #include <Vector.h>
00044
00045 class UmfpackGenLinSolver;
00046
00047 class UmfpackGenLinSOE : public LinearSOE
00048 {
00049 public:
00050 UmfpackGenLinSOE(UmfpackGenLinSolver &theSolver);
00051 UmfpackGenLinSOE(int N, int NNZ, int *rowStartA, int *colA,
00052 UmfpackGenLinSolver &theSolver);
00053
00054 ~UmfpackGenLinSOE();
00055
00056 int getNumEqn(void) const;
00057 int setSize(Graph &theGraph);
00058 int addA(const Matrix &, const ID &, double fact = 1.0);
00059 int addB(const Vector &, const ID &, double fact = 1.0);
00060 int setB(const Vector &, double fact = 1.0);
00061
00062 void zeroA(void);
00063 void zeroB(void);
00064
00065 const Vector &getX(void);
00066 const Vector &getB(void);
00067 double normRHS(void);
00068
00069 void setX(int loc, double value);
00070 void setX(const Vector &x);
00071 int setUmfpackGenLinSolver(UmfpackGenLinSolver &newSolver);
00072
00073 int sendSelf(int commitTag, Channel &theChannel);
00074 int recvSelf(int commitTag, Channel &theChannel,
00075 FEM_ObjectBroker &theBroker);
00076
00077 friend class UmfpackGenLinSolver;
00078
00079 protected:
00080
00081 private:
00082 int size;
00083 int nnz;
00084 double *A, *B, *X;
00085 int *colA, *rowStartA;
00086 int lValue;
00087 int *index;
00088
00089 Vector *vectX;
00090 Vector *vectB;
00091 int Asize, Bsize;
00092 bool factored;
00093
00094 };
00095
00096
00097 #endif
00098