Coin-All: /home/coin/SVN-release/CoinAll-1.1.0/Clp/src/ClpCholeskyDense.hpp Source File

00001 // Copyright (C) 2003, International Business Machines
00002 // Corporation and others.  All Rights Reserved.
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   // Copy
00120   ClpCholeskyDense(const ClpCholeskyDense&);
00121   // Assignment
00122   ClpCholeskyDense& operator=(const ClpCholeskyDense&);
00124   virtual ClpCholeskyBase * clone() const ;
00126    
00127     
00128 private:
00131 
00132   bool borrowSpace_;
00134 };
00135 
00136 #endif