hyper_reduction_comparison.h

#ifndef __HYPER_REDUCTION_COMPARISON_H__
#define __HYPER_REDUCTION_COMPARISON_H__

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

namespace PHiLiP {
namespace Tests {

using Eigen::MatrixXd;

template <int dim, int nstate>
class HyperReductionComparison: public TestsBase
{
public:
    HyperReductionComparison(const Parameters::AllParameters *const parameters_input,
                 const dealii::ParameterHandler &parameter_handler_input);
    
    Parameters::AllParameters reinit_params(const int max_iter) 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 std::shared_ptr<Epetra_Vector> ptr_weights;
};
} // End of Tests namespace
} // End of PHiLiP namespace

#endif