constrained.hh

Go to the documentation of this file.

 
#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::unique_ptr<SparseMatrix<Real> > CCt_;
    const Vector<Real>& g_;
    Vector<Real> SRg_;
    std::unique_ptr<Operator<Real> > CCtinv_, solver_;
    std::unique_ptr<concepts::Compose<Real> > R_, Qmin_, CtC_, QS_, QSQ_;
    std::unique_ptr<LiCoI<Real> > Q_;
    std::unique_ptr<LiCo<Real> > Sp_;
    SolverFabric<Real>& fabric_;
  };
 
} // namespace concepts
 
#endif // constrainedSP_hh