Commit 2871d6f7 authored by GILLES Sebastien's avatar GILLES Sebastien
Browse files

#9 FSI: add some comments (for data attributes introduced in Poromechanics model in the same time).

parent a04f0617
......@@ -123,13 +123,13 @@ namespace HappyHeart
const auto& inlet_border_felt_space = god_of_dof.GetFEltSpace(EnumUnderlyingType(FEltSpaceIndex::inlet_border));
inlet_pressure_source_ =
inlet_pressure_source_operator_ =
std::make_unique<GVO::TransientSource>(inlet_border_felt_space,
fluid_velocity,
mesh_dimension,
default_quadrature_rule_set,
DoComputeProcessorWiseLocal2Global::yes,
*inlet_pressure_); // \todo #608 Make this one a parameter rather than an hardcoded one!
GetInletPressure());
}
}
......@@ -213,7 +213,7 @@ namespace HappyHeart
const auto time = GetTimeManager().GetTime();
GlobalVectorWithCoefficient vector(rhs, 1.);
GetInletPressureOperator().Assemble(std::make_tuple(std::ref(vector)), time);
GetInletPressureSourceOperator().Assemble(std::make_tuple(std::ref(vector)), time);
}
break;
}
......
......@@ -157,6 +157,12 @@ namespace HappyHeart
///@}
private:
/// \name Accessors
///@{
const ExplicitStepVariationalFormulation& GetExplicitStepVariationalFormulation() const noexcept;
......@@ -184,11 +190,16 @@ namespace HappyHeart
GlobalVector& GetNonCstResidual() noexcept;
GlobalVector& GetNonCstMixedResidual() noexcept;
//! Constant accessor to the inlet pressure parameter.
const Parameter<ParameterNS::Type::vector>& GetInletPressure() const noexcept;
// \todo #608 Find a less crappy name!
const ConformInterpolatorNS::SubsetOrSuperset& GetTowardVelocityPressureInterpolator() const noexcept;
///@}
public:
// \todo #608 Find a less crappy name! It is actually the revert operation of GetTowardVelocityPressureInterpolator().
......@@ -199,7 +210,8 @@ namespace HappyHeart
private:
const GlobalVariationalOperatorNS::TransientSource& GetInletPressureOperator() const noexcept;
//! Constant accessor to the inlet pressure (acts as a source).
const GlobalVariationalOperatorNS::TransientSource& GetInletPressureSourceOperator() const noexcept;
//! Set the residual matrix and rhs.
void SetResidualLinearAlgebra();
......@@ -215,7 +227,8 @@ namespace HappyHeart
GlobalVariationalOperatorNS::ScalarDivVectorial::const_unique_ptr scalar_div_vectorial_ = nullptr;
GlobalVariationalOperatorNS::TransientSource::const_unique_ptr inlet_pressure_source_ = nullptr;
//! Inlet pressure (acts as a source).
GlobalVariationalOperatorNS::TransientSource::const_unique_ptr inlet_pressure_source_operator_ = nullptr;
///@}
......@@ -276,7 +289,8 @@ namespace HappyHeart
const ExplicitStepVariationalFormulation& explicit_variational_formulation_;
Parameter<ParameterNS::Type::vector>::unique_ptr inlet_pressure_;
//! Inlet pressure parameter.
Parameter<ParameterNS::Type::vector>::const_unique_ptr inlet_pressure_ = nullptr;
};
......
......@@ -106,12 +106,12 @@ namespace HappyHeart
inline const GlobalVariationalOperatorNS::TransientSource& ImplicitStepFluidVariationalFormulation
::GetInletPressureOperator() const noexcept
::GetInletPressureSourceOperator() const noexcept
{
assert(!(!inlet_pressure_source_));
return *inlet_pressure_source_;
assert(!(!inlet_pressure_source_operator_));
return *inlet_pressure_source_operator_;
}
inline const GlobalMatrix& ImplicitStepFluidVariationalFormulation::GetResidualMatrix() const noexcept
{
......@@ -158,8 +158,16 @@ namespace HappyHeart
return *mixed_residual_;
}
inline const Parameter<ParameterNS::Type::vector>& ImplicitStepFluidVariationalFormulation
::GetInletPressure() const noexcept
{
assert(!(!inlet_pressure_));
return *inlet_pressure_;
}
} // namespace FSI_EI_NS
......
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