Get dofs/node for SolidLine element from parent

Author: amgebauer

src/solid_3D_ele/4C_solid_3D_ele_line.cpp

@@ -35,6 +35,13 @@ std::shared_ptr<Core::Elements::Element> Discret::Elements::SolidLineType<dim>::
   return nullptr;
 }
 
+
+template <unsigned dim>
+Discret::Elements::SolidLine<dim>::SolidLine(const Discret::Elements::SolidLine<dim>& old) noexcept
+    : Core::Elements::FaceElement(old), num_dof_per_node_(old.num_dof_per_node_)
+{
+}
+
 template <unsigned dim>
 Discret::Elements::SolidLine<dim>::SolidLine(int id, int owner, int nnode, const int* nodeids,
     Core::Nodes::Node** nodes, Core::Elements::Element* parent, const int lline)
@@ -43,6 +50,15 @@ Discret::Elements::SolidLine<dim>::SolidLine(int id, int owner, int nnode, const
   set_node_ids(nnode, nodeids);
   build_nodal_pointers(nodes);
   set_parent_master_element(parent, lline);
+
+  num_dof_per_node_ = parent_element()->num_dof_per_node(*SolidLine::nodes()[0]);
+  // Safety check if all nodes have the same number of dofs!
+  for (int nlid = 1; nlid < num_node(); ++nlid)
+  {
+    if (num_dof_per_node_ != parent_master_element()->num_dof_per_node(*SolidLine::nodes()[nlid]))
+      FOUR_C_THROW("You need different NumDofPerNode for each node on this solid line? ({} != {})",
+          num_dof_per_node_, parent_master_element()->num_dof_per_node(*SolidLine::nodes()[nlid]));
+  }
 }
 
 template <unsigned dim>
@@ -72,6 +88,18 @@ Core::FE::CellType Discret::Elements::SolidLine<dim>::shape() const
       });
 }
 
+template <unsigned dim>
+void Discret::Elements::SolidLine<dim>::pack(Core::Communication::PackBuffer& data) const
+{
+  data.add_to_pack(num_dof_per_node_);
+}
+
+template <unsigned dim>
+void Discret::Elements::SolidLine<dim>::unpack(Core::Communication::UnpackBuffer& buffer)
+{
+  buffer.extract_from_pack(num_dof_per_node_);
+}
+
 template <unsigned dim>
 void Discret::Elements::SolidLine<dim>::print(std::ostream& os) const
 {

src/solid_3D_ele/4C_solid_3D_ele_line.hpp

@@ -68,22 +68,24 @@ namespace Discret::Elements
     SolidLine(int id, int owner, int nnode, const int* nodeids, Core::Nodes::Node** nodes,
         Core::Elements::Element* parent, const int lline);
 
+    explicit SolidLine(const SolidLine<dim>& old) noexcept;
+
     [[nodiscard]] Core::Elements::Element* clone() const override;
 
     [[nodiscard]] inline int unique_par_object_id() const override
     {
       return SolidLineType<dim>::instance().unique_par_object_id();
     }
 
-    void pack(Core::Communication::PackBuffer& data) const override {};
+    void pack(Core::Communication::PackBuffer& data) const override;
 
-    void unpack(Core::Communication::UnpackBuffer& buffer) override {};
+    void unpack(Core::Communication::UnpackBuffer& buffer) override;
 
     [[nodiscard]] Core::FE::CellType shape() const override;
 
     [[nodiscard]] inline int num_dof_per_node(const Core::Nodes::Node& node) const override
     {
-      return dim;
+      return num_dof_per_node_;
     }
 
     [[nodiscard]] inline int num_dof_per_element() const override { return 0; }
@@ -99,6 +101,9 @@ namespace Discret::Elements
         const Core::Conditions::Condition& condition, std::vector<int>& lm,
         Core::LinAlg::SerialDenseVector& elevec1,
         Core::LinAlg::SerialDenseMatrix* elemat1 = nullptr) override;
+
+   private:
+    int num_dof_per_node_ = 0;
   };
 }  // namespace Discret::Elements