operator/spooles.hh

Go to the documentation of this file.

// interface to SPOOLES

#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