pod_basis_online.h

#ifndef __POD_BASIS_ONLINE__
#define __POD_BASIS_ONLINE__

#include <deal.II/base/conditional_ostream.h>
#include <deal.II/lac/full_matrix.h>
#include <deal.II/lac/la_parallel_vector.h>
#include <deal.II/lac/trilinos_sparse_matrix.h>
#include <deal.II/lac/vector_operation.h>
#include <deal.II/numerics/vector_tools.h>

#include <eigen/Eigen/Dense>

#include "dg/dg_base.hpp"
#include "parameters/all_parameters.h"
#include "pod_basis_base.h"

namespace PHiLiP {
namespace ProperOrthogonalDecomposition {
using Eigen::MatrixXd;
using Eigen::VectorXd;

template <int dim>
class OnlinePOD: public PODBase<dim>
{
public:
    explicit OnlinePOD(std::shared_ptr<dealii::TrilinosWrappers::SparseMatrix> _system_matrix);

    std::shared_ptr<dealii::TrilinosWrappers::SparseMatrix> getPODBasis() override;

    dealii::LinearAlgebra::ReadWriteVector<double> getReferenceState() override;

    MatrixXd getSnapshotMatrix() override;

    void addSnapshot(dealii::LinearAlgebra::distributed::Vector<double> snapshot);

    void computeBasis();

    std::shared_ptr<dealii::TrilinosWrappers::SparseMatrix> basis;

    dealii::LinearAlgebra::ReadWriteVector<double> referenceState;

    std::shared_ptr<dealii::TrilinosWrappers::SparseMatrix> system_matrix;

    dealii::LAPACKFullMatrix<double> dealiiSnapshotMatrix;

    MatrixXd snapshotMatrix;

    const MPI_Comm mpi_communicator; 
    const int mpi_rank; 


    dealii::ConditionalOStream pcout;


};

}
}

#endif