use link instead of raw pointers

Author: fredroy

src/BeamAdapter/component/BaseBeamInterpolation.h

@@ -84,7 +84,7 @@ class BaseBeamInterpolation : public virtual sofa::core::objectmodel::BaseObject
     using VecEdgeID = type::vector<BaseMeshTopology::EdgeID>;
     using VecEdges = type::vector<BaseMeshTopology::Edge>;
 
-    BaseBeamInterpolation(/*sofa::component::engine::WireRestShape<DataTypes> *_restShape = nullptr*/);
+    BaseBeamInterpolation();
 
     virtual ~BaseBeamInterpolation() = default;
 
@@ -212,7 +212,6 @@ class BaseBeamInterpolation : public virtual sofa::core::objectmodel::BaseObject
     typename MechanicalObject<StateDataTypes>::SPtr m_StateNodes;
 
     Data< VecEdgeID > d_edgeList;
-    const VecEdges* m_topologyEdges{ nullptr };
 
     ///2. Vector of length of each beam. Same size as @sa d_edgeList
     Data< type::vector< Real > >    d_lengthList;
@@ -230,13 +229,12 @@ class BaseBeamInterpolation : public virtual sofa::core::objectmodel::BaseObject
     Data< sofa::type::vector<EdgeID> > d_beamCollision;
 
     Data<bool>          d_dofsAndBeamsAligned;
+    
+    /// link to the (edge) topology
+    SingleLink<BaseBeamInterpolation<DataTypes>, BaseMeshTopology, BaseLink::FLAG_STOREPATH|BaseLink::FLAG_STRONGLINK> l_topology;
 
-
-    /// pointer to the topology
-    BaseMeshTopology* m_topology{ nullptr };
-
-    /// pointer on mechanical state
-    MechanicalState<DataTypes>* m_mstate{ nullptr };
+    /// link to the mechanical state
+    SingleLink<BaseBeamInterpolation<DataTypes>, MechanicalState<DataTypes>, BaseLink::FLAG_STOREPATH|BaseLink::FLAG_STRONGLINK> l_mstate;
 };
 
 

src/BeamAdapter/component/BaseBeamInterpolation.inl

@@ -96,7 +96,7 @@ void BaseBeamInterpolation<DataTypes>::RotateFrameForAlignNormalizedX(const Quat
 
 
 template <class DataTypes>
-BaseBeamInterpolation<DataTypes>::BaseBeamInterpolation(/*sofa::component::engine::WireRestShape<DataTypes> *_restShape*/)
+BaseBeamInterpolation<DataTypes>::BaseBeamInterpolation()
     : m_StateNodes(sofa::core::objectmodel::New< sofa::component::statecontainer::MechanicalObject<StateDataTypes> >())
     , d_edgeList(initData(&d_edgeList, "edgeList", "list of the edge in the topology that are concerned by the Interpolation"))
     , d_lengthList(initData(&d_lengthList, "lengthList", "list of the length of each beam"))
@@ -106,8 +106,8 @@ BaseBeamInterpolation<DataTypes>::BaseBeamInterpolation(/*sofa::component::engin
     , d_beamCollision(initData(&d_beamCollision, "beamCollision", "list of beam (in edgeList) that needs to be considered for collision"))
     , d_dofsAndBeamsAligned(initData(&d_dofsAndBeamsAligned, true, "dofsAndBeamsAligned",
         "if false, a transformation for each beam is computed between the DOF and the beam nodes"))
-    , m_topology(nullptr)
-    , m_mstate(nullptr)
+    , l_topology(initLink("topology", "link to the topology (must contain edges)"))
+    , l_mstate(initLink("mstate", "link to the mechanical container"))
 {
 
 
@@ -122,33 +122,47 @@ void BaseBeamInterpolation<DataTypes>::init()
     Inherit::init();
 
     BaseContext* context = getContext();
-
-    m_mstate = dynamic_cast<sofa::core::behavior::MechanicalState<DataTypes> *> (context->getMechanicalState());
-    if (m_mstate == nullptr)
+    
+    if (!l_mstate)
     {
-        msg_error() << "No MechanicalState found. Component is de-activated.";
-        d_componentState.setValue(ComponentState::Invalid);
+        auto mstate = dynamic_cast<sofa::core::behavior::MechanicalState<DataTypes> *> (context->getMechanicalState());
+        l_mstate.set(mstate);
+    }
+    if (l_mstate)
+    {
+        msg_info() << "Found mechanical object named "<< l_mstate->getName() ;
+    }
+    else
+    {
+        msg_error() << "Cannot find the mechanical object. Please specify the link to the topology or insert one in the same node.";
+        this->d_componentState.setValue(sofa::core::objectmodel::ComponentState::Invalid);
         return;
     }
+    
+    if (!l_topology)
+    {
+        l_topology.set(this->getContext()->getMeshTopologyLink());
+    }
 
-    /// Get the topology from context and check if it is valid.
-    m_topology = context->getMeshTopology();
-    if (!m_topology)
+    if (l_topology)
     {
-        msg_error() << "No Topology found. Component is de-activated.";
-        d_componentState.setValue(ComponentState::Invalid);
+        msg_info() << "Found topology named "<< l_topology->getName() ;
+    }
+    else
+    {
+        msg_error() << "Cannot find a topology container. Please specify the link to the topology or insert one in the same node.";
+        this->d_componentState.setValue(sofa::core::objectmodel::ComponentState::Invalid);
         return;
     }
-
-    /// Check the topology properties
-    if (m_topology->getNbEdges() == 0)
+    
+    if(l_topology->getNbEdges() == 0)
     {
-        msg_error() << "Found a topology but it is empty. Component is de-activated";
+        msg_error() << "Found a topology but it is empty (no edges). Component is de-activated";
         d_componentState.setValue(ComponentState::Invalid);
         return;
     }
-
-    m_topologyEdges = &m_topology->getEdges();
+    
+    this->d_componentState.setValue(sofa::core::objectmodel::ComponentState::Valid);
 }
 
 
@@ -402,15 +416,15 @@ int BaseBeamInterpolation<DataTypes>::getNodeIndices(const EdgeID edgeInList,
     unsigned int& node0Idx,
     unsigned int& node1Idx)
 {
-    if (m_topologyEdges == nullptr)
+    if(!this->isComponentStateValid())
     {
-        msg_error() << "This object does not have edge topology defined (computation halted). ";
+        msg_error() << "(getNodeIndices) This component is invalid, check the other error messages. ";
         return -1;
     }
 
     /// 1. Get the indices of element and nodes
     const EdgeID& e = d_edgeList.getValue()[edgeInList];
-    const BaseMeshTopology::Edge& edge = (*m_topologyEdges)[e];
+    const BaseMeshTopology::Edge& edge = l_topology->getEdge(e);
     node0Idx = edge[0];
     node1Idx = edge[1];
 
@@ -437,14 +451,14 @@ void BaseBeamInterpolation<DataTypes>::getSplinePoints(const EdgeID edgeInList,
 template<class DataTypes>
 unsigned int BaseBeamInterpolation<DataTypes>::getStateSize() const
 {
-    if (m_mstate == nullptr)
+    if(!this->isComponentStateValid())
     {
-        msg_error() << "No _mstate found (Aborting)";
+        msg_error() << "(getStateSize) This component is invalid, check the other error messages. ";
         return 0;
     }
     else
     {
-        return m_mstate->getSize();
+        return l_mstate->getSize();
     }
 }
 

src/BeamAdapter/component/BeamInterpolation.inl

@@ -139,7 +139,7 @@ void BeamInterpolation<DataTypes>::computeCrossSectionInertiaMatrix()
     }
     else
     {
-        Size nbEdges = this->m_topology->getNbEdges();
+        Size nbEdges = this->l_topology->getNbEdges();
         m_section.resize(nbEdges);
         if ( crossSectionShape.getValue().getSelectedItem() == "elliptic")
         {
@@ -197,7 +197,7 @@ void BeamInterpolation<DataTypes>::init()
     if (this->d_componentState.getValue() == ComponentState::Invalid)
         return;
 
-    Size nbEdges = this->m_topology->getNbEdges();
+    Size nbEdges = this->l_topology->getNbEdges();
     checkDataSize(m_defaultRadius, d_radius, nbEdges);
     checkDataSize(m_defaultInnerRadius, d_innerRadius, nbEdges);
     checkDataSize(m_defaultLengthY, d_lengthY, nbEdges);
@@ -211,7 +211,7 @@ void BeamInterpolation<DataTypes>::init()
         auto edgeList = sofa::helper::getWriteOnlyAccessor(this->d_edgeList);
         edgeList.clear();
 
-        for (unsigned int i=0; i<this->m_topology->getNbEdges(); i++)
+        for (unsigned int i=0; i<this->l_topology->getNbEdges(); i++)
         {
             edgeList.push_back(i);
         }
@@ -225,7 +225,7 @@ void BeamInterpolation<DataTypes>::init()
             DOF1TransformNode1.resize(edgeList.size());
         }
 
-        ReadAccessor<Data<VecCoord> > statePos = this->m_mstate->read(sofa::core::vec_id::read_access::position) ;
+        ReadAccessor<Data<VecCoord> > statePos = this->l_mstate->read(sofa::core::vec_id::read_access::position) ;
 
         auto lengthList = sofa::helper::getWriteOnlyAccessor(this->d_lengthList);
         lengthList.clear();
@@ -335,17 +335,18 @@ bool BeamInterpolation<DataTypes>::interpolationIsAlreadyInitialized()
 template<class DataTypes>
 bool BeamInterpolation<DataTypes>::verifyTopology()
 {
+    const auto nbEdges = this->l_topology->getNbEdges();
+    
     //TODO(dmarchal) This contains "code" specific slang that cannot be understood by user.
-    dmsg_info() << "The vector _topologyEdges is now set with " << this->m_topologyEdges->size() << " edges" ;
-
+    dmsg_info() << "The vector _topologyEdges is now set with " << nbEdges << " edges" ;
 
     const VecElementID &edgeList = this->d_edgeList.getValue();
     for (unsigned int j = 0; j < edgeList.size(); j++)
     {
-        if(edgeList[j] > this->m_topologyEdges->size())
+        if(edgeList[j] > nbEdges)
         {
             msg_error() << "The provided edge index '" << edgeList[j] << "'is larger than '"
-                        << this->m_topologyEdges->size() << "' the amount of edges in the topology. " ;
+                        << nbEdges << "' the amount of edges in the topology. " ;
             return false;
         }
     }
@@ -420,7 +421,7 @@ void BeamInterpolation<DataTypes>::getMechanicalSampling(Real& dx, const Real x_
 {
     SOFA_UNUSED(x_localcurv_abs);
 
-    const auto numLines = this->m_topologyEdges->size();
+    const auto numLines = this->l_topology->getNbEdges();
     dx = getRestTotalLength()/numLines;
 }
 
@@ -429,14 +430,14 @@ void BeamInterpolation<DataTypes>::getCollisionSampling(Real &dx, const Real x_l
 {
     SOFA_UNUSED(x_localcurv_abs);
 
-    const auto numLines = this->m_topologyEdges->size();
+    const auto numLines = this->l_topology->getNbEdges();
     dx = getRestTotalLength()/numLines;
 }
 
 template <class DataTypes>
 void BeamInterpolation<DataTypes>::getNumberOfCollisionSegment(Real &dx, unsigned int &numLines)
 {
-    numLines = static_cast<unsigned int>(this->m_topologyEdges->size());
+    numLines = static_cast<unsigned int>(this->l_topology->getNbEdges());
     dx = getRestTotalLength()/numLines;
 }
 
@@ -599,16 +600,16 @@ void BeamInterpolation<DataTypes>::updateInterpolation(){
     if(d_vecID.getValue().getSelectedItem() == "current")
     {
         dmsg_info() <<" position " << msgendl
-                   << "      ="<< this->m_mstate->read( sofa::core::vec_id::read_access::position )->getValue( ) ;
-        x=this->m_mstate->read( sofa::core::vec_id::read_access::position );
+                   << "      ="<< this->l_mstate->read( sofa::core::vec_id::read_access::position )->getValue( ) ;
+        x=this->l_mstate->read( sofa::core::vec_id::read_access::position );
     }
     else if(d_vecID.getValue().getSelectedItem() == "free")
     {
-        x=this->m_mstate->read( sofa::core::vec_id::read_access::freePosition ) ;
+        x=this->l_mstate->read( sofa::core::vec_id::read_access::freePosition ) ;
     }
     else /// rest position
     {
-        x=this->m_mstate->read( sofa::core::vec_id::read_access::restPosition ) ;
+        x=this->l_mstate->read( sofa::core::vec_id::read_access::restPosition ) ;
         computeVel = false;
     }