SimpleTesselatedTetraTopologicalMapping

Special case of mapping where TetrahedronSetTopology is converted into a finer TetrahedronSetTopology

Target: Sofa.Component.Mapping.Linear

namespace: sofa::component::mapping::linear

parents:

• TopologicalMapping

Data

Name	Description	Default value
name	object name	unnamed
printLog	if true, emits extra messages at runtime.	0
tags	list of the subsets the objet belongs to
bbox	this object bounding box
componentState	The state of the component among (Dirty, Valid, Undefined, Loading, Invalid).	Undefined
listening	if true, handle the events, otherwise ignore the events	0
tetrahedraMappedFromTetra	Each Tetrahedron of the input topology is mapped to the 8 tetrahedrons in which it can be divided
tetraSource	Which tetra from the input topology map to a given tetra in the output topology (sofa::InvalidID if none)
pointMappedFromPoint	Each point of the input topology is mapped to the same point
pointMappedFromEdge	Each edge of the input topology is mapped to his midpoint
pointSource	Which input topology element map to a given point in the output topology : 0 -> none, > 0 -> point index + 1, < 0 , - edge index -1

Links

Name	Description	Destination type name
context	Graph Node containing this object (or BaseContext::getDefault() if no graph is used)	BaseContext
slaves	Sub-objects used internally by this object	BaseObject
master	nullptr for regular objects, or master object for which this object is one sub-objects	BaseObject
input	Input topology to map	BaseMeshTopology
output	Output topology to map	BaseMeshTopology

Examples

SimpleTesselatedTetraTopologicalMapping.scn

XMLPython

<!-- -->
<Node name="root" dt="0.05" showBoundingTree="0" gravity="0 0 0">
    <RequiredPlugin name="Sofa.Component.Collision.Detection.Algorithm"/> <!-- Needed to use components [BVHNarrowPhase BruteForceBroadPhase CollisionPipeline] -->
    <RequiredPlugin name="Sofa.Component.Collision.Detection.Intersection"/> <!-- Needed to use components [DiscreteIntersection] -->
    <RequiredPlugin name="Sofa.Component.Collision.Geometry"/> <!-- Needed to use components [TriangleCollisionModel] -->
    <RequiredPlugin name="Sofa.Component.Collision.Response.Contact"/> <!-- Needed to use components [CollisionResponse] -->
    <RequiredPlugin name="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedProjectiveConstraint] -->
    <RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshGmshLoader] -->
    <RequiredPlugin name="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
    <RequiredPlugin name="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [SimpleTesselatedTetraMechanicalMapping SimpleTesselatedTetraTopologicalMapping] -->
    <RequiredPlugin name="Sofa.Component.Mass"/> <!-- Needed to use components [DiagonalMass] -->
    <RequiredPlugin name="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->
    <RequiredPlugin name="Sofa.Component.SolidMechanics.FEM.Elastic"/> <!-- Needed to use components [TetrahedralCorotationalFEMForceField] -->
    <RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
    <RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TetrahedronSetGeometryAlgorithms TetrahedronSetTopologyContainer TetrahedronSetTopologyModifier] -->
    <RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->
    <VisualStyle displayFlags="showBehaviorModels showForceFields showCollisionModels showVisual" />
    <CollisionPipeline verbose="0" name="CollisionPipeline" />
    <BruteForceBroadPhase/>
    <BVHNarrowPhase/>
    <CollisionResponse response="PenalityContactForceField" name="collision response" />
    <DiscreteIntersection/>
    <DefaultAnimationLoop/>

    <Node name="TetraTopology1">
        <EulerImplicitSolver name="cg_odesolver" printLog="false"  rayleighStiffness="0.1" rayleighMass="0.1" />
        <CGLinearSolver iterations="25" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" />
        <MeshGmshLoader name="meshLoader0" filename="mesh/liver.msh" />
        <TetrahedronSetTopologyContainer name="Container1" src="@meshLoader0" />
        <TetrahedronSetTopologyModifier />
        <TetrahedronSetGeometryAlgorithms template="Vec3" drawEdges="1" drawColorEdges="0 1 0" />
        <MechanicalObject name="dofs" />
        <FixedProjectiveConstraint name="FixedProjectiveConstraint" indices="3 39 64" />
        <DiagonalMass massDensity="1" name="computed using mass density" />
        <TetrahedralCorotationalFEMForceField name="FEM" youngModulus="3000" poissonRatio="0.3" computeGlobalMatrix="false" method="large" />
        <Node name="TetraTopology2">
            <TetrahedronSetTopologyContainer name="Container2" />
            <TetrahedronSetTopologyModifier />
            <TetrahedronSetGeometryAlgorithms template="Vec3" drawTetrahedra="1" drawColorTetrahedra="1 0 0 1" />
            <SimpleTesselatedTetraTopologicalMapping input="@Container1" output="@Container2" />
            <MechanicalObject />
            <SimpleTesselatedTetraMechanicalMapping />
            <TriangleCollisionModel />
        </Node>
    </Node>
</Node>

def createScene(root_node):

   root = root_node.addChild('root', dt="0.05", showBoundingTree="0", gravity="0 0 0")

   root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Detection.Algorithm")
   root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Detection.Intersection")
   root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Geometry")
   root.addObject('RequiredPlugin', name="Sofa.Component.Collision.Response.Contact")
   root.addObject('RequiredPlugin', name="Sofa.Component.Constraint.Projective")
   root.addObject('RequiredPlugin', name="Sofa.Component.IO.Mesh")
   root.addObject('RequiredPlugin', name="Sofa.Component.LinearSolver.Iterative")
   root.addObject('RequiredPlugin', name="Sofa.Component.Mapping.Linear")
   root.addObject('RequiredPlugin', name="Sofa.Component.Mass")
   root.addObject('RequiredPlugin', name="Sofa.Component.ODESolver.Backward")
   root.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.FEM.Elastic")
   root.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
   root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Dynamic")
   root.addObject('RequiredPlugin', name="Sofa.Component.Visual")
   root.addObject('VisualStyle', displayFlags="showBehaviorModels showForceFields showCollisionModels showVisual")
   root.addObject('CollisionPipeline', verbose="0", name="CollisionPipeline")
   root.addObject('BruteForceBroadPhase', )
   root.addObject('BVHNarrowPhase', )
   root.addObject('CollisionResponse', response="PenalityContactForceField", name="collision response")
   root.addObject('DiscreteIntersection', )
   root.addObject('DefaultAnimationLoop', )

   tetra_topology1 = root.addChild('TetraTopology1')

   tetra_topology1.addObject('EulerImplicitSolver', name="cg_odesolver", printLog="false", rayleighStiffness="0.1", rayleighMass="0.1")
   tetra_topology1.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
   tetra_topology1.addObject('MeshGmshLoader', name="meshLoader0", filename="mesh/liver.msh")
   tetra_topology1.addObject('TetrahedronSetTopologyContainer', name="Container1", src="@meshLoader0")
   tetra_topology1.addObject('TetrahedronSetTopologyModifier', )
   tetra_topology1.addObject('TetrahedronSetGeometryAlgorithms', template="Vec3", drawEdges="1", drawColorEdges="0 1 0")
   tetra_topology1.addObject('MechanicalObject', name="dofs")
   tetra_topology1.addObject('FixedProjectiveConstraint', name="FixedProjectiveConstraint", indices="3 39 64")
   tetra_topology1.addObject('DiagonalMass', massDensity="1", name="computed using mass density")
   tetra_topology1.addObject('TetrahedralCorotationalFEMForceField', name="FEM", youngModulus="3000", poissonRatio="0.3", computeGlobalMatrix="false", method="large")

   tetra_topology2 = TetraTopology1.addChild('TetraTopology2')

   tetra_topology2.addObject('TetrahedronSetTopologyContainer', name="Container2")
   tetra_topology2.addObject('TetrahedronSetTopologyModifier', )
   tetra_topology2.addObject('TetrahedronSetGeometryAlgorithms', template="Vec3", drawTetrahedra="1", drawColorTetrahedra="1 0 0 1")
   tetra_topology2.addObject('SimpleTesselatedTetraTopologicalMapping', input="@Container1", output="@Container2")
   tetra_topology2.addObject('MechanicalObject', )
   tetra_topology2.addObject('SimpleTesselatedTetraMechanicalMapping', )
   tetra_topology2.addObject('TriangleCollisionModel', )

SimpleTesselatedTetraTopologicalMapping.scn

XMLPython

<Node name="root" dt="0.05" showBehaviorModels="1" showCollisionModels="1" showMappings="0" showForceFields="1" showBoundingTree="0" gravity="0 0 0">
    <Node name="RequiredPlugins">
        <RequiredPlugin name="Sofa.Component.Collision.Detection.Algorithm"/> <!-- Needed to use components [BVHNarrowPhase, BruteForceBroadPhase, CollisionPipeline] -->  
        <RequiredPlugin name="Sofa.Component.Collision.Geometry"/> <!-- Needed to use components [TriangleCollisionModel] -->  
        <RequiredPlugin name="Sofa.Component.Collision.Response.Contact"/> <!-- Needed to use components [CollisionResponse] -->  
        <RequiredPlugin name="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedConstraint] -->  
        <RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshGmshLoader] -->  
        <RequiredPlugin name="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->  
        <RequiredPlugin name="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [SimpleTesselatedTetraMechanicalMapping, SimpleTesselatedTetraTopologicalMapping] -->  
        <RequiredPlugin name="Sofa.Component.Mass"/> <!-- Needed to use components [DiagonalMass] -->  
        <RequiredPlugin name="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->  
        <RequiredPlugin name="Sofa.Component.SolidMechanics.FEM.Elastic"/> <!-- Needed to use components [TetrahedralCorotationalFEMForceField] -->  
        <RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->  
        <RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TetrahedronSetGeometryAlgorithms, TetrahedronSetTopologyContainer, TetrahedronSetTopologyModifier] -->  
        <RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->  
    </Node>
    <CollisionPipeline verbose="0" name="CollisionPipeline"/>
    <BruteForceBroadPhase/>
    <BVHNarrowPhase/>
    <CollisionResponse response="PenalityContactForceField" name="collision response"/>
    <Node name="TetraTopology1">
        <EulerImplicitSolver name="cg_odesolver" printLog="false" rayleighStiffness="0.1" rayleighMass="0.1" />
        <CGLinearSolver iterations="25" name="linear solver" tolerance="1.0e-9" threshold="1.0e-9" />
        <MeshGmshLoader name="meshLoader0" filename="mesh/liver.msh" />
        <TetrahedronSetTopologyContainer name="Container1" src="@meshLoader0" />
        <TetrahedronSetTopologyModifier/>
        <TetrahedronSetGeometryAlgorithms template="Vec3d" />
        <MechanicalObject name="dofs" />
        <FixedConstraint name="FixedConstraint" indices="3 39 64" />
        <DiagonalMass massDensity="1" name="computed using mass density"/>
        <TetrahedralCorotationalFEMForceField name="FEM" youngModulus="500" poissonRatio="0.3" computeGlobalMatrix="false" method="large"/> 
        <Node name="TetraTopology2">
            <TetrahedronSetTopologyContainer name="Container2"/>
            <TetrahedronSetTopologyModifier/>
            <TetrahedronSetGeometryAlgorithms template="Vec3d" />
            <SimpleTesselatedTetraTopologicalMapping input="@Container1" output="@Container2"/>
            <MechanicalObject/>
            <SimpleTesselatedTetraMechanicalMapping/>
            <TriangleCollisionModel />
        </Node>
    </Node>
</Node>

def createScene(root_node):

   root = root_node.addChild('root', dt="0.05", showBehaviorModels="1", showCollisionModels="1", showMappings="0", showForceFields="1", showBoundingTree="0", gravity="0 0 0")

   required_plugins = root.addChild('RequiredPlugins')

   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Collision.Detection.Algorithm")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Collision.Geometry")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Collision.Response.Contact")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Constraint.Projective")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.IO.Mesh")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.LinearSolver.Iterative")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Mapping.Linear")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Mass")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.ODESolver.Backward")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.FEM.Elastic")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Dynamic")
   required_plugins.addObject('RequiredPlugin', name="Sofa.Component.Visual")

   root.addObject('CollisionPipeline', verbose="0", name="CollisionPipeline")
   root.addObject('BruteForceBroadPhase', )
   root.addObject('BVHNarrowPhase', )
   root.addObject('CollisionResponse', response="PenalityContactForceField", name="collision response")

   tetra_topology1 = root.addChild('TetraTopology1')

   tetra_topology1.addObject('EulerImplicitSolver', name="cg_odesolver", printLog="false", rayleighStiffness="0.1", rayleighMass="0.1")
   tetra_topology1.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
   tetra_topology1.addObject('MeshGmshLoader', name="meshLoader0", filename="mesh/liver.msh")
   tetra_topology1.addObject('TetrahedronSetTopologyContainer', name="Container1", src="@meshLoader0")
   tetra_topology1.addObject('TetrahedronSetTopologyModifier', )
   tetra_topology1.addObject('TetrahedronSetGeometryAlgorithms', template="Vec3d")
   tetra_topology1.addObject('MechanicalObject', name="dofs")
   tetra_topology1.addObject('FixedConstraint', name="FixedConstraint", indices="3 39 64")
   tetra_topology1.addObject('DiagonalMass', massDensity="1", name="computed using mass density")
   tetra_topology1.addObject('TetrahedralCorotationalFEMForceField', name="FEM", youngModulus="500", poissonRatio="0.3", computeGlobalMatrix="false", method="large")

   tetra_topology2 = TetraTopology1.addChild('TetraTopology2')

   tetra_topology2.addObject('TetrahedronSetTopologyContainer', name="Container2")
   tetra_topology2.addObject('TetrahedronSetTopologyModifier', )
   tetra_topology2.addObject('TetrahedronSetGeometryAlgorithms', template="Vec3d")
   tetra_topology2.addObject('SimpleTesselatedTetraTopologicalMapping', input="@Container1", output="@Container2")
   tetra_topology2.addObject('MechanicalObject', )
   tetra_topology2.addObject('SimpleTesselatedTetraMechanicalMapping', )
   tetra_topology2.addObject('TriangleCollisionModel', )