@jscad/modeling 2.9.1 → 2.9.4

package/src/operations/modifiers/generalize.test.js

@@ -104,38 +104,6 @@ test('generalize: generalize of a geom3 produces an expected geom3', (t) => {
   ]
   t.notThrows(() => geom3.validate(result))
   t.true(comparePolygonsAsPoints(pts, exp))
-
-  // apply repairs only (triangles)
-  result = generalize({ repair: true }, geometry2)
-  pts = geom3.toPoints(result)
-  exp = [
-    [[-1.5707963267948966, -0.7853981633974483, -3.141592653589793], [-1.5707963267948966, -0.7853981633974483, 3.141592653589793],
-      [-1.5707963267948966, 0.7853981633974483, 3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, -3.141592653589793], [-1.5707963267948966, 0.7853981633974483, 3.141592653589793],
-      [-1.5707963267948966, 0.7853981633974483, -3.141592653589793]],
-    [[1.5707963267948966, -0.7853981633974483, -3.141592653589793], [1.5707963267948966, 0.7853981633974483, -3.141592653589793],
-      [1.5707963267948966, 0.7853981633974483, 3.141592653589793]],
-    [[1.5707963267948966, -0.7853981633974483, -3.141592653589793], [1.5707963267948966, 0.7853981633974483, 3.141592653589793],
-      [1.5707963267948966, -0.7853981633974483, 3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, -3.141592653589793], [1.5707963267948966, -0.7853981633974483, -3.141592653589793],
-      [1.5707963267948966, -0.7853981633974483, 3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, -3.141592653589793], [1.5707963267948966, -0.7853981633974483, 3.141592653589793],
-      [-1.5707963267948966, -0.7853981633974483, 3.141592653589793]],
-    [[-1.5707963267948966, 0.7853981633974483, -3.141592653589793], [-1.5707963267948966, 0.7853981633974483, 3.141592653589793],
-      [1.5707963267948966, 0.7853981633974483, 3.141592653589793]],
-    [[-1.5707963267948966, 0.7853981633974483, -3.141592653589793], [1.5707963267948966, 0.7853981633974483, 3.141592653589793],
-      [1.5707963267948966, 0.7853981633974483, -3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, -3.141592653589793], [-1.5707963267948966, 0.7853981633974483, -3.141592653589793],
-      [1.5707963267948966, 0.7853981633974483, -3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, -3.141592653589793], [1.5707963267948966, 0.7853981633974483, -3.141592653589793],
-      [1.5707963267948966, -0.7853981633974483, -3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, 3.141592653589793], [1.5707963267948966, -0.7853981633974483, 3.141592653589793],
-      [1.5707963267948966, 0.7853981633974483, 3.141592653589793]],
-    [[-1.5707963267948966, -0.7853981633974483, 3.141592653589793], [1.5707963267948966, 0.7853981633974483, 3.141592653589793],
-      [-1.5707963267948966, 0.7853981633974483, 3.141592653589793]]
-  ]
-  t.notThrows(() => geom3.validate(result))
-  t.true(comparePolygonsAsPoints(pts, exp))
 })
 
 test('generalize: generalize of a geom3 with T junctions produces an expected geom3', (t) => {

package/src/operations/modifiers/insertTjunctions.js

@@ -259,7 +259,7 @@ const insertTjunctions = (polygons) => {
                     // split the side by inserting the vertex:
                     const newvertices = polygon.vertices.slice(0)
                     newvertices.splice(insertionvertextagindex, 0, endvertex)
-                    const newpolygon = poly3.fromPoints(newvertices)
+                    const newpolygon = poly3.create(newvertices)
 
                     newpolygons[polygonindex] = newpolygon
 

package/src/operations/modifiers/insertTjunctions.test.js

@@ -56,33 +56,33 @@ test('insertTjunctions: insertTjunctions produces expected polygons', (t) => {
 
   const result3 = insertTjunctions(geom3.toPolygons(geometry3))
   let exp = [
-    poly3.fromPoints([[-1, -1, -1], [-1, -1, 1], [-1, 1, 1], [-1, 1, -1]]),
-    poly3.fromPoints([[1, -1, -1], [1, 1, -1], [1, 1, 1], [1, -1, 1]]),
-    poly3.fromPoints([[-1, -1, -1], [1, -1, -1], [1, -1, 1], [-1, -1, 1]]),
-    poly3.fromPoints([[-1, 1, -1], [-1, 1, 1], [1, 1, 1], [1, 1, -1]]),
-    poly3.fromPoints([[-1, -1, -1], [-1, 1, -1], [1, 1, -1], [1, -1, -1]]),
-    poly3.fromPoints([[0, 0, 1], [-1, -1, 1], [1, -1, 1], [1, 1, 1]]),
-    poly3.fromPoints([[1, 1, 1], [-1, 1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[-1, 1, 1], [-1, -1, 1], [0, 0, 1]])
+    poly3.create([[-1, -1, -1], [-1, -1, 1], [-1, 1, 1], [-1, 1, -1]]),
+    poly3.create([[1, -1, -1], [1, 1, -1], [1, 1, 1], [1, -1, 1]]),
+    poly3.create([[-1, -1, -1], [1, -1, -1], [1, -1, 1], [-1, -1, 1]]),
+    poly3.create([[-1, 1, -1], [-1, 1, 1], [1, 1, 1], [1, 1, -1]]),
+    poly3.create([[-1, -1, -1], [-1, 1, -1], [1, 1, -1], [1, -1, -1]]),
+    poly3.create([[0, 0, 1], [-1, -1, 1], [1, -1, 1], [1, 1, 1]]),
+    poly3.create([[1, 1, 1], [-1, 1, 1], [0, 0, 1]]),
+    poly3.create([[-1, 1, 1], [-1, -1, 1], [0, 0, 1]])
   ]
   t.not(result3, geom3.toPolygons(geometry3))
   t.true(comparePolygonLists(result3, exp))
 
   const result4 = insertTjunctions(geom3.toPolygons(geometry4))
   exp = [
-    poly3.fromPoints([[-1, -1, -1], [-1, -1, 1], [-1, 0, 1], [-1, 1, 1], [-1, 1, -1]]),
-    poly3.fromPoints([[1, -1, -1], [1, 1, -1], [1, 1, 1], [1, 0, 1], [1, -1, 1]]),
-    poly3.fromPoints([[-1, -1, -1], [1, -1, -1], [1, -1, 1], [0, -1, 1], [-1, -1, 1]]),
-    poly3.fromPoints([[-1, 1, -1], [-1, 1, 1], [0, 1, 1], [1, 1, 1], [1, 1, -1]]),
-    poly3.fromPoints([[-1, -1, -1], [-1, 1, -1], [1, 1, -1], [1, -1, -1]]),
-    poly3.fromPoints([[-1, -1, 1], [0, -1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[-1, 0, 1], [-1, -1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[0, -1, 1], [1, -1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[1, -1, 1], [1, 0, 1], [0, 0, 1]]),
-    poly3.fromPoints([[1, 0, 1], [1, 1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[1, 1, 1], [0, 1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[0, 1, 1], [-1, 1, 1], [0, 0, 1]]),
-    poly3.fromPoints([[-1, 1, 1], [-1, 0, 1], [0, 0, 1]])
+    poly3.create([[-1, -1, -1], [-1, -1, 1], [-1, 0, 1], [-1, 1, 1], [-1, 1, -1]]),
+    poly3.create([[1, -1, -1], [1, 1, -1], [1, 1, 1], [1, 0, 1], [1, -1, 1]]),
+    poly3.create([[-1, -1, -1], [1, -1, -1], [1, -1, 1], [0, -1, 1], [-1, -1, 1]]),
+    poly3.create([[-1, 1, -1], [-1, 1, 1], [0, 1, 1], [1, 1, 1], [1, 1, -1]]),
+    poly3.create([[-1, -1, -1], [-1, 1, -1], [1, 1, -1], [1, -1, -1]]),
+    poly3.create([[-1, -1, 1], [0, -1, 1], [0, 0, 1]]),
+    poly3.create([[-1, 0, 1], [-1, -1, 1], [0, 0, 1]]),
+    poly3.create([[0, -1, 1], [1, -1, 1], [0, 0, 1]]),
+    poly3.create([[1, -1, 1], [1, 0, 1], [0, 0, 1]]),
+    poly3.create([[1, 0, 1], [1, 1, 1], [0, 0, 1]]),
+    poly3.create([[1, 1, 1], [0, 1, 1], [0, 0, 1]]),
+    poly3.create([[0, 1, 1], [-1, 1, 1], [0, 0, 1]]),
+    poly3.create([[-1, 1, 1], [-1, 0, 1], [0, 0, 1]])
   ]
   t.not(result4, geometry4)
   t.true(comparePolygonLists(result4, exp))

package/src/operations/modifiers/mergePolygons.js

@@ -1,3 +1,4 @@
+const aboutEqualNormals = require('../../maths/utils/aboutEqualNormals')
 const vec3 = require('../../maths/vec3')
 
 const poly3 = require('../../geometries/poly3')
@@ -53,9 +54,12 @@ const calculateAnglesBetween = (current, opposite, normal) => {
   return [angle1, angle2]
 }
 
+const v1 = vec3.create()
+const v2 = vec3.create()
+
 const calculateAngle = (prevpoint, point, nextpoint, normal) => {
-  const d0 = vec3.subtract(vec3.create(), point, prevpoint)
-  const d1 = vec3.subtract(vec3.create(), nextpoint, point)
+  const d0 = vec3.subtract(v1, point, prevpoint)
+  const d1 = vec3.subtract(v2, nextpoint, point)
   vec3.cross(d0, d0, d1)
   return vec3.dot(d0, normal)
 }
@@ -76,7 +80,7 @@ const createPolygonAnd = (edge) => {
 
     edge = next
   }
-  if (points.length > 0) polygon = poly3.fromPoints(points)
+  if (points.length > 0) polygon = poly3.create(points)
   return polygon
 }
 
@@ -85,7 +89,7 @@ const createPolygonAnd = (edge) => {
  * @param {poly3[]} sourcepolygons - list of polygons
  * @returns {poly3[]} new set of polygons
  */
-const mergeCoplanarPolygons = (epsilon, sourcepolygons) => {
+const mergeCoplanarPolygons = (sourcepolygons) => {
   if (sourcepolygons.length < 2) return sourcepolygons
 
   const normal = sourcepolygons[0].plane
@@ -167,18 +171,11 @@ const mergeCoplanarPolygons = (epsilon, sourcepolygons) => {
     if (polygon) destpolygons.push(polygon)
   })
 
+  edgeList.clear()
+
   return destpolygons
 }
 
-// Normals are directional vectors with component values from 0 to 1.0, requiring specialized comparision
-// This EPS is derived from a serieas of tests to determine the optimal precision for comparing coplanar polygons,
-// as provided by the sphere primitive at high segmentation
-// This EPS is for 64 bit Number values
-const NEPS = 1e-13
-
-// Compare two normals (unit vectors) for equality.
-const aboutEqualNormals = (a, b) => (Math.abs(a[0] - b[0]) <= NEPS && Math.abs(a[1] - b[1]) <= NEPS && Math.abs(a[2] - b[2]) <= NEPS)
-
 const coplanar = (plane1, plane2) => {
   // expect the same distance from the origin, within tolerance
   if (Math.abs(plane1[3] - plane2[3]) < 0.00000015) {
@@ -202,7 +199,7 @@ const mergePolygons = (epsilon, polygons) => {
   let destpolygons = []
   polygonsPerPlane.forEach((mapping) => {
     const sourcepolygons = mapping[1]
-    const retesselayedpolygons = mergeCoplanarPolygons(epsilon, sourcepolygons)
+    const retesselayedpolygons = mergeCoplanarPolygons(sourcepolygons)
     destpolygons = destpolygons.concat(retesselayedpolygons)
   })
   return destpolygons

package/src/operations/{booleans → modifiers}/reTesselateCoplanarPolygons.js

@@ -23,15 +23,14 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
   const orthobasis = new OrthoNormalBasis(plane)
   const polygonvertices2d = [] // array of array of Vector2D
   const polygontopvertexindexes = [] // array of indexes of topmost vertex per polygon
-  const topy2polygonindexes = {}
-  const ycoordinatetopolygonindexes = {}
-
-  const ycoordinatebins = {}
+  const topy2polygonindexes = new Map()
+  const ycoordinatetopolygonindexes = new Map()
 
   // convert all polygon vertices to 2D
   // Make a list of all encountered y coordinates
   // And build a map of all polygons that have a vertex at a certain y coordinate:
-  const ycoordinateBinningFactor = 1.0 / EPS * 10
+  const ycoordinatebins = new Map()
+  const ycoordinateBinningFactor = 10 / EPS
   for (let polygonindex = 0; polygonindex < numpolygons; polygonindex++) {
     const poly3d = sourcepolygons[polygonindex]
     let vertices2d = []
@@ -46,15 +45,15 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
         // close to each other, give them the same y coordinate:
         const ycoordinatebin = Math.floor(pos2d[1] * ycoordinateBinningFactor)
         let newy
-        if (ycoordinatebin in ycoordinatebins) {
-          newy = ycoordinatebins[ycoordinatebin]
-        } else if (ycoordinatebin + 1 in ycoordinatebins) {
-          newy = ycoordinatebins[ycoordinatebin + 1]
-        } else if (ycoordinatebin - 1 in ycoordinatebins) {
-          newy = ycoordinatebins[ycoordinatebin - 1]
+        if (ycoordinatebins.has(ycoordinatebin)) {
+          newy = ycoordinatebins.get(ycoordinatebin)
+        } else if (ycoordinatebins.has(ycoordinatebin + 1)) {
+          newy = ycoordinatebins.get(ycoordinatebin + 1)
+        } else if (ycoordinatebins.has(ycoordinatebin - 1)) {
+          newy = ycoordinatebins.get(ycoordinatebin - 1)
         } else {
           newy = pos2d[1]
-          ycoordinatebins[ycoordinatebin] = pos2d[1]
+          ycoordinatebins.set(ycoordinatebin, pos2d[1])
         }
         pos2d = vec2.fromValues(pos2d[0], newy)
         vertices2d.push(pos2d)
@@ -66,10 +65,12 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
         if ((i === 0) || (y > maxy)) {
           maxy = y
         }
-        if (!(y in ycoordinatetopolygonindexes)) {
-          ycoordinatetopolygonindexes[y] = {}
+        let polygonindexes = ycoordinatetopolygonindexes.get(y)
+        if (!polygonindexes) {
+          polygonindexes = {} // PERF
+          ycoordinatetopolygonindexes.set(y, polygonindexes)
         }
-        ycoordinatetopolygonindexes[y][polygonindex] = true
+        polygonindexes[polygonindex] = true
       }
       if (miny >= maxy) {
         // degenerate polygon, all vertices have same y coordinate. Just ignore it from now:
@@ -77,10 +78,12 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
         numvertices = 0
         minindex = -1
       } else {
-        if (!(miny in topy2polygonindexes)) {
-          topy2polygonindexes[miny] = []
+        let polygonindexes = topy2polygonindexes.get(miny)
+        if (!polygonindexes) {
+          polygonindexes = []
+          topy2polygonindexes.set(miny, polygonindexes)
         }
-        topy2polygonindexes[miny].push(polygonindex)
+        polygonindexes.push(polygonindex)
       }
     } // if(numvertices > 0)
     // reverse the vertex order:
@@ -89,8 +92,9 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
     polygonvertices2d.push(vertices2d)
     polygontopvertexindexes.push(minindex)
   }
+
   const ycoordinates = []
-  for (const ycoordinate in ycoordinatetopolygonindexes) ycoordinates.push(ycoordinate)
+  ycoordinatetopolygonindexes.forEach((polylist, y) => ycoordinates.push(y))
   ycoordinates.sort(fnNumberSort)
 
   // Now we will iterate over all y coordinates, from lowest to highest y coordinate
@@ -108,15 +112,14 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
   let prevoutpolygonrow = []
   for (let yindex = 0; yindex < ycoordinates.length; yindex++) {
     const newoutpolygonrow = []
-    const ycoordinateasstring = ycoordinates[yindex]
-    const ycoordinate = Number(ycoordinateasstring)
+    const ycoordinate = ycoordinates[yindex]
 
     // update activepolygons for this y coordinate:
     // - Remove any polygons that end at this y coordinate
     // - update leftvertexindex and rightvertexindex (which point to the current vertex index
     //   at the the left and right side of the polygon
     // Iterate over all polygons that have a corner at this y coordinate:
-    const polygonindexeswithcorner = ycoordinatetopolygonindexes[ycoordinateasstring]
+    const polygonindexeswithcorner = ycoordinatetopolygonindexes.get(ycoordinate)
     for (let activepolygonindex = 0; activepolygonindex < activepolygons.length; ++activepolygonindex) {
       const activepolygon = activepolygons[activepolygonindex]
       const polygonindex = activepolygon.polygonindex
@@ -166,7 +169,7 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
       nextycoordinate = Number(ycoordinates[yindex + 1])
       const middleycoordinate = 0.5 * (ycoordinate + nextycoordinate)
       // update activepolygons by adding any polygons that start here:
-      const startingpolygonindexes = topy2polygonindexes[ycoordinateasstring]
+      const startingpolygonindexes = topy2polygonindexes.get(ycoordinate)
       for (const polygonindexKey in startingpolygonindexes) {
         const polygonindex = startingpolygonindexes[polygonindexKey]
         const vertices2d = polygonvertices2d[polygonindex]
@@ -212,8 +215,6 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
         })
       } // for(let polygonindex in startingpolygonindexes)
     } //  yindex < ycoordinates.length-1
-    // if( (yindex === ycoordinates.length-1) || (nextycoordinate - ycoordinate > EPS) )
-    // FIXME : what ???
 
     // Now activepolygons is up to date
     // Build the output polygons for the next row in newoutpolygonrow:
@@ -252,19 +253,19 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
     } // for(activepolygon in activepolygons)
     if (yindex > 0) {
       // try to match the new polygons against the previous row:
-      const prevcontinuedindexes = {}
-      const matchedindexes = {}
+      const prevcontinuedindexes = new Set()
+      const matchedindexes = new Set()
       for (let i = 0; i < newoutpolygonrow.length; i++) {
         const thispolygon = newoutpolygonrow[i]
         for (let ii = 0; ii < prevoutpolygonrow.length; ii++) {
-          if (!matchedindexes[ii]) { // not already processed?
+          if (!matchedindexes.has(ii)) { // not already processed?
             // We have a match if the sidelines are equal or if the top coordinates
             // are on the sidelines of the previous polygon
             const prevpolygon = prevoutpolygonrow[ii]
             if (vec2.distance(prevpolygon.bottomleft, thispolygon.topleft) < EPS) {
               if (vec2.distance(prevpolygon.bottomright, thispolygon.topright) < EPS) {
                 // Yes, the top of this polygon matches the bottom of the previous:
-                matchedindexes[ii] = true
+                matchedindexes.add(ii)
                 // Now check if the joined polygon would remain convex:
                 const v1 = line2.direction(thispolygon.leftline)
                 const v2 = line2.direction(prevpolygon.leftline)
@@ -284,16 +285,16 @@ const reTesselateCoplanarPolygons = (sourcepolygons) => {
                   thispolygon.outpolygon = prevpolygon.outpolygon
                   thispolygon.leftlinecontinues = leftlinecontinues
                   thispolygon.rightlinecontinues = rightlinecontinues
-                  prevcontinuedindexes[ii] = true
+                  prevcontinuedindexes.add(ii)
                 }
                 break
               }
             }
-          } // if(!prevcontinuedindexes[ii])
+          } // if(!prevcontinuedindexes.has(ii))
         } // for ii
       } // for i
       for (let ii = 0; ii < prevoutpolygonrow.length; ii++) {
-        if (!prevcontinuedindexes[ii]) {
+        if (!prevcontinuedindexes.has(ii)) {
           // polygon ends here
           // Finish the polygon with the last point(s):
           const prevpolygon = prevoutpolygonrow[ii]

package/src/operations/{booleans → modifiers}/reTesselateCoplanarPolygons.test.js

@@ -17,11 +17,11 @@ const rotatePoly3 = (angles, polygon) => {
 }
 
 test.only('retessellateCoplanarPolygons: should merge coplanar polygons', (t) => {
-  const polyA = poly3.fromPoints([[-5, -5, 0], [5, -5, 0], [5, 5, 0], [-5, 5, 0]])
-  const polyB = poly3.fromPoints([[5, -5, 0], [8, 0, 0], [5, 5, 0]])
-  const polyC = poly3.fromPoints([[-5, 5, 0], [-8, 0, 0], [-5, -5, 0]])
-  const polyD = poly3.fromPoints([[-5, 5, 0], [5, 5, 0], [0, 8, 0]])
-  const polyE = poly3.fromPoints([[5, -5, 0], [-5, -5, 0], [0, -8, 0]])
+  const polyA = poly3.create([[-5, -5, 0], [5, -5, 0], [5, 5, 0], [-5, 5, 0]])
+  const polyB = poly3.create([[5, -5, 0], [8, 0, 0], [5, 5, 0]])
+  const polyC = poly3.create([[-5, 5, 0], [-8, 0, 0], [-5, -5, 0]])
+  const polyD = poly3.create([[-5, 5, 0], [5, 5, 0], [0, 8, 0]])
+  const polyE = poly3.create([[5, -5, 0], [-5, -5, 0], [0, -8, 0]])
 
   // combine polygons in each direction
   let obs = reTesselateCoplanarPolygons([polyA, polyB])

package/src/operations/{booleans → modifiers}/retessellate.js

@@ -1,16 +1,9 @@
 const geom3 = require('../../geometries/geom3')
 const poly3 = require('../../geometries/poly3')
 
-const reTesselateCoplanarPolygons = require('./reTesselateCoplanarPolygons')
-
-// Normals are directional vectors with component values from 0 to 1.0, requiring specialized comparison
-// This EPS is derived from a series of tests to determine the optimal precision for comparing coplanar polygons,
-// as provided by the sphere primitive at high segmentation
-// This EPS is for 64 bit Number values
-const NEPS = 1e-13
+const aboutEqualNormals = require('../../maths/utils/aboutEqualNormals')
 
-// Compare two normals (unit vectors) for equality.
-const aboutEqualNormals = (a, b) => (Math.abs(a[0] - b[0]) <= NEPS && Math.abs(a[1] - b[1]) <= NEPS && Math.abs(a[2] - b[2]) <= NEPS)
+const reTesselateCoplanarPolygons = require('./reTesselateCoplanarPolygons')
 
 const coplanar = (plane1, plane2) => {
   // expect the same distance from the origin, within tolerance

package/src/operations/modifiers/snapPolygons.test.js

@@ -7,49 +7,49 @@ const snapPolygons = require('./snapPolygons')
 test('snapPolygons: snap of polygons produces expected results', (t) => {
   const polygons = [
     // valid polygons
-    poly3.fromPoints([[0, 0, 0], [0, 10, 0], [0, 10, 10]]), // OK
-    poly3.fromPoints([[0, 0, 0], [0, 10, 0], [0, 10, 10], [0, 0, 10]]), // OK
+    poly3.create([[0, 0, 0], [0, 10, 0], [0, 10, 10]]), // OK
+    poly3.create([[0, 0, 0], [0, 10, 0], [0, 10, 10], [0, 0, 10]]), // OK
     // invalid polygons
     poly3.create(),
-    poly3.fromPoints([[0, 0, 0]]),
-    poly3.fromPoints([[0, 0, 0], [0, 10, 0]]),
+    poly3.create([[0, 0, 0]]),
+    poly3.create([[0, 0, 0], [0, 10, 0]]),
     // duplicated vertices
-    poly3.fromPoints([
+    poly3.create([
       [-24.445112000000115, 19.346837333333426, 46.47572533333356],
       [-24.44446933333345, 19.346837333333426, 46.47508266666689],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019]
     ]), // OK
-    poly3.fromPoints([
+    poly3.create([
       [-24.445112000000115, 19.346837333333426, 46.47572533333356],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019]
     ]),
-    poly3.fromPoints([
+    poly3.create([
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019]
     ]),
     // duplicate vertices after snap
-    poly3.fromPoints([
+    poly3.create([
       [-24.445112000000115, 19.346837333333426, 46.47572533333356],
       [-24.44446933333345, 19.346837333333426, 46.47508266666689],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678 - 0.00001234, 18.79864266666676 + 0.000001234, 39.56448800000019 + 0.00001234]
     ]), // OK
-    poly3.fromPoints([
+    poly3.create([
       [-24.445112000000115, 19.346837333333426, 46.47572533333356],
       [-23.70540266666678 - 0.00001234, 18.79864266666676 + 0.000001234, 39.56448800000019 + 0.00001234],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678 - 0.00001234, 18.79864266666676 + 0.000001234, 39.56448800000019 + 0.00001234]
     ]),
-    poly3.fromPoints([
+    poly3.create([
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678 - 0.00001234, 18.79864266666676 + 0.000001234, 39.56448800000019 + 0.00001234],
       [-23.70540266666678, 18.79864266666676, 39.56448800000019],
       [-23.70540266666678 - 0.00001234, 18.79864266666676 + 0.000001234, 39.56448800000019 + 0.00001234]
     ]),
     // inverted polygon
-    poly3.fromPoints([
+    poly3.create([
       [20.109133333333336, -4.894033333333335, -1.0001266666666668],
       [20.021120000000003, -5.1802133333333344, -1.0001266666666668],
       [20.020300000000002, -5.182946666666668, -1.0001266666666668],
@@ -61,7 +61,7 @@ test('snapPolygons: snap of polygons produces expected results', (t) => {
   const results = snapPolygons(0.0001, polygons)
   t.is(results.length, 5)
 
-  const exp3 = poly3.fromPoints([
+  const exp3 = poly3.create([
     [-24.4451, 19.3468, 46.4757],
     [-24.4445, 19.3468, 46.475100000000005],
     [-23.7054, 18.7986, 39.5645]

package/src/operations/modifiers/triangulatePolygons.js

@@ -10,15 +10,15 @@ const triangulatePolygon = (epsilon, polygon, triangles) => {
       polygon.vertices.forEach((vertice) => vec3.add(midpoint, midpoint, vertice))
       vec3.snap(midpoint, vec3.divide(midpoint, midpoint, [nv, nv, nv]), epsilon)
       for (let i = 0; i < nv; i++) {
-        const poly = poly3.fromPoints([midpoint, polygon.vertices[i], polygon.vertices[(i + 1) % nv]])
+        const poly = poly3.create([midpoint, polygon.vertices[i], polygon.vertices[(i + 1) % nv]])
         if (polygon.color) poly.color = polygon.color
         triangles.push(poly)
       }
       return
     }
     // exactly 4 vertices, use simple triangulation
-    const poly0 = poly3.fromPoints([polygon.vertices[0], polygon.vertices[1], polygon.vertices[2]])
-    const poly1 = poly3.fromPoints([polygon.vertices[0], polygon.vertices[2], polygon.vertices[3]])
+    const poly0 = poly3.create([polygon.vertices[0], polygon.vertices[1], polygon.vertices[2]])
+    const poly1 = poly3.create([polygon.vertices[0], polygon.vertices[2], polygon.vertices[3]])
     if (polygon.color) {
       poly0.color = polygon.color
       poly1.color = polygon.color

package/src/operations/transforms/center.js

@@ -39,7 +39,7 @@ const center = (options, ...objects) => {
   const defaults = {
     axes: [true, true, true],
     relativeTo: [0, 0, 0]
-  // TODO : Add addition 'methods' of centering; midpoint, centeriod
+  // TODO: Add additional 'methods' of centering: midpoint, centroid
   }
   const { axes, relativeTo } = Object.assign({}, defaults, options)
 

package/src/primitives/cuboid.js

@@ -43,7 +43,7 @@ const cuboid = (options) => {
         ]
         return pos
       })
-      return poly3.fromPoints(points)
+      return poly3.create(points)
     })
   )
   return result

package/src/primitives/cylinderElliptic.js

@@ -90,7 +90,7 @@ const cylinderElliptic = (options) => {
   // adjust the points to center
   const fromPoints = (...points) => {
     const newpoints = points.map((point) => vec3.add(vec3.create(), point, center))
-    return poly3.fromPoints(newpoints)
+    return poly3.create(newpoints)
   }
 
   const polygons = []

package/src/primitives/ellipsoid.js

@@ -66,7 +66,7 @@ const ellipsoid = (options) => {
           point = vec3.subtract(vec3.create(), vec3.scale(p1, prevcylinderpoint, cospitch), vec3.scale(p2, zvector, sinpitch))
           points.push(vec3.add(point, point, center))
 
-          polygons.push(poly3.fromPoints(points))
+          polygons.push(poly3.create(points))
 
           points = []
           point = vec3.add(vec3.create(), vec3.scale(p1, prevcylinderpoint, prevcospitch), vec3.scale(p2, zvector, prevsinpitch))
@@ -81,7 +81,7 @@ const ellipsoid = (options) => {
           points.push(vec3.add(vec3.create(), center, point))
           points.reverse()
 
-          polygons.push(poly3.fromPoints(points))
+          polygons.push(poly3.create(points))
         }
         prevcospitch = cospitch
         prevsinpitch = sinpitch

package/src/primitives/polyhedron.js

@@ -60,7 +60,7 @@ const polyhedron = (options) => {
   }
 
   const polygons = faces.map((face, findex) => {
-    const polygon = poly3.fromPoints(face.map((pindex) => points[pindex]))
+    const polygon = poly3.create(face.map((pindex) => points[pindex]))
     if (colors && colors[findex]) polygon.color = colors[findex]
     return polygon
   })

package/src/primitives/roundedCuboid.js

@@ -57,10 +57,10 @@ const stitchCorners = (previousCorners, currentCorners) => {
     const previous = previousCorners[i]
     const current = currentCorners[i]
     for (let j = 0; j < (previous.length - 1); j++) {
-      polygons.push(poly3.fromPoints([previous[j], previous[j + 1], current[j]]))
+      polygons.push(poly3.create([previous[j], previous[j + 1], current[j]]))
 
       if (j < (current.length - 1)) {
-        polygons.push(poly3.fromPoints([current[j], previous[j + 1], current[j + 1]]))
+        polygons.push(poly3.create([current[j], previous[j + 1], current[j + 1]]))
       }
     }
   }
@@ -81,7 +81,7 @@ const stitchWalls = (previousCorners, currentCorners) => {
     const p1 = previous[0]
     const c1 = current[0]
 
-    polygons.push(poly3.fromPoints([p0, p1, c1, c0]))
+    polygons.push(poly3.create([p0, p1, c1, c0]))
   }
   return polygons
 }
@@ -104,7 +104,7 @@ const stitchSides = (bottomCorners, topCorners) => {
   const polygons = []
   for (let i = 0; i < topPoints.length; i++) {
     const j = (i + 1) % topPoints.length
-    polygons.push(poly3.fromPoints([bottomPoints[i], bottomPoints[j], topPoints[j], topPoints[i]]))
+    polygons.push(poly3.create([bottomPoints[i], bottomPoints[j], topPoints[j], topPoints[i]]))
   }
   return polygons
 }
@@ -168,10 +168,10 @@ const roundedCuboid = (options) => {
     if (slice === segments) {
       // add the top
       let points = cornersPos.map((corner) => corner[0])
-      polygons.push(poly3.fromPoints(points))
+      polygons.push(poly3.create(points))
       // add the bottom
       points = cornersNeg.map((corner) => corner[0])
-      polygons.push(poly3.fromPoints(points))
+      polygons.push(poly3.create(points))
     }
 
     prevCornersPos = cornersPos

package/src/primitives/roundedCylinder.js

@@ -66,7 +66,7 @@ const roundedCylinder = (options) => {
   const fromPoints = (points) => {
     // adjust the points to center
     const newpoints = points.map((point) => vec3.add(point, point, center))
-    return poly3.fromPoints(newpoints)
+    return poly3.create(newpoints)
   }
 
   const polygons = []

package/src/primitives/torus.test.js

@@ -10,24 +10,24 @@ test('torus (defaults)', (t) => {
   const obs = torus()
   const pts = geom3.toPoints(obs)
 
-  t.notThrows.skip(() => geom3.validate(obs))
+  t.notThrows(() => geom3.validate(obs))
   t.is(pts.length, 2048) // 32 * 32 * 2 (polys/segment) = 2048
 
   const bounds = measureBoundingBox(obs)
   const expectedBounds = [[-5, -5, -1], [5, 5, 1]]
-  t.notThrows.skip(() => geom3.validate(obs))
+  t.notThrows(() => geom3.validate(obs))
   t.true(comparePoints(bounds, expectedBounds), 'Bounding box was not as expected: ' + JSON.stringify(bounds))
 })
 
-test('torus (Simple options)', (t) => {
+test('torus (simple options)', (t) => {
   const obs = torus({ innerRadius: 0.5, innerSegments: 4, outerRadius: 5, outerSegments: 8 })
   const pts = geom3.toPoints(obs)
-  t.notThrows.skip(() => geom3.validate(obs))
+  t.notThrows(() => geom3.validate(obs))
   t.is(pts.length, 64) // 4 * 8 * 2 (polys/segment) = 64
 
   const bounds = measureBoundingBox(obs)
   const expectedBounds = [[-5.5, -5.5, -0.5], [5.5, 5.5, 0.5]]
-  t.notThrows.skip(() => geom3.validate(obs))
+  t.notThrows(() => geom3.validate(obs))
   t.true(comparePoints(bounds, expectedBounds), 'Bounding box was not as expected: ' + JSON.stringify(bounds))
 })
 
@@ -48,6 +48,6 @@ test('torus (square by square)', (t) => {
 
   const bounds = measureBoundingBox(obs)
   const expectedBounds = [[-5, -5, -1], [5, 5, 1]]
-  t.notThrows.skip(() => geom3.validate(obs))
+  t.notThrows(() => geom3.validate(obs))
   t.true(comparePoints(bounds, expectedBounds), 'Bounding box was not as expected: ' + JSON.stringify(bounds))
 })

package/src/primitives/triangle.js

@@ -1,11 +1,10 @@
+const { NEPS } = require('../maths/constants')
 const vec2 = require('../maths/vec2')
 
 const geom2 = require('../geometries/geom2')
 
 const { isNumberArray } = require('./commonChecks')
 
-const NEPS = 1e-13
-
 // returns angle C
 const solveAngleFromSSS = (a, b, c) => Math.acos(((a * a) + (b * b) - (c * c)) / (2 * a * b))
 

package/src/utils/trigonometry.js

@@ -1,5 +1,4 @@
-
-const NEPS = 1e-13
+const { NEPS } = require('../maths/constants')
 
 /*
  * Returns zero if n is within epsilon of zero, otherwise return n