[PL/BC] Set coordinates in spatial positions.

ProcessLib/BoundaryCondition/DirichletBoundaryConditionAuxiliaryFunctions.cpp

@@ -94,6 +94,7 @@ void getEssentialBCValuesLocal(
         if (global_index >= 0)
         {
             pos.setNodeID(id);
+            pos.setCoordinates(*node);
             bc_values.ids.emplace_back(global_index);
             bc_values.values.emplace_back(parameter(t, pos).front());
         }

ProcessLib/BoundaryCondition/PrimaryVariableConstraintDirichletBoundaryCondition.cpp

@@ -87,6 +87,7 @@ void PrimaryVariableConstraintDirichletBoundaryCondition::getEssentialBCValues(
             // fetch the value of the primary variable
             auto const local_x = x.get(std::vector{global_index});
             pos.setNodeID(id);
+            pos.setCoordinates(*node);
             if (_less && local_x[0] < _threshold_parameter(t, pos).front())
             {
                 bc_values.ids.emplace_back(global_index);

ProcessLib/BoundaryCondition/RobinBoundaryConditionLocalAssembler.h

@@ -31,6 +31,7 @@ class RobinBoundaryConditionLocalAssembler final
 {
     using Base = GenericNaturalBoundaryConditionLocalAssembler<
         ShapeFunction, IntegrationMethod, GlobalDim>;
+    using ShapeMatricesType = ShapeMatrixPolicyType<ShapeFunction, GlobalDim>;
 
 public:
     RobinBoundaryConditionLocalAssembler(MeshLib::Element const& e,
@@ -65,16 +66,18 @@ public:
             _data.u_0.getNodalValuesOnElement(Base::_element, t)
                 .template topRows<ShapeFunction::MeshElement::n_all_nodes>();
 
-        ParameterLib::SpatialPosition position;
-        position.setElementID(Base::_element.getID());
-
         for (unsigned ip = 0; ip < n_integration_points; ++ip)
         {
-            position.setIntegrationPoint(ip);
             auto const& ip_data = Base::_ns_and_weights[ip];
             auto const& N = ip_data.N;
             auto const& w = ip_data.weight;
 
+            ParameterLib::SpatialPosition const position{
+                boost::none, Base::_element.getID(), ip,
+                MathLib::Point3d(
+                    interpolateCoordinates<ShapeFunction, ShapeMatricesType>(
+                        Base::_element, N))};
+
             double integral_measure = 1.0;
             if (_data.integral_measure)
             {