1 /* $Id: ClpCholeskyPardiso.hpp 1665 2011-01-04 17:55:54Z lou $ */ 2 // Copyright (C) 2003, International Business Machines 3 // Corporation and others. All Rights Reserved. 4 // This code is licensed under the terms of the Eclipse Public License (EPL). 5 6 #ifndef ClpCholeskyPardiso_H 7 #define ClpCholeskyPardiso_H 8 #ifdef PARDISO_BARRIER 9 10 #include "ClpCholeskyBase.hpp" 11 #include "mkl_pardiso.h" 12 #include "mkl_types.h" 13 14 class ClpMatrixBase; 15 class ClpCholeskyDense; 16 17 /** Pardiso class for Clp Cholesky factorization 18 19 */ 20 class ClpCholeskyPardiso : public ClpCholeskyBase { 21 22 public: 23 /**@name Virtual methods that the derived classes provides */ 24 //@{ 25 /** Orders rows and saves pointer to matrix.and model. 26 Returns non-zero if not enough memory */ 27 virtual int order(ClpInterior *model); 28 /** Does Symbolic factorization given permutation. 29 This is called immediately after order. If user provides this then 30 user must provide factorize and solve. Otherwise the default factorization is used 31 returns non-zero if not enough memory */ 32 virtual int symbolic(); 33 /** Factorize - filling in rowsDropped and returning number dropped. 34 If return code negative then out of memory */ 35 virtual int factorize(const double *diagonal, int *rowsDropped); 36 /** Uses factorization to solve. */ 37 virtual void solve(double *region); 38 //@} 39 40 /**@name Constructors, destructor */ 41 //@{ 42 /** Constructor which has dense columns activated. 43 Default is off. */ 44 ClpCholeskyPardiso(int denseThreshold = -1); 45 /** Destructor */ 46 virtual ~ClpCholeskyPardiso(); 47 // Copy 48 ClpCholeskyPardiso(const ClpCholeskyPardiso &); 49 // Assignment 50 ClpCholeskyPardiso &operator=(const ClpCholeskyPardiso &); 51 /// Clone 52 virtual ClpCholeskyBase *clone() const; 53 //@} 54 55 private: 56 /**@name Data members */ 57 //@{ 58 int lastNumberDropped_; 59 //MKL_INT iparm_[64]; 60 //@} 61 }; 62 63 #endif 64 #endif 65 66 /* vi: softtabstop=2 shiftwidth=2 expandtab tabstop=2 67 */ 68