[SOLVED] Extracting tetrahedron forcefield indices

Hi everyone,

I have a forcefield which I have written in the framework of the native TatrahedronHyperelasticityFEMForceField. My addForce and addDForce implementations are quite different from the native ones.

What I would love to do is to apply internal pressurization to each tetrahedron index from my main SOFA C++ loop. I have the scene file similar to the following gist

https://gist.github.com/lakehanne/eb13b29b17e833112b0b83f40529851b

In my C++ file, I am accessing the dome ring node as follows:

    auto dome_ring = root->getChild("DomeRingNode");
    // auto dome_ring_tetras = dome_ring->getObject("dofs");
    using MechanicalObject = sofa::component::container::MechanicalObject<Vec3Types>;
    typename MechanicalObject::SPtr dome_ring_dofs = dome_ring->get<MechanicalObject>(groot->SearchDown);
    auto dome_ring_pos_ = dome_ring_dofs.get()->read(sofa::core::ConstVecCoordId::position());
 sofa::helper::ReadAccessor<Data<Vec3Types::VecCoord>>dr_pos_vecs(dome_ring_dofs.get()->read(sofa::core::ConstVecCoordId::position()));

What I am not sure of is if dr_pos_vecs contain the dofs for the tetras only or for just the Vec3d templated mechanical object. It appears that the size of this vector is different from that of the ones I created with CGAL. Would appreciate your help.