spooles.hh

Go to the documentation of this file.

 
#ifndef SPOOLES_HH
#define SPOOLES_HH
 
#include "basics/typedefs.hh"
#include "compositions.hh"
 
namespace concepts {
 
  using concepts::Real;
 
  // forwards
  template<typename F>
  class SparseMatrix;
 
  // *************************************************************** Spooles **
 
  class Spooles : public Operator<Real> {
  public:
    Spooles(const SparseMatrix<Real>& A, bool symmetric = false);
    virtual ~Spooles();
    virtual void operator()(const Function<Real>& fncY, Function<Real>& fncX);
    void operator()(const Vector<Real>& fncY, Vector<Real>& fncX);
    virtual const Space<Real>& spaceX() const { return spcX_; }
    virtual const Space<Real>& spaceY() const { return spcY_; }
    uint iterations() const { return 0; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    const uint n_;
    const Space<Real>& spcX_;
    const Space<Real>& spcY_;
    bool symmetric_;
    bool factored_;
 
    void* mtxA, * frontmtx, * mtxmanager, * chvmanager, * rootchv,
      * frontETree, * graph, * newToOldIV, * oldToNewIV, * adjIVL,
      * symbfacIVL, * bridge;
    int* oldToNew, * newToOld;
    int rc;
  };
 
} // namespace concepts
 
#endif // SPOOLES_HH