<?xml version="1.0" ?>

<Node name="root" gravity="-1.8 0 100"  dt="0.1" >
    <VisualStyle displayFlags="showVisual showForceFields" />
    
    <Node name="Springs">

		<EulerImplicitSolver name="cg_odesolver"/>
		<CGLinearSolver iterations="100"/>
        	<MeshGmshLoader name="meshLoader" filename="truthcylinder1.msh" />
		<MechanicalObject dx="15" src="@meshLoader"/>
		<TetrahedronSetTopologyContainer name="topo" src="@meshLoader" />
		<TetrahedronSetTopologyModifier name="modif" />
		<TetrahedronSetTopologyAlgorithms name="algo" />
		<TetrahedronSetGeometryAlgorithms name="tetraGeo"/>
		<UniformMass totalmass="15" />
		<FixedConstraint indices="0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 268 269 270 271 343 345" />
		<MeshSpringForceField name="Spring" tetrasStiffness="1870" tetrasDamping="0" />
		
    </Node>
    
    <Node name="CoFEM">
		<EulerImplicitSolver name="cg_odesolver"/>
		<CGLinearSolver iterations="100" />
        	<MeshGmshLoader name="meshLoader" filename="truthcylinder1.msh"/>
		<MechanicalObject dx="30" src="@meshLoader"/>
		<TetrahedronSetTopologyContainer name="topo" src="@meshLoader" />
		<TetrahedronSetTopologyModifier name="modif" />
		<TetrahedronSetTopologyAlgorithms name="algo" />
		<TetrahedronSetGeometryAlgorithms name="tetraGeo"/>
		<UniformMass totalmass="15" />
		<FixedConstraint indices="0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 268 269 270 271 343 345" />
		<TetrahedronFEMForceField name="FEM" youngModulus="1116" poissonRatio="0.49" method="polar" />
		
    </Node>

    <Node name="CoFEM_firstOrder">
		<EulerImplicitSolver name="cg_odesolver" firstOrder="1"/>
		<CGLinearSolver iterations="100" />
        	<MeshGmshLoader name="meshLoader" filename="truthcylinder1.msh"/>
		<MechanicalObject dx="45" src="@meshLoader"/>
		<TetrahedronSetTopologyContainer name="topo" src="@meshLoader" />
		<TetrahedronSetTopologyModifier name="modif" />
		<TetrahedronSetTopologyAlgorithms name="algo" />
		<TetrahedronSetGeometryAlgorithms name="tetraGeo"/>
		<UniformMass totalmass="15" />
		<FixedConstraint indices="0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 268 269 270 271 343 345" />
		<TetrahedronFEMForceField name="FEM" youngModulus="1116" poissonRatio="0.49" method="polar" />
		
    </Node>
	
    <Node name="LinearFEM">
		<EulerImplicitSolver name="cg_odesolver" />
		<CGLinearSolver iterations="100" />
        	<MeshGmshLoader name="meshLoader" filename="truthcylinder1.msh"/>
		<MechanicalObject dx="60" src="@meshLoader"/>
		<TetrahedronSetTopologyContainer name="topo" src="@meshLoader" />
		<TetrahedronSetTopologyModifier name="modif" />
		<TetrahedronSetTopologyAlgorithms name="algo" />
		<TetrahedronSetGeometryAlgorithms name="tetraGeo"/>
		<UniformMass totalmass="15" />
		<FixedConstraint indices="0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 268 269 270 271 343 345" />
		<TetrahedronFEMForceField name="FEM" youngModulus="1116" poissonRatio="0.49" method="small" />
		
    </Node>
    
</Node>