test_location_base.h

#ifndef __TEST_LOCATION_BASE__
#define __TEST_LOCATION_BASE__

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

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

template <int dim, int nstate>
class TestLocationBase
{
public:
    TestLocationBase(const RowVectorXd& parameter, std::unique_ptr<ROMSolution < dim, nstate>> rom_solution);

    virtual ~TestLocationBase () {};

    void compute_FOM_to_initial_ROM_error();

    virtual void compute_initial_rom_to_final_rom_error(std::shared_ptr<ProperOrthogonalDecomposition::PODBase<dim>> pod_updated) = 0;

    void compute_total_error();

    RowVectorXd parameter;

    std::unique_ptr<ROMSolution<dim, nstate>> rom_solution;

    double fom_to_initial_rom_error;

    double initial_rom_to_final_rom_error;

    double total_error;

    const MPI_Comm mpi_communicator; 
    const int mpi_rank; 


    dealii::ConditionalOStream pcout;

};

}
}


#endif