Reset Jacobian matrix to zero every residualImpl call with wantJac

Author: jbreue16

src/libcadet/model/ColumnModel1D.cpp

@@ -923,10 +923,11 @@ int ColumnModel1D::residualImpl(double t, unsigned int secIdx, StateType const*
 {
 	if (wantRes)
 	{
-		double* const resPtr = reinterpret_cast<double* const>(res);
-		Eigen::Map<Eigen::VectorXd> resi(resPtr, numDofs());
+		Eigen::Map<Eigen::Vector<ResidualType, Dynamic>> resi(res, numDofs());
 		resi.setZero();
 	}
+	if (wantJac)
+		_globalJac.coeffs().setZero();
 
 	LinearBufferAllocator tlmAlloc = threadLocalMem.get();
 	Indexer idxr(_disc);

src/libcadet/model/ColumnModel2D.cpp

@@ -1081,9 +1081,7 @@ int ColumnModel2D::residualImpl(double t, unsigned int secIdx, StateType const*
 		resi.setZero();
 	}
 	if (wantJac)
-	{
-		std::fill_n(_globalJac.valuePtr(), _globalJac.nonZeros(), 0.0);
-	}
+		_globalJac.coeffs().setZero();
 
 #ifdef CADET_PARALLELIZE
 	tbb::parallel_for(std::size_t(0), static_cast<std::size_t>(_disc.axNPoints * _disc.radNPoints * _disc.nParType + 1), [&](std::size_t pblk)