def test_meshStatistics():
mesh = loadLargeBunny()
area = pgf.read(pgf.area(mesh))
volume = pgf.read(pgf.volume(mesh))
nfaces = pgf.read(pgf.numFaces(mesh))
nverts = pgf.read(pgf.numVertices(mesh))
assert tu.equalf(area, 23.27821922302246)
assert tu.equalf(volume, 6.039211750030518)
assert nfaces == 69630
assert nverts == 34817
def test_meshPlaneClipping():
mesh = loadLargeBunny()
half = pgf.var_float(0.5)
plane = pgf.plane(pgf.vec3(half, half, half), pgf.vec3(half, half, half))
clipped = pgf.clipMesh(mesh, plane)
area = pgf.read(pgf.area(clipped))
nfaces = pgf.read(pgf.numFaces(clipped))
nverts = pgf.read(pgf.numVertices(clipped))
assert tu.equalf(area, 19.79486274)
assert nfaces == 59387
assert nverts == 30013
def test_rectangleMesh():
plane = pgf.plane(pgf.var_vec3((0., 0., 0.)), pgf.var_vec3((1., 1., 0.)))
box = pgf.box2(pgf.var_vec2((0., 0.)), pgf.var_vec2((15., 12.)))
edgeLn = pgf.var_float(1.)
rect = pgf.rectangleMesh(plane, box, edgeLn)
area = pgf.area(rect)
nfaces = pgf.numFaces(rect)
nverts = pgf.numVertices(rect)
assert 360 == pgf.read(nfaces)
assert 208 == pgf.read(nverts)
assert tu.equalf(180., pgf.read(area))
def test_scaleMesh():
mesh = loadLargeBunny()
factor = 10.0
scaled = pgf.scale(mesh, pgf.var_float(factor))
meshes = [mesh, scaled]
areas = [pgf.read(pgf.area(m)) for m in meshes]
volumes = [pgf.read(pgf.volume(m)) for m in meshes]
vcounts = [pgf.read(pgf.numVertices(m)) for m in meshes]
fcounts = [pgf.read(pgf.numFaces(m)) for m in meshes]
assert (abs(areas[1] - pow(factor, 2) * areas[0]) / areas[1]) < 1e-5
assert (abs(volumes[1] - pow(factor, 3) * volumes[0]) / volumes[1]) < 1e-5
assert vcounts[0] == vcounts[1]
assert fcounts[0] == fcounts[1]
