NewmarkSolverTest Class Reference

Inheritance diagram for NewmarkSolverTest:

Public Member Functions
	CPPUNIT_TEST_SUITE (NewmarkSolverTest)
	CPPUNIT_TEST (testNewmarkSolverConstantSecondOrderODESecondOrderStyle)
	CPPUNIT_TEST (testNewmarkSolverLinearTimeSecondOrderODESecondOrderStyle)
	CPPUNIT_TEST (testNewmarkSolverConstantSecondOrderODEFirstOrderStyle)
	CPPUNIT_TEST (testNewmarkSolverLinearTimeSecondOrderODEFirstOrderStyle)
	CPPUNIT_TEST_SUITE_END ()
void	testNewmarkSolverConstantSecondOrderODESecondOrderStyle ()
void	testNewmarkSolverLinearTimeSecondOrderODESecondOrderStyle ()
void	testNewmarkSolverConstantSecondOrderODEFirstOrderStyle ()
void	testNewmarkSolverLinearTimeSecondOrderODEFirstOrderStyle ()

Protected Member Functions
virtual void	aux_time_solver_init (NewmarkSolver &time_solver) override
void	set_beta (Real beta)
void	run_test_with_exact_soln (Real deltat, unsigned int n_timesteps)

Protected Attributes
Real	_beta

Detailed Description

Definition at line 82 of file second_order_unsteady_solver_test.C.

Member Function Documentation

aux_time_solver_init()

virtual void NewmarkSolverTestBase::aux_time_solver_init ( NewmarkSolver &  time_solver )

inlineoverrideprotectedvirtualinherited

Reimplemented from TimeSolverTestImplementation< NewmarkSolver >.

Definition at line 72 of file second_order_unsteady_solver_test.C.

  { time_solver.set_beta(_beta);
    time_solver.compute_initial_accel(); }

References NewmarkSolverTestBase::_beta.

CPPUNIT_TEST() [1/4]

NewmarkSolverTest::CPPUNIT_TEST

(

testNewmarkSolverConstantSecondOrderODEFirstOrderStyle

)

CPPUNIT_TEST() [2/4]

NewmarkSolverTest::CPPUNIT_TEST

(

testNewmarkSolverConstantSecondOrderODESecondOrderStyle

)

CPPUNIT_TEST() [3/4]

NewmarkSolverTest::CPPUNIT_TEST

(

testNewmarkSolverLinearTimeSecondOrderODEFirstOrderStyle

)

CPPUNIT_TEST() [4/4]

NewmarkSolverTest::CPPUNIT_TEST

(

testNewmarkSolverLinearTimeSecondOrderODESecondOrderStyle

)

CPPUNIT_TEST_SUITE()

NewmarkSolverTest::CPPUNIT_TEST_SUITE

(

NewmarkSolverTest

)

CPPUNIT_TEST_SUITE_END()

NewmarkSolverTest::CPPUNIT_TEST_SUITE_END

(

)

run_test_with_exact_soln()

void TimeSolverTestImplementation< NewmarkSolver >::run_test_with_exact_soln ( Real  deltat, unsigned int  n_timesteps  )

inlineprotectedinherited

Definition at line 27 of file time_solver_test_common.h.

  {
    Mesh mesh(*TestCommWorld);
    MeshTools::Generation::build_point(mesh);
    EquationSystems es(mesh);
    SystemType & system = es.add_system<SystemType>("ScalarSystem");
 
    system.time_solver = libmesh_make_unique<TimeSolverType>(system);
 
    es.init();
 
    DiffSolver & solver = *(system.time_solver->diff_solver().get());
    solver.relative_step_tolerance = std::numeric_limits<Real>::epsilon()*10;
    solver.relative_residual_tolerance = std::numeric_limits<Real>::epsilon()*10;
    solver.absolute_residual_tolerance = std::numeric_limits<Real>::epsilon()*10;
 
    NewtonSolver & newton = cast_ref<NewtonSolver &>(solver);
 
    // LASPACK GMRES + ILU defaults don't like these problems, so
    // we'll use a sophisticated "just divide the scalars" solver instead.
    newton.get_linear_solver().set_solver_type(JACOBI);
    newton.get_linear_solver().set_preconditioner_type(IDENTITY_PRECOND);
 
    system.deltat = deltat;
 
    TimeSolverType * time_solver = cast_ptr<TimeSolverType *>(system.time_solver.get());
    this->aux_time_solver_init(*time_solver);
 
    // We're going to want to check our solution, and when we run
    // "make check" with LIBMESH_RUN='mpirun -np N" for N>1 then we'll
    // need to keep that check in sync with the processors that are just
    // twiddling their thumbs, not owning our mesh point.
    std::vector<dof_id_type> solution_index;
    solution_index.push_back(0);
    const bool has_solution = system.get_dof_map().all_semilocal_indices(solution_index);
 
    for (unsigned int t_step=0; t_step != n_timesteps; ++t_step)
      {
        system.solve();
        system.time_solver->advance_timestep();
 
        Real rel_error = 0;
 
        if (has_solution)
          {
            Number exact_soln = system.u(system.time);
            rel_error =  std::abs((exact_soln - (*system.solution)(0))/exact_soln);
          }
        system.comm().max(rel_error);
 
        // Using relative error for comparison, so "exact" is 0
        LIBMESH_ASSERT_FP_EQUAL( 0.0,
                                 rel_error,
                                 std::numeric_limits<Real>::epsilon()*10 );
      }
  }

set_beta()

void NewmarkSolverTestBase::set_beta ( Real  beta )

inlineprotectedinherited

Definition at line 76 of file second_order_unsteady_solver_test.C.

  { _beta = beta; }

References NewmarkSolverTestBase::_beta.

Referenced by testNewmarkSolverLinearTimeSecondOrderODEFirstOrderStyle(), and testNewmarkSolverLinearTimeSecondOrderODESecondOrderStyle().

testNewmarkSolverConstantSecondOrderODEFirstOrderStyle()

void NewmarkSolverTest::testNewmarkSolverConstantSecondOrderODEFirstOrderStyle ( )

inline

Definition at line 113 of file second_order_unsteady_solver_test.C.

  {
    this->run_test_with_exact_soln<ConstantSecondOrderODE<SecondOrderScalarSystemFirstOrderTimeSolverBase>>(0.5,10);
  }

testNewmarkSolverConstantSecondOrderODESecondOrderStyle()

void NewmarkSolverTest::testNewmarkSolverConstantSecondOrderODESecondOrderStyle ( )

inline

Definition at line 99 of file second_order_unsteady_solver_test.C.

  {
    this->run_test_with_exact_soln<ConstantSecondOrderODE<SecondOrderScalarSystemSecondOrderTimeSolverBase>>(0.5,10);
  }

testNewmarkSolverLinearTimeSecondOrderODEFirstOrderStyle()

void NewmarkSolverTest::testNewmarkSolverLinearTimeSecondOrderODEFirstOrderStyle ( )

inline

Definition at line 118 of file second_order_unsteady_solver_test.C.

  {
    // For \beta = 1/6, we have the "linear acceleration method" for which
    // we should be able to exactly integrate linear (in time) acceleration
    // functions.
    this->set_beta(Real(1)/Real(6));
    this->run_test_with_exact_soln<LinearTimeSecondOrderODE<SecondOrderScalarSystemFirstOrderTimeSolverBase>>(0.5,10);
  }

References libMesh::Real, and NewmarkSolverTestBase::set_beta().

testNewmarkSolverLinearTimeSecondOrderODESecondOrderStyle()

void NewmarkSolverTest::testNewmarkSolverLinearTimeSecondOrderODESecondOrderStyle ( )

inline

Definition at line 104 of file second_order_unsteady_solver_test.C.

  {
    // For \beta = 1/6, we have the "linear acceleration method" for which
    // we should be able to exactly integrate linear (in time) acceleration
    // functions.
    this->set_beta(Real(1)/Real(6));
    this->run_test_with_exact_soln<LinearTimeSecondOrderODE<SecondOrderScalarSystemSecondOrderTimeSolverBase>>(0.5,10);
  }

References libMesh::Real, and NewmarkSolverTestBase::set_beta().

Member Data Documentation

_beta

Real NewmarkSolverTestBase::_beta

protectedinherited

Definition at line 79 of file second_order_unsteady_solver_test.C.

Referenced by NewmarkSolverTestBase::aux_time_solver_init(), and NewmarkSolverTestBase::set_beta().

---

The documentation for this class was generated from the following file:

• tests/solvers/second_order_unsteady_solver_test.C