def test_clipper():
vertices, faces = mmio.load_VTP('meshmagick/tests/data/SEAREV.vtp')
searev = Mesh(vertices, faces)
plane = Plane()
clipper = mc.MeshClipper(searev,
plane,
assert_closed_boundaries=True,
verbose=False)
for iter in range(50):
thetax, thetay = np.random.rand(2) * 2 * math.pi
plane.rotate_normal(thetax, thetay)
clipper.plane = plane
#!/usr/bin/env python
# -*- coding: utf-8 -*-
import meshmagick.mmio as mmio
from meshmagick.mesh import Mesh
from math import pi, fabs
import pytest
vertices, faces = mmio.load_VTP('meshmagick/tests/data/Cylinder.vtp')
cylinder = Mesh(vertices, faces)
xmin, xmax, ymin, ymax, zmin, zmax = cylinder.axis_aligned_bbox
R = (xmax - xmin) / 2.
h = zmax - zmin
V = pi * R**2 * h
Ixx = Iyy = V * (R**2 + h**2 / 3) / 4
Izz = V * R**2 / 2
def test_cylinder_inertia():
inertia = cylinder.eval_plain_mesh_inertias(rho_medium=1.)
assert pytest.approx(Ixx, rel=1e-1) == inertia.xx
assert pytest.approx(Izz, rel=1e-1) == inertia.zz
def test_move_cog():
# parallel axis theorem to find MOI at base of cylinder
Ixx_base = Ixx + V * (h / 2)**2
# use meshmagick to find moments of inertia
inertia = cylinder.eval_plain_mesh_inertias(rho_medium=1.)