hrom_test_location.h

#ifndef __HROM_TEST_LOCATION__
#define __HROM_TEST_LOCATION__

#include "parameters/all_parameters.h"
#include "pod_basis_base.h"
#include "reduced_order_solution.h"
#include <eigen/Eigen/Dense>
#include "test_location_base.h"

namespace PHiLiP {
namespace ProperOrthogonalDecomposition {
using Eigen::RowVectorXd;

// Based very closely on the ROMTestLocation class

/*
Based on the work in Donovan Blais' thesis:
Goal-Oriented Adaptive Sampling for Projection-Based Reduced-Order Models, 2022

Details on the ROM points/errors can be found in sections 5 and 6

Derivation of the new error indicator will likely be detailed in Calista Biondic's thesis
*/
template <int dim, int nstate>
class HROMTestLocation: public TestLocationBase<dim,nstate>
{
public:
    HROMTestLocation(const RowVectorXd& parameter, std::unique_ptr<ROMSolution < dim, nstate>> rom_solution, std::shared_ptr< DGBase<dim, double> > dg_input, Epetra_Vector weights);

    void compute_initial_rom_to_final_rom_error(std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim>> pod_updated) override;

    std::shared_ptr<Epetra_CrsMatrix> generate_hyper_reduced_jacobian(const Epetra_CrsMatrix &system_matrix);

    std::shared_ptr<Epetra_CrsMatrix> generate_test_basis(Epetra_CrsMatrix &epetra_system_matrix, const Epetra_CrsMatrix &pod_basis);

    std::shared_ptr<DGBase<dim,double>> dg;

    Epetra_Vector ECSW_weights;

};

}
}


#endif