Tempus Version of the Day
Time Integration
Loading...
Searching...
No Matches
Tempus_UnitTest_SDIRK_2Stage3rdOrder.cpp
Go to the documentation of this file.
1// @HEADER
2// ****************************************************************************
3// Tempus: Copyright (2017) Sandia Corporation
4//
5// Distributed under BSD 3-clause license (See accompanying file Copyright.txt)
6// ****************************************************************************
7// @HEADER
8
10
11
12namespace Tempus_Unit_Test {
13
14using Teuchos::RCP;
15using Teuchos::rcp;
16using Teuchos::rcp_const_cast;
17using Teuchos::rcp_dynamic_cast;
18using Teuchos::ParameterList;
19using Teuchos::sublist;
20
21
22// ************************************************************
23// ************************************************************
24TEUCHOS_UNIT_TEST(SDIRK_2Stage3rdOrder, Default_Construction)
25{
26 auto stepper = rcp(new Tempus::StepperSDIRK_2Stage3rdOrder<double>());
28
29 // Test stepper properties.
30 TEUCHOS_ASSERT(stepper->getOrder() == 3);
31 std::string gammaType = "3rd Order A-stable";
32 TEUCHOS_ASSERT(stepper->getGammaType() == gammaType);
33 double gamma = 0.7886751345948128;
34 TEUCHOS_ASSERT(stepper->getGamma() == gamma);
35 stepper->setGammaType(gammaType); stepper->initialize(); TEUCHOS_TEST_FOR_EXCEPT(!stepper->isInitialized());
36 stepper->setGamma(gamma); stepper->initialize(); TEUCHOS_TEST_FOR_EXCEPT(!stepper->isInitialized());
37
38}
39
40
41// ************************************************************
42// ************************************************************
43TEUCHOS_UNIT_TEST(SDIRK_2Stage3rdOrder, StepperFactory_Construction)
44{
45 auto model = rcp(new Tempus_Test::SinCosModel<double>());
46 testFactoryConstruction("SDIRK 2 Stage 3rd order", model);
47}
48
49
50// ************************************************************
51// ************************************************************
52TEUCHOS_UNIT_TEST(SDIRK_2Stage3rdOrder, AppAction)
53{
54 auto stepper = rcp(new Tempus::StepperSDIRK_2Stage3rdOrder<double>());
55 auto model = rcp(new Tempus_Test::SinCosModel<double>());
56 testRKAppAction(stepper, model, out, success);
57}
58
59
60} // namespace Tempus_Test
Sine-Cosine model problem from Rythmos. This is a canonical Sine-Cosine differential equation.
void testDIRKAccessorsFullConstruction(const RCP< Tempus::StepperDIRK< double > > &stepper)
Unit test utility for ExplicitRK Stepper construction and accessors.
void testRKAppAction(const Teuchos::RCP< Tempus::StepperRKBase< double > > &stepper, const Teuchos::RCP< const Thyra::ModelEvaluator< double > > &model, Teuchos::FancyOStream &out, bool &success)
Unit test utility for Stepper RK AppAction.
TEUCHOS_UNIT_TEST(BackwardEuler, Default_Construction)
void testFactoryConstruction(std::string stepperType, const Teuchos::RCP< const Thyra::ModelEvaluator< double > > &model)
Unit test utility for Stepper construction through StepperFactory.