About improving sofa run frame rate

[yangwxb]

hello!
I mainly learned sofa through the demo provided by the official. In the process of writing scripts by myself, I found that the more models I loaded, the lower the overall frame rate was. I tried to use CUDA to solve this problem, but many components do not support the data type CudaVec3f, I would like to ask what method can be used to improve the frame rate? Dear developers, please help to see if there is a problem with the script I wrote.I will upload the model if needed.
Uploading scene.zip…

[hugtalbot]

Dear @yangwxb

Have you checked the related documentation page: https://www.sofa-framework.org/community/doc/using-sofa/performances/improve-performances/ ?

[yangwxb]

<RequiredPlugin name="Sofa.Component.Collision.Detection.Algorithm"/> <!-- Needed to use components [BVHNarrowPhase, BruteForceBroadPhase, DefaultPipeline] -->  
<RequiredPlugin name="Sofa.Component.Collision.Detection.Intersection"/> <!-- Needed to use components [MinProximityIntersection] -->  
<RequiredPlugin name="Sofa.Component.Collision.Geometry"/> <!-- Needed to use components [PointCollisionModel, SphereCollisionModel, TriangleCollisionModel] -->  
<RequiredPlugin name="Sofa.Component.Collision.Response.Contact"/> <!-- Needed to use components [DefaultContactManager] -->  
<RequiredPlugin name="Sofa.Component.Constraint.Projective"/> <!-- Needed to use components [FixedConstraint] -->  
<RequiredPlugin name="Sofa.Component.Engine.Select"/> <!-- Needed to use components [BoxROI] -->  
<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.Mapping.NonLinear"/> <!-- Needed to use components [RigidMapping] -->


<RequiredPlugin name="Sofa.Component.Mass"/> <!-- Needed to use components [DiagonalMass, UniformMass] -->  
<RequiredPlugin name="Sofa.Component.ODESolver.Backward"/> <!-- Needed to use components [EulerImplicitSolver] -->  
<RequiredPlugin name="Sofa.Component.SceneUtility"/> <!-- Needed to use components [InfoComponent] -->  
<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.Constant"/> <!-- Needed to use components [MeshTopology] -->

<RequiredPlugin name="Sofa.Component.Topology.Container.Dynamic"/> <!-- Needed to use components [TetrahedronSetGeometryAlgorithms, TetrahedronSetTopologyContainer, TetrahedronSetTopologyModifier, TriangleSetGeometryAlgorithms, TriangleSetTopologyContainer, TriangleSetTopologyModifier] -->  
<RequiredPlugin name="Sofa.Component.Topology.Container.Grid"/> <!-- Needed to use components [RegularGridTopology] -->  
<RequiredPlugin name="Sofa.Component.Topology.Mapping"/> <!-- Needed to use components [Hexa2TetraTopologicalMapping, Tetra2TriangleTopologicalMapping] -->  
<RequiredPlugin name="Sofa.Component.Visual"/> <!-- Needed to use components [VisualStyle] -->  
<RequiredPlugin name="Sofa.GL.Component.Rendering3D"/> <!-- Needed to use components [OglModel] --> 
<RequiredPlugin name="SofaCarving"/>
<RequiredPlugin name="Sofa.Component.Constraint.Lagrangian.Solver"/> <!-- Needed to use components [LCPConstraintSolver] -->
<RequiredPlugin name="Sofa.Component.SolidMechanics.Spring"/> <!-- Needed to use components [RegularGridSpringForceField] -->
<RequiredPlugin name="Sofa.Component.IO.Mesh"/> <!-- Needed to use components [MeshOBJLoader] --> 
<RequiredPlugin name="Sofa.Component.Controller"/>
<RequiredPlugin name="Sofa.Component.Topology.Utility"/> <!-- Needed to use components [TopologicalChangeProcessor] -->
<RequiredPlugin name="Sofa.Component.Constraint.Lagrangian.Correction"/> <!-- Needed to use components [UncoupledConstraintCorrection] -->
<RequiredPlugin name="Sofa.Component.MechanicalLoad"/> <!-- Needed to use components [ConstantForceField,PlaneForceField] -->  
<RequiredPlugin name="SofaCUDA"/>
<RequiredPlugin name="Sofa.Component.AnimationLoop"/> <!-- Needed to use components [FreeMotionAnimationLoop] -->

<VisualStyle displayFlags="showVisual" />
<DefaultVisualManagerLoop />
<DefaultAnimationLoop/>
<DefaultPipeline verbose="0" depth="15" draw="0" />
<!-- <FreeMotionAnimationLoop/> -->
<BruteForceBroadPhase/>
<BVHNarrowPhase/>

<DefaultContactManager response="PenalityContactForceField" />
<!-- <DefaultContactManager response="StickContactForceField" /> -->

<MinProximityIntersection name="Proximity" alarmDistance="0.0003" contactDistance="0.0002" />

<LCPConstraintSolver tolerance="0.001" maxIt="1000"/>




<Node name="nangmo">
    
    <EulerImplicitSolver name="EulerImplicit"  rayleighStiffness="0.05" rayleighMass="0.0515" />
    <CGLinearSolver name="CG Solver3" iterations="15" tolerance="1e-9" threshold="1e-25"/>
    <MeshGmshLoader name="loader" filename="3D_YuanTeng_1500.msh" scale3d="1 1 1" rotation="-90 0 0" />
    <MechanicalObject src="@loader" name="Volume" translation ="0 0.001 0" template="CudaVec3f"/>
    
    <TetrahedronSetTopologyContainer  name="Tetra_topo" position="@Volume.position" src="@loader"/>
    <TetrahedronSetTopologyModifier   name="Modifier" />
    <TetrahedronSetGeometryAlgorithms name="GeomAlgo" template="CudaVec3f" drawTetrahedra="0" drawScaleTetrahedra="0.8"/>
    
   
    <!-- <TetrahedralCorotationalFEMForceField name="CFEM" youngModulus="20000" poissonRatio="0.38" method="large" /> -->
    <TetrahedronFEMForceField name="CFEM" youngModulus="20000" poissonRatio="0.38" method="large" />
    <ConstantForceField indices="0 1" forces="0 0 0  0 0 0" showArrowSize="0.5" printLog="1"/>
    <BoxROI name="box" box="-0.005  0.0015 -0.005 0.005 0.0018 0.005" drawBoxes="true" />
    <FixedConstraint  name="FixedConstraint" indices="@box.indices" activate_projectVelocity="1"/>
	<!-- <UncoupledConstraintCorrection useOdeSolverIntegrationFactors="0"/> -->
	
    <Node name="Surface">
       
        <TriangleSetTopologyContainer  name="Container"  />
        <TriangleSetTopologyModifier   name="Modifier"/>
        <TriangleSetGeometryAlgorithms name="GeomAlgo" template="CudaVec3f" drawTriangles="0" />
        <Tetra2TriangleTopologicalMapping input="@../Tetra_topo" output="@Container" />
        
		
		<!-- <TriangularFEMForceField name="FEM" youngModulus="20000" poissonRatio="0.38" method="large" template="CudaVec3f"/> -->
        <TriangularBendingSprings name="FEM-Bend" stiffness="1000" damping="80"  topology="@Container" template="CudaVec3f"/>
        
        
        <UniformMass totalMass="0.001"/>
        
        <TriangleCollisionModel name="triangleCol" bothSide="true" group="2" contactStiffness="0.05"/>
        <LineCollisionModel name='edgeCol' group="2" contactStiffness="0.05" />
        <PointCollisionModel name="pointCol" group="2" contactStiffness="0.05"/>
        
        <Node name="Visua1">
            <OglModel name="Visualmode1" color = "red"/>
            <IdentityMapping input="@../../Volume" output="@Visualmode1" />
        </Node>
    </Node>
</Node>

[yangwxb]

dear @hugtalbot
I modified my document according to the writing method of the example, but it is invalid, may I ask what is wrong with my writing method?

[hugtalbot]

I suggest we close this topic since you have several topics opened regarding SofaCUDA, correct @yangwxb ?

[yangwxb]

Yes, I'm sorry, but I can't find the setting to turn off this topic.

[hugtalbot]

Done 👍