operator/constrained.hh

Go to the documentation of this file.

// solves a constrained source problem

#ifndef constrainedSP_hh
#define constrainedSP_hh

#include <memory>
#include "compositions.hh"
#include "sparseMatrix.hh"

namespace concepts {

  // forward declaration
  template<typename F>
  class SolverFabric;

  // *********************************************************** Constrained **

  class Constrained : public Operator<Real> {
  public:
    Constrained(SparseMatrix<Real>& S, SparseMatrix<Real>& C,
                const Vector<Real>& g, SolverFabric<Real>& fabric);
    virtual void operator()(const Function<Real>& fncY, Function<Real>& fncX);
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    SparseMatrix<Real>& C_, &S_;
    SparseMatrix<Real> Ct_;
    std::auto_ptr<SparseMatrix<Real> > CCt_;
    const Vector<Real>& g_;
    Vector<Real> SRg_;
    std::auto_ptr<Operator<Real> > CCtinv_, solver_;
    std::auto_ptr<concepts::Compose<Real> > R_, Qmin_, CtC_, QS_, QSQ_;
    std::auto_ptr<LiCoI<Real> > Q_;
    std::auto_ptr<LiCo<Real> > Sp_;
    SolverFabric<Real>& fabric_;
  };

} // namespace concepts

#endif // constrainedSP_hh