constrained.hh

Go to the documentation of this file.

 
#ifndef constrainedEV_hh
#define constrainedEV_hh
 
#include <memory>
#include <cstdlib>
#include "eigens.hh"
#include "solverFabric.hh"
#include "space/space.hh"
#include "operator/compositions.hh"
#include "operator/sparseMatrix.hh"
#include "operator/permutation.hh"
#include "sparseqr/givensRotations.hh"
 
namespace eigensolver {
 
  // *********************************************************** Constrained **
 
  class Constrained : public EigenSolver<Real> {
  public:
    Constrained(concepts::Operator<Real>& A, concepts::Operator<Real>& M,
    const concepts::SparseMatrix<Real>& C, SolverFabric<Real>& solver) :
      solver_(solver), computed_(false), A_(A), M_(M), C_(C),
      lambda_(0), eigenvectors_(0), ev_(0), kmax_(0),
      Pr(0), Prt(nullptr), Qgt(0), Qg(0), ext(0), restr(0), PrQgt(nullptr),
      QgPrt(nullptr), restrPrQgt(nullptr), QgPrtext(nullptr), restrPrQgtA(nullptr),
      restrPrQgtAQgPrtext(nullptr), restrPrQgtM(nullptr), restrPrQgtMQgPrtext(nullptr),
      it_(0), k_conv_(0) {}
    virtual ~Constrained();
    virtual const concepts::Array<Real>& getEV();
    virtual const concepts::Array<concepts::Vector<Real>*>& getEF();
    virtual uint iterations() const { return it_; }
    virtual uint converged() const { return k_conv_; }
  protected:
    virtual std::ostream& info(std::ostream& os) const;
  private:
    SolverFabric<Real>& solver_;
 
    void compute_();
    bool computed_;
 
    concepts::Operator<Real>& A_;
    concepts::Operator<Real>& M_;
    const concepts::SparseMatrix<Real>& C_;
 
    concepts::Array<Real> lambda_;
    concepts::Array<Real> eigenvectors_;
    concepts::Array<concepts::Vector<Real>*> ev_;
 
    uint kmax_;
 
    concepts::Permutation<Real>* Pr;
    std::unique_ptr<concepts::Permutation<Real> > Prt;
    sparseqr::GivensRotations<Real>* Qgt, * Qg;
    concepts::TrivExtendRestrict<Real>* ext, * restr;
    std::unique_ptr<concepts::Compose<Real> > PrQgt, QgPrt, restrPrQgt,
      QgPrtext, restrPrQgtA, restrPrQgtAQgPrtext, restrPrQgtM,
      restrPrQgtMQgPrtext;
 
    uint it_, k_conv_;
  };
 
} // namespace eigensolver
 
#endif // constrainedEV_hh