SubsetTopologicalMapping

TopologicalMapping where the destination topology is a subset of the source topology. The implementation currently assumes that both topologies have been initialized correctly.

Target: Sofa.Component.Topology.Mapping

namespace: sofa::component::topology::mapping

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 object 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
samePoints	True if the same set of points is used in both topologies	0
handleEdges	True if edges events and mapping should be handled	0
handleTriangles	True if triangles events and mapping should be handled	0
handleQuads	True if quads events and mapping should be handled	0
handleTetrahedra	True if tetrahedra events and mapping should be handled	0
handleHexahedra	True if hexahedra events and mapping should be handled	0
pointS2D	Internal source -> destination topology points map
pointD2S	Internal destination -> source topology points map (link to SubsetMapping::indices to handle the mechanical-side of the mapping
edgeS2D	Internal source -> destination topology edges map
edgeD2S	Internal destination -> source topology edges map
triangleS2D	Internal source -> destination topology triangles map
triangleD2S	Internal destination -> source topology triangles map
quadS2D	Internal source -> destination topology quads map
quadD2S	Internal destination -> source topology quads map
tetrahedronS2D	Internal source -> destination topology tetrahedra map
tetrahedronD2S	Internal destination -> source topology tetrahedra map
hexahedronS2D	Internal source -> destination topology hexahedra map
hexahedronD2S	Internal destination -> source topology hexahedra map

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

SubsetTopologicalMapping2.scn

XMLPython

<?xml version="1.0" ?>
<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.Engine.Select"/> <!-- Needed to use components [BoxROI] -->
    <RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshOBJLoader] -->
    <RequiredPlugin name="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
    <RequiredPlugin name="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [IdentityMapping] -->
    <RequiredPlugin name="Sofa.Component.Mass"/> <!-- Needed to use components [DiagonalMass] -->
    <RequiredPlugin name="Sofa.Component.MechanicalLoad"/> <!-- Needed to use components [TaitSurfacePressureForceField] -->
    <RequiredPlugin name="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->
    <RequiredPlugin name="Sofa.Component.SolidMechanics.FEM.Elastic"/> <!-- Needed to use components [TriangularFEMForceFieldOptim] -->
    <RequiredPlugin name="Sofa.Component.SolidMechanics.Spring"/> <!-- Needed to use components [FastTriangularBendingSprings] -->
    <RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
    <RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TriangleSetGeometryAlgorithms TriangleSetTopologyContainer TriangleSetTopologyModifier] -->
    <RequiredPlugin name="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [SubsetTopologicalMapping] -->
    <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/>

<!--    <SplitAndAttachBodyButtonSetting button="Middle" stiffness="5000" arrowSize="0.2" snapDistance="2" />-->
    <Node name="FullTopology">
        <MeshOBJLoader name="loader" filename="mesh/sphere_05.obj" />
        <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" />
        <TriangleSetTopologyContainer name="Container1" src="@loader" />
        <TriangleSetTopologyModifier />
        <TriangleSetGeometryAlgorithms name="GeomAlgo" />
        <MechanicalObject name="dofs" />
        <BoxROI name="roi1" box="-55 -55 -55 55 55 -25" />
        <FixedProjectiveConstraint name="FixedProjectiveConstraint" indices="@roi1.indices" />
        <DiagonalMass massDensity="1" name="computed using mass density" />
        <TriangleCollisionModel group="1" printLog="1" />
        <TaitSurfacePressureForceField name="Pressure" gamma="5" B="10000" injectedVolume="0" printLog="1" />
        <Node name="SubsetTopology1">
            <BoxROI name="subsetROI" position="@../loader.position" triangles="@../loader.triangles" computeEdges="0" computeTetrahedra="0" box="-55 -55 25 55 55 55" />
            <TriangleSetTopologyContainer name="Container2" position="@../loader.position" triangles="@subsetROI.trianglesInROI" />
            <TriangleSetTopologyModifier />
            <TriangleSetGeometryAlgorithms name="GeomAlgo" />
            <SubsetTopologicalMapping input="@Container1" output="@Container2" samePoints="true" handleTriangles="true" printLog="1" />
            <MechanicalObject />
            <IdentityMapping />
            <TriangularFEMForceFieldOptim name="FEM1" youngModulus="1000" poissonRatio="0.3" />
            <FastTriangularBendingSprings name="Bending" bendingStiffness="100" />
            <!--<TriangleCollisionModel color="1 0 0 1" group="1" printLog="1" />-->
        </Node>
        <Node name="SubsetTopology2">
            <BoxROI name="subsetROI" position="@../loader.position" triangles="@../loader.triangles" computeEdges="0" computeTetrahedra="0" box="-55 -55 -55 55 55 25" />
            <TriangleSetTopologyContainer name="Container2" position="@../loader.position" triangles="@subsetROI.trianglesInROI" />
            <TriangleSetTopologyModifier />
            <TriangleSetGeometryAlgorithms name="GeomAlgo" />
            <SubsetTopologicalMapping input="@Container1" output="@Container2" samePoints="true" handleTriangles="true" printLog="1" />
            <MechanicalObject />
            <IdentityMapping />
            <TriangularFEMForceFieldOptim name="FEM2" youngModulus="10000" poissonRatio="0.4" />
            <FastTriangularBendingSprings name="Bending" bendingStiffness="1000" />
            <!--<TriangleCollisionModel group="1" printLog="1" />-->
        </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.Engine.Select")
   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.MechanicalLoad")
   root.addObject('RequiredPlugin', name="Sofa.Component.ODESolver.Backward")
   root.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.FEM.Elastic")
   root.addObject('RequiredPlugin', name="Sofa.Component.SolidMechanics.Spring")
   root.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
   root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Dynamic")
   root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Mapping")
   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', )

   full_topology = root.addChild('FullTopology')

   full_topology.addObject('MeshOBJLoader', name="loader", filename="mesh/sphere_05.obj")
   full_topology.addObject('EulerImplicitSolver', name="cg_odesolver", printLog="false", rayleighStiffness="0.1", rayleighMass="0.1")
   full_topology.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
   full_topology.addObject('TriangleSetTopologyContainer', name="Container1", src="@loader")
   full_topology.addObject('TriangleSetTopologyModifier', )
   full_topology.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo")
   full_topology.addObject('MechanicalObject', name="dofs")
   full_topology.addObject('BoxROI', name="roi1", box="-55 -55 -55 55 55 -25")
   full_topology.addObject('FixedProjectiveConstraint', name="FixedProjectiveConstraint", indices="@roi1.indices")
   full_topology.addObject('DiagonalMass', massDensity="1", name="computed using mass density")
   full_topology.addObject('TriangleCollisionModel', group="1", printLog="1")
   full_topology.addObject('TaitSurfacePressureForceField', name="Pressure", gamma="5", B="10000", injectedVolume="0", printLog="1")

   subset_topology1 = FullTopology.addChild('SubsetTopology1')

   subset_topology1.addObject('BoxROI', name="subsetROI", position="@../loader.position", triangles="@../loader.triangles", computeEdges="0", computeTetrahedra="0", box="-55 -55 25 55 55 55")
   subset_topology1.addObject('TriangleSetTopologyContainer', name="Container2", position="@../loader.position", triangles="@subsetROI.trianglesInROI")
   subset_topology1.addObject('TriangleSetTopologyModifier', )
   subset_topology1.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo")
   subset_topology1.addObject('SubsetTopologicalMapping', input="@Container1", output="@Container2", samePoints="true", handleTriangles="true", printLog="1")
   subset_topology1.addObject('MechanicalObject', )
   subset_topology1.addObject('IdentityMapping', )
   subset_topology1.addObject('TriangularFEMForceFieldOptim', name="FEM1", youngModulus="1000", poissonRatio="0.3")
   subset_topology1.addObject('FastTriangularBendingSprings', name="Bending", bendingStiffness="100")

   subset_topology2 = FullTopology.addChild('SubsetTopology2')

   subset_topology2.addObject('BoxROI', name="subsetROI", position="@../loader.position", triangles="@../loader.triangles", computeEdges="0", computeTetrahedra="0", box="-55 -55 -55 55 55 25")
   subset_topology2.addObject('TriangleSetTopologyContainer', name="Container2", position="@../loader.position", triangles="@subsetROI.trianglesInROI")
   subset_topology2.addObject('TriangleSetTopologyModifier', )
   subset_topology2.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo")
   subset_topology2.addObject('SubsetTopologicalMapping', input="@Container1", output="@Container2", samePoints="true", handleTriangles="true", printLog="1")
   subset_topology2.addObject('MechanicalObject', )
   subset_topology2.addObject('IdentityMapping', )
   subset_topology2.addObject('TriangularFEMForceFieldOptim', name="FEM2", youngModulus="10000", poissonRatio="0.4")
   subset_topology2.addObject('FastTriangularBendingSprings', name="Bending", bendingStiffness="1000")

SubsetTopologicalMapping.scn

XMLPython

<?xml version="1.0" ?>
<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.Engine.Select"/> <!-- Needed to use components [BoxROI] -->
    <RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshOBJLoader] -->
    <RequiredPlugin name="Sofa.Component.LinearSolver.Iterative"/> <!-- Needed to use components [CGLinearSolver] -->
    <RequiredPlugin name="Sofa.Component.Mapping.Linear"/> <!-- Needed to use components [IdentityMapping] -->
    <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 [TriangularFEMForceField] -->
    <RequiredPlugin name="Sofa.Component.StateContainer"/> <!-- Needed to use components [MechanicalObject] -->
    <RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TriangleSetGeometryAlgorithms TriangleSetTopologyContainer TriangleSetTopologyModifier] -->
    <RequiredPlugin name="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [SubsetTopologicalMapping] -->
    <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="FullTopology">
        <MeshOBJLoader name="loader" filename="mesh/sphere_05.obj" />
        <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" />
        <TriangleSetTopologyContainer name="Container1" src="@loader" />
        <TriangleSetTopologyModifier />
        <TriangleSetGeometryAlgorithms name="GeomAlgo" />
        <MechanicalObject name="dofs" />
        <BoxROI name="roi1" box="-55 -55 -55 55 55 -25" />
        <FixedProjectiveConstraint name="FixedProjectiveConstraint" indices="@roi1.indices" />
        <DiagonalMass massDensity="1" name="computed using mass density" />
        <TriangularFEMForceField name="FEM" youngModulus="3000" poissonRatio="0.4" />
        <Node name="SubsetTopology1">
            <BoxROI name="subsetROI" position="@../loader.position" triangles="@../loader.triangles" computeEdges="0" computeTetrahedra="0" box="-55 -55 25 55 55 55" />
            <TriangleSetTopologyContainer name="Container2" position="@../loader.position" triangles="@subsetROI.trianglesInROI" />
            <TriangleSetTopologyModifier />
            <TriangleSetGeometryAlgorithms name="GeomAlgo" />
            <SubsetTopologicalMapping input="@Container1" output="@Container2" samePoints="true" />
            <MechanicalObject />
            <IdentityMapping />
            <TriangleCollisionModel color="1 0 0 1" group="1" printLog="1" />
        </Node>
        <Node name="SubsetTopology2">
            <BoxROI name="subsetROI" position="@../loader.position" triangles="@../loader.triangles" computeEdges="0" computeTetrahedra="0" box="-55 -55 -55 55 55 25" />
            <TriangleSetTopologyContainer name="Container2" position="@../loader.position" triangles="@subsetROI.trianglesInROI" />
            <TriangleSetTopologyModifier />
            <TriangleSetGeometryAlgorithms name="GeomAlgo" />
            <SubsetTopologicalMapping input="@Container1" output="@Container2" samePoints="true" />
            <MechanicalObject />
            <IdentityMapping />
            <TriangleCollisionModel group="1" printLog="1" />
        </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.Engine.Select")
   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.Topology.Mapping")
   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', )

   full_topology = root.addChild('FullTopology')

   full_topology.addObject('MeshOBJLoader', name="loader", filename="mesh/sphere_05.obj")
   full_topology.addObject('EulerImplicitSolver', name="cg_odesolver", printLog="false", rayleighStiffness="0.1", rayleighMass="0.1")
   full_topology.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
   full_topology.addObject('TriangleSetTopologyContainer', name="Container1", src="@loader")
   full_topology.addObject('TriangleSetTopologyModifier', )
   full_topology.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo")
   full_topology.addObject('MechanicalObject', name="dofs")
   full_topology.addObject('BoxROI', name="roi1", box="-55 -55 -55 55 55 -25")
   full_topology.addObject('FixedProjectiveConstraint', name="FixedProjectiveConstraint", indices="@roi1.indices")
   full_topology.addObject('DiagonalMass', massDensity="1", name="computed using mass density")
   full_topology.addObject('TriangularFEMForceField', name="FEM", youngModulus="3000", poissonRatio="0.4")

   subset_topology1 = FullTopology.addChild('SubsetTopology1')

   subset_topology1.addObject('BoxROI', name="subsetROI", position="@../loader.position", triangles="@../loader.triangles", computeEdges="0", computeTetrahedra="0", box="-55 -55 25 55 55 55")
   subset_topology1.addObject('TriangleSetTopologyContainer', name="Container2", position="@../loader.position", triangles="@subsetROI.trianglesInROI")
   subset_topology1.addObject('TriangleSetTopologyModifier', )
   subset_topology1.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo")
   subset_topology1.addObject('SubsetTopologicalMapping', input="@Container1", output="@Container2", samePoints="true")
   subset_topology1.addObject('MechanicalObject', )
   subset_topology1.addObject('IdentityMapping', )
   subset_topology1.addObject('TriangleCollisionModel', color="1 0 0 1", group="1", printLog="1")

   subset_topology2 = FullTopology.addChild('SubsetTopology2')

   subset_topology2.addObject('BoxROI', name="subsetROI", position="@../loader.position", triangles="@../loader.triangles", computeEdges="0", computeTetrahedra="0", box="-55 -55 -55 55 55 25")
   subset_topology2.addObject('TriangleSetTopologyContainer', name="Container2", position="@../loader.position", triangles="@subsetROI.trianglesInROI")
   subset_topology2.addObject('TriangleSetTopologyModifier', )
   subset_topology2.addObject('TriangleSetGeometryAlgorithms', name="GeomAlgo")
   subset_topology2.addObject('SubsetTopologicalMapping', input="@Container1", output="@Container2", samePoints="true")
   subset_topology2.addObject('MechanicalObject', )
   subset_topology2.addObject('IdentityMapping', )
   subset_topology2.addObject('TriangleCollisionModel', group="1", printLog="1")