antioch-0.4.0
List of all members | Public Member Functions | Protected Member Functions | Protected Attributes
AntiochTesting::ArrheniusRateEigenDoubleTest Class Reference
Inheritance diagram for AntiochTesting::ArrheniusRateEigenDoubleTest:
AntiochTesting::ArrheniusRateEigenTest< double > AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > > AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >, Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >

Public Member Functions

 CPPUNIT_TEST_SUITE (ArrheniusRateEigenDoubleTest)
 
 CPPUNIT_TEST (test_standard_rate)
 
 CPPUNIT_TEST (test_standard_deriv)
 
 CPPUNIT_TEST (test_standard_rate_and_deriv)
 
 CPPUNIT_TEST_SUITE_END ()
 
virtual void setUp ()
 
virtual void tearDown ()
 
virtual void init ()
 
virtual void clear ()
 
void test_standard_rate ()
 
void test_standard_deriv ()
 
void test_standard_rate_and_deriv ()
 
void test_rate (const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &reaction_rate, const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &T, typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typetol)
 
void test_deriv (const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &reaction_rate, const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &T, typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typetol)
 
void test_rate_and_deriv (const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &reaction_rate, const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &T, typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typetol)
 

Protected Member Functions

virtual Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 > 
exact_rate (Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 >T)
 
virtual Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 > 
exact_deriv (Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 >T)
 
void reset_params (Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeCf, Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeEa)
 
void reset_params (Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeCf, Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeEa, Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeR)
 
Antioch::value_type
< Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 >
>::type 
value (Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeT)
 
Antioch::value_type
< Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 >
>::type 
deriv (Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeT)
 
Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 > 
setup_T (const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &example)
 

Protected Attributes

Antioch::ArrheniusRate
< typename Antioch::value_type
< Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 >
>::type > * 
_rate
 
Antioch::value_type
< Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 >
>::type 
_Cf
 
Antioch::value_type
< Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 >
>::type 
_Ea
 
Antioch::value_type
< Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 >
>::type 
_R
 
Eigen::Array< double,
2 *ANTIOCH_N_TUPLES, 1 > * 
_example
 

Detailed Description

Definition at line 77 of file arrhenius_rate_eigen_test.C.

Member Function Documentation

virtual void AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::clear ( )
inlinevirtualinherited

Definition at line 64 of file arrhenius_rate_vector_test_base.h.

65  {
66  delete _rate;
67  }
Antioch::ArrheniusRate< typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > * _rate
AntiochTesting::ArrheniusRateEigenDoubleTest::CPPUNIT_TEST ( test_standard_rate  )
AntiochTesting::ArrheniusRateEigenDoubleTest::CPPUNIT_TEST ( test_standard_deriv  )
AntiochTesting::ArrheniusRateEigenDoubleTest::CPPUNIT_TEST ( test_standard_rate_and_deriv  )
AntiochTesting::ArrheniusRateEigenDoubleTest::CPPUNIT_TEST_SUITE ( ArrheniusRateEigenDoubleTest  )
AntiochTesting::ArrheniusRateEigenDoubleTest::CPPUNIT_TEST_SUITE_END ( )
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::deriv ( Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  T)
inlineprotectedinherited

Definition at line 65 of file arrhenius_rate_test_helper.h.

66  {
67  using std::exp;
68  return _Ea/(_R*T*T)*_Cf *exp(-_Ea/(_R*T));
69  }
virtual Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::exact_deriv ( Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 >  T)
inlineprotectedvirtualinherited

Implements AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >, Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >.

Definition at line 111 of file arrhenius_rate_vector_test_base.h.

112  {
113  PairScalars e_deriv = *(this->_example);
114  for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
115  {
116  e_deriv[2*tuple] = this->deriv(T[2*tuple]);
117  e_deriv[2*tuple+1] = this->deriv(T[2*tuple+1]);
118  }
119  return e_deriv;
120  }
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type deriv(Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeT)
virtual Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::exact_rate ( Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 >  T)
inlineprotectedvirtualinherited

Implements AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >, Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >.

Definition at line 100 of file arrhenius_rate_vector_test_base.h.

101  {
102  PairScalars e_rate = *(this->_example);
103  for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
104  {
105  e_rate[2*tuple] = this->value(T[2*tuple]);
106  e_rate[2*tuple+1] = this->value(T[2*tuple+1]);
107  }
108  return e_rate;
109  }
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type value(Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeT)
virtual void AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::init ( )
inlinevirtualinherited

Definition at line 52 of file arrhenius_rate_vector_test_base.h.

53  {
54  typedef typename Antioch::value_type<PairScalars>::type Scalar;
55 
56  Scalar Cf = 1.4L;
57  Scalar Ea = 5.0L;
58  Scalar R = 1.0L; // Ea in K
59 
60  this->reset_params(Cf,Ea,R);
62  }
Antioch::ArrheniusRate< typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > * _rate
void reset_params(Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeCf, Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typeEa)
void AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::reset_params ( Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  Cf,
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  Ea 
)
inlineprotectedinherited

Definition at line 46 of file arrhenius_rate_test_helper.h.

47  {
48  _Cf = Cf;
49  _Ea = Ea;
50  }
void AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::reset_params ( Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  Cf,
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  Ea,
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  R 
)
inlineprotectedinherited

Definition at line 52 of file arrhenius_rate_test_helper.h.

53  {
54  _Cf = Cf;
55  _Ea = Ea;
56  _R = R;
57  }
virtual void AntiochTesting::ArrheniusRateEigenTest< double >::setUp ( )
inlinevirtualinherited
Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > , Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::setup_T ( const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &  example)
inlineprotectedinherited

Definition at line 123 of file reaction_rate_vector_test_base.h.

124  {
125  // Construct from example to avoid resizing issues
126  PairScalars T = example;
127  for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
128  {
129  T[2*tuple] = 1500.1;
130  T[2*tuple+1] = 1600.1;
131  }
132  return T;
133  }
virtual void AntiochTesting::ArrheniusRateEigenTest< double >::tearDown ( )
inlinevirtualinherited
void AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > , Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::test_deriv ( const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &  reaction_rate,
const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &  T,
typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  tol 
)
inlineinherited

Definition at line 64 of file reaction_rate_vector_test_base.h.

67  {
68  const PairScalars deriv = reaction_rate.derivative(T);
69  const PairScalars exact_deriv = this->exact_deriv(T);
70 
71  for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
72  {
73  CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple],
74  exact_deriv[0],
75  tol );
76 
77  CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple+1],
78  exact_deriv[1],
79  tol );
80 
81  }
82  }
_Cf VectorStateType VectorStateType derivative(const KineticsConditions< StateType, VectorStateType > &T) const ANTIOCH_AUTOFUNC(StateType
void AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > , Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::test_rate ( const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &  reaction_rate,
const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &  T,
typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  tol 
)
inlineinherited

Definition at line 44 of file reaction_rate_vector_test_base.h.

47  {
48  const PairScalars rate = reaction_rate(T);
49  const PairScalars exact_rate = this->exact_rate(T);
50 
51  for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
52  {
53  CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple],
54  exact_rate[0],
55  tol );
56 
57  CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple+1],
58  exact_rate[1],
59  tol );
60 
61  }
62  }
void AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > , Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::test_rate_and_deriv ( const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &  reaction_rate,
const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &  T,
typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  tol 
)
inlineinherited

Definition at line 84 of file reaction_rate_vector_test_base.h.

87  {
88  // Init using T as the "example"
89  PairScalars rate = T;
90  PairScalars deriv = T;
91 
92  reaction_rate.rate_and_derivative(T,rate,deriv);
93 
94  const PairScalars exact_rate = this->exact_rate(T);
95  const PairScalars exact_deriv = this->exact_deriv(T);
96 
97  for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
98  {
99  CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple],
100  exact_rate[0],
101  tol );
102 
103  CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple+1],
104  exact_rate[1],
105  tol );
106 
107  CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple],
108  exact_deriv[0],
109  tol );
110 
111  CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple+1],
112  exact_deriv[1],
113  tol );
114  }
115  }
_Cf VectorStateType VectorStateType(*) VectorStateType voi rate_and_derivative)(const KineticsConditions< StateType, VectorStateType > &T, StateType &rate, StateType &drate_dT) const
void AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::test_standard_deriv ( )
inlineinherited

Definition at line 78 of file arrhenius_rate_vector_test_base.h.

79  {
80  typedef typename Antioch::value_type<PairScalars>::type Scalar;
81 
82  const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 100;
83  PairScalars T = this->setup_T(*(this->_example));
84  this->test_deriv( *_rate, T, tol );
85  }
void test_deriv(const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &reaction_rate, const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &T, typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typetol)
Antioch::ArrheniusRate< typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > * _rate
void AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::test_standard_rate ( )
inlineinherited

Definition at line 69 of file arrhenius_rate_vector_test_base.h.

70  {
71  typedef typename Antioch::value_type<PairScalars>::type Scalar;
72 
73  const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 100;
74  PairScalars T = this->setup_T(*(this->_example));
75  this->test_rate( *_rate, T, tol );
76  }
Antioch::ArrheniusRate< typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > * _rate
void test_rate(const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &reaction_rate, const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &T, typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typetol)
void AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::test_standard_rate_and_deriv ( )
inlineinherited

Definition at line 87 of file arrhenius_rate_vector_test_base.h.

88  {
89  typedef typename Antioch::value_type<PairScalars>::type Scalar;
90 
91  const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 100;
92  PairScalars T = this->setup_T(*(this->_example));
93  this->test_rate_and_deriv( *_rate, T, tol );
94  }
void test_rate_and_deriv(const Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > &reaction_rate, const Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > &T, typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::typetol)
Antioch::ArrheniusRate< typename Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > * _rate
Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::value ( Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type  T)
inlineprotectedinherited

Definition at line 59 of file arrhenius_rate_test_helper.h.

60  {
61  using std::exp;
62  return _Cf*exp(-_Ea/(_R*T));
63  }

Member Data Documentation

Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::_Cf
protectedinherited

Definition at line 44 of file arrhenius_rate_test_helper.h.

Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::_Ea
protectedinherited

Definition at line 44 of file arrhenius_rate_test_helper.h.

Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > * AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type > , Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::_example
protectedinherited

Definition at line 121 of file reaction_rate_vector_test_base.h.

Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type >::_R
protectedinherited

Definition at line 44 of file arrhenius_rate_test_helper.h.

Antioch::ArrheniusRate<typename Antioch::value_type<Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::type>* AntiochTesting::ArrheniusRateVectorTestBase< Eigen::Array< double, 2 *ANTIOCH_N_TUPLES, 1 > >::_rate
protectedinherited

Definition at line 98 of file arrhenius_rate_vector_test_base.h.


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

Generated on Thu Jul 7 2016 11:09:51 for antioch-0.4.0 by  doxygen 1.8.8