first_order_unsteady_solver_test.C

Go to the documentation of this file.

#include <libmesh/equation_systems.h>
#include <libmesh/mesh.h>
#include <libmesh/mesh_generation.h>
#include <libmesh/quadrature.h>
#include <libmesh/diff_solver.h>
#include <libmesh/euler_solver.h>
#include <libmesh/euler2_solver.h>
 
#include "solvers/time_solver_test_common.h"
 
 
class ConstantFirstOrderODE : public FirstOrderScalarSystemBase
{
public:
  ConstantFirstOrderODE(EquationSystems & es,
                        const std::string & name_in,
                        const unsigned int number_in)
    : FirstOrderScalarSystemBase(es, name_in, number_in)
  {}
 
  virtual Number F( FEMContext & /*context*/, unsigned int /*qp*/ )
  { return 5.0; }
 
  virtual Number M( FEMContext & /*context*/, unsigned int /*qp*/ )
  { return Real(21)/10; }
 
  virtual Number u( Real t )
  { return Real(50)/21*t; }
};
 
class LinearTimeFirstOrderODE : public FirstOrderScalarSystemBase
{
public:
  LinearTimeFirstOrderODE(EquationSystems & es,
                          const std::string & name_in,
                          const unsigned int number_in)
    : FirstOrderScalarSystemBase(es, name_in, number_in)
  {}
 
  virtual Number F( FEMContext & context, unsigned int /*qp*/ )
  { return 2.0*context.get_time(); }
 
  virtual Number M( FEMContext & /*context*/, unsigned int /*qp*/ )
  { return 5.0; }
 
  virtual Number u( Real t )
  { return 1/Real(5)*t*t; }
};
 
template<typename TimeSolverType>
class ThetaSolverTestBase : public TimeSolverTestImplementation<TimeSolverType>
{
public:
  ThetaSolverTestBase()
    : TimeSolverTestImplementation<TimeSolverType>(),
    _theta(1.0)
  {}
 
protected:
 
  virtual void aux_time_solver_init( TimeSolverType & time_solver )
  { time_solver.theta = _theta; }
 
  void set_theta( Real theta )
  { _theta = theta; }
 
  Real _theta;
};
 
class EulerSolverTest : public CppUnit::TestCase,
                        public ThetaSolverTestBase<EulerSolver>
{
public:
  LIBMESH_CPPUNIT_TEST_SUITE( EulerSolverTest );
 
#ifdef LIBMESH_HAVE_SOLVER
  CPPUNIT_TEST( testEulerSolverConstantFirstOrderODE );
  CPPUNIT_TEST( testEulerSolverLinearTimeFirstOrderODE );
#endif
 
  CPPUNIT_TEST_SUITE_END();
 
public:
 
  void testEulerSolverConstantFirstOrderODE()
  {
    LOG_UNIT_TEST;
 
    this->set_theta(1.0);
    this->run_test_with_exact_soln<ConstantFirstOrderODE>(0.5,10);
 
    this->set_theta(0.5);
    this->run_test_with_exact_soln<ConstantFirstOrderODE>(0.5,10);
  }
 
  void testEulerSolverLinearTimeFirstOrderODE()
  {
    LOG_UNIT_TEST;
 
    // Need \theta = 0.5 since this has t in F.
    this->set_theta(0.5);
    this->run_test_with_exact_soln<LinearTimeFirstOrderODE>(0.5,10);
  }
 
};
 
class Euler2SolverTest : public CppUnit::TestCase,
                         public ThetaSolverTestBase<Euler2Solver>
{
public:
  LIBMESH_CPPUNIT_TEST_SUITE( Euler2SolverTest );
 
#ifdef LIBMESH_HAVE_SOLVER
  CPPUNIT_TEST( testEuler2SolverConstantFirstOrderODE );
  CPPUNIT_TEST( testEuler2SolverLinearTimeFirstOrderODE );
#endif
 
  CPPUNIT_TEST_SUITE_END();
 
public:
  void testEuler2SolverConstantFirstOrderODE()
  {
    LOG_UNIT_TEST;
 
    this->set_theta(1.0);
    this->run_test_with_exact_soln<ConstantFirstOrderODE>(0.5,10);
 
    this->set_theta(0.5);
    this->run_test_with_exact_soln<ConstantFirstOrderODE>(0.5,10);
  }
 
  void testEuler2SolverLinearTimeFirstOrderODE()
  {
    LOG_UNIT_TEST;
 
    // Need \theta = 0.5 since this has t in F.
    this->set_theta(0.5);
    this->run_test_with_exact_soln<LinearTimeFirstOrderODE>(0.5,10);
  }
};
 
CPPUNIT_TEST_SUITE_REGISTRATION( EulerSolverTest );
CPPUNIT_TEST_SUITE_REGISTRATION( Euler2SolverTest );