Merge branch 'feature/get-source-deriv-from-problem' into 'master'

[assembly][analytic] Get source derivative per default from problem

See merge request !757

dumux/common/fvproblem.hh

@@ -63,6 +63,7 @@ class FVProblem
     using PointSourceHelper = typename GET_PROP_TYPE(TypeTag, PointSourceHelper);
     using SolutionVector = typename GET_PROP_TYPE(TypeTag, SolutionVector);
     using ElementSolutionVector = typename GET_PROP_TYPE(TypeTag, ElementSolutionVector);
+    using VolumeVariables = typename GET_PROP_TYPE(TypeTag, VolumeVariables);
     using ElementVolumeVariables = typename GET_PROP_TYPE(TypeTag, ElementVolumeVariables);
 
     enum {
@@ -382,6 +383,17 @@ public:
     void pointSourceAtPos(PointSource& pointSource,
                           const GlobalPosition &globalPos) const {}
 
+    /*!
+     * \brief Add source term derivative to the Jacobian
+     * \note Only needed in case of analytic differentiation and solution dependent sources
+     */
+    template<class MatrixBlock>
+    void addSourceDerivatives(MatrixBlock& block,
+                              const Element& element,
+                              const FVElementGeometry& fvGeometry,
+                              const VolumeVariables& volVars,
+                              const SubControlVolume& scv) const {}
+
     /*!
      * \brief Adds contribution of point sources for a specific sub control volume
      *        to the values.

dumux/porousmediumflow/1p/incompressiblelocalresidual.hh

@@ -74,7 +74,10 @@ public:
                               const Element& element,
                               const FVElementGeometry& fvGeometry,
                               const VolumeVariables& curVolVars,
-                              const SubControlVolume& scv) const {}
+                              const SubControlVolume& scv) const
+    {
+        problem.addSourceDerivatives(partialDerivatives, element, fvGeometry, curVolVars, scv);
+    }
 
     //! flux derivatives for the cell-centered tpfa scheme
     template<class PartialDerivativeMatrices, class T = TypeTag>