superLU.hh

Go to the documentation of this file.

 
#ifndef SuperLU_hh
#define SuperLU_hh
 
#include "compositions.hh"
#include "basics/typedefs.hh"
 
namespace concepts {
 
  // *************************************************************** SuperLU **
 
  template<class F>
  class SuperLU : public VecOperator<F> {
  public:
    enum order { NATURAL = 0, MMD_PROD = 1, MMD_SUM = 2, COLAMD = 3 };
    SuperLU(const Operator<F>& A, enum order ispec = COLAMD);
    virtual ~SuperLU();
 
    uint iterations() const { return 0; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    void* A_;
 
    uint n_;
 
    bool factorized_;
    enum order ispec_;
 
    void* L_, * U_;
    int *perm_r_, *perm_c_;
 
    virtual void apply_(const Vector<F>& fncY, Vector<F>& fncX);
    virtual void apply_();
  };
 
} // namespace concepts
 
#endif // SuperLU_hh