[implicit][localjac] specialize updateGlobalJacobian routine to box or cc

dumux/implicit/localjacobian.hh

@@ -226,13 +226,15 @@ protected:
     /*!
      * \brief Updates the current global Jacobian matrix with the
      *        partial derivatives of all equations in regard to the
-     *        primary variable 'pvIdx' at dof 'col' .
+     *        primary variable 'pvIdx' at dof 'col'. Specialization for cc methods.
      */
-    void updateGlobalJacobian_(JacobianMatrix& matrix,
-                               const int globalI,
-                               const int globalJ,
-                               const int pvIdx,
-                               const PrimaryVariables &partialDeriv)
+    template<class T = TypeTag>
+    typename std::enable_if<!GET_PROP_VALUE(T, ImplicitIsBox)>::type
+    updateGlobalJacobian_(JacobianMatrix& matrix,
+                          const int globalI,
+                          const int globalJ,
+                          const int pvIdx,
+                          const PrimaryVariables &partialDeriv)
     {
         for (int eqIdx = 0; eqIdx < numEq; eqIdx++)
         {
@@ -245,6 +247,30 @@ protected:
         }
     }
 
+    /*!
+     * \brief Updates the current global Jacobian matrix with the
+     *        partial derivatives of all equations in regard to the
+     *        primary variable 'pvIdx' at dof 'col'. Specialization for the box method.
+     */
+    template<class T = TypeTag>
+    typename std::enable_if<GET_PROP_VALUE(T, ImplicitIsBox)>::type
+    updateGlobalJacobian_(JacobianMatrix& matrix,
+                          const int globalI,
+                          const int globalJ,
+                          const int pvIdx,
+                          const PrimaryVariables &partialDeriv)
+    {
+        for (int eqIdx = 0; eqIdx < numEq; eqIdx++)
+        {
+            // A[i][col][eqIdx][pvIdx] is the rate of change of
+            // the residual of equation 'eqIdx' at dof 'i'
+            // depending on the primary variable 'pvIdx' at dof
+            // 'col'.
+            matrix[globalI][globalJ][eqIdx][pvIdx] += partialDeriv[eqIdx];
+            Valgrind::CheckDefined(matrix[globalI][globalJ][eqIdx][pvIdx]);
+        }
+    }
+
     // The problem we would like to solve
     Problem *problemPtr_;