improve code comments regarding reaction jacobian

Author: jbreue16

src/libcadet/model/ColumnModel1D.hpp

@@ -501,7 +501,7 @@ class ColumnModel1D : public UnitOperationBase
 			}
 		}
 
-		// particle jacobian (including film diffusion, isotherm and time derivative)
+		// particle jacobian (including film diffusion, isotherm, reaction and time derivative)
 		for (int colNode = 0; colNode < _disc.nPoints; colNode++)
 		{
 			for (int type = 0; type < _disc.nParType; type++)
@@ -519,7 +519,7 @@ class ColumnModel1D : public UnitOperationBase
 		// inlet and bulk jacobian
 		_convDispOp.calcTransportJacobian(_globalJac, _jacInlet, idxr.offsetC());
 
-		// particle jacobian (without isotherm, which is handled in residualKernel)
+		// particle transport jacobian
 		for (int colNode = 0; colNode < _disc.nPoints; colNode++)
 		{
 			for (int parType = 0; parType < _disc.nParType; parType++)

src/libcadet/model/parts/ParticleDiffusionOperatorDG.cpp

@@ -912,11 +912,12 @@ namespace parts
 	// ========================================						DG particle Jacobian							=============================================  //
 	// ==========================================================================================================================================================  //
 
-	void ParticleDiffusionOperatorDG::parBindingPattern(std::vector<Eigen::Triplet<double>>& tripletList, const int offset, const unsigned int colNode)
+	void ParticleDiffusionOperatorDG::parBindingAndReactionPattern(std::vector<Eigen::Triplet<double>>& tripletList, const int offset, const unsigned int colNode)
 	{
-		// every bound state might depend on every bound and liquid state
 		for (int parNode = 0; parNode < _nParPoints; parNode++)
 		{
+			// binding: every bound state might depend on every bound and liquid state
+			// reaction: every liquid and bound state might depend on every liquid and bound state
 			for (int state = 0; state < strideParPoint(); state++)
 			{
 				for (int stateDep = 0; stateDep < strideParPoint(); stateDep++) {
@@ -928,14 +929,14 @@ namespace parts
 			}
 		}
 	}
+
 	/**
 	 * @brief calculates the particle dispersion jacobian Pattern, including entries for the dependence of particle entries on bulk entries through film diffusion boundary condition
 	 * @detail Does NOT add film diffusion entries for the dependence of bulk conc. on particle conc.
 	*/
 	void ParticleDiffusionOperatorDG::setParticleJacobianPattern(std::vector<ParticleDiffusionOperatorDG::T>& tripletList, unsigned int offsetPar, unsigned int offsetBulk, unsigned int colNode, unsigned int secIdx)
 	{
-
-		parBindingPattern(tripletList, offsetPar, colNode);
+		parBindingAndReactionPattern(tripletList, offsetPar, colNode);
 
 		// Ordering of particle surface diffusion:
 		// bnd0comp0, bnd1comp0, bnd0comp1, bnd1comp1, bnd0comp2, bnd1comp2

src/libcadet/model/parts/ParticleDiffusionOperatorDG.hpp

@@ -143,7 +143,7 @@ namespace parts
 
 	protected:
 
-		void parBindingPattern(std::vector<Eigen::Triplet<double>>& tripletList, const int offset, const unsigned int colNode);
+		void parBindingAndReactionPattern(std::vector<Eigen::Triplet<double>>& tripletList, const int offset, const unsigned int colNode);
 
 		template <typename StateType, typename ResidualType, typename ParamType, bool wantJac, bool wantRes>
 		int residualImpl(double t, unsigned int secIdx, StateType const* yPar, StateType const* yBulk, double const* yDotPar, ResidualType* resPar, linalg::BandedEigenSparseRowIterator& jacBase);