pod_basis_base.h

#ifndef __POD_BASIS_INTERFACE__
#define __POD_BASIS_INTERFACE__

#include <deal.II/numerics/vector_tools.h>
#include <deal.II/lac/la_parallel_vector.h>
#include <deal.II/lac/trilinos_sparse_matrix.h>
#include <eigen/Eigen/Dense>

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

template <int dim>
class PODBase
{
public:
    virtual ~PODBase() = default;

    virtual std::shared_ptr<dealii::TrilinosWrappers::SparseMatrix> getPODBasis() = 0;

    virtual dealii::LinearAlgebra::ReadWriteVector<double> getReferenceState() = 0;

    virtual MatrixXd getSnapshotMatrix() = 0;
};

}
}


#endif