File reorganization to move toward standard python package layout

1 files changed, 0 insertions, 123 deletions

diff --git a/make.py b/make.py
deleted file mode 100644
index 0623b2b..0000000
--- a/make.py
+++ /dev/null
@@ -1,123 +0,0 @@
-import numpy as np
-from chroma.geometry import Mesh
-from chroma.transform import rotate
-from chroma.itertoolset import *
-
-def mesh_grid(grid):
-    return np.vstack(zip(grid[:-1].flatten(),grid[1:].flatten(),np.roll(grid[1:],-1,1).flatten()) + zip(grid[:-1].flatten(),np.roll(grid[1:],-1,1).flatten(),np.roll(grid[:-1],-1,1).flatten()))
-
-def linear_extrude(x1, y1, height, x2=None, y2=None, center=None):
-    """
-    Return the solid mesh formed by linearly extruding the polygon formed by
-    the x and y points `x1` and `y1` by a distance `height`. If `x2` and `y2`
-    are given extrude by connecting the points `x1` and `y1` to `x2` and `y2`;
-    this allows the creation of tapered solids.
-
-    .. note::
-        The path traced by the points `x` and `y` should go counter-clockwise,
-        otherwise the mesh will be inside out.
-
-    Example:
-        >>> # create a hexagon prism
-        >>> angles = np.linspace(0, 2*np.pi, 6, endpoint=False)
-        >>> m = linear_extrude(np.cos(angles), np.sin(angles), 2.0)
-    """
-    if len(x1) != len(y1):
-        raise Exception('`x` and `y` arrays must have the same length.')
-
-    if x2 is None:
-        x2 = x1
-
-    if y2 is None:
-        y2 = y1
-
-    if len(x2) != len(y2) or len(x2) != len(x1):
-        raise Exception('`x` and `y` arrays must have the same length.')
-
-    n = len(x1)
-
-    vertex_iterators = [izip(repeat(0,n),repeat(0,n),repeat(-height/2.0,n)),izip(x1,y1,repeat(-height/2.0,n)),izip(x2,y2,repeat(height/2.0,n)),izip(repeat(0,n),repeat(0,n),repeat(height/2.0,n))]
-
-    vertices = np.fromiter(flatten(roundrobin(*vertex_iterators)), float)
-    vertices = vertices.reshape((len(vertices)//3,3))
-
-    if center is not None:
-        vertices += center
-
-    triangles = mesh_grid(np.arange(len(vertices)).reshape((len(x1),len(vertices)//len(x1))).transpose()[::-1])
-
-    return Mesh(vertices, triangles, remove_duplicate_vertices=True)
-
-def rotate_extrude(x, y, nsteps=64):
-    """
-    Return the solid mesh formed by extruding the profile defined by the x and
-    y points `x` and `y` around the y axis.
-
-    .. note::
-        The path traced by the points `x` and `y` should go counter-clockwise,
-        otherwise the mesh will be inside out.
-
-    Example:
-        >>> # create a bipyramid
-        >>> m = rotate_extrude([0,1,0], [-1,0,1], nsteps=4)
-    """
-    if len(x) != len(y):
-        raise Exception('`x` and `y` arrays must have the same length.')
-
-    points = np.array([x,y,np.zeros(len(x))]).transpose()
-
-    steps = np.linspace(0, 2*np.pi, nsteps, endpoint=False)
-    vertices = np.vstack([rotate(points,angle,(0,-1,0)) for angle in steps])
-    triangles = mesh_grid(np.arange(len(vertices)).reshape((len(steps),len(points))).transpose()[::-1])
-
-    return Mesh(vertices, triangles, remove_duplicate_vertices=True)
-
-def box(dx, dy, dz, center=(0,0,0)):
-    "Return a box with linear dimensions `dx`, `dy`, and `dz`."
-    return linear_extrude([-dx/2.0,dx/2.0,dx/2.0,-dx/2.0],[-dy/2.0,-dy/2.0,dy/2.0,dy/2.0],height=dz,center=center)
-
-def cube(size=1, height=None, center=(0,0,0)):
-    "Return a cube mesh whose sides have length `size`."
-    if height is None:
-        height = size
-
-    return linear_extrude([-size/2.0,size/2.0,size/2.0,-size/2.0],[-size/2.0,-size/2.0,size/2.0,size/2.0], height=size, center=center)
-
-def cylinder(radius=1, height=2, radius2=None, nsteps=64):
-    """
-    Return a cylinder mesh with a radius of length `radius`, and a height of
-    length `height`. If `radius2` is specified, return a cone shaped cylinder
-    with bottom radius `radius`, and top radius `radius2`.
-    """
-    if radius2 is None:
-        radius2 = radius
-
-    return rotate_extrude([0,radius,radius2,0], [-height/2.0, -height/2.0, height/2.0, height/2.0], nsteps)
-
-def segmented_cylinder(radius, height=2, nsteps=64, nsegments=100):
-    """
-    Return a cylinder mesh segmented into `nsegments` points along its profile.
-    """
-    nsegments_radius = int((nsegments*radius/(2*radius+height))/2)
-    nsegments_height = int((nsegments*height/(2*radius+height))/2)
-    x = np.concatenate([np.linspace(0,radius,nsegments_radius,endpoint=False),[radius]*nsegments_height,np.linspace(radius,0,nsegments_radius,endpoint=False),[0]])
-    y = np.concatenate([[-height/2.0]*nsegments_radius,np.linspace(-height/2.0,height/2.0,nsegments_height,endpoint=False),[height/2.0]*(nsegments_radius+1)])
-    return rotate_extrude(x, y, nsteps)
-
-def sphere(radius=1, nsteps=64):
-    "Return a sphere mesh."
-    profile_angles = np.linspace(-np.pi/2, np.pi/2, nsteps)
-    return rotate_extrude(radius*np.cos(profile_angles), radius*np.sin(profile_angles), nsteps)
-
-def torus(radius=1, offset=3, nsteps=64, circle_steps=None):
-    """
-    Return a torus mesh. `offset` is the distance from the center of the torus
-    to the center of the torus barrel. `radius` is the radius of the torus
-    barrel. `nsteps` is the number of steps in the rotational extrusion of the
-    circle. `circle_steps` if specified is the number of steps around the
-    circumference of the torus barrel, else it defaults to `nsteps`.
-    """
-    if circle_steps is None:
-        circle_steps = nsteps
-    profile_angles = np.linspace(0, 2*np.pi, circle_steps)
-    return rotate_extrude(np.cos(profile_angles) + offset, np.sin(profile_angles), nsteps)