umfpack.hh

Go to the documentation of this file.

 
#ifndef UMFPACK_HH
#define UMFPACK_HH
 
#include "basics/typedefs.hh"
#include "compositions.hh"
 
namespace concepts {
 
  // forwards
  template<typename F>
  class SparseMatrix;
 
  // *************************************************************** Umfpack **
 
  class Umfpack : public Operator<Real> {
  public:
    Umfpack(const SparseMatrix<Real>& A, bool symmetric = false);
    virtual ~Umfpack();
    virtual void operator()(const Function<Real>& fncY, Function<Real>& fncX);
    void operator()(const Vector<Real>& fncY, Vector<Real>& fncX);
    uint iterations() const { return 0; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    const uint n_;
 
    concepts::Array<long int> Ap_, Ai_;
    concepts::Array<Real> Ax_;
    void* symbolic_, * numeric_;
 
    bool factored_;
  };
 
} // namespace concepts
 
#endif // UMFPACK_HH