homogeneous_isotropic_turbulence_initialization_check.h

#ifndef __homogeneous_isotropic_turbulence_initialization_check__
#define __homogeneous_isotropic_turbulence_initialization_check__

#include "dg/dg_base.hpp"
#include "parameters/all_parameters.h"
#include "tests.h"

namespace PHiLiP {
namespace Tests {

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

    const dealii::ParameterHandler &parameter_handler;
    
    const double kinetic_energy_expected;

    int run_test () const override;
};

} // End of Tests namespace
} // End of PHiLiP namespace

#endif