diff --git a/ProcessLib/HydroMechanics/HydroMechanicsProcess-impl.h b/ProcessLib/HydroMechanics/HydroMechanicsProcess-impl.h index 29c79cd51255e8fcb869f8c632da0df7845e09fd..d9804f7ce7e467155d7fb95c79399fb8c6362ee2 100644 --- a/ProcessLib/HydroMechanics/HydroMechanicsProcess-impl.h +++ b/ProcessLib/HydroMechanics/HydroMechanicsProcess-impl.h @@ -298,31 +298,26 @@ void HydroMechanicsProcess<DisplacementDim>:: } // For the staggered scheme - // For the equations of displacement - if (_coupled_solutions->process_id == 1) + if (_coupled_solutions->process_id == 0) + { + DBUG( + "Assemble the Jacobian equations of liquid fluid process in " + "HydroMechanics for the staggered scheme."); + } + else { DBUG( "Assemble the Jacobian equations of mechanical process in " "HydroMechanics for the staggered scheme."); - - GlobalExecutor::executeMemberDereferenced( - _global_assembler, &VectorMatrixAssembler::assembleWithJacobian, - _local_assemblers, *_local_to_global_index_map, t, x, xdot, - dxdot_dx, dx_dx, M, K, b, Jac, _coupled_solutions, - _local_to_global_index_map_with_base_nodes.get()); - return; } - // For the equations of pressure - DBUG( - "Assemble the Jacobian equations of liquid fluid process in " - "HydroMechanics for the staggered scheme."); GlobalExecutor::executeMemberDereferenced( _global_assembler, &VectorMatrixAssembler::assembleWithJacobian, _local_assemblers, *_local_to_global_index_map, t, x, xdot, dxdot_dx, dx_dx, M, K, b, Jac, _coupled_solutions, _local_to_global_index_map_with_base_nodes.get()); } + template <int DisplacementDim> void HydroMechanicsProcess<DisplacementDim>::preTimestepConcreteProcess( GlobalVector const& x, double const t, double const dt,