yade-dev team mailing list archive
-
yade-dev team
-
Mailing list archive
-
Message #12484
[Branch ~yade-pkg/yade/git-trunk] Rev 3761: Examples for creating pfacets
------------------------------------------------------------
revno: 3761
committer: Klaus Thoeni <klaus.thoeni@xxxxxxxxx>
timestamp: Tue 2015-12-15 23:59:52 +1100
message:
Examples for creating pfacets
added:
examples/pfacet/
examples/pfacet/box.gts
examples/pfacet/gts-pfacet.py
examples/pfacet/mesh-pfacet.py
examples/pfacet/octahedron.gts
examples/pfacet/pfacetcreators.py
examples/pfacet/sphere.gts
examples/pfacet/sphere.mesh
--
lp:yade
https://code.launchpad.net/~yade-pkg/yade/git-trunk
Your team Yade developers is subscribed to branch lp:yade.
To unsubscribe from this branch go to https://code.launchpad.net/~yade-pkg/yade/git-trunk/+edit-subscription
=== added directory 'examples/pfacet'
=== added file 'examples/pfacet/box.gts'
--- examples/pfacet/box.gts 1970-01-01 00:00:00 +0000
+++ examples/pfacet/box.gts 2015-12-15 12:59:52 +0000
@@ -0,0 +1,39 @@
+8 18 12 GtsSurface GtsFace GtsEdge GtsVertex
+-0.5 0.5 -0.5
+-0.5 0.5 0.5
+-0.5 -0.5 0.5
+0.5 -0.5 -0.5
+0.5 0.5 -0.5
+0.5 0.5 0.5
+0.5 -0.5 0.5
+-0.5 -0.5 -0.5
+1 2
+3 1
+2 3
+4 5
+1 4
+5 1
+5 6
+5 2
+6 2
+7 5
+7 4
+4 8
+4 3
+3 8
+7 6
+3 7
+8 1
+2 7
+1 2 3
+4 5 6
+7 8 9
+4 10 11
+12 13 14
+15 10 7
+16 13 11
+17 5 12
+1 8 6
+16 18 3
+14 2 17
+9 18 15
=== added file 'examples/pfacet/gts-pfacet.py'
--- examples/pfacet/gts-pfacet.py 1970-01-01 00:00:00 +0000
+++ examples/pfacet/gts-pfacet.py 2015-12-15 12:59:52 +0000
@@ -0,0 +1,88 @@
+from yade import qt
+from yade.gridpfacet import *
+import gts, os.path, locale
+
+locale.setlocale(locale.LC_ALL, 'en_US.UTF-8') # Note: gts is locale-dependend. If, for example, german locale is used, gts.read()-function does not import floats normally
+
+'''
+if you get "Error: unsupported locale setting"
+-> type as root: "dpkg-reconfigure locales"
+-> choose "en_US.UTF-8" (press space to choose)
+'''
+
+################
+### ENGINES ###
+################
+
+O.engines=[
+ ForceResetter(),
+ InsertionSortCollider([
+ Bo1_Wall_Aabb(),
+ Bo1_PFacet_Aabb(),
+ ],sortThenCollide=True),
+ InteractionLoop(
+ [
+ Ig2_GridNode_GridNode_GridNodeGeom6D(),
+ Ig2_Wall_PFacet_ScGeom(),Ig2_Wall_Sphere_ScGeom()
+ ],
+ [
+ Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True),
+ Ip2_FrictMat_FrictMat_FrictPhys()],
+ [
+ Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),
+ Law2_ScGeom_FrictPhys_CundallStrack(),
+ Law2_ScGridCoGeom_FrictPhys_CundallStrack(),
+ Law2_GridCoGridCoGeom_FrictPhys_CundallStrack()
+ ]),
+ GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.5,label='ts'),
+ NewtonIntegrator(gravity=(0,-9.81,0),damping=0.1,label='newton')
+]
+
+
+#################
+### MATERIAL ###
+#################
+
+O.materials.append(CohFrictMat(young=1e8,poisson=0.3,density=2650,frictionAngle=radians(20),normalCohesion=1e100,shearCohesion=1e100,momentRotationLaw=True,label='gridNodeMat'))
+O.materials.append(FrictMat(young=1e8,poisson=0.3,density=2650,frictionAngle=radians(20),label='pFacetMat'))
+
+
+###################
+### IMPORT MESH ###
+###################
+radius=1e-02
+wire=False
+fixed=False
+
+z=-1.2
+color=[0,0,1]
+
+nodesIds0,cylIds0,pfIds0 = gtsPFacet('octahedron.gts',shift=(0,0,0),scale=1.,radius=radius,wire=wire,fixed=fixed,materialNodes='gridNodeMat',material='pFacetMat',color=color)
+
+nodesIds1,cylIds1,pfIds1 = gtsPFacet('box.gts',shift=(3.,0.,0.),scale=2,radius=radius,wire=wire,fixed=fixed,materialNodes='gridNodeMat',material='pFacetMat',color=color)
+
+nodesIds2,cylIds2,pfIds2 = gtsPFacet('sphere.gts',shift=(6.,0.,0.),scale=1.0,radius=radius,wire=wire,fixed=fixed,materialNodes='gridNodeMat',material='pFacetMat',color=color)
+
+
+
+#####################
+##### Wall ###
+#####################
+
+O.bodies.append(utils.wall(position=z,sense=0, axis=1,color=Vector3(1,0,0),material='pFacetMat'))
+
+
+##########
+## VIEW ##
+##########
+
+qt.Controller()
+qtv = qt.View()
+qtr = qt.Renderer()
+qtr.light2=True
+qtr.lightPos=Vector3(1200,1500,500)
+qtr.bgColor=[1,1,1]
+qtv.ortho=True
+
+
+O.saveTmp()
=== added file 'examples/pfacet/mesh-pfacet.py'
--- examples/pfacet/mesh-pfacet.py 1970-01-01 00:00:00 +0000
+++ examples/pfacet/mesh-pfacet.py 2015-12-15 12:59:52 +0000
@@ -0,0 +1,104 @@
+# -*- coding: utf-8
+from yade import qt,plot
+from yade.gridpfacet import *
+
+##################
+### PARAMETERS ###
+##################
+phi=20.
+E=3.*1e8
+color=[255./255.,102./255.,0./255.]
+r=0.005
+
+# position of imported mesh
+xpafet=0.22
+ypafet=0.05
+
+
+################
+### ENGINES ###
+################
+O.engines=[
+ ForceResetter(),
+ InsertionSortCollider([
+ Bo1_PFacet_Aabb(),
+ ]),
+ InteractionLoop([
+ Ig2_GridNode_GridNode_GridNodeGeom6D(),
+ Ig2_GridConnection_GridConnection_GridCoGridCoGeom(),
+ Ig2_Sphere_PFacet_ScGridCoGeom(),
+ Ig2_PFacet_PFacet_ScGeom(),
+ ],
+ [
+ Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True),
+ Ip2_FrictMat_FrictMat_FrictPhys()
+ ],
+ [
+ Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),
+ Law2_ScGeom_FrictPhys_CundallStrack(),
+ Law2_ScGridCoGeom_FrictPhys_CundallStrack(),
+ Law2_GridCoGridCoGeom_FrictPhys_CundallStrack()
+ ]
+ ),
+ GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.5,label='ts'),
+ NewtonIntegrator(gravity=(0,-9.81,0),damping=0.0,label='newton'),
+ PyRunner(iterPeriod=200,command='history()'),
+]
+
+
+################
+### MATERIAL ###
+################
+O.materials.append( CohFrictMat( young=E,poisson=0.3,density=2650,frictionAngle=radians(phi),normalCohesion=3e100,shearCohesion=3e100,momentRotationLaw=True,label='gridNodeMat' ) ) # material to create the gridConnections
+O.materials.append( FrictMat( young=E,poisson=0.3,density=2650,frictionAngle=radians(phi),label='pFacetMat' ) ) # material for general interactions
+
+
+########################################
+### GENERATE THE WALL AND THE SPHERE ###
+########################################
+
+# FIXED WALL
+O.bodies.append( gridNode([-0.25,-0.22,-.25],r,wire=False,fixed=True,material='gridNodeMat',color=color) )
+O.bodies.append( gridNode([.5,0.,-.25],r,wire=False,fixed=True,material='gridNodeMat',color=color) )
+O.bodies.append( gridNode([-.25,-0.22,.25],r,wire=False,fixed=True,material='gridNodeMat',color=color) )
+O.bodies.append( gridNode([.5,0.,.25],r,wire=False,fixed=True,material='gridNodeMat',color=color) )
+O.bodies.append( gridConnection(0,1,r,color=color,material='gridNodeMat') )
+O.bodies.append( gridConnection(2,3,r,color=color,material='gridNodeMat') )
+O.bodies.append( gridConnection(2,1,r,color=color,material='gridNodeMat') )
+O.bodies.append( gridConnection(2,0,r,color=color,material='gridNodeMat') )
+O.bodies.append( gridConnection(3,1,r,color=color,material='gridNodeMat') )
+O.bodies.append( pfacet(2,1,0,wire=False,material='pFacetMat',color=color) )
+O.bodies.append( pfacet(2,3,1,wire=False,material='pFacetMat',color=color) )
+
+# IMPORT MESH
+oriBody = Quaternion(Vector3(0,0,1),pi/2.)
+nodesIds,cylIds,pfIds = gmshPFacet( meshfile='sphere.mesh', shift=Vector3(xpafet,ypafet,0.), orientation=oriBody, radius=r, wire=False, fixed=False, materialNodes='gridNodeMat', material='pFacetMat', color=[1,0,0] )
+
+
+############
+### PLOT ###
+############
+def history():
+ xyz=[]
+ for k in [0,1,2]:
+ ksum=0
+ for i in nodesIds:
+ ksum+=O.bodies[i].state.pos[k]
+ xyz.append(ksum/len(nodesIds)) # take average value as reference
+ plot.addData(i=O.iter,t=O.time,x=xyz[0],y=xyz[1],z=xyz[2])
+plot.plots={'x':'y'}
+plot.plot()
+
+
+##########
+## VIEW ##
+##########
+
+qt.Controller()
+qtv = qt.View()
+qtr = qt.Renderer()
+qtr.light2=True
+qtr.lightPos=Vector3(1200,1500,500)
+qtr.bgColor=[1,1,1]
+
+O.saveTmp()
=== added file 'examples/pfacet/octahedron.gts'
--- examples/pfacet/octahedron.gts 1970-01-01 00:00:00 +0000
+++ examples/pfacet/octahedron.gts 2015-12-15 12:59:52 +0000
@@ -0,0 +1,1539 @@
+258 768 512 GtsSurface GtsFace GtsEdge GtsVertex
+0 0.625 0.625
+0.1909179688 0.517578125 0.6655273438
+0 0.4765625 0.7685546875
+-0.3125 -0.890625 0
+-0.1826171875 -0.91015625 0.1826171875
+-0.3530273438 -0.8041992188 0.18359375
+0.1826171875 0.91015625 0.1826171875
+0 0.890625 0.3125
+0.18359375 0.8041992188 0.3530273438
+0 0.1494140625 0.97265625
+-0.1826171875 0.1826171875 0.91015625
+0 0.3125 0.890625
+-0.3671875 -0.3671875 0.671875
+-0.18359375 -0.3530273438 0.8041992188
+-0.3530273438 -0.18359375 0.8041992188
+0.625 0 0.625
+0.7685546875 0 0.4765625
+0.6655273438 -0.1909179688 0.517578125
+-0.97265625 0 -0.1494140625
+-1 0 0
+-0.97265625 -0.1494140625 0
+0 -0.4765625 -0.7685546875
+-0.1909179688 -0.517578125 -0.6655273438
+0 -0.625 -0.625
+0 0.97265625 0.1494140625
+-0.1494140625 0.97265625 0
+0 1 0
+-0.18359375 0.8041992188 0.3530273438
+-0.1826171875 0.91015625 0.1826171875
+-0.373046875 0.529296875 -0.529296875
+-0.3671875 0.3671875 -0.671875
+-0.529296875 0.373046875 -0.529296875
+-0.890625 -0.3125 0
+-0.91015625 -0.1826171875 0.1826171875
+0.625 0.625 0
+0.517578125 0.6655273438 0.1909179688
+0.6655273438 0.517578125 0.1909179688
+-0.8041992188 -0.18359375 -0.3530273438
+-0.671875 -0.3671875 -0.3671875
+-0.8041992188 -0.3530273438 -0.18359375
+0.1909179688 0.6655273438 -0.517578125
+0.373046875 0.529296875 -0.529296875
+0.3671875 0.671875 -0.3671875
+0.1909179688 -0.517578125 -0.6655273438
+0.373046875 -0.529296875 -0.529296875
+0.3671875 -0.3671875 -0.671875
+0 -0.97265625 0.1494140625
+0 -0.890625 0.3125
+-0.373046875 -0.529296875 0.529296875
+-0.529296875 -0.373046875 0.529296875
+0.517578125 -0.1909179688 -0.6655273438
+0.625 0 -0.625
+0.4765625 0 -0.7685546875
+0.97265625 0 -0.1494140625
+1 0 0
+0.97265625 0.1494140625 0
+0.517578125 0.6655273438 -0.1909179688
+0.3530273438 0.8041992188 -0.18359375
+0.4765625 0.7685546875 0
+-0.1909179688 0.517578125 0.6655273438
+-0.1909179688 0.6655273438 0.517578125
+-0.4765625 -0.7685546875 0
+-0.517578125 -0.6655273438 0.1909179688
+0 0.7685546875 0.4765625
+0.1909179688 0.6655273438 0.517578125
+0.3671875 0.3671875 -0.671875
+0.1909179688 0.517578125 -0.6655273438
+0.373046875 0.529296875 0.529296875
+0.3671875 0.3671875 0.671875
+0 0.890625 -0.3125
+0.1826171875 0.91015625 -0.1826171875
+0 0.97265625 -0.1494140625
+0.91015625 0.1826171875 0.1826171875
+0.97265625 0 0.1494140625
+0.890625 0 0.3125
+-0.18359375 -0.3530273438 -0.8041992188
+-0.3530273438 -0.18359375 -0.8041992188
+-0.3671875 -0.3671875 -0.671875
+-0.8041992188 0.18359375 -0.3530273438
+-0.91015625 0.1826171875 -0.1826171875
+-0.890625 0 -0.3125
+-0.4765625 0 0.7685546875
+-0.3530273438 0.18359375 0.8041992188
+-0.517578125 0.1909179688 0.6655273438
+-0.890625 0.3125 0
+-0.8041992188 0.3530273438 -0.18359375
+0.1494140625 0.97265625 0
+0.529296875 -0.373046875 0.529296875
+0.671875 -0.3671875 0.3671875
+0.8041992188 0.3530273438 -0.18359375
+0.91015625 0.1826171875 -0.1826171875
+0.890625 0.3125 0
+0 -0.3125 0.890625
+-0.1826171875 -0.1826171875 0.91015625
+0 -0.1494140625 0.97265625
+-0.529296875 0.529296875 0.373046875
+-0.529296875 0.373046875 0.529296875
+-0.671875 0.3671875 0.3671875
+-0.18359375 0.3530273438 0.8041992188
+-0.3671875 0.3671875 0.671875
+-0.91015625 -0.1826171875 -0.1826171875
+0.8041992188 -0.18359375 0.3530273438
+-0.529296875 0.529296875 -0.373046875
+-0.671875 0.3671875 -0.3671875
+0.18359375 -0.3530273438 0.8041992188
+0 -0.4765625 0.7685546875
+0.529296875 -0.373046875 -0.529296875
+0.1826171875 -0.91015625 0.1826171875
+0.8041992188 0.18359375 0.3530273438
+0.6655273438 0.1909179688 0.517578125
+-0.18359375 -0.8041992188 0.3530273438
+-0.3671875 -0.671875 0.3671875
+-0.1909179688 -0.6655273438 0.517578125
+-0.1909179688 0.6655273438 -0.517578125
+0 0.625 -0.625
+-0.1909179688 0.517578125 -0.6655273438
+-0.373046875 0.529296875 0.529296875
+-0.1909179688 -0.517578125 0.6655273438
+-0.529296875 -0.529296875 0.373046875
+-0.6655273438 -0.517578125 0.1909179688
+0.6655273438 -0.1909179688 -0.517578125
+0.7685546875 0 -0.4765625
+-0.517578125 0.1909179688 -0.6655273438
+-0.6655273438 0.1909179688 -0.517578125
+0.3530273438 -0.18359375 -0.8041992188
+0.517578125 -0.6655273438 -0.1909179688
+0.529296875 -0.529296875 -0.373046875
+0.6655273438 -0.517578125 -0.1909179688
+0.890625 0 -0.3125
+-0.517578125 -0.1909179688 0.6655273438
+0.529296875 0.373046875 0.529296875
+0.517578125 0.1909179688 0.6655273438
+-0.8041992188 -0.3530273438 0.18359375
+-0.8041992188 -0.18359375 0.3530273438
+-0.18359375 0.3530273438 -0.8041992188
+-0.1826171875 0.1826171875 -0.91015625
+-0.3530273438 0.18359375 -0.8041992188
+0.1826171875 0.1826171875 -0.91015625
+0.18359375 0.3530273438 -0.8041992188
+0.3530273438 0.18359375 -0.8041992188
+-0.1826171875 0.91015625 -0.1826171875
+-0.18359375 0.8041992188 -0.3530273438
+0.8041992188 -0.3530273438 -0.18359375
+0.8041992188 -0.18359375 -0.3530273438
+0.91015625 -0.1826171875 -0.1826171875
+0.3125 0.890625 0
+0.18359375 0.3530273438 0.8041992188
+0.1826171875 0.1826171875 0.91015625
+0.3530273438 0.18359375 0.8041992188
+0 0.1494140625 -0.97265625
+-0.1494140625 0 -0.97265625
+-0.517578125 -0.6655273438 -0.1909179688
+-0.625 -0.625 0
+-0.97265625 0 0.1494140625
+-0.890625 0 0.3125
+0.3671875 0.671875 0.3671875
+0.529296875 0.529296875 0.373046875
+0.529296875 0.529296875 -0.373046875
+-0.4765625 0 -0.7685546875
+0.671875 0.3671875 -0.3671875
+0.6655273438 0.517578125 -0.1909179688
+0 -0.7685546875 0.4765625
+0.671875 -0.3671875 -0.3671875
+0.4765625 -0.7685546875 0
+0.625 -0.625 0
+0.3530273438 -0.18359375 0.8041992188
+0.3671875 -0.3671875 0.671875
+0.3125 0 0.890625
+0.671875 0.3671875 0.3671875
+-0.3671875 0.671875 0.3671875
+0.7685546875 0.4765625 0
+0.8041992188 0.3530273438 0.18359375
+0.1909179688 -0.517578125 0.6655273438
+0.373046875 -0.529296875 0.529296875
+-0.6655273438 -0.517578125 -0.1909179688
+-0.7685546875 -0.4765625 0
+-0.1494140625 0 0.97265625
+0 0 1
+-0.3671875 -0.671875 -0.3671875
+-0.3530273438 -0.8041992188 -0.18359375
+0 -0.625 0.625
+0 -0.1494140625 -0.97265625
+-0.1826171875 -0.1826171875 -0.91015625
+0.3671875 -0.671875 -0.3671875
+-0.7685546875 0 0.4765625
+-0.6655273438 0.1909179688 0.517578125
+-0.8041992188 0.18359375 0.3530273438
+0.3125 -0.890625 0
+0.3530273438 -0.8041992188 -0.18359375
+0.1494140625 0 0.97265625
+-0.517578125 -0.1909179688 -0.6655273438
+-0.529296875 -0.373046875 -0.529296875
+-0.373046875 -0.529296875 -0.529296875
+0.3125 0 -0.890625
+0.8041992188 0.18359375 -0.3530273438
+0.1909179688 -0.6655273438 0.517578125
+0.18359375 -0.8041992188 0.3530273438
+0 0.3125 -0.890625
+0 0.4765625 -0.7685546875
+-0.1826171875 -0.91015625 -0.1826171875
+-0.18359375 -0.8041992188 -0.3530273438
+-0.7685546875 0 -0.4765625
+-0.625 0 -0.625
+-0.3125 0.890625 0
+-0.3530273438 0.8041992188 -0.18359375
+0 0.7685546875 -0.4765625
+-0.97265625 0.1494140625 0
+-0.3125 0 0.890625
+0.1494140625 -0.97265625 0
+0 -1 0
+0.18359375 -0.3530273438 -0.8041992188
+0.6655273438 0.1909179688 -0.517578125
+0.517578125 0.1909179688 -0.6655273438
+-0.6655273438 0.517578125 0.1909179688
+-0.671875 -0.3671875 0.3671875
+-0.3671875 0.671875 -0.3671875
+-0.91015625 0.1826171875 0.1826171875
+-0.6655273438 -0.1909179688 -0.517578125
+0.517578125 -0.6655273438 0.1909179688
+0.1826171875 -0.91015625 -0.1826171875
+-0.529296875 -0.529296875 -0.373046875
+0.529296875 0.373046875 -0.529296875
+0.529296875 -0.529296875 0.373046875
+0.6655273438 -0.517578125 0.1909179688
+0.91015625 -0.1826171875 0.1826171875
+0.7685546875 -0.4765625 0
+-0.8041992188 0.3530273438 0.18359375
+0.3530273438 -0.8041992188 0.18359375
+-0.625 0.625 0
+-0.4765625 0.7685546875 0
+-0.517578125 0.6655273438 -0.1909179688
+0.3530273438 0.8041992188 0.18359375
+-0.3125 0 -0.890625
+0 -0.97265625 -0.1494140625
+-0.1494140625 -0.97265625 0
+0.8041992188 -0.3530273438 0.18359375
+0.97265625 -0.1494140625 0
+-0.517578125 0.6655273438 0.1909179688
+-0.3530273438 0.8041992188 0.18359375
+-0.1909179688 -0.6655273438 -0.517578125
+0.3671875 -0.671875 0.3671875
+0 -0.7685546875 -0.4765625
+0 -0.890625 -0.3125
+0.18359375 -0.8041992188 -0.3530273438
+0.890625 -0.3125 0
+0.4765625 0 0.7685546875
+0.517578125 -0.1909179688 0.6655273438
+0.1909179688 -0.6655273438 -0.517578125
+-0.625 0 0.625
+-0.6655273438 0.517578125 -0.1909179688
+-0.6655273438 -0.1909179688 0.517578125
+0 -0.3125 -0.890625
+0.1826171875 -0.1826171875 -0.91015625
+0.1826171875 -0.1826171875 0.91015625
+0.18359375 0.8041992188 -0.3530273438
+-0.7685546875 0.4765625 0
+0.1494140625 0 -0.97265625
+0 0 -1
+1 2
+1 3
+2 3
+4 5
+6 5
+6 4
+7 8
+8 9
+7 9
+10 11
+12 10
+11 12
+13 14
+15 14
+15 13
+16 17
+18 16
+17 18
+19 20
+21 19
+21 20
+22 23
+24 22
+23 24
+25 26
+27 25
+26 27
+28 8
+8 29
+28 29
+30 31
+31 32
+30 32
+33 21
+34 33
+21 34
+35 36
+37 36
+37 35
+38 39
+39 40
+38 40
+41 42
+41 43
+42 43
+44 45
+44 46
+45 46
+47 5
+48 47
+5 48
+49 13
+13 50
+49 50
+51 52
+53 51
+53 52
+54 55
+56 54
+55 56
+57 58
+58 59
+59 57
+60 1
+1 61
+60 61
+62 6
+6 63
+63 62
+9 64
+64 65
+65 9
+66 67
+42 67
+66 42
+2 68
+2 69
+68 69
+70 71
+70 72
+71 72
+73 74
+75 73
+75 74
+76 77
+76 78
+78 77
+79 80
+79 81
+81 80
+82 83
+83 84
+84 82
+85 86
+80 86
+80 85
+87 27
+87 25
+88 18
+89 88
+89 18
+90 91
+90 92
+92 91
+93 94
+95 93
+94 95
+96 97
+96 98
+98 97
+99 100
+60 99
+100 60
+81 19
+101 81
+19 101
+18 102
+102 89
+32 103
+32 104
+104 103
+93 105
+93 106
+105 106
+46 107
+46 51
+107 51
+48 108
+108 47
+109 110
+110 17
+17 109
+111 112
+113 111
+112 113
+114 115
+115 116
+114 116
+61 117
+117 60
+118 106
+106 14
+14 118
+63 119
+119 120
+120 63
+121 122
+52 121
+52 122
+32 123
+123 124
+124 32
+125 46
+51 125
+67 41
+126 127
+127 128
+128 126
+129 54
+91 129
+54 91
+13 130
+50 130
+131 132
+132 110
+110 131
+11 99
+99 83
+83 11
+133 134
+134 34
+34 133
+135 136
+136 137
+137 135
+138 139
+139 140
+140 138
+141 70
+70 142
+141 142
+143 144
+144 145
+145 143
+146 87
+7 87
+146 7
+147 148
+148 149
+149 147
+50 119
+119 49
+3 147
+3 12
+147 12
+45 107
+107 127
+127 45
+136 150
+150 151
+151 136
+152 153
+62 152
+62 153
+34 154
+155 34
+155 154
+36 156
+157 156
+36 157
+43 158
+43 57
+158 57
+123 137
+137 159
+159 123
+158 160
+161 158
+161 160
+48 111
+162 48
+111 162
+107 163
+121 107
+121 163
+164 165
+126 164
+165 126
+105 166
+105 167
+167 166
+153 63
+120 153
+91 56
+56 92
+148 168
+168 149
+131 157
+131 169
+169 157
+61 170
+117 170
+110 169
+72 87
+72 27
+76 22
+23 76
+92 171
+172 171
+92 172
+167 173
+174 173
+167 174
+40 175
+175 176
+176 40
+25 7
+25 8
+95 177
+178 95
+178 177
+101 21
+33 101
+14 93
+14 94
+179 152
+180 179
+152 180
+61 64
+64 28
+28 61
+106 173
+106 181
+173 181
+176 33
+133 176
+33 133
+151 182
+182 183
+183 151
+40 101
+40 33
+126 184
+127 184
+185 186
+186 187
+187 185
+162 113
+181 162
+113 181
+116 30
+116 31
+188 164
+189 188
+164 189
+148 10
+10 190
+190 148
+191 78
+192 78
+191 192
+78 23
+193 23
+78 193
+125 53
+194 125
+194 53
+129 195
+91 195
+196 162
+162 197
+197 196
+122 129
+195 122
+184 45
+198 139
+198 199
+139 199
+200 180
+200 4
+4 180
+179 201
+180 201
+124 202
+203 124
+203 202
+204 141
+205 141
+205 204
+206 70
+142 206
+9 156
+156 65
+31 123
+137 31
+194 138
+140 194
+154 20
+207 154
+20 207
+73 172
+73 92
+202 81
+38 202
+81 38
+192 193
+147 2
+208 11
+83 208
+51 121
+153 176
+176 120
+47 209
+47 210
+210 209
+177 10
+11 177
+189 126
+184 189
+211 44
+46 211
+63 112
+119 112
+42 158
+212 52
+122 212
+150 138
+150 198
+138 198
+53 140
+140 213
+213 53
+214 96
+214 98
+128 165
+172 169
+169 109
+172 109
+171 35
+171 37
+163 127
+50 215
+215 119
+202 79
+142 205
+142 216
+216 205
+120 133
+206 114
+115 206
+217 155
+154 217
+218 203
+202 218
+180 62
+4 62
+219 165
+164 219
+144 121
+122 144
+57 161
+57 35
+35 161
+215 120
+133 215
+163 128
+37 172
+124 104
+84 97
+97 186
+186 84
+209 188
+220 209
+188 220
+192 218
+39 192
+39 218
+221 152
+152 175
+175 221
+30 114
+213 222
+222 212
+212 213
+223 219
+219 224
+224 223
+212 195
+160 212
+195 160
+65 2
+68 65
+75 102
+225 102
+225 75
+226 128
+128 143
+143 226
+187 227
+227 217
+217 187
+139 66
+67 139
+201 200
+199 67
+199 115
+67 115
+19 207
+195 90
+169 37
+104 79
+79 124
+197 108
+108 228
+228 197
+181 196
+173 196
+229 230
+231 230
+229 231
+9 232
+232 7
+8 64
+193 221
+221 192
+233 137
+233 159
+174 223
+223 88
+88 174
+200 234
+234 235
+235 200
+217 207
+85 217
+207 85
+207 80
+80 19
+89 236
+102 236
+94 15
+94 208
+208 15
+218 38
+12 99
+99 3
+237 145
+145 54
+54 237
+238 239
+239 230
+230 238
+240 179
+201 240
+241 219
+228 241
+219 228
+177 94
+74 225
+225 237
+237 74
+237 55
+24 240
+23 240
+24 44
+44 22
+160 90
+171 90
+240 242
+242 201
+220 243
+243 244
+220 244
+35 59
+170 96
+170 238
+96 238
+69 131
+69 132
+201 243
+243 200
+66 140
+123 203
+159 203
+135 116
+31 135
+185 155
+155 187
+143 163
+163 144
+114 142
+3 60
+236 245
+226 236
+226 245
+246 16
+247 246
+16 247
+242 24
+52 213
+184 248
+45 248
+245 237
+145 245
+134 155
+185 134
+210 234
+209 234
+158 222
+160 222
+90 161
+241 197
+249 84
+186 249
+58 43
+246 166
+166 247
+109 73
+109 75
+222 42
+222 66
+100 117
+97 117
+97 100
+100 83
+84 100
+165 224
+230 204
+230 205
+64 1
+65 1
+116 199
+190 168
+249 185
+224 89
+236 224
+238 229
+191 159
+203 191
+168 246
+166 168
+174 241
+241 223
+108 188
+188 228
+243 242
+228 164
+250 231
+250 229
+215 134
+251 215
+134 251
+15 82
+208 82
+240 193
+193 179
+179 221
+22 211
+22 252
+211 252
+49 113
+112 49
+161 171
+130 251
+130 249
+249 251
+173 105
+146 58
+59 146
+209 108
+156 232
+177 208
+244 242
+242 248
+248 244
+159 77
+77 191
+199 135
+205 231
+231 216
+135 198
+98 186
+187 98
+5 111
+10 178
+178 190
+236 225
+245 225
+175 39
+239 170
+77 233
+247 167
+102 17
+196 241
+213 66
+129 145
+56 73
+74 56
+82 130
+82 249
+88 167
+251 50
+118 13
+234 220
+234 243
+231 103
+103 250
+196 174
+87 71
+157 37
+125 211
+244 189
+244 184
+183 76
+77 183
+29 25
+26 29
+238 214
+229 214
+86 79
+211 253
+253 125
+132 16
+16 110
+254 105
+166 254
+43 255
+58 255
+255 41
+246 132
+233 183
+103 30
+227 85
+256 227
+85 256
+12 148
+144 129
+224 226
+165 226
+254 190
+168 254
+68 131
+157 68
+204 239
+235 4
+251 185
+245 143
+103 216
+130 15
+181 118
+235 47
+5 235
+221 39
+253 194
+257 253
+257 194
+256 250
+250 86
+86 256
+197 48
+198 136
+190 95
+95 254
+253 182
+182 257
+204 26
+141 26
+255 70
+255 71
+149 69
+132 149
+136 233
+151 233
+98 227
+227 214
+74 55
+153 175
+170 28
+223 89
+113 118
+69 147
+115 41
+41 206
+156 68
+17 75
+214 256
+256 229
+111 6
+112 6
+258 151
+182 258
+104 250
+26 72
+72 141
+252 253
+252 182
+138 257
+18 247
+36 59
+149 246
+183 252
+252 76
+232 146
+189 220
+239 28
+216 30
+210 235
+101 38
+29 239
+29 204
+59 232
+232 36
+218 191
+118 49
+248 44
+247 88
+206 255
+71 58
+86 104
+117 96
+71 146
+154 21
+258 150
+254 93
+248 24
+216 114
+258 257
+257 150
+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 48
+49 50 51
+52 53 54
+55 56 57
+58 59 60
+61 62 63
+64 65 66
+67 68 69
+70 71 72
+73 74 75
+76 77 78
+79 80 81
+82 83 84
+85 86 87
+88 89 90
+91 92 93
+94 95 96
+97 26 98
+99 100 101
+102 103 104
+105 106 107
+108 109 110
+111 112 113
+114 115 116
+117 101 118
+119 120 121
+122 123 124
+125 126 127
+128 50 129
+130 131 132
+133 134 135
+136 137 138
+66 139 140
+141 142 143
+144 145 146
+147 148 149
+150 151 152
+126 153 154
+155 43 74
+156 157 158
+159 160 161
+53 162 163
+164 165 166
+167 168 169
+170 171 172
+173 174 175
+176 177 178
+179 180 181
+182 183 184
+185 186 187
+188 189 190
+54 191 192
+193 194 195
+196 197 198
+199 200 201
+202 203 204
+205 206 207
+208 209 210
+211 212 213
+214 215 216
+217 218 219
+220 221 222
+223 224 225
+226 227 228
+229 230 231
+232 146 233
+234 104 235
+189 236 237
+238 239 240
+139 241 242
+239 166 243
+244 245 97
+246 22 247
+248 249 250
+251 252 253
+254 255 256
+257 258 7
+259 260 261
+262 263 34
+264 105 265
+266 267 268
+269 270 271
+272 273 274
+275 276 277
+278 279 280
+281 282 263
+283 284 156
+285 286 287
+288 289 290
+291 292 31
+293 294 295
+296 297 298
+299 300 301
+302 303 304
+305 306 307
+308 309 160
+310 311 312
+313 314 308
+315 198 284
+316 317 318
+319 320 321
+322 323 267
+324 325 326
+327 328 329
+180 330 331
+332 72 333
+334 335 214
+336 178 337
+338 339 340
+250 341 342
+343 344 345
+304 346 300
+3 193 347
+348 169 349
+350 55 148
+351 233 352
+353 354 355
+356 10 357
+358 359 283
+47 360 361
+362 363 144
+45 211 364
+149 365 366
+367 368 369
+370 371 372
+109 373 374
+228 158 375
+376 377 378
+379 39 380
+197 223 381
+191 382 383
+89 384 343
+385 386 387
+352 388 276
+389 390 136
+207 391 392
+326 393 394
+395 321 396
+397 398 226
+399 147 400
+401 402 403
+388 404 405
+157 381 406
+380 407 249
+120 152 408
+409 410 411
+412 413 414
+415 416 417
+418 419 420
+138 291 421
+422 423 424
+425 426 427
+428 429 430
+431 76 432
+433 434 435
+436 437 438
+439 440 441
+442 443 73
+444 319 323
+445 446 447
+181 385 328
+340 448 19
+449 102 309
+407 450 376
+451 408 452
+453 454 455
+274 456 457
+458 459 460
+9 461 462
+295 358 227
+8 463 70
+464 465 346
+215 466 467
+468 469 470
+471 472 473
+474 475 476
+477 478 448
+118 479 480
+481 482 483
+484 417 40
+485 194 486
+487 488 489
+490 491 492
+493 494 322
+495 496 497
+107 259 498
+499 500 501
+502 489 58
+24 503 504
+505 23 506
+48 125 196
+430 507 449
+145 383 404
+103 508 248
+509 510 494
+511 512 513
+402 63 514
+161 234 59
+515 516 517
+277 172 35
+518 519 164
+268 395 203
+520 521 444
+177 442 522
+523 216 524
+369 316 176
+292 525 526
+527 287 528
+529 530 182
+331 389 531
+532 2 64
+533 534 535
+536 537 538
+12 485 167
+503 539 509
+540 424 365
+541 542 315
+543 544 487
+545 546 527
+32 334 150
+222 288 134
+355 547 548
+549 217 550
+507 219 551
+552 455 496
+553 411 554
+526 175 335
+212 555 61
+556 557 537
+558 559 83
+476 96 477
+75 560 561
+562 563 564
+92 565 566
+426 397 567
+568 329 569
+71 570 571
+137 446 572
+236 298 573
+574 554 285
+282 256 275
+479 575 576
+577 492 458
+524 578 579
+580 581 556
+582 583 468
+454 584 585
+510 586 520
+587 585 293
+460 588 589
+590 591 592
+593 483 594
+595 493 596
+596 597 464
+598 599 600
+135 601 602
+603 403 379
+604 605 606
+124 272 607
+551 603 508
+62 608 609
+129 353 610
+461 332 611
+497 587 398
+452 324 384
+612 357 348
+613 614 615
+616 617 578
+572 618 525
+619 387 620
+618 317 621
+286 622 623
+51 220 624
+297 625 626
+627 533 628
+41 629 254
+528 441 391
+516 630 490
+467 631 616
+557 231 632
+18 117 633
+634 312 552
+635 561 422
+98 257 186
+488 636 159
+637 82 638
+605 639 640
+253 470 641
+163 604 642
+643 143 13
+644 645 511
+504 595 303
+646 647 588
+127 350 224
+457 648 252
+57 372 540
+81 244 649
+210 650 38
+361 651 153
+566 564 409
+265 481 14
+652 653 359
+654 85 655
+650 240 450
+656 25 657
+658 577 659
+660 88 95
+478 90 114
+661 662 651
+165 663 664
+665 229 666
+667 668 555
+378 558 341
+648 634 582
+44 669 667
+663 670 536
+405 590 170
+671 655 631
+33 119 672
+486 532 112
+673 674 675
+676 11 296
+183 677 636
+678 567 679
+617 87 299
+573 680 681
+682 238 683
+168 111 565
+491 684 568
+320 473 685
+686 606 574
+535 438 687
+620 688 646
+162 15 689
+151 523 325
+690 273 141
+653 615 541
+691 49 692
+693 420 629
+694 695 696
+697 698 699
+307 337 370
+700 128 453
+621 701 173
+702 703 680
+704 705 695
+131 664 16
+706 707 327
+708 79 709
+519 710 711
+712 201 713
+714 374 715
+716 501 502
+255 717 351
+241 271 718
+713 280 671
+469 719 100
+687 184 544
+521 645 471
+456 289 310
+423 550 429
+204 69 232
+679 375 436
+142 123 264
+465 693 416
+290 690 720
+371 522 635
+318 445 443
+410 110 622
+366 428 314
+721 190 710
+722 390 723
+333 432 724
+675 699 94
+559 132 725
+592 686 546
+719 427 575
+394 484 344
+726 659 727
+701 368 199
+84 435 499
+728 133 729
+623 714 439
+65 570 269
+730 731 278
+626 260 702
+396 6 67
+437 406 529
+647 121 732
+681 666 581
+512 586 613
+733 734 707
+597 266 418
+735 736 704
+696 737 336
+600 735 661
+724 683 209
+60 638 716
+154 305 56
+377 243 130
+538 738 17
+514 739 37
+27 245 733
+740 237 580
+741 742 654
+187 462 743
+171 545 206
+640 93 553
+440 673 475
+414 744 294
+270 463 28
+279 736 741
+718 745 630
+746 672 688
+472 547 747
+311 221 700
+345 748 115
+749 750 684
+751 752 739
+68 729 362
+392 474 339
+517 658 373
+689 593 639
+579 753 393
+720 754 601
+602 192 363
+301 415 753
+213 401 218
+755 46 542
+754 643 52
+756 99 738
+723 757 669
+30 749 745
+113 140 562
+624 728 5
+174 712 466
+42 281 748
+757 330 708
+709 758 668
+727 589 697
+513 652 744
+759 451 660
+480 627 434
+734 80 179
+760 108 563
+583 495 425
+758 761 608
+77 347 721
+364 549 560
+242 515 760
+36 205 762
+200 763 730
+500 628 543
+594 349 91
+86 247 302
+703 106 764
+750 657 706
+447 722 155
+765 505 755
+677 400 313
+766 421 746
+685 692 4
+614 539 765
+767 763 768
+569 619 459
+195 676 188
+571 1 431
+235 342 637
+761 649 185
+698 732 759
+632 641 756
+576 678 534
+609 743 751
+386 531 766
+261 625 356
+715 726 674
+725 633 433
+742 599 246
+29 258 656
+230 607 251
+21 762 338
+584 610 412
+116 262 20
+747 354 691
+530 225 399
+764 122 665
+711 740 670
+419 202 717
+705 731 767
+382 642 591
+662 694 306
+78 518 682
+768 367 737
+506 598 360
+752 611 208
+548 644 413
+482 498 612
=== added file 'examples/pfacet/pfacetcreators.py'
--- examples/pfacet/pfacetcreators.py 1970-01-01 00:00:00 +0000
+++ examples/pfacet/pfacetcreators.py 2015-12-15 12:59:52 +0000
@@ -0,0 +1,119 @@
+# encoding: utf-8
+from yade import qt
+from yade.gridpfacet import *
+
+
+###########################
+##### ENGINES #####
+###########################
+
+O.engines=[
+ ForceResetter(),
+ InsertionSortCollider([
+ Bo1_Sphere_Aabb(),
+ Bo1_Wall_Aabb(),
+ Bo1_PFacet_Aabb(),
+ ],sortThenCollide=True),
+ InteractionLoop(
+ [
+ Ig2_GridNode_GridNode_GridNodeGeom6D(),
+ Ig2_Sphere_PFacet_ScGridCoGeom(),
+ Ig2_Wall_PFacet_ScGeom(),Ig2_Wall_Sphere_ScGeom()
+ ],
+ [
+ Ip2_CohFrictMat_CohFrictMat_CohFrictPhys(setCohesionNow=True,setCohesionOnNewContacts=True),
+ Ip2_FrictMat_FrictMat_FrictPhys()],
+ [
+ Law2_ScGeom6D_CohFrictPhys_CohesionMoment(),
+ Law2_ScGeom_FrictPhys_CundallStrack(),
+ Law2_ScGridCoGeom_FrictPhys_CundallStrack(),
+ Law2_GridCoGridCoGeom_FrictPhys_CundallStrack()
+ ]),
+ GlobalStiffnessTimeStepper(timestepSafetyCoefficient=0.8,label='ts'),
+ NewtonIntegrator(gravity=(0,-9.81,0),damping=0.1,label='newton')
+]
+
+
+O.materials.append(CohFrictMat(young=1e7,poisson=1,density=1e2,frictionAngle=radians(30),normalCohesion=3e7,shearCohesion=3e7,momentRotationLaw=True,label='gridNodeMat'))
+O.materials.append(FrictMat(young=1e7,poisson=1,density=1e2,frictionAngle=radians(30),label='gridConnectionMat'))
+
+
+###################################
+##### PFacet creators #####
+###################################
+
+fixed = False
+color=[255./255.,102./255.,0./255.]
+
+nodesIds=[]
+cylIds=[]
+
+#position of the node in the middle
+a=0.00
+r=0.03
+
+
+## Option 1: pfacet(id1,id2,id3) -> based on 3 gridNodes already connected via 3 gridConnections
+
+nodesIds.append( O.bodies.append(gridNode([0,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color)) )
+nodesIds.append( O.bodies.append(gridNode([1,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color)) )
+nodesIds.append( O.bodies.append(gridNode([0.5,1,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color)) )
+
+cylIds.append(O.bodies.append( gridConnection(0,1,r,color=color) ))
+cylIds.append(O.bodies.append( gridConnection(1,2,r,color=color) ))
+cylIds.append(O.bodies.append( gridConnection(2,0,r,color=color) ))
+
+O.bodies.append( pfacet(nodesIds[0],nodesIds[1],nodesIds[2],wire=False,color=color,highlight=False,material=O.materials[1]) )
+
+
+## Option 2: pfacetCreator1(vertices) -> based on 3 vertices
+
+v1=Vector3(2,0,0)
+v2=Vector3(3,0,0)
+v3=Vector3(2.5,1,0)
+vertices=[v1,v2,v3]
+pfacetCreator1(vertices,r,nodesIds=[],cylIds=[],pfIds=[],wire=False,color=color,fixed=fixed,materialNodes='gridNodeMat',material='gridConnectionMat')
+
+
+## Option 3: pfacetCreator2(id1,id2,vertex) -> based on 2 gridNodes connected via a gridConnection and a vertex
+
+nodesIds.append( O.bodies.append(gridNode([4,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color)) )
+nodesIds.append( O.bodies.append(gridNode([5,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color)) )
+vertex=Vector3(4.5,1,0)
+cylIds.append(O.bodies.append( gridConnection(nodesIds[3],nodesIds[4],r,color=color) ) )
+
+
+pfacetCreator2(nodesIds[3],nodesIds[4],vertex,r,nodesIds=nodesIds,wire=True,materialNodes='gridNodeMat',material='gridConnectionMat',color=color,fixed=fixed)
+
+
+## Option 4: pfacetCreator3(id1,id2,id3) -> based on 3 gridNodes
+
+a = O.bodies.append(gridNode([6,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color))
+b = O.bodies.append(gridNode([7,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color))
+c = O.bodies.append(gridNode([6.5,1,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color))
+
+pfacetCreator3(a,b,c,cylIds=[],pfIds=[],wire=False,material=-1,color=color)
+
+
+## Option 4: pfacetCreator4(id1,id2,id3) -> based on 3 gridConnections
+
+a = O.bodies.append(gridNode([8,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color))
+b = O.bodies.append(gridNode([9,0,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color))
+c = O.bodies.append(gridNode([8.5,1,0],r,wire=False,fixed=fixed,material='gridNodeMat',color=color))
+n=len(cylIds)
+cylIds.append(O.bodies.append( gridConnection(a,b,r,color=color) ) )
+cylIds.append(O.bodies.append( gridConnection(b,c,r,color=color) ) )
+cylIds.append(O.bodies.append( gridConnection(c,a,r,color=color) ) )
+
+pfacetCreator4(cylIds[n],cylIds[n+1],cylIds[n+2],pfIds=[],wire=False,material=-1,color=color)
+
+
+#####################
+##### Wall ###
+#####################
+
+O.bodies.append(utils.wall(position=-1,sense=0, axis=1,color=Vector3(1,0,0),material='gridConnectionMat'))
+
+
+qt.View()
+O.saveTmp()
=== added file 'examples/pfacet/sphere.gts'
--- examples/pfacet/sphere.gts 1970-01-01 00:00:00 +0000
+++ examples/pfacet/sphere.gts 2015-12-15 12:59:52 +0000
@@ -0,0 +1,243 @@
+42 120 80 GtsSurface GtsFace GtsEdge GtsVertex
+0 0 -1
+-0.3090169728 0.5 -0.8090170026
+0.3090169728 0.5 -0.8090170026
+-1 0 0
+-0.850650847 -0.5257311463 0
+-0.8090170026 -0.3090169728 0.5
+0.5 0.8090170026 0.3090169728
+0.5 0.8090170026 -0.3090169728
+0.850650847 0.5257311463 0
+0 1 0
+0 0.850650847 0.5257311463
+-0.5 0.8090170026 0.3090169728
+0 0.850650847 -0.5257311463
+-0.5 0.8090170026 -0.3090169728
+0 -0.850650847 0.5257311463
+0 -1 0
+0.5 -0.8090170026 0.3090169728
+0.5257311463 0 -0.850650847
+0.3090169728 -0.5 -0.8090170026
+0.8090170026 -0.3090169728 -0.5
+-0.5 -0.8090170026 -0.3090169728
+0 -0.850650847 -0.5257311463
+-0.3090169728 -0.5 -0.8090170026
+-0.8090170026 -0.3090169728 -0.5
+-0.5257311463 0 -0.850650847
+-0.8090170026 0.3090169728 0.5
+-0.8090170026 0.3090169728 -0.5
+0.5 -0.8090170026 -0.3090169728
+0.850650847 -0.5257311463 0
+0 0 1
+-0.5257311463 0 0.850650847
+-0.3090169728 -0.5 0.8090170026
+-0.850650847 0.5257311463 0
+0.3090169728 0.5 0.8090170026
+-0.3090169728 0.5 0.8090170026
+-0.5 -0.8090170026 0.3090169728
+0.5257311463 0 0.850650847
+0.8090170026 0.3090169728 0.5
+0.8090170026 -0.3090169728 0.5
+0.3090169728 -0.5 0.8090170026
+0.8090170026 0.3090169728 -0.5
+1 0 0
+1 2
+2 3
+3 1
+4 5
+5 6
+4 6
+7 8
+7 9
+9 8
+10 11
+12 10
+12 11
+13 10
+13 14
+10 14
+15 16
+17 16
+17 15
+18 19
+20 18
+19 20
+21 22
+23 21
+22 23
+24 23
+25 24
+23 25
+6 26
+26 4
+27 24
+27 25
+28 17
+28 29
+29 17
+30 31
+31 32
+30 32
+33 4
+26 33
+34 35
+34 11
+11 35
+1 25
+25 2
+36 21
+36 5
+5 21
+21 24
+11 7
+10 7
+16 22
+21 16
+2 27
+27 14
+14 2
+19 23
+19 22
+37 34
+38 37
+34 38
+38 7
+7 34
+35 12
+17 39
+39 40
+40 17
+15 36
+16 36
+8 13
+3 8
+13 3
+41 3
+18 41
+3 18
+9 42
+41 42
+41 9
+32 15
+36 32
+6 36
+38 39
+37 39
+39 29
+6 31
+31 26
+42 20
+20 41
+14 12
+14 33
+33 12
+24 4
+24 5
+42 29
+39 42
+16 28
+4 27
+8 10
+20 28
+29 20
+28 19
+22 28
+33 27
+37 30
+40 30
+40 37
+32 40
+15 40
+35 30
+35 31
+30 34
+42 38
+1 19
+23 1
+18 1
+2 13
+9 38
+32 6
+26 35
+12 26
+8 41
+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
+6 28 29
+30 31 26
+32 33 34
+35 36 37
+38 29 39
+40 41 42
+43 44 1
+45 46 47
+48 25 23
+10 49 50
+51 52 22
+53 54 55
+56 57 24
+58 59 60
+60 61 62
+12 63 42
+64 65 66
+16 67 68
+69 70 71
+72 73 74
+44 31 53
+75 76 77
+78 79 67
+80 5 46
+81 59 82
+34 83 64
+28 84 85
+86 87 76
+88 89 90
+68 45 52
+4 91 92
+93 94 83
+95 32 17
+15 88 11
+96 30 91
+50 7 97
+98 99 33
+57 100 101
+47 92 48
+38 102 96
+51 101 95
+103 104 105
+106 78 107
+35 108 109
+93 99 86
+103 58 110
+87 20 73
+111 81 94
+112 56 113
+18 66 107
+114 3 74
+37 106 104
+2 115 71
+110 40 108
+75 116 111
+114 19 112
+43 113 27
+117 80 79
+21 98 100
+13 97 69
+118 85 109
+41 62 49
+61 116 8
+90 39 119
+119 118 63
+115 55 14
+120 72 70
+9 77 120
+36 84 117
+54 102 89
+65 82 105
=== added file 'examples/pfacet/sphere.mesh'
--- examples/pfacet/sphere.mesh 1970-01-01 00:00:00 +0000
+++ examples/pfacet/sphere.mesh 2015-12-15 12:59:52 +0000
@@ -0,0 +1,253 @@
+ MeshVersionFormatted 2
+ Dimension
+ 3
+ Vertices
+ 67
+ 0 0 0 1
+ 0.1 0 0 2
+ 0 0.1 0 3
+ -0.1 0 0 4
+ 0 -0.1 0 5
+ 0 0 0.1 6
+ 0 0 -0.1 7
+ 0.092387953251167 0.038268343236416 0 1
+ 0.070710678118792 0.070710678118518 0 1
+ 0.038268343236596 0.092387953251093 0 1
+ -0.038268343236416 0.092387953251167 0 2
+ -0.070710678118518 0.070710678118792 0 2
+ -0.092387953251093 0.038268343236596 0 2
+ -0.092387953251167 -0.038268343236416 0 3
+ -0.070710678118792 -0.070710678118518 0 3
+ -0.038268343236596 -0.092387953251093 0 3
+ 0.038268343236416 -0.092387953251167 0 4
+ 0.070710678118518 -0.070710678118792 0 4
+ 0.092387953251093 -0.038268343236596 0 4
+ 0 0.038268343236416 0.092387953251167 6
+ 0 0.070710678118518 0.070710678118792 6
+ 0 0.092387953251093 0.038268343236596 6
+ 0 -0.038268343236416 0.092387953251167 7
+ 0 -0.070710678118518 0.070710678118792 7
+ 0 -0.092387953251093 0.038268343236596 7
+ 0.038268343236416 0 0.092387953251167 8
+ 0.070710678118518 0 0.070710678118792 8
+ 0.092387953251093 0 0.038268343236596 8
+ -0.038268343236416 0 0.092387953251167 9
+ -0.070710678118518 0 0.070710678118792 9
+ -0.092387953251093 0 0.038268343236596 9
+ 0 0.038268343236416 -0.092387953251167 10
+ 0 0.070710678118518 -0.070710678118792 10
+ 0 0.092387953251093 -0.038268343236596 10
+ 0.038268343236416 0 -0.092387953251167 11
+ 0.070710678118518 0 -0.070710678118792 11
+ 0.092387953251093 0 -0.038268343236596 11
+ 0 -0.038268343236416 -0.092387953251167 12
+ 0 -0.070710678118518 -0.070710678118792 12
+ 0 -0.092387953251093 -0.038268343236596 12
+ -0.038268343236416 0 -0.092387953251167 13
+ -0.070710678118518 0 -0.070710678118792 13
+ -0.092387953251093 0 -0.038268343236596 13
+ 0.077742384372375 0.044475395850299 -0.044475395850482 15
+ 0.044475395850283 0.044475395850283 -0.077742384372498 15
+ 0.044475395850498 0.077742384372253 -0.044475395850498 15
+ 0.044475395850299 -0.077742384372375 -0.044475395850482 17
+ 0.044475395850283 -0.044475395850283 -0.077742384372498 17
+ 0.077742384372253 -0.044475395850498 -0.044475395850498 17
+ -0.044475395850283 -0.044475395850283 -0.077742384372498 19
+ -0.077742384372375 -0.044475395850299 -0.044475395850482 19
+ -0.044475395850498 -0.077742384372253 -0.044475395850498 19
+ -0.044475395850283 0.044475395850283 -0.077742384372498 21
+ -0.044475395850299 0.077742384372375 -0.044475395850482 21
+ -0.077742384372253 0.044475395850498 -0.044475395850498 21
+ 0.044475395850299 -0.077742384372375 0.044475395850482 23
+ 0.044475395850283 -0.044475395850283 0.077742384372498 23
+ 0.077742384372253 -0.044475395850498 0.044475395850498 23
+ 0.044475395850283 0.044475395850283 0.077742384372498 25
+ 0.044475395850498 0.077742384372253 0.044475395850498 25
+ 0.077742384372375 0.044475395850299 0.044475395850482 25
+ -0.044475395850283 -0.044475395850283 0.077742384372498 27
+ -0.077742384372375 -0.044475395850299 0.044475395850482 27
+ -0.044475395850498 -0.077742384372253 0.044475395850498 27
+ -0.044475395850283 0.044475395850283 0.077742384372498 29
+ -0.044475395850299 0.077742384372375 0.044475395850482 29
+ -0.077742384372253 0.044475395850498 0.044475395850498 29
+ Edges
+ 48
+ 2 8 1
+ 8 9 1
+ 9 10 1
+ 10 3 1
+ 3 11 2
+ 11 12 2
+ 12 13 2
+ 13 4 2
+ 4 14 3
+ 14 15 3
+ 15 16 3
+ 16 5 3
+ 5 17 4
+ 17 18 4
+ 18 19 4
+ 19 2 4
+ 6 20 6
+ 20 21 6
+ 21 22 6
+ 22 3 6
+ 6 23 7
+ 23 24 7
+ 24 25 7
+ 25 5 7
+ 6 26 8
+ 26 27 8
+ 27 28 8
+ 28 2 8
+ 6 29 9
+ 29 30 9
+ 30 31 9
+ 31 4 9
+ 7 32 10
+ 32 33 10
+ 33 34 10
+ 34 3 10
+ 7 35 11
+ 35 36 11
+ 36 37 11
+ 37 2 11
+ 7 38 12
+ 38 39 12
+ 39 40 12
+ 40 5 12
+ 7 41 13
+ 41 42 13
+ 42 43 13
+ 43 4 13
+ Triangles
+ 128
+ 2 8 37 15
+ 8 44 37 15
+ 8 9 44 15
+ 37 44 36 15
+ 36 45 35 15
+ 45 32 35 15
+ 45 33 32 15
+ 35 32 7 15
+ 36 44 45 15
+ 44 46 45 15
+ 44 9 46 15
+ 45 46 33 15
+ 9 10 46 15
+ 10 34 46 15
+ 10 3 34 15
+ 46 34 33 15
+ 5 17 40 17
+ 17 47 40 17
+ 17 18 47 17
+ 40 47 39 17
+ 39 48 38 17
+ 48 35 38 17
+ 48 36 35 17
+ 38 35 7 17
+ 39 47 48 17
+ 47 49 48 17
+ 47 18 49 17
+ 48 49 36 17
+ 18 19 49 17
+ 19 37 49 17
+ 19 2 37 17
+ 49 37 36 17
+ 7 38 41 19
+ 38 50 41 19
+ 38 39 50 19
+ 41 50 42 19
+ 42 51 43 19
+ 51 14 43 19
+ 51 15 14 19
+ 43 14 4 19
+ 39 52 50 19
+ 52 51 50 19
+ 52 15 51 19
+ 50 51 42 19
+ 5 16 40 19
+ 16 52 40 19
+ 16 15 52 19
+ 40 52 39 19
+ 7 41 32 21
+ 41 53 32 21
+ 41 42 53 21
+ 32 53 33 21
+ 33 54 34 21
+ 54 11 34 21
+ 54 12 11 21
+ 34 11 3 21
+ 12 55 13 21
+ 55 43 13 21
+ 55 42 43 21
+ 13 43 4 21
+ 12 54 55 21
+ 54 53 55 21
+ 54 33 53 21
+ 55 53 42 21
+ 5 17 25 23
+ 17 56 25 23
+ 17 18 56 23
+ 25 56 24 23
+ 24 57 23 23
+ 57 26 23 23
+ 57 27 26 23
+ 23 26 6 23
+ 27 58 28 23
+ 58 19 28 23
+ 58 18 19 23
+ 28 19 2 23
+ 27 57 58 23
+ 57 56 58 23
+ 57 24 56 23
+ 58 56 18 23
+ 6 26 20 25
+ 26 59 20 25
+ 26 27 59 25
+ 20 59 21 25
+ 21 60 22 25
+ 60 10 22 25
+ 60 9 10 25
+ 22 10 3 25
+ 9 61 8 25
+ 61 28 8 25
+ 61 27 28 25
+ 8 28 2 25
+ 21 59 60 25
+ 59 61 60 25
+ 59 27 61 25
+ 60 61 9 25
+ 6 23 29 27
+ 23 62 29 27
+ 23 24 62 27
+ 29 62 30 27
+ 30 63 31 27
+ 63 14 31 27
+ 63 15 14 27
+ 31 14 4 27
+ 30 62 63 27
+ 62 64 63 27
+ 62 24 64 27
+ 63 64 15 27
+ 24 25 64 27
+ 25 16 64 27
+ 25 5 16 27
+ 64 16 15 27
+ 6 29 20 29
+ 29 65 20 29
+ 29 30 65 29
+ 20 65 21 29
+ 21 66 22 29
+ 66 11 22 29
+ 66 12 11 29
+ 22 11 3 29
+ 21 65 66 29
+ 65 67 66 29
+ 65 30 67 29
+ 66 67 12 29
+ 30 31 67 29
+ 31 13 67 29
+ 31 4 13 29
+ 67 13 12 29
+ End
