@woosh/meep-engine 2.109.1 → 2.109.2

package/src/core/geom/3d/tetrahedra/visualize_tetrahedral_mesh.d.ts

@@ -0,0 +1,11 @@
+/**
+ *
+ * @param {TetrahedralMesh} mesh
+ * @param {number[]} positions
+ * @param {string|number} color CSS color definition
+ * @param {number} opacity
+ * @returns {Entity}
+ */
+export function visualize_tetrahedral_mesh({ mesh, positions, color, opacity }: TetrahedralMesh): Entity;
+import Entity from "../../../../engine/ecs/Entity.js";
+//# sourceMappingURL=visualize_tetrahedral_mesh.d.ts.map

package/src/core/geom/3d/tetrahedra/visualize_tetrahedral_mesh.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"visualize_tetrahedral_mesh.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/tetrahedra/visualize_tetrahedral_mesh.js"],"names":[],"mappings":"AAQA;;;;;;;GAOG;AACH,kGAFa,MAAM,CA2BlB;mBAxCkB,kCAAkC"}

package/src/core/geom/3d/tetrahedra/visualize_tetrahedral_mesh.js

@@ -0,0 +1,42 @@
+import { LineBasicMaterial } from "three";
+import Entity from "../../../../engine/ecs/Entity.js";
+import { Transform } from "../../../../engine/ecs/transform/Transform.js";
+import { DrawMode } from "../../../../engine/graphics/ecs/mesh-v2/DrawMode.js";
+import { ShadedGeometry } from "../../../../engine/graphics/ecs/mesh-v2/ShadedGeometry.js";
+import { ShadedGeometryFlags } from "../../../../engine/graphics/ecs/mesh-v2/ShadedGeometryFlags.js";
+import { build_tetrahedral_mesh_buffer_geometry } from "./build_tetrahedral_mesh_buffer_geometry.js";
+
+/**
+ *
+ * @param {TetrahedralMesh} mesh
+ * @param {number[]} positions
+ * @param {string|number} color CSS color definition
+ * @param {number} opacity
+ * @returns {Entity}
+ */
+export function visualize_tetrahedral_mesh({
+                                               mesh, positions,
+                                               color = '#FF00FF',
+                                               opacity = 0.35
+                                           }) {
+
+    const geometry = build_tetrahedral_mesh_buffer_geometry(mesh, positions);
+
+    const sg = ShadedGeometry.from(geometry, new LineBasicMaterial({
+        color: color,
+        transparent: true,
+        opacity,
+        depthWrite: false
+    }), DrawMode.LineSegments);
+
+    sg.clearFlag(ShadedGeometryFlags.CastShadow);
+    sg.clearFlag(ShadedGeometryFlags.ReceiveShadow);
+
+    const entity = new Entity();
+
+    entity
+        .add(sg)
+        .add(new Transform());
+
+    return entity;
+}

package/src/core/geom/3d/topology/samples/sampleFloodFill.js

@@ -21,7 +21,7 @@ import { seededRandom } from "../../../../math/random/seededRandom.js";
 import { v3_dot } from "../../../vec3/v3_dot.js";
 import { compute_aabb_from_points } from "../../aabb/compute_aabb_from_points.js";
 import { SurfacePoint3 } from "../../SurfacePoint3.js";
-import { computeTriangleSurfaceArea } from "../../triangle/computeTriangleSurfaceArea.js";
+import { compute_triangle_area_3d } from "../../triangle/compute_triangle_area_3d.js";
 import { expandConnectivityByLocality } from "../expandConnectivityByLocality.js";
 import { TopoMesh } from "../struct/TopoMesh.js";
 import { three_buffer_geometry_to_topo_mesh } from "../three_buffer_geometry_to_topo_mesh.js";
@@ -226,7 +226,7 @@ async function main(engine) {
                     filter(current, prev) {
                         current.computeNormal();
 
-                        const area = computeTriangleSurfaceArea(
+                        const area = compute_triangle_area_3d(
                             current.vertices[0].x,
                             current.vertices[0].z,
                             current.vertices[0].y,

package/src/core/geom/3d/topology/tm_face_area.js

@@ -1,4 +1,4 @@
-import { computeTriangleSurfaceArea } from "../triangle/computeTriangleSurfaceArea.js";
+import { compute_triangle_area_3d } from "../triangle/compute_triangle_area_3d.js";
 
 /**
  *
@@ -10,7 +10,7 @@ export function tm_face_area(face) {
     const v1 = face.vertices[1];
     const v2 = face.vertices[2];
 
-    return computeTriangleSurfaceArea(
+    return compute_triangle_area_3d(
         v0.x, v0.y, v0.z,
         v1.x, v1.y, v1.z,
         v2.x, v2.y, v2.z

package/src/core/geom/3d/triangle/{computeTriangleSurfaceArea.d.ts → compute_triangle_area_3d.d.ts}

@@ -11,5 +11,5 @@
  * @param {number} z2
  * @returns {number}
  */
-export function computeTriangleSurfaceArea(x0: number, y0: number, z0: number, x1: number, y1: number, z1: number, x2: number, y2: number, z2: number): number;
-//# sourceMappingURL=computeTriangleSurfaceArea.d.ts.map
+export function compute_triangle_area_3d(x0: number, y0: number, z0: number, x1: number, y1: number, z1: number, x2: number, y2: number, z2: number): number;
+//# sourceMappingURL=compute_triangle_area_3d.d.ts.map

package/src/core/geom/3d/triangle/compute_triangle_area_3d.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"compute_triangle_area_3d.d.ts","sourceRoot":"","sources":["../../../../../../src/core/geom/3d/triangle/compute_triangle_area_3d.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;;GAYG;AACH,6CAXW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACJ,MAAM,CAyBlB"}

package/src/core/geom/3d/triangle/{computeTriangleSurfaceArea.js → compute_triangle_area_3d.js}

@@ -11,7 +11,7 @@
  * @param {number} z2
  * @returns {number}
  */
-export function computeTriangleSurfaceArea(
+export function compute_triangle_area_3d(
     x0, y0, z0,
     x1, y1, z1,
     x2, y2, z2
@@ -33,7 +33,5 @@ export function computeTriangleSurfaceArea(
     //area is equal to half-magnitude of the cross-product
     const magnitude = Math.sqrt(x * x + y * y + z * z);
 
-    const area = magnitude / 2;
-
-    return area;
+    return magnitude * 0.5;
 }

package/src/core/geom/Quaternion.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Quaternion.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Quaternion.js"],"names":[],"mappings":";AAqBA;IAk8CI;;;;OAIG;IACH,kCAFa,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,0BALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,6BALW,UAAU,QACV,UAAU,MACV,UAAU,aACV,MAAM,QAkBhB;IAp/CD;;;;;;;OAOG;IACH,gBANW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,EA+BhB;IA1BG;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV;;;OAGG;IACH,oBAFU,OAAO,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAEnD;IAqBjC,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IA+BD;;;;;;;;OAQG;IACH,kBAPW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,QAiDhB;IAED;;;;;OAKG;IACH,uBAJW,OAAO,QACP,OAAO,QAMjB;IAED;;;;OAIG;IACH,WAHW,UAAU,GACT,MAAM,CAQjB;IAED;;;OAGG;IACH,mBAFW,UAAU,QAKpB;IAED;;OAEG;IACH,eAqBC;IAED;;;;OAIG;IACH,eAHW,UAAU,GACT,MAAM,CAoBjB;IAED;;;;OAIG;IACH,oBAHW,OAAO,SACP,MAAM,QAIhB;IAED;;;;;;OAMG;IACH,mBALW,MAAM,MACN,MAAM,MACN,MAAM,SACN,MAAM,QAwBhB;IAED;;;;;OAKG;IACH,2BAJW,OAAO,SACP,UAAU,SACV,UAAU,QA+BpB;IAED;;;;OAIG;IACH,wBAHW,OAAO,GACL,MAAM,CAYlB;IAED;;;;OAIG;IACH,kBAHW,OAAO,GACL,MAAM,CAkBlB;IAED,kBAWC;IAED;;;;OAIG;IACH,oBAHW,MAAM,GACL,UAAU,CASrB;IAED;;OAEG;IACH,gBAFW,UAAU,QAIpB;IAED;;;;OAIG;IACH,2BAHW,UAAU,UACV,UAAU,QAcpB;IAED;;;;;;;;;;;OAWG;IACH,yBAVW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACL,UAAU,CAYrB;IAED;;;OAGG;IACH,UAFY,MAAM,CASjB;IAED;;;;OAIG;IACH,qBAHW,UAAU,aACV,MAAM,QAIhB;IAED;;;;;OAKG;IACH,4BAJW,OAAO,OACP,OAAO,YACP,OAAO,QAuBjB;IAED;;;;OAIG;IACH,eAHW,OAAO,UACP,OAAO,QAUjB;IAED;;;OAGG;IACH,wBAFsB,MAAM,QAM3B;IAED;;;;;;;OAOG;IACH,kBANW,MAAM,KACN,MAAM,KACN,MAAM,mBAEJ,UAAU,CAmCtB;IAED;;;;OAIG;IACH,yBAFW,OAAO,QA4BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA2BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA0BjB;IAED;;;;;;OAMG;IACH,mBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAIhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;OAIG;IACH,yBAHW,OAAO,MACP,OAAO,QA+CjB;IAED;;;;OAIG;IACH,uBAHW,OAAO,MACP,OAAO,QAiEjB;IAED;;;;OAIG;IACH,sBAHW,OAAO,MACP,OAAO,QAwCjB;IAED;;;OAGG;IACH,wCASC;IAED;;;;;;;;;;;;;;OAcG;IACH,6BAXW,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,GACJ,UAAU,CA4EtB;IAED;;;;OAIG;IACH,YAHW,UAAU,KACV,MAAM,QAMhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,UACV,UAAU,KACV,MAAM,QAgBhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,MACV,UAAU,KACV,MAAM,QAiDhB;IAGD;;;;OAIG;IACH,aAHW,UAAU,KACV,MAAM,QAIhB;IAED;;;;;;OAMG;IACH,qCAJW,UAAU,KACV,UAAU,KACV,MAAM,QAIhB;IAED;;;OAGG;IACH,iCAIC;IAED;;;;OAIG;IACH,YAHW,UAAU,GACR,UAAU,CAItB;IAED;;;OAGG;IACH,SAFa,UAAU,CAQtB;IAED;;;;;;;OAOG;IACH,OANW,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CA+BtB;IAED;;;OAGG;IACH,aAFa,UAAU,CAItB;IAED;;;;;MAOC;IAED,yBAEC;IAED;;;OAGG;IACH,2CAKC;IAED;;;OAGG;IACH,6CAOC;IAED;;;OAGG;IACH,kDAKC;IAED;;;OAGG;IACH,oDAOC;IAED;;;;OAIG;IACH,wBAFW,MAAM,QAIhB;IAED;;;;OAIG;IACH,kBAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAShB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAOhB;IAED,oBAEC;IAED;;;;OAIG;IACH,cAHW,UAAU,GACR,OAAO,CAQnB;IAED;;;OAGG;IACH,QAFa,MAAM,CAQlB;IAED;;;;;OAKG;IACH,qBAJW,UAAU,cACV,MAAM,GACL,OAAO,CAIlB;IAED;;;;;;;;;OASG;IACH,uBAKC;IAED;;;;OAIG;IACH,sBAHsB,MAAM,GAChB,UAAU,CAkBrB;IAED,mBAEC;IAwDL,mBAjKe,MAAM,EAAE,WACR,MAAM,UAgKS;IAC9B,iBApJe,MAAM,EAAE,WACR,MAAM,UAmJO;IAh7CxB;;OAEG;IACH,sDAOC;CAm6CJ;;kBAOS,UAAU;;mBAjhDD,4BAA4B;oBAU3B,cAAc"}
+{"version":3,"file":"Quaternion.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Quaternion.js"],"names":[],"mappings":";AAsBA;IA+2CI;;;;OAIG;IACH,kCAFa,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,0BALW,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CAQtB;IAED;;;;;;OAMG;IACH,6BALW,UAAU,QACV,UAAU,MACV,UAAU,aACV,MAAM,QAkBhB;IAj6CD;;;;;;;OAOG;IACH,gBANW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,EA+BhB;IA1BG;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IACV;;;OAGG;IACH,GAFU,MAAM,CAEN;IAEV;;;OAGG;IACH,oBAFU,OAAO,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,EAAE,MAAM,CAAC,CAEnD;IAqBjC,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IAkBD,mBAEC;IAlBD,gBAEC;IA+BD;;;;;;;;OAQG;IACH,kBAPW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,QAiDhB;IAED;;;;;OAKG;IACH,uBAJW,OAAO,QACP,OAAO,QAMjB;IAED;;;;OAIG;IACH,WAHW,UAAU,GACT,MAAM,CAQjB;IAED;;;OAGG;IACH,mBAFW,UAAU,QAKpB;IAED;;OAEG;IACH,eAqBC;IAED;;;;OAIG;IACH,eAHW,UAAU,GACT,MAAM,CAoBjB;IAED;;;;OAIG;IACH,oBAHW,OAAO,SACP,MAAM,QAIhB;IAED;;;;;;OAMG;IACH,mBALW,MAAM,MACN,MAAM,MACN,MAAM,SACN,MAAM,QAwBhB;IAED;;;;;OAKG;IACH,2BAJW,OAAO,SACP,UAAU,SACV,UAAU,QA+BpB;IAED;;;;OAIG;IACH,wBAHW,OAAO,GACL,MAAM,CAYlB;IAED;;;;OAIG;IACH,kBAHW,OAAO,GACL,MAAM,CAkBlB;IAED,kBAWC;IAED;;;;OAIG;IACH,oBAHW,MAAM,GACL,UAAU,CASrB;IAED;;OAEG;IACH,gBAFW,UAAU,QAIpB;IAED;;;;OAIG;IACH,2BAHW,UAAU,UACV,UAAU,QAcpB;IAED;;;;;;;;;;;OAWG;IACH,yBAVW,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,MACN,MAAM,GACL,UAAU,CAYrB;IAED;;;OAGG;IACH,UAFY,MAAM,CASjB;IAED;;;;OAIG;IACH,qBAHW,UAAU,aACV,MAAM,QAIhB;IAED;;;;;OAKG;IACH,eAJW,OAAO,UACP,OAAO,OACP,OAAO,QAUjB;IAED;;;OAGG;IACH,wBAFsB,MAAM,QAM3B;IAED;;;;;;;OAOG;IACH,kBANW,MAAM,KACN,MAAM,KACN,MAAM,mBAEJ,UAAU,CAmCtB;IAED;;;;OAIG;IACH,yBAFW,OAAO,QA4BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA2BjB;IAED;;;OAGG;IACH,yBAFW,OAAO,QA0BjB;IAED;;;;;;OAMG;IACH,mBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAIhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;;;;;OAQG;IACH,sBAJW,MAAM,KACN,MAAM,KACN,MAAM,QAqBhB;IAED;;;;OAIG;IACH,yBAHW,OAAO,MACP,OAAO,QA+CjB;IAED;;;;;OAKG;IACH,sBAHW,OAAO,MACP,OAAO,QAiEjB;IAED;;;OAGG;IACH,wCASC;IAED;;;;;;;;;;;;;;OAcG;IACH,6BAXW,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,OACN,MAAM,GACJ,UAAU,CA4EtB;IAED;;;;OAIG;IACH,YAHW,UAAU,KACV,MAAM,QAMhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,UACV,UAAU,KACV,MAAM,QAgBhB;IAED;;;;;OAKG;IACH,uBAJW,UAAU,MACV,UAAU,KACV,MAAM,QAiDhB;IAGD;;;;OAIG;IACH,aAHW,UAAU,KACV,MAAM,QAIhB;IAED;;;OAGG;IACH,iCAIC;IAED;;;;OAIG;IACH,YAHW,UAAU,GACR,UAAU,CAItB;IAED;;;OAGG;IACH,SAFa,UAAU,CAQtB;IAED;;;;;;;OAOG;IACH,OANW,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,GACJ,UAAU,CA+BtB;IAED;;;OAGG;IACH,aAFa,UAAU,CAItB;IAED;;;;;MAOC;IAED,yBAEC;IAED;;;OAGG;IACH,2CAKC;IAED;;;OAGG;IACH,6CAOC;IAED;;;OAGG;IACH,kDAKC;IAED;;;OAGG;IACH,oDAOC;IAED;;;;OAIG;IACH,wBAFW,MAAM,QAIhB;IAED;;;;OAIG;IACH,kBAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAShB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAOhB;IAED,oBAEC;IAED;;;;OAIG;IACH,cAHW,UAAU,GACR,OAAO,CAQnB;IAED;;;OAGG;IACH,QAFa,MAAM,CAQlB;IAED;;;;;OAKG;IACH,qBAJW,UAAU,cACV,MAAM,GACL,OAAO,CAIlB;IAED;;;;;;;;;OASG;IACH,uBAKC;IAED;;;;OAIG;IACH,sBAHsB,MAAM,GAChB,UAAU,CAkBrB;IAED,mBAEC;IAwDL,mBAjKe,MAAM,EAAE,WACR,MAAM,UAgKS;IAC9B,iBApJe,MAAM,EAAE,WACR,MAAM,UAmJO;IA71CxB;;OAEG;IACH,sDAOC;CAg1CJ;;kBAOS,UAAU;;mBA/7CD,4BAA4B;oBAW3B,cAAc"}

package/src/core/geom/Quaternion.js

@@ -9,6 +9,7 @@ import { computeHashFloat } from "../primitives/numbers/computeHashFloat.js";
 import { quat_decode_from_uint32 } from "./3d/quaternion/quat_decode_from_uint32.js";
 import { quat_encode_to_uint32 } from "./3d/quaternion/quat_encode_to_uint32.js";
 import { v3_dot } from "./vec3/v3_dot.js";
+import { v4_length } from "./vec4/v4_length.js";
 import Vector3 from "./Vector3.js";
 
 
@@ -451,7 +452,7 @@ class Quaternion {
         const z = this.z;
         const w = this.w;
 
-        return Math.sqrt(x * x + y * y + z * z + w * w);
+        return v4_length(x, y, z, w);
     }
 
     /**
@@ -463,48 +464,20 @@ class Quaternion {
         Quaternion.rotateTowards(this, this, other, max_delta);
     }
 
-    /**
-     *
-     * @param {Vector3} direction
-     * @param {Vector3} [up]
-     * @param {Vector3} [forward]
-     */
-    alignToDirection(direction, up = Vector3.up, forward = Vector3.forward) {
-
-        const dot = up.dot(direction);
-
-
-        if (Math.abs(dot + 1) < 0.000001) {
-            // vector a and b point exactly in the opposite direction,
-            // so it is a 180 degrees turn around the up-axis
-            this.set(forward.x, forward.y, forward.z, Math.PI);
-        } else if (Math.abs(dot - 1) < 0.000001) {
-            // vector a and b point exactly in the same direction
-            // so we return the identity quaternion
-            this.set(0, 0, 0, 1);
-        }
-
-        const angle = Math.acos(dot);
-
-        axis.crossVectors(up, direction);
-        axis.normalize();
-
-        this.fromAxisAngle(axis, angle);
-    }
-
     /**
      *
      * @param {Vector3} source
      * @param {Vector3} target
+     * @param {Vector3} [up]
      */
-    lookAt(source, target) {
+    lookAt(source, target, up = Vector3.up) {
         const forward = new Vector3();
 
         forward.subVectors(target, source);
 
         forward.normalize();
 
-        this.alignToDirection(forward);
+        this.lookRotation(forward, up);
     }
 
     /**
@@ -887,22 +860,23 @@ class Quaternion {
         const qW_1 = qW + l0;
 
         //normalize result
-        const l1 = Math.sqrt(p0 + qW_1 * qW_1);
+        const l1 = 1 / Math.sqrt(p0 + qW_1 * qW_1);
 
-        const x = qX / l1;
-        const y = qY / l1;
-        const z = qZ / l1;
-        const w = qW / l1;
+        const x = qX * l1;
+        const y = qY * l1;
+        const z = qZ * l1;
+        const w = qW * l1;
 
         this.set(x, y, z, w);
     }
 
     /**
+     * NOTE: Vectors need to be normalized
      * Based on blog post: http://lolengine.net/blog/2013/09/18/beautiful-maths-quaternion-from-vectors
      * @param {Vector3} from
      * @param {Vector3} to
      */
-    fromUnitVectors2(from, to) {
+    fromUnitVectors(from, to) {
         assert.ok(from.isNormalized(), `from vector is not normalized, length = ${from.length()}`);
         assert.ok(to.isNormalized(), `to vector is not normalized, length = ${to.length()}`);
 
@@ -923,7 +897,7 @@ class Quaternion {
 
         const uv_dot = v3_dot(ax, ay, az, bx, by, bz);
 
-        if (uv_dot === -1) {
+        if (uv_dot < -0.9999999) {
             //to vector is opposite, produce a reversal quaternion
 
             tempvec3.crossVectors(Vector3.left, from);
@@ -967,51 +941,6 @@ class Quaternion {
         );
     }
 
-    /**
-     * Input vectors must be normalized
-     * @param {Vector3} from
-     * @param {Vector3} to
-     */
-    fromUnitVectors(from, to) {
-        assert.ok(from.isNormalized(), `from vector is not normalized, length = ${from.length()}`);
-        assert.ok(to.isNormalized(), `to vector is not normalized, length = ${to.length()}`);
-
-        const dot = from.dot(to);
-
-        let x, y, z, w;
-
-        if (dot < -0.9999999) {
-            const tempvec3 = new Vector3();
-
-            tempvec3.crossVectors(Vector3.left, from);
-
-            if (tempvec3.lengthSqr() < 0.00001) {
-                tempvec3.crossVectors(Vector3.up, from);
-            }
-
-            tempvec3.normalize();
-
-            this.fromAxisAngle(tempvec3, Math.PI);
-        } else if (dot > 0.9999999) {
-            x = 0;
-            y = 0;
-            z = 0;
-            w = 1;
-
-            this.set(x, y, z, w);
-        } else {
-            //cross product from x to
-
-            x = from.y * to.z - from.z * to.y;
-            y = from.z * to.x - from.x * to.z;
-            z = from.x * to.y - from.y * to.x;
-            w = 1 + dot;
-
-            this.set(x, y, z, w);
-            this.normalize();
-        }
-    }
-
     /**
      * @deprecated
      * @param {Matrix4} m
@@ -1216,17 +1145,6 @@ class Quaternion {
         this.slerpQuaternions(this, other, t);
     }
 
-    /**
-     * @see https://github.com/gareth-cross/quat/blob/master/include/quaternion.hpp
-     * TODO implement
-     * @param {Quaternion} a
-     * @param {Quaternion} b
-     * @param {number} t
-     */
-    quaternionsInterpolateRungeKutta4(a, b, t) {
-        throw new Error('Not Implemented');
-    }
-
     /**
      *
      * @param {function} handler
@@ -1557,13 +1475,6 @@ Quaternion.prototype.toArray = Quaternion.prototype.writeToArray;
 Quaternion.identity = Object.freeze(new Quaternion(0, 0, 0, 1));
 
 
-/**
- *
- * @type {Vector3}
- */
-const axis = new Vector3();
-
-
 const tempvec3 = new Vector3();
 
 export default Quaternion;

package/src/core/geom/Quaternion.spec.js

@@ -70,14 +70,14 @@ test('multiply quaternions', () => {
     });
 });
 
-test('fromUnitVectors2', () => {
+test('fromUnitVectors', () => {
     const q = new Quaternion();
 
     const v = new Vector3();
 
     function check(v0, v1) {
 
-        q.fromUnitVectors2(v0, v1);
+        q.fromUnitVectors(v0, v1);
 
         v.copy(v0);
 

package/src/core/geom/Vector1.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"Vector1.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Vector1.js"],"names":[],"mappings":";AAWA;IA2PI;;;;;OAKG;IACH,kBAJW,OAAO,KACP,OAAO,GACN,MAAM,CAIjB;IAlQD;;;;;;OAMG;IACH,gBALW,MAAM,EAchB;IAHG,UAAU;IAEV,0DAA6B;IAWjC;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,iBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;OAIG;IACH,OAHW,MAAM,GACJ,OAAO,CAgBnB;IAED;;;OAGG;IACH,aAFW,MAAM,QAOhB;IAED;;;OAGG;IACH,UAFY,OAAO,CAIlB;IAED,kBAEC;IAED,kBAEC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;OAGG;IACH,WAFW,OAAO,8BAAwB,WAIzC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,gBAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,kBAFW,MAAM,QAIhB;IAED;;;;;OAKG;IACH,WAJW,MAAM,QACN,MAAM,GACJ,OAAO,CAInB;IAED;;;OAGG;IACH,YAFW,OAAO,8BAAwB,QAIzC;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;OAEG;IACH,QAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,4BAGC;IAED,iBAEC;IAED,yBAEC;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED,oBAEC;IAED;;;OAGG;IACH,2CAEC;IAED;;;OAGG;IACH,6CAIC;IAgBD,mBAEC;IAND,gBAEC;IAaL;;;OAGG;IACH,oBAFU,OAAO,CAEU;IAXvB,sDAIC;CACJ;;;cAYS,OAAO;aAMP,OAAO;;mBAzSE,4BAA4B"}
+{"version":3,"file":"Vector1.d.ts","sourceRoot":"","sources":["../../../../src/core/geom/Vector1.js"],"names":[],"mappings":";AAWA;IAmQI;;;;;OAKG;IACH,kBAJW,OAAO,KACP,OAAO,GACN,MAAM,CAIjB;IA1QD;;;;;;OAMG;IACH,gBALW,MAAM,EAchB;IAHG,UAAU;IAEV,0DAA6B;IAWjC;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,YAFa,MAAM,CAIlB;IAED;;;;OAIG;IACH,iBAHW,OAAO,GACL,MAAM,CAIlB;IAED;;;;OAIG;IACH,OAHW,MAAM,GACJ,OAAO,CAgBnB;IAED;;;OAGG;IACH,aAFW,MAAM,QAOhB;IAED;;;OAGG;IACH,UAFY,OAAO,CAIlB;IAED,kBAEC;IAED,kBAEC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;OAGG;IACH,WAFW,OAAO,8BAAwB,WAIzC;IAED;;;;OAIG;IACH,iBAFY,OAAO,CAIlB;IAED;;;;OAIG;IACH,WAHW,OAAO,GACL,OAAO,CAInB;IAED;;;;OAIG;IACH,gBAHW,OAAO,GACL,OAAO,CAInB;IAED;;;OAGG;IACH,kBAFW,MAAM,QAIhB;IAED;;;;;OAKG;IACH,WAJW,MAAM,QACN,MAAM,GACJ,OAAO,CAInB;IAED;;;OAGG;IACH,eAEC;IAED;;;OAGG;IACH,YAFW,OAAO,8BAAwB,QAIzC;IAED;;;OAGG;IACH,SAFa,OAAO,CAInB;IAED;;;;OAIG;IACH,cAHW,OAAO,GACL,OAAO,CAInB;IAED;;OAEG;IACH,QAFa,MAAM,CAIlB;IAED;;;OAGG;IACH,4BAGC;IAED,iBAEC;IAED,yBAEC;IAED;;;;OAIG;IACH,qBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED;;;;OAIG;IACH,oBAHW,MAAM,EAAE,WACR,MAAM,QAIhB;IAED,oBAEC;IAED;;;OAGG;IACH,2CAEC;IAED;;;OAGG;IACH,6CAIC;IAgBD,mBAEC;IAND,gBAEC;IAaL;;;OAGG;IACH,oBAFU,OAAO,CAEU;IAXvB,sDAIC;CACJ;;;cAYS,OAAO;aAMP,OAAO;;mBAjTE,4BAA4B"}

package/src/core/geom/Vector1.js

@@ -20,7 +20,7 @@ class Vector1 extends Number {
     constructor(x = 0) {
         super();
 
-        assert.isNumber( x, 'x');
+        assert.isNumber(x, 'x');
         assert.notNaN(x, 'x');
 
         this.x = x;
@@ -67,7 +67,7 @@ class Vector1 extends Number {
      * @returns {Vector1}
      */
     set(x) {
-        assert.isNumber( x,'x');
+        assert.isNumber(x, 'x');
         assert.notNaN(x, 'x');
 
         const oldValue = this.x;
@@ -87,7 +87,7 @@ class Vector1 extends Number {
      * @param {number} x
      */
     setSilent(x) {
-        assert.isNumber( x,'x');
+        assert.isNumber(x, 'x');
         assert.notNaN(x, 'x');
 
         this.x = x;
@@ -171,6 +171,14 @@ class Vector1 extends Number {
         return this.set(clamp(this.x, low, high));
     }
 
+    /**
+     * Negate sign of the value, if value is 7, it will become -7 and vice versa
+     * Zero is unaffected
+     */
+    negate() {
+        this.set(-this.x);
+    }
+
     /**
      *
      * @param {Vector1|Vector2|Vector3|Vector4} other
@@ -225,7 +233,7 @@ class Vector1 extends Number {
      * @param {number[]} array
      * @param {number} offset
      */
-    readFromArray(array, offset=0) {
+    readFromArray(array, offset = 0) {
         this.set(array[offset]);
     }
 
@@ -234,11 +242,11 @@ class Vector1 extends Number {
      * @param {number[]} array
      * @param {number} offset
      */
-    writeToArray(array, offset=0) {
+    writeToArray(array, offset = 0) {
         array[offset] = this.x;
     }
 
-    asArray(){
+    asArray() {
         return [this.x];
     }
 

package/src/core/geom/vec3/v3_array_displace_in_direction.d.ts

@@ -0,0 +1,14 @@
+/**
+ *
+ * @param {number[]|Float32Array} output
+ * @param {number} output_offset
+ * @param {number} origin_x
+ * @param {number} origin_y
+ * @param {number} origin_z
+ * @param {number} direction_x
+ * @param {number} direction_y
+ * @param {number} direction_z
+ * @param {number} distance
+ */
+export function v3_array_displace_in_direction(output: number[] | Float32Array, output_offset: number, origin_x: number, origin_y: number, origin_z: number, direction_x: number, direction_y: number, direction_z: number, distance: number): void;
+//# sourceMappingURL=v3_array_displace_in_direction.d.ts.map

package/src/core/geom/vec3/v3_array_displace_in_direction.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v3_array_displace_in_direction.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_array_displace_in_direction.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,uDAVW,MAAM,EAAE,GAAC,YAAY,iBACrB,MAAM,YACN,MAAM,YACN,MAAM,YACN,MAAM,eACN,MAAM,eACN,MAAM,eACN,MAAM,YACN,MAAM,QAYhB"}

package/src/core/geom/vec3/{v3_shift_along_direction.js → v3_array_displace_in_direction.js}

@@ -10,7 +10,7 @@
  * @param {number} direction_z
  * @param {number} distance
  */
-export function v3_shift_along_direction(
+export function v3_array_displace_in_direction(
     output, output_offset,
     origin_x, origin_y, origin_z,
     direction_x, direction_y, direction_z,

package/src/core/geom/vec3/v3_displace_in_direction.d.ts

@@ -2,13 +2,13 @@
  *
  * @param {Vector3} result
  * @param {number} distance
- * @param {number} sourceX
- * @param {number} sourceY
- * @param {number} sourceZ
- * @param {number} directionX
- * @param {number} directionY
- * @param {number} directionZ
+ * @param {number} origin_x
+ * @param {number} source_y
+ * @param {number} source_z
+ * @param {number} direction_x
+ * @param {number} direction_y
+ * @param {number} direction_z
  * @returns {boolean} true if displaced successfully, false if direction had 0 magnitude
  */
-export function v3_displace_in_direction(result: Vector3, distance: number, sourceX: number, sourceY: number, sourceZ: number, directionX: number, directionY: number, directionZ: number): boolean;
+export function v3_displace_in_direction(result: Vector3, distance: number, origin_x: number, source_y: number, source_z: number, direction_x: number, direction_y: number, direction_z: number): boolean;
 //# sourceMappingURL=v3_displace_in_direction.d.ts.map

package/src/core/geom/vec3/v3_displace_in_direction.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"v3_displace_in_direction.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_displace_in_direction.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,oEATW,MAAM,WACN,MAAM,WACN,MAAM,WACN,MAAM,cACN,MAAM,cACN,MAAM,cACN,MAAM,GACJ,OAAO,CAgCnB"}
+{"version":3,"file":"v3_displace_in_direction.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_displace_in_direction.js"],"names":[],"mappings":"AAAA;;;;;;;;;;;GAWG;AACH,oEATW,MAAM,YACN,MAAM,YACN,MAAM,YACN,MAAM,eACN,MAAM,eACN,MAAM,eACN,MAAM,GACJ,OAAO,CAmCnB"}

package/src/core/geom/vec3/v3_displace_in_direction.js

@@ -2,40 +2,43 @@
  *
  * @param {Vector3} result
  * @param {number} distance
- * @param {number} sourceX
- * @param {number} sourceY
- * @param {number} sourceZ
- * @param {number} directionX
- * @param {number} directionY
- * @param {number} directionZ
+ * @param {number} origin_x
+ * @param {number} source_y
+ * @param {number} source_z
+ * @param {number} direction_x
+ * @param {number} direction_y
+ * @param {number} direction_z
  * @returns {boolean} true if displaced successfully, false if direction had 0 magnitude
  */
 export function v3_displace_in_direction(
     result,
     distance,
-    sourceX, sourceY, sourceZ,
-    directionX, directionY, directionZ
+    origin_x, source_y, source_z,
+    direction_x, direction_y, direction_z
 ) {
+    
     if (distance === 0) {
-        result.set(sourceX, sourceY, sourceZ);
+        result.set(origin_x, source_y, source_z);
         return true;
     }
 
-    //scale direction to have length equal to distance
+    // scale direction to have length equal to distance
 
-    const directionLength = Math.sqrt(directionX * directionX + directionY * directionY + directionZ * directionZ);
+    const direction_length_sqr = direction_x * direction_x + direction_y * direction_y + direction_z * direction_z;
 
-    if(directionLength === 0){
+    if (direction_length_sqr === 0) {
         // direction vector is zero length, avoid division by zero, clamp to source instead
-        result.set(sourceX, sourceY, sourceZ);
+        result.set(origin_x, source_y, source_z);
         return false;
     }
 
-    const m = distance / directionLength;
+    const direction_length = Math.sqrt(direction_length_sqr);
 
-    const x = sourceX + directionX * m;
-    const y = sourceY + directionY * m;
-    const z = sourceZ + directionZ * m;
+    const m = distance / direction_length;
+
+    const x = origin_x + direction_x * m;
+    const y = source_y + direction_y * m;
+    const z = source_z + direction_z * m;
 
     result.set(x, y, z);
 

package/src/core/geom/vec3/v3_displace_in_direction.spec.d.ts

@@ -0,0 +1,2 @@
+export {};
+//# sourceMappingURL=v3_displace_in_direction.spec.d.ts.map

package/src/core/geom/vec3/v3_displace_in_direction.spec.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v3_displace_in_direction.spec.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec3/v3_displace_in_direction.spec.js"],"names":[],"mappings":""}

package/src/core/geom/vec3/v3_displace_in_direction.spec.js

@@ -0,0 +1,80 @@
+import Vector3 from "../Vector3.js";
+import { v3_displace_in_direction } from "./v3_displace_in_direction.js";
+
+test("zero distance", () => {
+
+    const r = new Vector3();
+
+    expect(v3_displace_in_direction(
+        r, 0,
+        0, 0, 0,
+        0, 0, 1
+    )).toBe(true);
+
+});
+
+test("direction magnitude is 0", () => {
+
+    const r = new Vector3();
+
+    expect(v3_displace_in_direction(
+        r, 1,
+        0, 0, 0,
+        0, 0, 0
+    )).toBe(false);
+
+});
+
+test("cardinal directions", () => {
+
+    const r = new Vector3();
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, 0, 1
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, 2, 0]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, 0, -1
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, 2, 6]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, 1, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, -1, 3]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        0, -1, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([1, 5, 3]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        1, 0, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([-2, 2, 3]);
+
+    expect(v3_displace_in_direction(
+        r, -3,
+        1, 2, 3,
+        -1, 0, 0
+    )).toBe(true);
+
+    expect(r.asArray()).toEqual([4, 2, 3]);
+
+});

package/src/core/geom/vec4/v4_length.d.ts

@@ -0,0 +1,10 @@
+/**
+ *
+ * @param {number} x
+ * @param {number} y
+ * @param {number} z
+ * @param {number} w
+ * @returns {number}
+ */
+export function v4_length(x: number, y: number, z: number, w: number): number;
+//# sourceMappingURL=v4_length.d.ts.map

package/src/core/geom/vec4/v4_length.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"v4_length.d.ts","sourceRoot":"","sources":["../../../../../src/core/geom/vec4/v4_length.js"],"names":[],"mappings":"AAEA;;;;;;;GAOG;AACH,6BANW,MAAM,KACN,MAAM,KACN,MAAM,KACN,MAAM,GACJ,MAAM,CAIlB"}

package/src/core/geom/vec4/v4_length.js

@@ -0,0 +1,13 @@
+import { v4_length_sqr } from "./v4_length_sqr.js";
+
+/**
+ *
+ * @param {number} x
+ * @param {number} y
+ * @param {number} z
+ * @param {number} w
+ * @returns {number}
+ */
+export function v4_length(x, y, z, w) {
+    return Math.sqrt(v4_length_sqr(x, y, z, w));
+}