#1049 Fix Stokes model.

Data/Lua/demo_input_stokes.lua

@@ -188,6 +188,31 @@ Domain3 = {
 }
 
 
+Domain4 = {
+    -- Index of the geometric mesh upon which the domain is defined (as defined in the present file). Might be
+    -- left empty if domain not limited to one mesh; at most one value is expected here.
+    -- Expected format: {VALUE1}
+    -- Expected format: VALUE
+    mesh_index = { 1 },
+    
+    -- List of dimensions encompassed by the domain. Might be left empty if no restriction at all upon
+    -- dimensions.
+    -- Expected format: {VALUE1, VALUE2, ...}
+    -- Constraint: ops_in(v, {0, 1, 2, 3})
+    dimension_list = { },
+    
+    -- List of mesh labels encompassed by the domain. Might be left empty if no restriction at all upon mesh
+    -- labels. This parameter does not make sense if no mesh is defined for the domain.
+    -- Expected format: {VALUE1, VALUE2, ...}
+    mesh_label_list = { },
+    
+    -- List of geometric element types considered in the domain. Might be left empty if no restriction upon
+    -- these.
+    -- Expected format: {"VALUE1", "VALUE2", ...}
+    geometric_element_type_list = { }
+}
+
+
 EssentialBoundaryCondition1 = {
     
     -- Name of the boundary condition (must be unique).

Sources/ModelInstances/Stokes/InputParameterList.hpp

@@ -42,7 +42,8 @@ namespace HappyHeart
         {
             highest_dimension = 1,
             neumann = 2,
-            dirichlet = 3
+            dirichlet = 3,
+            full_mesh
         };
         
         
@@ -104,6 +105,7 @@ namespace HappyHeart
             InputParameter::Domain<EnumUnderlyingType(DomainIndex::highest_dimension)>,
             InputParameter::Domain<EnumUnderlyingType(DomainIndex::neumann)>,
             InputParameter::Domain<EnumUnderlyingType(DomainIndex::dirichlet)>,
+            InputParameter::Domain<EnumUnderlyingType(DomainIndex::full_mesh)>,
         
             InputParameter::FEltSpace<EnumUnderlyingType(FEltSpaceIndex::highest_dimension)>,
             InputParameter::FEltSpace<EnumUnderlyingType(FEltSpaceIndex::neumann)>,

Sources/ModelInstances/Stokes/VariationalFormulation.cpp

@@ -12,6 +12,8 @@
 
 #include "HappyHeart/Utilities/TimeKeep/TimeKeep.hpp"
 
+#include "HappyHeart/Geometry/Domain/DomainManager.hpp"
+
 #include "HappyHeart/ModelInstances/Stokes/VariationalFormulation.hpp"
 
 
@@ -175,11 +177,12 @@ namespace HappyHeart
         void VariationalFormulation::DefineOperators(const InputParameterList& input_parameter_data)
         {
             const auto& god_of_dof = this->GetGodOfDof();
-            const auto& mesh = god_of_dof.GetGeometricMeshRegion();
+
+            decltype(auto) domain = DomainManager::GetInstance().GetDomain(EnumUnderlyingType(DomainIndex::full_mesh));
             
             fluid_viscosity_ = InitScalarParameterFromInputData<InputParameter::Fluid::Viscosity>("Viscosity",
-                                                                                     mesh,
-                                                                                     input_parameter_data);
+                                                                                                  domain,
+                                                                                                  input_parameter_data);
             
             const auto& felt_space_highest_dimension = god_of_dof.GetFEltSpace(EnumUnderlyingType(FEltSpaceIndex::highest_dimension));
             const auto& felt_space_neumann = god_of_dof.GetFEltSpace(EnumUnderlyingType(FEltSpaceIndex::neumann));