HROM_error_post_sampling.h

#ifndef __HROM_ERROR_POST_SAMPLING_H__
#define __HROM_ERROR_POST_SAMPLING_H__

#include "tests.h"
#include "parameters/all_parameters.h"
#include "dg/dg_base.hpp"
#include <eigen/Eigen/Dense>
#include <Epetra_CrsMatrix.h>
#include <Epetra_Map.h>
#include <Epetra_Vector.h>

namespace PHiLiP {
namespace Tests {

using Eigen::MatrixXd;
using Eigen::VectorXd;

template <int dim, int nstate>
class HROMErrorPostSampling: public TestsBase
{
public:
    HROMErrorPostSampling(const Parameters::AllParameters *const parameters_input,
                 const dealii::ParameterHandler &parameter_handler_input);
    
    Parameters::AllParameters reinit_params(std::string path) const;

    bool getWeightsFromFile(std::shared_ptr<DGBase<dim,double>> &dg) const;

    int run_test () const override;

    const dealii::ParameterHandler &parameter_handler;

    mutable MatrixXd snapshot_parameters;

    mutable MatrixXd rom_points;

    mutable std::shared_ptr<Epetra_Vector> ptr_weights;
};
} // End of Tests namespace
} // End of PHiLiP namespace

#endif