#820 Apply boundary condition to delta residual.

Sources/ModelInstances/UnderDevelopment/Poromechanics/ImplicitStepFluid/InnerLoop/VariationalFormulation.hxx

@@ -377,6 +377,26 @@ namespace HappyHeart
                         GlobalVectorWithCoefficient with_coeff(delta_residual, 1.);
                         GetSolidDeltaResidualOperator().Assemble(std::make_tuple(std::ref(with_coeff)),
                                                                  this->GetVariableHolder().GetMidpointPosition());
+                        
+                        decltype(auto) bc_manager = DirichletBoundaryConditionManager::GetInstance();
+                        
+                        decltype(auto) god_of_dof_manager = GodOfDofManager::GetInstance();
+                        
+                        decltype(auto) solid_god_of_dof =
+                        god_of_dof_manager.GetGodOfDof(EnumUnderlyingType(MeshIndex::solid));
+
+                        auto& bc1 =
+                            bc_manager.GetDirichletBoundaryCondition(EnumUnderlyingType(BoundaryConditionIndex::solid_radial));
+
+                        solid_god_of_dof.ApplyBoundaryCondition<BoundaryConditionMethod::penalization>(bc1,
+                                                                                                       delta_residual);
+
+                        auto& bc2 =
+                            bc_manager.GetDirichletBoundaryCondition(EnumUnderlyingType(BoundaryConditionIndex::solid_vertical));
+                        
+                        solid_god_of_dof.ApplyBoundaryCondition<BoundaryConditionMethod::penalization>(bc2,
+                                                                                                       delta_residual);
+                        delta_residual.Assembly(__FILE__, __LINE__);
 
                         const auto filename = parent::GetOutputDirectory(numbering_subset)
                         + "/"

Sources/ModelInstances/UnderDevelopment/Poromechanics/LocalVariationalOperatorInstances/SolidDeltaResidual.hxx

@@ -45,7 +45,6 @@ namespace HappyHeart
             delta_fluid_mass_(delta_fluid_mass)
             {
                 const auto& elementary_data = parent::GetElementaryData();
-                const unsigned int Ncomponent = elementary_data.GetGeomEltDimension();
                 const unsigned int Nnode = elementary_data.Nnode();
                 
                 const auto mesh_dimension = elementary_data.GetGeomMeshRegionDimension();