lapack.hh

Go to the documentation of this file.

 
#pragma once
 
#include "basics/typedefs.hh"
#include "operator/compositions.hh"
 
namespace concepts {
  
  template<class F>
  class DenseMatrix;
 
  template<class F>
  class DiagonalMatrix;
 
  // ************************************************************ LapackChol **
 
  class LapackChol : public VecOperator<Real> {
  public:
    // Flag, if the upper or the lower triangular part is stored.
    enum UpLo {U = 0, L = 1};
 
    LapackChol(const ElementMatrix<Real> A, enum UpLo uplo = U);
 
 
    // Ax= b now solve  D*Ax = b where D is a diagonal matrix that is a variable at each apply step
    // this leads to a change of the rhs in the routine
    //@param Dinv is the inverse of the diagonal matrix D
    virtual void apply(DiagonalMatrix<Real>& Dinv, const Vector<Real>& fncY, Vector<Real>& fncX);
 
 
  protected:
    virtual std::ostream& info(std::ostream& os) const;
    virtual void apply_(const Vector<Real>& fncY, Vector<Real>& fncX);
    virtual void apply_();
  private:
    mutable ElementMatrix<Real> A_;
    char uplo_;
    mutable bool factorized_;
  };
} // namespace concepts