ClpCholeskyDense.hpp

Go to the documentation of this file.

/* $Id: ClpCholeskyDense.hpp 2385 2019-01-06 19:43:06Z unxusr $ */
/*
  Copyright (C) 2003, International Business Machines Corporation
  and others.  All Rights Reserved.

  This code is licensed under the terms of the Eclipse Public License (EPL).
*/
#ifndef ClpCholeskyDense_H
#define ClpCholeskyDense_H

#include "ClpCholeskyBase.hpp"
class ClpMatrixBase;

class ClpCholeskyDense : public ClpCholeskyBase {

public:
  virtual int order(ClpInterior *model);
  virtual int symbolic();
  virtual int factorize(const CoinWorkDouble *diagonal, int *rowsDropped);
  virtual void solve(CoinWorkDouble *region);
  int reserveSpace(const ClpCholeskyBase *factor, int numberRows);
  int space(int numberRows) const;
  void factorizePart2(int *rowsDropped);
  void factorizePart3(int *rowsDropped);
  void solveF1(longDouble *a, int n, CoinWorkDouble *region);
  void solveF2(longDouble *a, int n, CoinWorkDouble *region, CoinWorkDouble *region2);
  void solveB1(longDouble *a, int n, CoinWorkDouble *region);
  void solveB2(longDouble *a, int n, CoinWorkDouble *region, CoinWorkDouble *region2);
  int bNumber(const longDouble *array, int &, int &);
  inline longDouble *aMatrix() const
  {
    return sparseFactor_;
  }
  inline longDouble *diagonal() const
  {
    return diagonal_;
  }
  ClpCholeskyDense();
  virtual ~ClpCholeskyDense();
  ClpCholeskyDense(const ClpCholeskyDense &);
  ClpCholeskyDense &operator=(const ClpCholeskyDense &);
  virtual ClpCholeskyBase *clone() const;
private:
  bool borrowSpace_;
};

/* structure for C */
typedef struct {
  longDouble *diagonal_;
  longDouble *a;
  longDouble *work;
  int *rowsDropped;
  double doubleParameters_[1]; /* corresponds to 10 */
  int integerParameters_[2]; /* corresponds to 34, nThreads */
  int n;
  int numberBlocks;
} ClpCholeskyDenseC;

extern "C" {
void ClpCholeskySpawn(void *);
}
void ClpCholeskyCfactor(ClpCholeskyDenseC *thisStruct,
  longDouble *a, int n, int numberBlocks,
  longDouble *diagonal, longDouble *work, int *rowsDropped);

void ClpCholeskyCtriRec(ClpCholeskyDenseC *thisStruct,
  longDouble *aTri, int nThis,
  longDouble *aUnder, longDouble *diagonal,
  longDouble *work,
  int nLeft, int iBlock, int jBlock,
  int numberBlocks);
void ClpCholeskyCrecTri(ClpCholeskyDenseC *thisStruct,
  longDouble *aUnder, int nTri, int nDo,
  int iBlock, int jBlock, longDouble *aTri,
  longDouble *diagonal, longDouble *work,
  int numberBlocks);
void ClpCholeskyCrecRec(ClpCholeskyDenseC *thisStruct,
  longDouble *above, int nUnder, int nUnderK,
  int nDo, longDouble *aUnder, longDouble *aOther,
  longDouble *work,
  int iBlock, int jBlock,
  int numberBlocks);
void ClpCholeskyCfactorLeaf(ClpCholeskyDenseC *thisStruct,
  longDouble *a, int n,
  longDouble *diagonal, longDouble *work,
  int *rowsDropped);
void ClpCholeskyCtriRecLeaf(/*ClpCholeskyDenseC * thisStruct,*/
  longDouble *aTri, longDouble *aUnder,
  longDouble *diagonal, longDouble *work,
  int nUnder);
void ClpCholeskyCrecTriLeaf(/*ClpCholeskyDenseC * thisStruct, */
  longDouble *aUnder, longDouble *aTri,
  /*longDouble * diagonal,*/ longDouble *work, int nUnder);
void ClpCholeskyCrecRecLeaf(/*ClpCholeskyDenseC * thisStruct, */
  const longDouble *COIN_RESTRICT above,
  const longDouble *COIN_RESTRICT aUnder,
  longDouble *COIN_RESTRICT aOther,
  const longDouble *COIN_RESTRICT work,
  int nUnder);
#endif

/* vi: softtabstop=2 shiftwidth=2 expandtab tabstop=2
*/