RBGen_AnasaziPOD.h

Go to the documentation of this file.

#ifndef RBGEN_ANASAZI_POD_H
#define RBGEN_ANASAZI_POD_H

#include "RBGen_PODMethod.hpp"
#include "RBGen_Method.hpp"
#include "Epetra_MultiVector.h"
#include "Epetra_Operator.h"
#include "Teuchos_ParameterList.hpp"

namespace RBGen {
  
  class AnasaziPOD : public virtual Method<Epetra_MultiVector,Epetra_Operator>, public virtual PODMethod<double> {
    
  public:


    AnasaziPOD();

    virtual ~AnasaziPOD() {};




    void computeBasis();
    

    void updateBasis( const Teuchos::RCP< Epetra_MultiVector >& update_ss ) {};



    
    Teuchos::RCP<const Epetra_MultiVector> getBasis() const { return basis_; }

    std::vector<double> getSingularValues() const { return sv_; }

    double getCompTime() const { return comp_time_; }
    

    
    void Initialize( const Teuchos::RCP< Teuchos::ParameterList >& params,
                     const Teuchos::RCP< const Epetra_MultiVector >& ss,
         const Teuchos::RCP< RBGen::FileIOHandler< Epetra_Operator > >& fileio = Teuchos::null );

    void Reset( const Teuchos::RCP<Epetra_MultiVector>& new_ss ) { ss_ = new_ss; }

    void ResetOp( const Teuchos::RCP<Epetra_Operator>& new_op ) { op_ = new_op; }

    void ResetInnerProdOp( const Teuchos::RCP<Epetra_Operator>& new_op ) { inner_prod_op_ = new_op; }
   



    bool isInitialized() { return isInitialized_; }
 

  private:

    // Is this object initialized.
    bool isInitialized_;

    // Is the inner (A^T*A) or outer (A*A^T) product being used for the SVD computation
    bool isInner_;

    // Size of the basis that this method will compute.
    int basis_size_;

    // Computational time (in seconds, using wall clock).
    double comp_time_;

    // Pointers to the snapshots and reduced basis.
    Teuchos::RCP<const Epetra_MultiVector> ss_;
    Teuchos::RCP<Epetra_MultiVector> basis_;

    // Pointer to the SVD weighted operator
    Teuchos::RCP<Epetra_Operator> op_;

    // Pointer to the inner product operator
    Teuchos::RCP<Epetra_Operator> inner_prod_op_;

    // Vector holding singular values.
    std::vector<double> sv_;
  };
  
} // end of RBGen namespace

#endif // RBGEN_ANASAZI_POD_H