Skip to content
Snippets Groups Projects
Commit 09a9897c authored by wenqing's avatar wenqing
Browse files

Changed the name of coupling_term to coupled_solutions

parent 628a8250
No related branches found
No related tags found
No related merge requests found
Showing
with 51 additions and 52 deletions
...@@ -40,7 +40,7 @@ public: ...@@ -40,7 +40,7 @@ public:
std::vector<double>& /*local_K_data*/, std::vector<double>& /*local_K_data*/,
std::vector<double>& /*local_b_data*/, std::vector<double>& /*local_b_data*/,
std::vector<double>& /*local_Jac_data*/, std::vector<double>& /*local_Jac_data*/,
LocalCoupledSolutions const& /*coupling_term*/) LocalCoupledSolutions const& /*coupled_solutions*/)
{ {
} }
......
...@@ -63,7 +63,7 @@ void ComponentTransportProcess::assembleConcreteProcess( ...@@ -63,7 +63,7 @@ void ComponentTransportProcess::assembleConcreteProcess(
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void ComponentTransportProcess::assembleWithJacobianConcreteProcess( void ComponentTransportProcess::assembleWithJacobianConcreteProcess(
...@@ -77,7 +77,7 @@ void ComponentTransportProcess::assembleWithJacobianConcreteProcess( ...@@ -77,7 +77,7 @@ void ComponentTransportProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
} // namespace ComponentTransport } // namespace ComponentTransport
......
...@@ -68,7 +68,7 @@ void GroundwaterFlowProcess::assembleConcreteProcess(const double t, ...@@ -68,7 +68,7 @@ void GroundwaterFlowProcess::assembleConcreteProcess(const double t,
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void GroundwaterFlowProcess::assembleWithJacobianConcreteProcess( void GroundwaterFlowProcess::assembleWithJacobianConcreteProcess(
...@@ -82,7 +82,7 @@ void GroundwaterFlowProcess::assembleWithJacobianConcreteProcess( ...@@ -82,7 +82,7 @@ void GroundwaterFlowProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
} // namespace GroundwaterFlow } // namespace GroundwaterFlow
......
...@@ -62,7 +62,7 @@ void HTProcess::assembleConcreteProcess(const double t, ...@@ -62,7 +62,7 @@ void HTProcess::assembleConcreteProcess(const double t,
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void HTProcess::assembleWithJacobianConcreteProcess( void HTProcess::assembleWithJacobianConcreteProcess(
...@@ -76,7 +76,7 @@ void HTProcess::assembleWithJacobianConcreteProcess( ...@@ -76,7 +76,7 @@ void HTProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
} // namespace HT } // namespace HT
......
...@@ -95,7 +95,7 @@ void HeatConductionProcess::assembleConcreteProcess(const double t, ...@@ -95,7 +95,7 @@ void HeatConductionProcess::assembleConcreteProcess(const double t,
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void HeatConductionProcess::assembleWithJacobianConcreteProcess( void HeatConductionProcess::assembleWithJacobianConcreteProcess(
...@@ -109,7 +109,7 @@ void HeatConductionProcess::assembleWithJacobianConcreteProcess( ...@@ -109,7 +109,7 @@ void HeatConductionProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
void HeatConductionProcess::computeSecondaryVariableConcrete( void HeatConductionProcess::computeSecondaryVariableConcrete(
...@@ -119,7 +119,7 @@ void HeatConductionProcess::computeSecondaryVariableConcrete( ...@@ -119,7 +119,7 @@ void HeatConductionProcess::computeSecondaryVariableConcrete(
GlobalExecutor::executeMemberOnDereferenced( GlobalExecutor::executeMemberOnDereferenced(
&HeatConductionLocalAssemblerInterface::computeSecondaryVariable, &HeatConductionLocalAssemblerInterface::computeSecondaryVariable,
_local_assemblers, *_local_to_global_index_map, t, x, _local_assemblers, *_local_to_global_index_map, t, x,
_coupling_term); _coupled_solutions);
} }
} // namespace HeatConduction } // namespace HeatConduction
......
...@@ -180,7 +180,7 @@ void HydroMechanicsProcess<DisplacementDim>::assembleConcreteProcess( ...@@ -180,7 +180,7 @@ void HydroMechanicsProcess<DisplacementDim>::assembleConcreteProcess(
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
template <int DisplacementDim> template <int DisplacementDim>
...@@ -196,7 +196,7 @@ void HydroMechanicsProcess<DisplacementDim>:: ...@@ -196,7 +196,7 @@ void HydroMechanicsProcess<DisplacementDim>::
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
template <int DisplacementDim> template <int DisplacementDim>
void HydroMechanicsProcess<DisplacementDim>::preTimestepConcreteProcess( void HydroMechanicsProcess<DisplacementDim>::preTimestepConcreteProcess(
......
...@@ -408,7 +408,7 @@ void HydroMechanicsProcess<GlobalDim>::computeSecondaryVariableConcrete( ...@@ -408,7 +408,7 @@ void HydroMechanicsProcess<GlobalDim>::computeSecondaryVariableConcrete(
GlobalExecutor::executeMemberOnDereferenced( GlobalExecutor::executeMemberOnDereferenced(
&HydroMechanicsLocalAssemblerInterface::computeSecondaryVariable, &HydroMechanicsLocalAssemblerInterface::computeSecondaryVariable,
_local_assemblers, *_local_to_global_index_map, t, x, _local_assemblers, *_local_to_global_index_map, t, x,
_coupling_term); _coupled_solutions);
// Copy displacement jumps in a solution vector to mesh property // Copy displacement jumps in a solution vector to mesh property
// Remark: the copy is required because mesh properties for primary // Remark: the copy is required because mesh properties for primary
......
...@@ -75,7 +75,7 @@ private: ...@@ -75,7 +75,7 @@ private:
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _global_assembler, &VectorMatrixAssembler::assemble,
_local_assemblers, *_local_to_global_index_map, t, x, M, K, b, _local_assemblers, *_local_to_global_index_map, t, x, M, K, b,
_coupling_term); _coupled_solutions);
} }
void assembleWithJacobianConcreteProcess( void assembleWithJacobianConcreteProcess(
...@@ -89,7 +89,7 @@ private: ...@@ -89,7 +89,7 @@ private:
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, _local_assemblers, *_local_to_global_index_map, t, x, xdot,
dxdot_dx, dx_dx, M, K, b, Jac, _coupling_term); dxdot_dx, dx_dx, M, K, b, Jac, _coupled_solutions);
} }
void preTimestepConcreteProcess(GlobalVector const& x, double const t, void preTimestepConcreteProcess(GlobalVector const& x, double const t,
......
...@@ -73,7 +73,7 @@ private: ...@@ -73,7 +73,7 @@ private:
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _global_assembler, &VectorMatrixAssembler::assemble,
_local_assemblers, *_local_to_global_index_map, t, x, M, K, b, _local_assemblers, *_local_to_global_index_map, t, x, M, K, b,
_coupling_term); _coupled_solutions);
} }
void assembleWithJacobianConcreteProcess( void assembleWithJacobianConcreteProcess(
...@@ -87,7 +87,7 @@ private: ...@@ -87,7 +87,7 @@ private:
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, _local_assemblers, *_local_to_global_index_map, t, x, xdot,
dxdot_dx, dx_dx, M, K, b, Jac, _coupling_term); dxdot_dx, dx_dx, M, K, b, Jac, _coupled_solutions);
} }
void preTimestepConcreteProcess(GlobalVector const& x, double const t, void preTimestepConcreteProcess(GlobalVector const& x, double const t,
......
...@@ -268,7 +268,7 @@ LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>:: ...@@ -268,7 +268,7 @@ LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::
// the assert must be changed to perm.rows() == _element->getDimension() // the assert must be changed to perm.rows() == _element->getDimension()
assert(permeability.rows() == GlobalDim || permeability.rows() == 1); assert(permeability.rows() == GlobalDim || permeability.rows() == 1);
if (!_coupling_term) if (!_coupled_solutions)
{ {
computeDarcyVelocity(permeability, local_x, veloctiy_cache_vectors); computeDarcyVelocity(permeability, local_x, veloctiy_cache_vectors);
} }
...@@ -276,7 +276,7 @@ LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>:: ...@@ -276,7 +276,7 @@ LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::
{ {
auto const local_coupled_xs = auto const local_coupled_xs =
getCurrentLocalSolutionsOfCoupledProcesses( getCurrentLocalSolutionsOfCoupledProcesses(
_coupling_term->coupled_xs, indices); _coupled_solutions->coupled_xs, indices);
if (local_coupled_xs.empty()) if (local_coupled_xs.empty())
computeDarcyVelocity(permeability, local_x, veloctiy_cache_vectors); computeDarcyVelocity(permeability, local_x, veloctiy_cache_vectors);
else else
......
...@@ -43,15 +43,16 @@ class LiquidFlowLocalAssemblerInterface ...@@ -43,15 +43,16 @@ class LiquidFlowLocalAssemblerInterface
{ {
public: public:
LiquidFlowLocalAssemblerInterface( LiquidFlowLocalAssemblerInterface(
CoupledSolutionsForStaggeredScheme* const coupling_term) CoupledSolutionsForStaggeredScheme* const coupled_solutions)
: _coupling_term(coupling_term) : _coupled_solutions(coupled_solutions)
{ {
} }
void setCoupledSolutionsForStaggeredScheme(std::size_t const /*mesh_item_id*/, void setCoupledSolutionsForStaggeredScheme(
CoupledSolutionsForStaggeredScheme* const coupling_term) std::size_t const /*mesh_item_id*/,
CoupledSolutionsForStaggeredScheme* const coupled_solutions)
{ {
_coupling_term = coupling_term; _coupled_solutions = coupled_solutions;
} }
virtual std::vector<double> const& getIntPtDarcyVelocity( virtual std::vector<double> const& getIntPtDarcyVelocity(
...@@ -61,11 +62,8 @@ public: ...@@ -61,11 +62,8 @@ public:
std::vector<double>& /*cache*/) const = 0; std::vector<double>& /*cache*/) const = 0;
protected: protected:
// TODO: remove _coupling_term or move integration point data from local
// assembler class to a new class to make local assembler unique for each
//process.
/// Pointer that is set from a Process class. /// Pointer that is set from a Process class.
CoupledSolutionsForStaggeredScheme* _coupling_term; CoupledSolutionsForStaggeredScheme* _coupled_solutions;
}; };
template <typename ShapeFunction, typename IntegrationMethod, template <typename ShapeFunction, typename IntegrationMethod,
...@@ -95,8 +93,8 @@ public: ...@@ -95,8 +93,8 @@ public:
double const gravitational_acceleration, double const gravitational_acceleration,
double const reference_temperature, double const reference_temperature,
LiquidFlowMaterialProperties const& material_propertries, LiquidFlowMaterialProperties const& material_propertries,
CoupledSolutionsForStaggeredScheme* coupling_term) CoupledSolutionsForStaggeredScheme* coupled_solutions)
: LiquidFlowLocalAssemblerInterface(coupling_term), : LiquidFlowLocalAssemblerInterface(coupled_solutions),
_element(element), _element(element),
_integration_method(integration_order), _integration_method(integration_order),
_shape_matrices(initShapeMatrices<ShapeFunction, ShapeMatricesType, _shape_matrices(initShapeMatrices<ShapeFunction, ShapeMatricesType,
......
...@@ -66,7 +66,7 @@ void LiquidFlowProcess::initializeConcreteProcess( ...@@ -66,7 +66,7 @@ void LiquidFlowProcess::initializeConcreteProcess(
pv.getShapeFunctionOrder(), _local_assemblers, pv.getShapeFunctionOrder(), _local_assemblers,
mesh.isAxiallySymmetric(), integration_order, _gravitational_axis_id, mesh.isAxiallySymmetric(), integration_order, _gravitational_axis_id,
_gravitational_acceleration, _reference_temperature, _gravitational_acceleration, _reference_temperature,
*_material_properties, _coupling_term); *_material_properties, _coupled_solutions);
_secondary_variables.addSecondaryVariable( _secondary_variables.addSecondaryVariable(
"darcy_velocity", "darcy_velocity",
...@@ -85,7 +85,7 @@ void LiquidFlowProcess::assembleConcreteProcess(const double t, ...@@ -85,7 +85,7 @@ void LiquidFlowProcess::assembleConcreteProcess(const double t,
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void LiquidFlowProcess::assembleWithJacobianConcreteProcess( void LiquidFlowProcess::assembleWithJacobianConcreteProcess(
...@@ -99,7 +99,7 @@ void LiquidFlowProcess::assembleWithJacobianConcreteProcess( ...@@ -99,7 +99,7 @@ void LiquidFlowProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
void LiquidFlowProcess::computeSecondaryVariableConcrete(const double t, void LiquidFlowProcess::computeSecondaryVariableConcrete(const double t,
...@@ -108,15 +108,16 @@ void LiquidFlowProcess::computeSecondaryVariableConcrete(const double t, ...@@ -108,15 +108,16 @@ void LiquidFlowProcess::computeSecondaryVariableConcrete(const double t,
DBUG("Compute the velocity for LiquidFlowProcess."); DBUG("Compute the velocity for LiquidFlowProcess.");
GlobalExecutor::executeMemberOnDereferenced( GlobalExecutor::executeMemberOnDereferenced(
&LiquidFlowLocalAssemblerInterface::computeSecondaryVariable, &LiquidFlowLocalAssemblerInterface::computeSecondaryVariable,
_local_assemblers, *_local_to_global_index_map, t, x, _coupling_term); _local_assemblers, *_local_to_global_index_map, t, x, _coupled_solutions);
} }
void LiquidFlowProcess::setCoupledSolutionsForStaggeredSchemeToLocalAssemblers() void LiquidFlowProcess::setCoupledSolutionsForStaggeredSchemeToLocalAssemblers()
{ {
DBUG("Compute the velocity for LiquidFlowProcess."); DBUG("Compute the velocity for LiquidFlowProcess.");
GlobalExecutor::executeMemberOnDereferenced( GlobalExecutor::executeMemberOnDereferenced(
&LiquidFlowLocalAssemblerInterface::setCoupledSolutionsForStaggeredScheme, &LiquidFlowLocalAssemblerInterface::
_local_assemblers, _coupling_term); setCoupledSolutionsForStaggeredScheme,
_local_assemblers, _coupled_solutions);
} }
} // end of namespace } // end of namespace
......
...@@ -20,7 +20,7 @@ void LocalAssemblerInterface::assembleWithCoupledTerm( ...@@ -20,7 +20,7 @@ void LocalAssemblerInterface::assembleWithCoupledTerm(
std::vector<double>& /*local_M_data*/, std::vector<double>& /*local_M_data*/,
std::vector<double>& /*local_K_data*/, std::vector<double>& /*local_K_data*/,
std::vector<double>& /*local_b_data*/, std::vector<double>& /*local_b_data*/,
LocalCoupledSolutions const& /*coupling_term*/) LocalCoupledSolutions const& /*coupled_solutions*/)
{ {
OGS_FATAL( OGS_FATAL(
"The assembleWithCoupledTerm() function is not implemented in the " "The assembleWithCoupledTerm() function is not implemented in the "
...@@ -47,7 +47,7 @@ void LocalAssemblerInterface::assembleWithJacobianAndCoupling( ...@@ -47,7 +47,7 @@ void LocalAssemblerInterface::assembleWithJacobianAndCoupling(
std::vector<double>& /*local_K_data*/, std::vector<double>& /*local_K_data*/,
std::vector<double>& /*local_b_data*/, std::vector<double>& /*local_b_data*/,
std::vector<double>& /*local_Jac_data*/, std::vector<double>& /*local_Jac_data*/,
LocalCoupledSolutions const& /*coupling_term*/) LocalCoupledSolutions const& /*coupled_solutions*/)
{ {
OGS_FATAL( OGS_FATAL(
"The assembleWithJacobianAndCoupling() function is not implemented in" "The assembleWithJacobianAndCoupling() function is not implemented in"
......
...@@ -46,7 +46,7 @@ public: ...@@ -46,7 +46,7 @@ public:
std::vector<double>& local_M_data, std::vector<double>& local_M_data,
std::vector<double>& local_K_data, std::vector<double>& local_K_data,
std::vector<double>& local_b_data, std::vector<double>& local_b_data,
LocalCoupledSolutions const& coupling_term); LocalCoupledSolutions const& coupled_solutions);
virtual void assembleWithJacobian(double const t, virtual void assembleWithJacobian(double const t,
std::vector<double> const& local_x, std::vector<double> const& local_x,
...@@ -63,7 +63,7 @@ public: ...@@ -63,7 +63,7 @@ public:
const double dx_dx, std::vector<double>& local_M_data, const double dx_dx, std::vector<double>& local_M_data,
std::vector<double>& local_K_data, std::vector<double>& local_b_data, std::vector<double>& local_K_data, std::vector<double>& local_b_data,
std::vector<double>& local_Jac_data, std::vector<double>& local_Jac_data,
LocalCoupledSolutions const& coupling_term); LocalCoupledSolutions const& coupled_solutions);
virtual void computeSecondaryVariable(std::size_t const mesh_item_id, virtual void computeSecondaryVariable(std::size_t const mesh_item_id,
NumLib::LocalToGlobalIndexMap const& dof_table, NumLib::LocalToGlobalIndexMap const& dof_table,
......
...@@ -142,7 +142,7 @@ void PhaseFieldProcess<DisplacementDim>::assembleConcreteProcess( ...@@ -142,7 +142,7 @@ void PhaseFieldProcess<DisplacementDim>::assembleConcreteProcess(
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
template <int DisplacementDim> template <int DisplacementDim>
...@@ -157,7 +157,7 @@ void PhaseFieldProcess<DisplacementDim>::assembleWithJacobianConcreteProcess( ...@@ -157,7 +157,7 @@ void PhaseFieldProcess<DisplacementDim>::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
template <int DisplacementDim> template <int DisplacementDim>
......
...@@ -31,7 +31,7 @@ Process::Process( ...@@ -31,7 +31,7 @@ Process::Process(
_secondary_variables(std::move(secondary_variables)), _secondary_variables(std::move(secondary_variables)),
_named_function_caller(std::move(named_function_caller)), _named_function_caller(std::move(named_function_caller)),
_global_assembler(std::move(jacobian_assembler)), _global_assembler(std::move(jacobian_assembler)),
_coupling_term(nullptr), _coupled_solutions(nullptr),
_integration_order(integration_order), _integration_order(integration_order),
_process_variables(std::move(process_variables)), _process_variables(std::move(process_variables)),
_boundary_conditions(parameters) _boundary_conditions(parameters)
......
...@@ -71,9 +71,9 @@ public: ...@@ -71,9 +71,9 @@ public:
MathLib::MatrixSpecifications getMatrixSpecifications() const final; MathLib::MatrixSpecifications getMatrixSpecifications() const final;
void setCoupledSolutionsForStaggeredScheme(CoupledSolutionsForStaggeredScheme* const coupling_term) void setCoupledSolutionsForStaggeredScheme(CoupledSolutionsForStaggeredScheme* const coupled_solutions)
{ {
_coupling_term = coupling_term; _coupled_solutions = coupled_solutions;
} }
virtual void setCoupledSolutionsForStaggeredSchemeToLocalAssemblers() {} virtual void setCoupledSolutionsForStaggeredSchemeToLocalAssemblers() {}
...@@ -205,7 +205,7 @@ protected: ...@@ -205,7 +205,7 @@ protected:
/// Pointer to CoupledSolutionsForStaggeredScheme, which contains the references to the /// Pointer to CoupledSolutionsForStaggeredScheme, which contains the references to the
/// coupled processes and the references to the solutions of the coupled /// coupled processes and the references to the solutions of the coupled
/// processes. /// processes.
CoupledSolutionsForStaggeredScheme* _coupling_term; CoupledSolutionsForStaggeredScheme* _coupled_solutions;
/// Order of the integration method for element-wise integration. /// Order of the integration method for element-wise integration.
/// The Gauss-Legendre integration method and available orders is /// The Gauss-Legendre integration method and available orders is
......
...@@ -70,7 +70,7 @@ void RichardsComponentTransportProcess::assembleConcreteProcess( ...@@ -70,7 +70,7 @@ void RichardsComponentTransportProcess::assembleConcreteProcess(
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void RichardsComponentTransportProcess::assembleWithJacobianConcreteProcess( void RichardsComponentTransportProcess::assembleWithJacobianConcreteProcess(
...@@ -84,7 +84,7 @@ void RichardsComponentTransportProcess::assembleWithJacobianConcreteProcess( ...@@ -84,7 +84,7 @@ void RichardsComponentTransportProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
} // namespace RichardsComponentTransport } // namespace RichardsComponentTransport
......
...@@ -73,7 +73,7 @@ void RichardsFlowProcess::assembleConcreteProcess( ...@@ -73,7 +73,7 @@ void RichardsFlowProcess::assembleConcreteProcess(
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
void RichardsFlowProcess::assembleWithJacobianConcreteProcess( void RichardsFlowProcess::assembleWithJacobianConcreteProcess(
...@@ -87,7 +87,7 @@ void RichardsFlowProcess::assembleWithJacobianConcreteProcess( ...@@ -87,7 +87,7 @@ void RichardsFlowProcess::assembleWithJacobianConcreteProcess(
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
} }
} // namespace RichardsFlow } // namespace RichardsFlow
......
...@@ -150,7 +150,7 @@ void SmallDeformationProcess<DisplacementDim>::assembleConcreteProcess( ...@@ -150,7 +150,7 @@ void SmallDeformationProcess<DisplacementDim>::assembleConcreteProcess(
// Call global assembler for each local assembly item. // Call global assembler for each local assembly item.
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers, _global_assembler, &VectorMatrixAssembler::assemble, _local_assemblers,
*_local_to_global_index_map, t, x, M, K, b, _coupling_term); *_local_to_global_index_map, t, x, M, K, b, _coupled_solutions);
} }
template <int DisplacementDim> template <int DisplacementDim>
...@@ -166,7 +166,7 @@ void SmallDeformationProcess<DisplacementDim>:: ...@@ -166,7 +166,7 @@ void SmallDeformationProcess<DisplacementDim>::
GlobalExecutor::executeMemberDereferenced( GlobalExecutor::executeMemberDereferenced(
_global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _global_assembler, &VectorMatrixAssembler::assembleWithJacobian,
_local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx,
dx_dx, M, K, b, Jac, _coupling_term); dx_dx, M, K, b, Jac, _coupled_solutions);
b.copyValues(*_nodal_forces); b.copyValues(*_nodal_forces);
std::transform(_nodal_forces->begin(), _nodal_forces->end(), std::transform(_nodal_forces->begin(), _nodal_forces->end(),
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment