6 #ifndef ClpCholeskyBase_H
7 #define ClpCholeskyBase_H
12 #ifndef CLP_LONG_CHOLESKY
13 #define CLP_LONG_CHOLESKY 0
20 #if COIN_LONG_WORK == 0
21 #if CLP_LONG_CHOLESKY > 0
22 #define CHOLESKY_BAD_COMBINATION
25 #if CLP_LONG_CHOLESKY == 0
26 #define CHOLESKY_BAD_COMBINATION
29 #ifdef CHOLESKY_BAD_COMBINATION
30 #warning("Bad combination of CLP_LONG_CHOLESKY and COIN_BIG_DOUBLE/COIN_LONG_WORK");
31 "Bad combination of CLP_LONG_CHOLESKY and COIN_LONG_WORK"
33 #if CLP_LONG_CHOLESKY > 1
35 #define CHOL_SMALL_VALUE 1.0e-15
36 #elif CLP_LONG_CHOLESKY == 1
38 #define CHOL_SMALL_VALUE 1.0e-11
41 #define CHOL_SMALL_VALUE 1.0e-11
230 int symbolic1(
const int *Astart,
const int *Arow);
234 void symbolic2(
const int *Astart,
const int *Arow);
244 int preOrder(
bool lowerTriangular,
bool includeDiagonal,
bool doKKT);
longDouble * sparseFactor() const
Return sparseFactor.
void symbolic2(const int *Astart, const int *Arow)
Symbolic2 - Fills in indices Uses lower triangular so can do cliques etc.
int type_
type (may be useful) if > 20 do KKT
int size() const
Return size.
int sizeFactor_
sizeFactor.
virtual int symbolic()
Does Symbolic factorization given permutation.
int numberRowsDropped() const
numberRowsDropped. Number of rows gone
Abstract base class for Clp Matrices.
int * permuteInverse_
permute inverse.
This solves LPs using interior point methods.
char * whichDense_
Dense indicators.
char * rowsDropped() const
rowsDropped - which rows are gone
double goDense_
Go dense at this fraction.
double getDoubleParameter(int i)
get double parameter
longDouble * sparseFactor_
sparseFactor.
void setIntegerParameter(int i, int value)
Set integer parameter.
char * rowsDropped_
rowsDropped
int * indexStart_
Index starts.
void setType(int type)
Sets type.
int preOrder(bool lowerTriangular, bool includeDiagonal, bool doKKT)
Forms ADAT - returns nonzero if not enough memory.
double choleskyCondition_
choleskyCondition.
ClpMatrixBase * rowCopy_
Row copy of matrix.
double goDense() const
goDense i.e. use dense factoriaztion if > this (default 0.7).
int * choleskyStart_
choleskyStart - element starts
longDouble * workDouble_
double work array
virtual void solveKKT(CoinWorkDouble *region1, CoinWorkDouble *region2, const CoinWorkDouble *diagonal, CoinWorkDouble diagonalScaleFactor)
Uses factorization to solve.
int rank() const
rank. Returns rank
int numberRows_
numberRows. Number of Rows in factorization
int sizeIndex_
Size of index array.
int orderAMD()
AMD ordering.
virtual void solve(CoinWorkDouble *region)
Uses factorization to solve.
virtual ~ClpCholeskyBase()
Destructor (has to be public)
Base class for Clp Cholesky factorization Will do better factorization.
longDouble * workDouble() const
Return workDouble.
bool kkt() const
If KKT on.
int numberRowsDropped_
numberRowsDropped. Number of rows gone
void factorizePart2(int *rowsDropped)
Factorize - filling in rowsDropped and returning number dropped in integerParam.
void resetRowsDropped()
reset numberRowsDropped and rowsDropped.
int getIntegerParameter(int i)
get integer parameter
int symbolic1(const int *Astart, const int *Arow)
Symbolic1 - works out size without clever stuff.
ClpInterior * model_
model.
int status() const
status. Returns status
int numberRows() const
Return number of rows.
longDouble * diagonal() const
Return diagonal.
longDouble * diagonal_
Diagonal.
void setModel(ClpInterior *model)
model.
double doubleParameters_[64]
doubleParameters;
int firstDense_
First dense row.
virtual int order(ClpInterior *model)
Orders rows and saves pointer to matrix.and model.
ClpCholeskyDense * dense_
Dense cholesky.
ClpCholeskyBase & operator=(const ClpCholeskyBase &)
Assignment.
void updateDense(longDouble *d, int *first)
Updates dense part (broken out for profiling)
int numberTrials_
numberTrials. Number of trials before rejection
void setKKT(bool yesNo)
Set KKT.
virtual int factorize(const CoinWorkDouble *diagonal, int *rowsDropped)
Factorize - filling in rowsDropped and returning number dropped.
double choleskyCondition() const
choleskyCondition.
int * choleskyRow_
choleskyRow (can be shorter than sparsefactor)
int * permute_
main permute.
void setDoubleParameter(int i, double value)
Set double parameter.
longDouble * denseColumn_
Dense columns (updated)
virtual ClpCholeskyBase * clone() const
int type() const
Returns type.
ClpCholeskyBase(int denseThreshold=-1)
Constructor which has dense columns activated.
int status_
status. Status of factorization
bool doKKT_
Doing full KKT (only used if default symbolic and factorization)
void setGoDense(double value)
goDense i.e. use dense factoriaztion if > this (default 0.7).
int denseThreshold_
Dense threshold (for taking out of Cholesky)
int integerParameters_[64]
integerParameters