@@ -29,11 +29,21 @@ REAL_TYPE tolerance( REAL_TYPE const a, REAL_TYPE const b )
2929 return std::numeric_limits< double >::epsilon () * std::max ( fabs ( a ), fabs ( b ) ) * 10 ;
3030}
3131
32+ /* *
33+ * POD struct for transferring data between host and device for computeResidualAndJacobianTest.
34+ * @tparam numReactions Number of reactions.
35+ * @tparam numSpecies Number of species.
36+ */
3237template < int numReactions, int numSpecies >
3338struct ComputeResidualAndJacobianTestData
3439{
35- CArrayWrapper< double , numReactions, numReactions > jacobian;
40+ // / The reaction residuals
3641 double residual[numReactions] = { 0.0 };
42+
43+ // / The residual derivatives wrt reactions
44+ CArrayWrapper< double , numReactions, numReactions > jacobian;
45+
46+ // / the species concentrations
3747 double speciesConcentration[numSpecies];
3848};
3949
@@ -95,10 +105,17 @@ void computeResidualAndJacobianTest( PARAMS_DATA const & params,
95105
96106// ******************************************************************************
97107
108+ /* *
109+ * POD struct for transferring data between host and device for testEnforceEquilibrium.
110+ * @tparam numSpecies Number of species.
111+ */
98112template < int numSpecies >
99113struct TestEnforceEquilibriumData
100114{
115+ // / The initial species concentrations
101116 double speciesConcentration0[numSpecies];
117+
118+ // / The final species concentrations
102119 double speciesConcentration[numSpecies];
103120};
104121
0 commit comments