tests.h

#ifndef __TESTS_H__
#define __TESTS_H__

#include "parameters/all_parameters.h"
#include <string>
#include <deal.II/grid/tria.h>
#include <deal.II/base/conditional_ostream.h>
#include <deal.II/base/parameter_handler.h>

//#include "dg/dg.h"
namespace PHiLiP {
namespace Tests {


class TestsBase
{
public:
    const Parameters::AllParameters *const all_parameters; 
    //const int ndim; ///< Number of dimensions. Run-time variable instead of compile-time constant.
    //const int nstate; ///< Number of state variables. Run-time variable instead of compile-time constant.

    TestsBase () = delete;

    explicit TestsBase(const Parameters::AllParameters *const parameters_input);

    virtual ~TestsBase() = default;


    virtual int run_test() const = 0;
protected:
    const MPI_Comm mpi_communicator; 
    const int mpi_rank; 
    const int n_mpi; 

    dealii::ConditionalOStream pcout;


    std::vector<int> get_number_1d_cells(const int ngrids) const;

    // /// Evaluates the number of cells to generate the grids for 1D grid based on input file.
    // void globally_refine_and_interpolate(DGBase<dim, double> &dg) const;

    std::string get_pde_string(const Parameters::AllParameters *const param) const;

    std::string get_conv_num_flux_string(const Parameters::AllParameters *const param) const;

    std::string get_diss_num_flux_string(const Parameters::AllParameters *const param) const;

    std::string get_manufactured_solution_string(const Parameters::AllParameters *const param) const;
};

template<int dim, int nstate, typename MeshType = dealii::Triangulation<dim>>
class TestsFactory
{
public:

    static std::unique_ptr< TestsBase > create_test(
        const Parameters::AllParameters *const parameters_input,
        dealii::ParameterHandler &parameter_handler_input);


    static std::unique_ptr< TestsBase > select_mesh(
        const Parameters::AllParameters *const parameters_input,
        dealii::ParameterHandler &parameter_handler_input);


    static std::unique_ptr< TestsBase > select_test(
        const Parameters::AllParameters *const parameters_input,
        dealii::ParameterHandler &parameter_handler_input);
};

} // Tests namespace
} // PHiLiP namespace
#endif