Commit 9eb9fa35 authored by Timo Koch's avatar Timo Koch
Browse files

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

parent 9b307303
......@@ -62,6 +62,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 {
......@@ -378,6 +379,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.
......
......@@ -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>
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment