00001
00002
00003 #ifndef ClpCholeskyDense_H
00004 #define ClpCholeskyDense_H
00005
00006 #include "ClpCholeskyBase.hpp"
00007 class ClpMatrixBase;
00008
00012 class ClpCholeskyDense : public ClpCholeskyBase {
00013
00014 public:
00019 virtual int order(ClpInterior * model) ;
00024 virtual int symbolic();
00027 virtual int factorize(const double * diagonal, int * rowsDropped) ;
00029 virtual void solve (double * region) ;
00031
00037 int reserveSpace(const ClpCholeskyBase * factor, int numberRows) ;
00039 CoinBigIndex space( int numberRows) const;
00041 void factorizePart2(int * rowsDropped) ;
00043 void factorizePart3(int * rowsDropped) ;
00045 void factor(longDouble * a, int n, int numberBlocks,
00046 longDouble * diagonal, longDouble * work, int * rowsDropped);
00048 void triRec(longDouble * aTri, int nThis, longDouble * aUnder, longDouble * diagonal, longDouble * work,
00049 int nLeft, int iBlock, int jBlock,
00050 int numberBlocks);
00052 void recTri(longDouble * aUnder, int nTri, int nDo,
00053 int iBlock, int jBlock,longDouble * aTri,
00054 longDouble * diagonal, longDouble * work,
00055 int numberBlocks);
00060 void recRec(longDouble * above, int nUnder, int nUnderK,
00061 int nDo, longDouble * aUnder, longDouble *aOther,longDouble * diagonal, longDouble * work,
00062 int kBlock,int iBlock, int jBlock,
00063 int numberBlocks);
00065 void factorLeaf(longDouble * a, int n,
00066 longDouble * diagonal, longDouble * work, int * rowsDropped);
00068 void triRecLeaf(longDouble * aTri, longDouble * aUnder,
00069 longDouble * diagonal, longDouble * work,
00070 int nUnder);
00072 void recTriLeaf(longDouble * aUnder, longDouble * aTri,
00073 longDouble * diagonal, longDouble * work, int nUnder);
00078 void recRecLeaf(longDouble * above,
00079 longDouble * aUnder, longDouble *aOther, longDouble * diagonal, longDouble * work,
00080 int nUnder);
00082 void solveF1(longDouble * a,int n,double * region);
00083 void solveF2(longDouble * a,int n,double * region,double * region2);
00085 void solveB1(longDouble * a,int n,double * region);
00086 void solveB2(longDouble * a,int n,double * region,double * region2);
00088 void solveLong (longDouble * region) ;
00090 void solveF1Long(longDouble * a,int n,longDouble * region);
00091 void solveF2Long(longDouble * a,int n,longDouble * region,longDouble * region2);
00093 void solveB1Long(longDouble * a,int n,longDouble * region);
00094 void solveB2Long(longDouble * a,int n,longDouble * region,longDouble * region2);
00096 void solveLongWork (longWork * region) ;
00098 void solveF1LongWork(longDouble * a,int n,longWork * region);
00099 void solveF2LongWork(longDouble * a,int n,longWork * region,longWork * region2);
00101 void solveB1LongWork(longDouble * a,int n,longWork * region);
00102 void solveB2LongWork(longDouble * a,int n,longWork * region,longWork * region2);
00103 int bNumber(const longDouble * array,int &, int&);
00105 inline longDouble * aMatrix() const
00106 { return sparseFactor_;}
00108 inline longDouble * diagonal() const
00109 { return diagonal_;}
00111
00112
00116 ClpCholeskyDense();
00118 virtual ~ClpCholeskyDense();
00119
00120 ClpCholeskyDense(const ClpCholeskyDense&);
00121
00122 ClpCholeskyDense& operator=(const ClpCholeskyDense&);
00124 virtual ClpCholeskyBase * clone() const ;
00126
00127
00128 private:
00131
00132 bool borrowSpace_;
00134 };
00135
00136 #endif