Commit 2190ab5d authored by Samuel Burbulla's avatar Samuel Burbulla
Browse files

[bugfix] Update some function calls.

parent c13bef75
......@@ -365,6 +365,12 @@ public:
const SolutionVector& curSol() const
{ return ParentType::curSol(); }
/*!
* \brief Returns the coupling context of the coupled problem
*/
const PoroMechanicsCouplingContext& poroMechanicsCouplingContext() const
{ return poroMechCouplingContext_; }
private:
/*!
* \brief Initializes the pm flow domain coupling map. Since the elements
......
......@@ -180,9 +180,9 @@ public:
const SolutionVector& curSol)
{
curSol_ = curSol;
BulkFacetFlowManager::init(matrixFlowProblem, facetFlowProblem, matrixGridVars, facetFlowGridVars, bulkFacetFlowMapper, curSol);
BulkFacetMechManager::init(mechProblem, facetFlowProblem, mechGridVars, facetFlowGridVars, bulkFacetMechMapper, curSol);
PoroMechManager::init(matrixFlowProblem, mechProblem, matrixGridVars, mechGridVars, curSol);
BulkFacetFlowManager::init(matrixFlowProblem, facetFlowProblem, bulkFacetFlowMapper, curSol);
BulkFacetMechManager::init(mechProblem, facetFlowProblem, bulkFacetMechMapper, curSol);
PoroMechManager::init(matrixFlowProblem, mechProblem, curSol);
// set up the map of bulk dofs coinciding with the lower-dimensional elements
// in order to be able to compute the deformation-dependent aperture
......@@ -231,6 +231,12 @@ public:
}
}
/*!
* \brief Returns the coupling context of the coupled problem
*/
const auto& poroMechanicsCouplingContext() const
{ return PoroMechManager::poroMechanicsCouplingContext(); }
/*!
* \brief Computes the aperture of a sub-control volume within
* a given lower-dimensional element as a function of the
......
......@@ -159,19 +159,22 @@ int main(int argc, char** argv) try
couplingMapperFlow->update(fvGridGeometry[bulkFlowId], fvGridGeometry[facetFlowId], gridManager.getEmbeddings());
couplingMapperMech->update(fvGridGeometry[bulkMechId], fvGridGeometry[facetFlowId], gridManager.getEmbeddings());
// the grid variables
using GridVariables = MultiDomainFVGridVariables<Traits>;
GridVariables gridVars(fvGridGeometry.getTuple(), problem.getTuple());
gridVars.init(x);
// initialize the coupling manager
couplingManager->init(problem.get(bulkFlowId),
problem.get(facetFlowId),
problem.get(bulkMechId),
gridVars.get(bulkFlowId),
gridVars.get(facetFlowId),
gridVars.get(bulkMechId),
couplingMapperFlow,
couplingMapperMech,
x);
// the grid variables
using GridVariables = MultiDomainFVGridVariables<Traits>;
GridVariables gridVars(fvGridGeometry.getTuple(), problem.getTuple());
gridVars.init(x);
// intialize the vtk output module
MultiDomainVtkOutputModule<Traits> vtkWriter;
......
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