AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest Class Reference

Inheritance diagram for AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest:

Public Member Functions
	CPPUNIT_TEST_SUITE (ArrheniusRateMetaPhysicLDoubleTest)
	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< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > &reaction_rate, const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &T, typename Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typetol)
void	test_deriv (const Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > &reaction_rate, const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &T, typename Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typetol)
void	test_rate_and_deriv (const Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > &reaction_rate, const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &T, typename Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typetol)

Protected Member Functions
virtual MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float >	exact_rate (MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float >T)
virtual MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float >	exact_deriv (MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float >T)
void	reset_params (Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeCf, Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeEa)
void	reset_params (Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeCf, Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeEa, Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeR)
Antioch::value_type < MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > >::type	value (Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeT)
Antioch::value_type < MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > >::type	deriv (Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::typeT)
MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float >	setup_T (const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &example)

Protected Attributes
Antioch::ArrheniusRate < typename Antioch::value_type < MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > >::type > *	_rate
Antioch::value_type < MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > >::type	_Cf
Antioch::value_type < MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > >::type	_Ea
Antioch::value_type < MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > >::type	_R
MetaPhysicL::NumberArray < 2 *ANTIOCH_N_TUPLES, float > *	_example

Detailed Description

Definition at line 77 of file arrhenius_rate_metaphysicl_test.C.

Member Function Documentation

virtual void AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::clear ( )

inlinevirtualinherited

Definition at line 64 of file arrhenius_rate_vector_test_base.h.

Referenced by AntiochTesting::ArrheniusRateMetaPhysicLTest< long double >::tearDown().

    {
      delete _rate;
    }

AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest::CPPUNIT_TEST

(

test_standard_rate

)

AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest::CPPUNIT_TEST

(

test_standard_deriv

)

AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest::CPPUNIT_TEST

(

test_standard_rate_and_deriv

)

AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest::CPPUNIT_TEST_SUITE

(

ArrheniusRateMetaPhysicLDoubleTest

)

AntiochTesting::ArrheniusRateMetaPhysicLDoubleTest::CPPUNIT_TEST_SUITE_END

(

)

Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::deriv ( Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  T )

inlineprotectedinherited

Definition at line 65 of file arrhenius_rate_test_helper.h.

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

virtual MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::exact_deriv ( MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float >  T )

inlineprotectedvirtualinherited

Implements AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >, MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >.

Definition at line 111 of file arrhenius_rate_vector_test_base.h.

    {
      PairScalars e_deriv = *(this->_example);
      for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
        {
          e_deriv[2*tuple]   =  this->deriv(T[2*tuple]);
          e_deriv[2*tuple+1] =  this->deriv(T[2*tuple+1]);
        }
      return e_deriv;
    }

virtual MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::exact_rate ( MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float >  T )

inlineprotectedvirtualinherited

Implements AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >, MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >.

Definition at line 100 of file arrhenius_rate_vector_test_base.h.

    {
      PairScalars e_rate = *(this->_example);
      for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
        {
          e_rate[2*tuple]   =  this->value(T[2*tuple]);
          e_rate[2*tuple+1] =  this->value(T[2*tuple+1]);
        }
      return e_rate;
    }

virtual void AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::init ( )

inlinevirtualinherited

Definition at line 52 of file arrhenius_rate_vector_test_base.h.

Referenced by AntiochTesting::ArrheniusRateMetaPhysicLTest< long double >::setUp().

    {
      typedef typename Antioch::value_type<PairScalars>::type Scalar;

      Scalar Cf = 1.4L;
      Scalar Ea = 5.0L;
      Scalar R = 1.0L; // Ea in K

      this->reset_params(Cf,Ea,R);
      _rate = new Antioch::ArrheniusRate<Scalar>(Cf,Ea,R);
    }

void AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::reset_params ( Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  Cf, Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  Ea  )

inlineprotectedinherited

Definition at line 46 of file arrhenius_rate_test_helper.h.

    {
      _Cf = Cf;
      _Ea = Ea;
    }

void AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::reset_params ( Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  Cf, Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  Ea, Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  R  )

inlineprotectedinherited

Definition at line 52 of file arrhenius_rate_test_helper.h.

    {
      _Cf = Cf;
      _Ea = Ea;
      _R = R;
    }

virtual void AntiochTesting::ArrheniusRateMetaPhysicLTest< double >::setUp ( )

inlinevirtualinherited

Definition at line 50 of file arrhenius_rate_metaphysicl_test.C.

    {
      this->init();
      this->_example = new MetaPhysicL::NumberArray<2*ANTIOCH_N_TUPLES, float>(0);
    }

MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > , MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::setup_T ( const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &  example )

inlineprotectedinherited

Definition at line 123 of file reaction_rate_vector_test_base.h.

    {
      // Construct from example to avoid resizing issues
      PairScalars T = example;
      for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
        {
          T[2*tuple]   = 1500.1;
          T[2*tuple+1] = 1600.1;
        }
      return T;
    }

virtual void AntiochTesting::ArrheniusRateMetaPhysicLTest< double >::tearDown ( )

inlinevirtualinherited

Definition at line 56 of file arrhenius_rate_metaphysicl_test.C.

    {
      this->clear();
      delete this->_example;
    }

void AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > , MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::test_deriv ( const Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > &  reaction_rate, const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &  T, typename Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  tol  )

inlineinherited

Definition at line 64 of file reaction_rate_vector_test_base.h.

    {
      const PairScalars deriv = reaction_rate.derivative(T);
      const PairScalars exact_deriv = this->exact_deriv(T);

      for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
        {
          CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple],
                                        exact_deriv[0],
                                        tol );

          CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple+1],
                                        exact_deriv[1],
                                        tol );

        }
    }

void AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > , MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::test_rate ( const Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > &  reaction_rate, const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &  T, typename Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  tol  )

inlineinherited

Definition at line 44 of file reaction_rate_vector_test_base.h.

    {
      const PairScalars rate = reaction_rate(T);
      const PairScalars exact_rate = this->exact_rate(T);

      for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
        {
          CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple],
                                        exact_rate[0],
                                        tol );

          CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple+1],
                                        exact_rate[1],
                                        tol );

        }
    }

void AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > , MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::test_rate_and_deriv ( const Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > &  reaction_rate, const MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > &  T, typename Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  tol  )

inlineinherited

Definition at line 84 of file reaction_rate_vector_test_base.h.

    {
      // Init using T as the "example"
      PairScalars rate = T;
      PairScalars deriv = T;

      reaction_rate.rate_and_derivative(T,rate,deriv);

      const PairScalars exact_rate = this->exact_rate(T);
      const PairScalars exact_deriv = this->exact_deriv(T);

      for (unsigned int tuple=0; tuple != ANTIOCH_N_TUPLES; ++tuple)
        {
          CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple],
                                        exact_rate[0],
                                        tol );

          CPPUNIT_ASSERT_DOUBLES_EQUAL( rate[2*tuple+1],
                                        exact_rate[1],
                                        tol );

          CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple],
                                        exact_deriv[0],
                                        tol );

          CPPUNIT_ASSERT_DOUBLES_EQUAL( deriv[2*tuple+1],
                                        exact_deriv[1],
                                        tol );
        }
    }

void AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::test_standard_deriv ( )

inlineinherited

Definition at line 78 of file arrhenius_rate_vector_test_base.h.

    {
      typedef typename Antioch::value_type<PairScalars>::type Scalar;

      const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 100;
      PairScalars T = this->setup_T(*(this->_example));
      this->test_deriv( *_rate, T, tol );
    }

void AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::test_standard_rate ( )

inlineinherited

Definition at line 69 of file arrhenius_rate_vector_test_base.h.

    {
      typedef typename Antioch::value_type<PairScalars>::type Scalar;

      const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 100;
      PairScalars T = this->setup_T(*(this->_example));
      this->test_rate( *_rate, T, tol );
    }

void AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::test_standard_rate_and_deriv ( )

inlineinherited

Definition at line 87 of file arrhenius_rate_vector_test_base.h.

    {
      typedef typename Antioch::value_type<PairScalars>::type Scalar;

      const Scalar tol = std::numeric_limits<Scalar>::epsilon() * 100;
      PairScalars T = this->setup_T(*(this->_example));
      this->test_rate_and_deriv( *_rate, T, tol );
    }

Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::value ( Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type  T )

inlineprotectedinherited

Definition at line 59 of file arrhenius_rate_test_helper.h.

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

Member Data Documentation

Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::_Cf

protectedinherited

Definition at line 44 of file arrhenius_rate_test_helper.h.

Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::_Ea

protectedinherited

Definition at line 44 of file arrhenius_rate_test_helper.h.

MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > * AntiochTesting::ReactionRateVectorTestBase< Antioch::ArrheniusRate< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type > , MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::_example

protectedinherited

Definition at line 121 of file reaction_rate_vector_test_base.h.

Referenced by AntiochTesting::ArrheniusRateMetaPhysicLTest< long double >::setUp(), and AntiochTesting::ArrheniusRateMetaPhysicLTest< long double >::tearDown().

Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type AntiochTesting::ArrheniusRateTestHelper< Antioch::value_type< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type >::_R

protectedinherited

Definition at line 44 of file arrhenius_rate_test_helper.h.

Antioch::ArrheniusRate<typename Antioch::value_type<MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::type>* AntiochTesting::ArrheniusRateVectorTestBase< MetaPhysicL::NumberArray< 2 *ANTIOCH_N_TUPLES, float > >::_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:

• test/metaphysicl_unit/arrhenius_rate_metaphysicl_test.C