Bugs in SofaPython3

[Hoangn06]

Hello,

I've testing the Sofapython3 by launching python examples in the SOFA Documentation. I've launched lots of examples but the simulations always couldn't show some objects I imported in or sometimes couldn't show any objects. The simulations couldn't detect the AnimationLoop even though the AnimationLoop has been called inside the simulation.

Here is a few example for the bug.

1. Python script of LinearForceField provided in https://sofa-framework.github.io/doc/components/mechanicalload/linearforcefield/#__tabbed_1_2
   In this example, there should be 2 objects: TorusRigidX and TorusRigidY but SOFA showed only the TorusRigidX and more over, SOFA can't detect AnimationLoop even though it's declared inside the script. And I received also
   
   

2. Python script for Gravity provided in https://sofa-framework.github.io/doc/components/mechanicalload/gravity/
   In this example, there should be 4 objects but the simulation couldn't detect any of the object. Also couldn't detect the AnimationLoop.
   
   

Can someone please help me fix this error? I've launched the src scripts and it worked perfectly fine but the python Scripts didn't work at all.

Thank you and have a great day.

[hugtalbot]

Hi @Hoangn06

1. You refer to the scene example associated with the LinearForceField while the error refers to the script 15-python3-addGUI.py. I am confused. In your script, your error shows a problem with something named TorusRigidX which is not defined.
2. This error is explicit : unterminated string literal at line 88

Cheers

[Hoangn06]

Hello @hugtalbot,

Sorry for the confusion. So in order to test scripts fast, I created a copy of the script 15-python-3-addGUI.py and copy the codes of another script inside. The script I launch is actually the LinearForceField script:
"
def createScene(root_node):

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

root.addObject('RequiredPlugin', name="Sofa.Component.IO.Mesh")
root.addObject('RequiredPlugin', name="Sofa.Component.LinearSolver.Iterative")
root.addObject('RequiredPlugin', name="Sofa.Component.Mapping.NonLinear")
root.addObject('RequiredPlugin', name="Sofa.Component.Mass")
root.addObject('RequiredPlugin', name="Sofa.Component.MechanicalLoad")
root.addObject('RequiredPlugin', name="Sofa.Component.ODESolver.Forward")
root.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
root.addObject('RequiredPlugin', name="Sofa.Component.Visual")
root.addObject('RequiredPlugin', name="Sofa.GL.Component.Rendering3D")
root.addObject('VisualStyle', displayFlags="showBehaviorModels showForceFields")
root.addObject('DefaultAnimationLoop', )

torus_rigid_x = root.addChild('TorusRigidX')

torus_rigid_x.addObject('EulerExplicitSolver', )
torus_rigid_x.addObject('CGLinearSolver', iterations="25", threshold="0.00000001", tolerance="1e-5")
torus_rigid_x.addObject('MechanicalObject', template="Rigid3", dx="2", dy="0", dz="0", rx="0", ry="0", rz="0", scale="1.0")
torus_rigid_x.addObject('UniformMass', totalMass="1.0")
torus_rigid_x.addObject('LinearForceField', points="0", forces="0 0 0 0 0 0 1 0 0 0 0 0 -1 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 0 0 0", force="2.0", times="0 4 8 10 12")

visu = TorusRigidX.addChild('Visu')

visu.addObject('MeshOBJLoader', name="meshLoader_0", filename="mesh/torus.obj", scale="0.3", handleSeams="1")
visu.addObject('OglModel', name="Visual", src="@meshLoader_0", color="gray")
visu.addObject('RigidMapping', input="@..", output="@visual")

torus_rigid_y = root.addChild('TorusRigidY')

torus_rigid_y.addObject('EulerExplicitSolver', )
torus_rigid_y.addObject('CGLinearSolver', iterations="25", threshold="0.00000001", tolerance="1e-5")
torus_rigid_y.addObject('MechanicalObject', template="Rigid3", dx="2", dy="2", dz="0", rx="0", ry="0", rz="0", scale="1.0")
torus_rigid_y.addObject('UniformMass', totalMass="1.0")
torus_rigid_y.addObject('LinearForceField', points="0", forces="0 0 0 0 0 0 0 1 0 0 0 0 0 -1 0 0 0 0 0 -1 0 0 0 0 0 0 0 0 0 0", force="2.0", times="0 4 8 10 12")

visu = TorusRigidY.addChild('Visu')

visu.addObject('MeshOBJLoader', name="meshLoader_1", filename="mesh/torus.obj", scale="0.3", handleSeams="1")
visu.addObject('OglModel', name="Visual", src="@meshLoader_1", color="gray")
visu.addObject('RigidMapping', input="@..", output="@visual")
"
I dont know what happened and why it couldn't detect/define the TorusRigidX.

Same problem with the script of Gravity. Even I delete the floor part and have the script like this:

def createScene(root_node):

root = root_node.addChild('root', dt="0.03333")

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.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.Spring")
root.addObject('RequiredPlugin', name="Sofa.Component.StateContainer")
root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Constant")
root.addObject('RequiredPlugin', name="Sofa.Component.Topology.Container.Grid")
root.addObject('RequiredPlugin', name="Sofa.GL.Component.Rendering3D")
root.addObject('RequiredPlugin', name="Sofa.Component.MechanicalLoad")
root.addObject('CollisionPipeline', verbose="0", depth="10", draw="0")
root.addObject('BruteForceBroadPhase', )
root.addObject('BVHNarrowPhase', )
root.addObject('MinProximityIntersection', name="Proximity", alarmDistance="0.75", contactDistance="0.5")
root.addObject('CollisionResponse', name="Response", response="PenalityContactForceField")
root.addObject('DefaultAnimationLoop', )

torus1 = root.addChild('Torus1')

torus1.addObject('Gravity', gravity="0 -10 0")
torus1.addObject('EulerImplicitSolver', name="cg_odesolver", printLog="false", rayleighStiffness="0.1", rayleighMass="0.1")
torus1.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
torus1.addObject('MechanicalObject', dx="0", dy="20", dz="0", ry="90")
torus1.addObject('UniformMass', totalMass="10")
torus1.addObject('RegularGridTopology', nx="6", ny="5", nz="2", xmin="-7.5", xmax="7.5", ymin="-6", ymax="6", zmin="-1.75", zmax="1.75")
torus1.addObject('RegularGridSpringForceField', name="Springs", stiffness="350", damping="1")

visu = Torus1.addChild('Visu')

visu.addObject('MeshOBJLoader', name="meshLoader_0", filename="mesh/torus2_scale3.obj", handleSeams="1")
visu.addObject('OglModel', name="Visual", src="@meshLoader_0", color="blue")
visu.addObject('BarycentricMapping', input="@..", output="@visual")

surf = Torus1.addChild('Surf')

surf.addObject('MeshOBJLoader', name="loader", filename="mesh/torus2_scale3.obj")
surf.addObject('MeshTopology', src="@loader")
surf.addObject('MechanicalObject', src="@loader")
surf.addObject('TriangleCollisionModel', )
surf.addObject('LineCollisionModel', )
surf.addObject('PointCollisionModel', )
surf.addObject('BarycentricMapping', )

torus2 = root.addChild('Torus2')

torus2.addObject('Gravity', gravity="0 10 0")
torus2.addObject('EulerImplicitSolver', name="cg_odesolver", printLog="false")
torus2.addObject('CGLinearSolver', iterations="25", name="linear solver", tolerance="1.0e-9", threshold="1.0e-9")
torus2.addObject('MechanicalObject', dx="0", dy="-20", dz="0")
torus2.addObject('UniformMass', totalMass="10")
torus2.addObject('RegularGridTopology', nx="6", ny="5", nz="2", xmin="-7.5", xmax="7.5", ymin="-6", ymax="6", zmin="-1.75", zmax="1.75")
torus2.addObject('RegularGridSpringForceField', name="Springs", stiffness="350", damping="1")

visu = Torus2.addChild('Visu')

visu.addObject('MeshOBJLoader', name="meshLoader_3", filename="mesh/torus2_scale3.obj", handleSeams="1")
visu.addObject('OglModel', name="Visual", src="@meshLoader_3", color="blue")
visu.addObject('BarycentricMapping', input="@..", output="@visual")

surf = Torus2.addChild('Surf')

surf.addObject('MeshOBJLoader', name="loader", filename="mesh/torus2_scale3.obj")
surf.addObject('MeshTopology', src="@loader")
surf.addObject('MechanicalObject', src="@loader")
surf.addObject('TriangleCollisionModel', )
surf.addObject('LineCollisionModel', )
surf.addObject('PointCollisionModel', )
surf.addObject('BarycentricMapping', )

I still always received a message like this.

Can you please help me @hugtalbot! Thank you in advance

[Hoangn06]

Hello,

Can you help me to solve this problem please?

I can't manage to use Python Script beacause SOFA always can not detect objects in my scene..

Thank you and have a nice day

[hugtalbot]

hey @Hoangn06

you are mixing 'Torus1' as name of the node and the name of the python variable torus1
See :

torus1 = root.addChild('Torus1')

Do not hesitate to use LeChat or ChatGPT it would have guided you straightaway

[Hoangn06]

@hugtalbot Yes it's working now. Thank you!

Can ChatGPT solve the problem in the SOFA script? I

[hugtalbot]

Yes it can!