Fix reaction config scope in particles

Author: jbreue16

src/libcadet/model/particle/GeneralRateParticle.cpp

@@ -162,6 +162,7 @@ namespace model
 				paramProvider.popScope();
 			}
 
+			MultiplexedScopeSelector scopeGuard(paramProvider, "reaction", _dynReaction->usesParamProviderInDiscretizationConfig());
 			reactionConfSuccess = _dynReaction->configureModelDiscretization(paramProvider, _nComp, _nBound.get(), _parDiffOp->offsetBoundComp()) && reactionConfSuccess;
 		}
 

src/libcadet/model/particle/HomogeneousParticle.cpp

@@ -176,6 +176,7 @@ namespace model
 				paramProvider.popScope();
 			}
 
+			MultiplexedScopeSelector scopeGuard(paramProvider, "reaction", _dynReaction->usesParamProviderInDiscretizationConfig());
 			reactionConfSuccess = _dynReaction->configureModelDiscretization(paramProvider, _nComp, _nBound.get(), _boundOffset) && reactionConfSuccess;
 		}
 
@@ -404,24 +405,22 @@ namespace model
 
 	void HomogeneousParticle::setParJacPattern(std::vector<Eigen::Triplet<double>>& tripletList, const unsigned int offsetPar, const unsigned int offsetBulk, unsigned int colNode, unsigned int secIdx) const
 	{
-		// binding Jacobian pattern
-		// add dense nComp * nBound blocks, since all solid and liquid entries can be coupled through binding.
+		// binding and reaction Jacobian pattern
+		// add dense nComp * nBound blocks, since all solid and liquid entries can be coupled through reaction and or binding.
 		for (unsigned int parState = 0; parState < _nComp + _strideBound; parState++) {
 			for (unsigned int toParState = 0; toParState < _nComp + _strideBound; toParState++) {
 				tripletList.push_back(Eigen::Triplet<double>(offsetPar + parState, offsetPar + toParState, 0.0));
 			}
 		}
 
 		// flux Jacobian pattern
-
 		for (unsigned int comp = 0; comp < _nComp; comp++)
 		{
 			tripletList.push_back(Eigen::Triplet<double>(offsetPar + comp, offsetBulk + comp, 0.0));
 			tripletList.push_back(Eigen::Triplet<double>(offsetBulk + comp, offsetPar + comp, 0.0));
 		}
 	}
 
-
 	unsigned int HomogeneousParticle::jacobianNNZperParticle() const
 	{
 		return (_nComp + _strideBound) * (_nComp + _strideBound) + _nComp * 4; // reaction, binding patter + film diffusion pattern for one particle