No commit message

No commit message

benchmarks/boolean/aorta_bound.lua

-
-
-options = gmshOptions()
-options:initOptions()
-
-options:numberSet('Mesh', 0, 'Algorithm3D', 4.0)
-options:numberSet('Mesh', 0, 'Optimize', 1.0)
---options:numberSet('Mesh', 0, 'CharacteristicLengthFactor', 0.2)
-options:numberSet('Mesh', 0, 'CharacteristicLengthExtendFromBoundary',  1)
-
-print'*** create GModel from stl ***'
-model = GModel()
-model:load ('aortaADAPT.stl')
-
-model:setFactory('Gmsh')
-
-print'*** create Topology ***'
-model:createTopology()
-
-e1 = model:getEdgeByTag(1)
-e2 = model:getEdgeByTag(2)
-e3 = model:getEdgeByTag(3)
-e4 = model:getEdgeByTag(4)
-e5 = model:getEdgeByTag(5)
-
-print'*** add Planar face ***'
-f0 = model:getFaceByTag(1) 
-
-f1 = model:addPlanarFace({{e1}})
-f2 = model:addPlanarFace({{e2}}) 
-f3 = model:addPlanarFace({{e3}}) 
-f4 = model:addPlanarFace({{e4}}) 
-f5 = model:addPlanarFace({{e5}}) 
-
-print'*** add Volume ***'
-r1 = model:addVolume({{f0,f1,f2,f3,f4,f5}});
-
-print'*** meshing ***'
-model:mesh(3);
-model:save("aorta_bound_lua.msh")
-
-
-

benchmarks/boolean/constraintSurface.lua

-
-g = GModel()
-v1 = g:addVertex(0, 0, 0, .1)
-v2 = g:addVertex(1, 0, 0, .1)
-v3 = g:addVertex(1, 1, 0, .1)
-v4 = g:addVertex(0, 3, 1, .1)
-v5 = g:addVertex(2, 3, 0, .1)
-v6 = g:addVertex(-2, 3, 0, .1)
-v7 = g:addVertex(-1, 3, 0, .1)
-
---v11 = g:addVertex(2.5, 2.5, 0, .1)
---v12 = g:addVertex(2.6, 2.5, 0, .1)
---v13 = g:addVertex(2.6, 2.6, 0, .1)
---v14 = g:addVertex(2.5, 2.6, 0, .1)
-
-
-e7 = g:addBezier (v2,v1, {{v3:x(),v3:y(),0},{v4:x(),v4:y(),0},{v5:x(),v5:y(),0},{v6:x(),v6:y(),0},{v7:x(),v7:y(),0}}); 							   
---e8 = g:addLine (v2,v1);
-
-v221 = g:addVertex(0, 0, 10, .1)
-v121 = g:addVertex(1, 0, 10, .1)
-e17 = g:addBezier (v121,v221, {{v3:x()+1,v3:y(),10},{v4:x(),v4:y()+2,10},{v5:x(),v5:y()-3,10},{v6:x(),v6:y(),10},{v7:x(),v7:y(),10}}); 							   
-
---e11 = g:addLine (v11,v12);
---e12 = g:addLine (v12,v13);
---e13 = g:addLine (v13,v14);
---e14 = g:addLine (v14,v11);
-
-
-g:addRuledFaces ({{e7},{e17}})
---g:addFace ({e7,e8},{{0,0,0}})
---g:addFace ({e11,e12,e13,e14})

benchmarks/boolean/cube1.lua

-
-g = GModel()
-
-v1 = g:addVertex(0, 0, 0, 1)
-v2 = g:addVertex(1, 0, 0, 1)
-e1 = g:addLine(v1, v2)
-f1 = g:extrude(e1, {0,0,0}, {0,10,0})
-r1 = g:extrude(f1, {0,0,0}, {0,0,1})
-r1 = g:revolve(f1, {-.1,0,0}, {-.1,1,0}, 3.1415 / 2)
-
---g:mesh(3)
-
-print('number of mesh vertices in volume ', r1, ' = ', r1:getNumMeshVertices())

benchmarks/boolean/nurbs.lua

-g = GModel();
-
-v1 = g:addVertex(0,0,0,.1); 
-v2 = g:addVertex(1,-.2,0,.1); 
-
-v3 = g:addVertex(1,0,0,.1); 
-v4 = g:addVertex(1,1,0,.1); 
-v5 = g:addVertex(2,1,0,.1); 
-v6 = g:addVertex(2.5,-.05,0,.1); 
-v7 = g:addVertex(1.5,-1,0,.1); 
-
---g:addLine(v1,v3);
---g:addLine(v3,v4);
---g:addLine(v4,v5);
---g:addLine(v5,v6);
---g:addLine(v6,v7);
---g:addLine(v7,v2);
-
---g:addNURBS (v1,v2, {{0.1,0.1,0},{0.2,0.1,0},{0.3,0.2,0},{0.4,0.6,0},{0.9,0.8,0}},
---	   {0,0,0,0,0.25,0.5,0.75,1,1,1,1},{1,1,1,1,1,1,1},{4,1,1,1,1,1,4}); 
-
---g:addBezier (v1,v2, {{v3:x(),v3:y(),0},{v4:x(),v4:y(),0},{v5:x(),v5:y(),0},{v6:x(),v6:y(),0},{v7:x(),v7:y(),0}}); 
-
-g:addNURBS (v1, v2, {{v3:x(),v3:y(),0},{v4:x(),v4:y(),0},{v5:x(),v5:y(),0},{v6:x(),v6:y(),0},{v7:x(),v7:y(),0}},          {0,0.25,0.5,0.75,1},{1,1,1,1,1,1,1},{4,1,1,1,4}); 							   
---g:addNURBS (v1,v2, {{v3:x(),v3:y(),0},{v4:x(),v4:y(),0},{v5:x(),v5:y(),0},{v6:x(),v6:y(),0},{v7:x(),v7:y(),0}},	   	   {0,0.5,1},{1,1,1,1,1,1,1},{4,3,4}); 							   
-
---e1 = g:addNURBS (v1,v2, {{v3:x(),v3:y(),0},{v4:x(),v4:y(),0},{v5:x(),v5:y(),0},{v6:x(),v6:y(),0},{v7:x(),v7:y(),0}},	   
---	   {0,1/3.,2./3.,1},{1,1,1,1,1,1,1},{5,1,1,5}); 							   
-
---f1 = g:extrude(e1, {0,0,0}, {0,0,1})

benchmarks/boolean/square2.lua

-
-g = GModel()
-v1 = g:addVertex(0, 0, 0, 1)
-v2 = g:addVertex(1, 0, 0, 1)
-e1 = g:addLine(v1, v2)
-f1 = g:extrude(e1, {0,0,0}, {0,1,0})
-

benchmarks/boolean/wikipedia.lua

-
--- from http://en.wikipedia.org/wiki/Constructive_solid_geometry
-
-R = 1.4;
-s = .7;
-t = 1.25;
-myModel = GModel();
-myModel:addBlock({-R,-R,-R},{R,R,R});
-
-myTool = GModel();
-myTool:addSphere(0,0,0,R*t);
-
-myModel:computeIntersection(myTool,0);
-
-myTool2 = GModel();
-myTool2:addCylinder({-2*R,0,0},{2*R,0,0},R*s);
-
-myTool3 = GModel();
-myTool3:addCylinder({0,-2*R,0},{0,2*R,0},R*s);
-
-myModel2 = GModel();
-myModel2:addCylinder({0,0,-2*R},{0,0,2*R},R*s);
-myModel2:computeUnion(myTool2,0);
-myModel2:computeUnion(myTool3,0);
-
-myModel:computeDifference(myModel2,0);
-
-myModel:setAsCurrent();