@flighthq/geometry 0.5.1-next.903.f434bc2 → 0.5.1-next.905.5fbf787

package/README.md

@@ -13,5 +13,5 @@ Import the supported application-facing API from `@flighthq/geometry`. The `@fli
 This package is part of the locked-version Flight SDK graph. Applications may instead install and import `@flighthq/sdk` when package-level tree shaking is sufficient.
 
 - [Flight project](https://github.com/flighthq/flight)
-- [Source for @flighthq/geometry@0.5.1-next.903.f434bc2](https://github.com/flighthq/flight/tree/f434bc233778d75aeea61995df3a8cfaec459a85/packages/geometry)
-- [License](https://github.com/flighthq/flight/blob/f434bc233778d75aeea61995df3a8cfaec459a85/LICENSE.md)
+- [Source for @flighthq/geometry@0.5.1-next.905.5fbf787](https://github.com/flighthq/flight/tree/5fbf7874064c384307542d68b4dcb6a895ff0dcf/packages/geometry)
+- [License](https://github.com/flighthq/flight/blob/5fbf7874064c384307542d68b4dcb6a895ff0dcf/LICENSE.md)

package/package.json

@@ -1,6 +1,6 @@
 {
   "name": "@flighthq/geometry",
-  "version": "0.5.1-next.903.f434bc2",
+  "version": "0.5.1-next.905.5fbf787",
   "author": "Joshua Granick and other contributors",
   "license": "MIT",
   "repository": {
@@ -38,10 +38,10 @@
     "clean:dist": "tsx ../../scripts/clean-package-dist.ts"
   },
   "dependencies": {
-    "@flighthq/entity": "0.5.1-next.903.f434bc2",
-    "@flighthq/log": "0.5.1-next.903.f434bc2",
-    "@flighthq/math": "0.5.1-next.903.f434bc2",
-    "@flighthq/types": "0.5.1-next.903.f434bc2"
+    "@flighthq/entity": "0.5.1-next.905.5fbf787",
+    "@flighthq/log": "0.5.1-next.905.5fbf787",
+    "@flighthq/math": "0.5.1-next.905.5fbf787",
+    "@flighthq/types": "0.5.1-next.905.5fbf787"
   },
   "devDependencies": {
     "typescript": "^5.3.0"

package/src/aabb.test.ts

@@ -22,6 +22,8 @@ import {
   unionAabb,
 } from '@flighthq/geometry/contract';
 
+import { initializeAabb } from './aabb';
+
 describe('cloneAabb', () => {
   it('creates an independent copy with independent corner vectors', () => {
     const a = createAabb(-1, -2, -3, 4, 5, 6);
@@ -198,6 +200,12 @@ describe('getClosestPointOnAabb', () => {
   });
 });
 
+describe('initializeAabb', () => {
+  it('is the construction initializer of createAabb', () => {
+    expect(typeof initializeAabb).toBe('function');
+  });
+});
+
 describe('intersectAabb', () => {
   it('produces the overlapping sub-box for two overlapping boxes', () => {
     const a = createAabb(0, 0, 0, 4, 4, 4);
@@ -361,7 +369,6 @@ describe('transformAabbByMatrix4', () => {
     expect(a.max.z).toBeCloseTo(6, 6);
   });
 });
-
 describe('unionAabb', () => {
   it('encloses both boxes', () => {
     const a = createAabb(0, 0, 0, 1, 1, 1);

package/src/boundingSphere.test.ts

@@ -15,6 +15,8 @@ import {
   transformBoundingSphereByMatrix4,
 } from '@flighthq/geometry/contract';
 
+import { initializeBoundingSphere } from './boundingSphere';
+
 describe('cloneBoundingSphere', () => {
   it('creates an independent copy with an independent center', () => {
     const s = createBoundingSphere(1, 2, 3, 4);
@@ -122,6 +124,12 @@ describe('getClosestPointOnBoundingSphere', () => {
   });
 });
 
+describe('initializeBoundingSphere', () => {
+  it('is the construction initializer of createBoundingSphere', () => {
+    expect(typeof initializeBoundingSphere).toBe('function');
+  });
+});
+
 describe('isBoundingSphereIntersectingBoundingSphere', () => {
   it('returns true for overlapping spheres', () => {
     const a = createBoundingSphere(0, 0, 0, 2);
@@ -242,7 +250,6 @@ describe('setBoundingSphereFromAabb', () => {
     expect(s.radius).toBe(-1);
   });
 });
-
 describe('transformBoundingSphereByMatrix4', () => {
   it('translates the center and preserves the radius', () => {
     const s = createBoundingSphere(0, 0, 0, 2);

package/src/capsule.test.ts

@@ -12,6 +12,8 @@ import {
   setCapsule,
 } from '@flighthq/geometry/contract';
 
+import { initializeCapsule } from './capsule';
+
 describe('createCapsule', () => {
   it('stores all fields', () => {
     const c = createCapsule(1, 2, 3, 4, 5, 6, 0.5);
@@ -109,6 +111,12 @@ describe('getClosestPointOnCapsule', () => {
   });
 });
 
+describe('initializeCapsule', () => {
+  it('is the construction initializer of createCapsule', () => {
+    expect(typeof initializeCapsule).toBe('function');
+  });
+});
+
 describe('intersectRay3DCapsule', () => {
   it('hits a unit capsule from the side', () => {
     // Capsule from (0,-1,0) to (0,1,0) with radius 1.
@@ -308,20 +316,6 @@ describe('isCapsuleIntersectingSphere', () => {
   });
 });
 
-describe('setCapsule', () => {
-  it('updates all fields in place', () => {
-    const c = createCapsule(0, 0, 0, 0, 0, 0, 0);
-    setCapsule(c, 1, 2, 3, 4, 5, 6, 7);
-    expect(c.startX).toBe(1);
-    expect(c.startY).toBe(2);
-    expect(c.startZ).toBe(3);
-    expect(c.endX).toBe(4);
-    expect(c.endY).toBe(5);
-    expect(c.endZ).toBe(6);
-    expect(c.radius).toBe(7);
-  });
-});
-
 // The six coordinates of a capsule axis: start x, y, z then end x, y, z.
 type SegmentArguments = [number, number, number, number, number, number];
 
@@ -347,3 +341,16 @@ function measureSegmentDistance(first: readonly number[], second: readonly numbe
   }
   return smallest;
 }
+describe('setCapsule', () => {
+  it('updates all fields in place', () => {
+    const c = createCapsule(0, 0, 0, 0, 0, 0, 0);
+    setCapsule(c, 1, 2, 3, 4, 5, 6, 7);
+    expect(c.startX).toBe(1);
+    expect(c.startY).toBe(2);
+    expect(c.startZ).toBe(3);
+    expect(c.endX).toBe(4);
+    expect(c.endY).toBe(5);
+    expect(c.endZ).toBe(6);
+    expect(c.radius).toBe(7);
+  });
+});

package/src/frustum.test.ts

@@ -18,6 +18,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Matrix4 } from '@flighthq/types/contract';
 
+import { initializeFrustum } from './frustum';
+
 // View-projection for a camera at (0,0,5) looking toward the origin down -z.
 function createTestViewProjection(): Matrix4 {
   const projection = createPerspectiveMatrix4(0.5, 1, 0.1, 100);
@@ -96,6 +98,12 @@ describe('getFrustumCorners', () => {
   });
 });
 
+describe('initializeFrustum', () => {
+  it('is the construction initializer of createFrustum', () => {
+    expect(typeof initializeFrustum).toBe('function');
+  });
+});
+
 describe('isFrustumContainingPoint', () => {
   it('contains a point in front of the camera', () => {
     const f = createFrustum();
@@ -196,7 +204,6 @@ describe('isFrustumIntersectingSphere', () => {
     expect(isFrustumIntersectingSphere(f, createBoundingSphere(0, 0, 5.5, 2))).toBe(true);
   });
 });
-
 describe('setFrustumFromMatrix4', () => {
   it('produces unit-length plane normals', () => {
     const f = createFrustum();

package/src/matrix.test.ts

@@ -43,6 +43,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Matrix, Matrix3Like, Matrix4Like } from '@flighthq/types/contract';
 
+import { initializeMatrix } from './matrix';
+
 describe('cloneMatrix', () => {
   it('should clone the matrix3x2 correctly', () => {
     const m1 = createMatrix(2, 3, 4, 5, 6, 7);
@@ -305,6 +307,12 @@ describe('equalsMatrix', () => {
   });
 });
 
+describe('initializeMatrix', () => {
+  it('is the construction initializer of createMatrix', () => {
+    expect(typeof initializeMatrix).toBe('function');
+  });
+});
+
 describe('inverseMatrix', () => {
   it('should invert the matrix3x2 correctly', () => {
     // Create a matrix3x2 with scaling of 2 and translation of (5, 3)
@@ -1180,7 +1188,6 @@ describe('translateMatrixByVectorXY', () => {
     expect(matrix).toMatchObject({ a: 2, b: 3, c: 4, d: 5, tx: 28, ty: 36 });
   });
 });
-
 describe('writeMatrixToFloat32Array', () => {
   it('writes 6 values at the offset', () => {
     const array = new Float32Array(6);

package/src/matrix3.test.ts

@@ -30,6 +30,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Matrix3 } from '@flighthq/types/contract';
 
+import { initializeMatrix3 } from './matrix3';
+
 // Matrix3 storage is column-major: element (row r, column c) lives at m[3 * c + r], matching
 // Matrix4 and the GL/GLSL/WGSL uniform ABI. A matrix-like {m} built from a raw Float32Array below
 // therefore lists its values column-by-column. getMatrix3Element(r, c) / column / row accessors are
@@ -363,6 +365,12 @@ describe('getMatrix3Element', () => {
   });
 });
 
+describe('initializeMatrix3', () => {
+  it('is the construction initializer of createMatrix3', () => {
+    expect(typeof initializeMatrix3).toBe('function');
+  });
+});
+
 describe('inverseMatrix3', () => {
   it('should invert the matrix correctly', () => {
     // A scale of 2 with translation of (5, 3): rows [2,0,5],[0,2,3],[0,0,1].
@@ -833,7 +841,6 @@ describe('transposeMatrix3', () => {
     expect(Array.from(out.m)).toEqual([1, 0, 0, 0, 1, 0, 0, 0, 1]);
   });
 });
-
 describe('writeMatrix3ToFloat32Array', () => {
   it('writes 9 values starting at the given offset', () => {
     const m = createMatrix3(1, 2, 3, 4, 5, 6, 7, 8, 9);

package/src/matrix4.test.ts

@@ -58,6 +58,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Matrix, Matrix4 } from '@flighthq/types/contract';
 
+import { initializeMatrix4 } from './matrix4';
+
 const X_AXIS = { x: 1, y: 0, z: 0, w: 0 };
 const Y_AXIS = { x: 0, y: 1, z: 0, w: 0 };
 const Z_AXIS = { x: 0, y: 0, z: 1, w: 0 };
@@ -837,6 +839,12 @@ describe('getMatrix4Position', () => {
   });
 });
 
+describe('initializeMatrix4', () => {
+  it('is the construction initializer of createMatrix4', () => {
+    expect(typeof initializeMatrix4).toBe('function');
+  });
+});
+
 describe('interpolateMatrix4', () => {
   it('returns a at t=0', () => {
     const a = createMatrix4();
@@ -2044,6 +2052,11 @@ describe('transposeMatrix4', () => {
   });
 });
 
+function expectMatrix4Close(actual: Matrix4, expected: Matrix4): void {
+  for (let i = 0; i < 16; i++) {
+    expect(actual.m[i]).toBeCloseTo(expected.m[i]);
+  }
+}
 describe('writeMatrix4ToFloat32Array', () => {
   it('writes 16 elements at offset 0', () => {
     const m = createMatrix4();
@@ -2068,9 +2081,3 @@ describe('writeMatrix4ToFloat32Array', () => {
     expect(arr[3]).toBe(0);
   });
 });
-
-function expectMatrix4Close(actual: Matrix4, expected: Matrix4): void {
-  for (let i = 0; i < 16; i++) {
-    expect(actual.m[i]).toBeCloseTo(expected.m[i]);
-  }
-}

package/src/obb.test.ts

@@ -20,6 +20,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Matrix4, Obb, Vector3 } from '@flighthq/types/contract';
 
+import { initializeObb } from './obb';
+
 describe('createObb', () => {
   it('stores center, half-extents, and orientation', () => {
     const o = createObb(1, 2, 3, 4, 5, 6, 0, 0, 0, 1);
@@ -72,6 +74,12 @@ describe('getClosestPointOnObb', () => {
   });
 });
 
+describe('initializeObb', () => {
+  it('is the construction initializer of createObb', () => {
+    expect(typeof initializeObb).toBe('function');
+  });
+});
+
 describe('intersectRay3DObb', () => {
   it('hits a unit OBB at the expected t from outside', () => {
     const o = createObb(0, 0, 0, 1, 1, 1, 0, 0, 0, 1);
@@ -293,91 +301,6 @@ describe('setObb', () => {
   });
 });
 
-describe('transformObbByMatrix4', () => {
-  it('translates the center by a translation matrix', () => {
-    const o = createObb(0, 0, 0, 1, 1, 1, 0, 0, 0, 1);
-    const m = createMatrix4();
-    setMatrix4Position(m, createVector3(5, 0, 0));
-    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
-    transformObbByMatrix4(out, o, m);
-    expect(out.centerX).toBeCloseTo(5, 6);
-    expect(out.centerY).toBeCloseTo(0, 6);
-    expect(out.halfExtentX).toBeCloseTo(1, 6);
-    expect(out.orientationW).toBeCloseTo(1, 5);
-  });
-
-  it('scales half-extents by a uniform scale matrix', () => {
-    const o = createObb(0, 0, 0, 1, 1, 1, 0, 0, 0, 1);
-    const m = createMatrix4(2, 0, 0, 0, 0, 2, 0, 0, 0, 0, 2, 0, 0, 0, 0, 1);
-    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
-    transformObbByMatrix4(out, o, m);
-    expect(out.halfExtentX).toBeCloseTo(2, 6);
-    expect(out.halfExtentY).toBeCloseTo(2, 6);
-    expect(out.halfExtentZ).toBeCloseTo(2, 6);
-  });
-
-  it('supports out === obb', () => {
-    const o = createObb(1, 0, 0, 1, 1, 1, 0, 0, 0, 1);
-    const m = createMatrix4();
-    setMatrix4Position(m, createVector3(0, 3, 0));
-    transformObbByMatrix4(o, o, m);
-    expect(o.centerX).toBeCloseTo(1, 6);
-    expect(o.centerY).toBeCloseTo(3, 6);
-    expect(o.halfExtentX).toBeCloseTo(1, 6);
-  });
-
-  it('rotates the box the same way the matrix rotates points', () => {
-    // The orientation extracted from the matrix must turn the box in the matrix's own direction.
-    // A conjugated extraction turns it by the same angle about the opposite axis, which is
-    // invisible for axis-aligned boxes and half turns but wrong for every general rotation.
-    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(1, 0, 0, 0.7), createTiltedObb());
-    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(1, 2, 3, 0.9), createTiltedObb());
-  });
-
-  it('rotates the box correctly near a half turn about each axis', () => {
-    // A half turn drives the three `trace <= 0` branches of the matrix-to-quaternion extraction,
-    // one per dominant diagonal term. Just short of pi so the conjugate is a different rotation:
-    // at exactly pi the quaternion equals its own conjugate up to sign and hides an error.
-    const nearHalfTurn = Math.PI - 0.01;
-    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(1, 0, 0, nearHalfTurn), createTiltedObb());
-    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(0, 1, 0, nearHalfTurn), createTiltedObb());
-    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(0, 0, 1, nearHalfTurn), createTiltedObb());
-  });
-
-  it('flattens the box along a collapsed matrix column and keeps the other extents', () => {
-    // A zero-length column has no direction to read, so that axis falls back to its identity
-    // direction rather than dividing by zero, and the half-extent along it goes to zero.
-    const o = createObb(0, 0, 0, 2, 3, 4, 0, 0, 0, 1);
-    const m = createMatrix4(0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1);
-    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
-    transformObbByMatrix4(out, o, m);
-    expect(out.halfExtentX).toBe(0);
-    expect(out.halfExtentY).toBeCloseTo(3, 6);
-    expect(out.halfExtentZ).toBeCloseTo(4, 6);
-    expect(out.orientationW).toBeCloseTo(1, 6);
-  });
-
-  it('collapses the box to a point at the matrix position when every column is zero', () => {
-    // All three columns collapse, so the fallback rotation is the identity and the box keeps the
-    // orientation it came in with — the transform contributes nothing but the translation.
-    const tilt = orientationFromAxisAngle(1, 1, 1, 0.6);
-    const o = createObb(1, 2, 3, 2, 3, 4, ...tilt);
-    const m = createMatrix4(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 6, 7, 1);
-    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
-    transformObbByMatrix4(out, o, m);
-    expect(out.centerX).toBeCloseTo(5, 6);
-    expect(out.centerY).toBeCloseTo(6, 6);
-    expect(out.centerZ).toBeCloseTo(7, 6);
-    expect(out.halfExtentX).toBe(0);
-    expect(out.halfExtentY).toBe(0);
-    expect(out.halfExtentZ).toBe(0);
-    expect(out.orientationX).toBeCloseTo(tilt[0], 6);
-    expect(out.orientationY).toBeCloseTo(tilt[1], 6);
-    expect(out.orientationZ).toBeCloseTo(tilt[2], 6);
-    expect(out.orientationW).toBeCloseTo(tilt[3], 6);
-  });
-});
-
 // A thin rod: long along one of its own local axes, barely thick along the other two, so the
 // dominant feature of the box is a single edge direction.
 function createRodObb(
@@ -486,3 +409,87 @@ const FACE_AXIS_SEPARATIONS: ReadonlyArray<{
   { axis: 1, clear: 4.7, extents: [0.5, 1.5, 2.5], ofTurnedBox: true },
   { axis: 2, clear: 3.7, extents: [2.5, 1.5, 0.5], ofTurnedBox: true },
 ];
+describe('transformObbByMatrix4', () => {
+  it('translates the center by a translation matrix', () => {
+    const o = createObb(0, 0, 0, 1, 1, 1, 0, 0, 0, 1);
+    const m = createMatrix4();
+    setMatrix4Position(m, createVector3(5, 0, 0));
+    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
+    transformObbByMatrix4(out, o, m);
+    expect(out.centerX).toBeCloseTo(5, 6);
+    expect(out.centerY).toBeCloseTo(0, 6);
+    expect(out.halfExtentX).toBeCloseTo(1, 6);
+    expect(out.orientationW).toBeCloseTo(1, 5);
+  });
+
+  it('scales half-extents by a uniform scale matrix', () => {
+    const o = createObb(0, 0, 0, 1, 1, 1, 0, 0, 0, 1);
+    const m = createMatrix4(2, 0, 0, 0, 0, 2, 0, 0, 0, 0, 2, 0, 0, 0, 0, 1);
+    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
+    transformObbByMatrix4(out, o, m);
+    expect(out.halfExtentX).toBeCloseTo(2, 6);
+    expect(out.halfExtentY).toBeCloseTo(2, 6);
+    expect(out.halfExtentZ).toBeCloseTo(2, 6);
+  });
+
+  it('supports out === obb', () => {
+    const o = createObb(1, 0, 0, 1, 1, 1, 0, 0, 0, 1);
+    const m = createMatrix4();
+    setMatrix4Position(m, createVector3(0, 3, 0));
+    transformObbByMatrix4(o, o, m);
+    expect(o.centerX).toBeCloseTo(1, 6);
+    expect(o.centerY).toBeCloseTo(3, 6);
+    expect(o.halfExtentX).toBeCloseTo(1, 6);
+  });
+
+  it('rotates the box the same way the matrix rotates points', () => {
+    // The orientation extracted from the matrix must turn the box in the matrix's own direction.
+    // A conjugated extraction turns it by the same angle about the opposite axis, which is
+    // invisible for axis-aligned boxes and half turns but wrong for every general rotation.
+    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(1, 0, 0, 0.7), createTiltedObb());
+    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(1, 2, 3, 0.9), createTiltedObb());
+  });
+
+  it('rotates the box correctly near a half turn about each axis', () => {
+    // A half turn drives the three `trace <= 0` branches of the matrix-to-quaternion extraction,
+    // one per dominant diagonal term. Just short of pi so the conjugate is a different rotation:
+    // at exactly pi the quaternion equals its own conjugate up to sign and hides an error.
+    const nearHalfTurn = Math.PI - 0.01;
+    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(1, 0, 0, nearHalfTurn), createTiltedObb());
+    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(0, 1, 0, nearHalfTurn), createTiltedObb());
+    expectTransformedObbMatchesRotatedPoints(rotationMatrix4(0, 0, 1, nearHalfTurn), createTiltedObb());
+  });
+
+  it('flattens the box along a collapsed matrix column and keeps the other extents', () => {
+    // A zero-length column has no direction to read, so that axis falls back to its identity
+    // direction rather than dividing by zero, and the half-extent along it goes to zero.
+    const o = createObb(0, 0, 0, 2, 3, 4, 0, 0, 0, 1);
+    const m = createMatrix4(0, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1, 0, 0, 0, 0, 1);
+    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
+    transformObbByMatrix4(out, o, m);
+    expect(out.halfExtentX).toBe(0);
+    expect(out.halfExtentY).toBeCloseTo(3, 6);
+    expect(out.halfExtentZ).toBeCloseTo(4, 6);
+    expect(out.orientationW).toBeCloseTo(1, 6);
+  });
+
+  it('collapses the box to a point at the matrix position when every column is zero', () => {
+    // All three columns collapse, so the fallback rotation is the identity and the box keeps the
+    // orientation it came in with — the transform contributes nothing but the translation.
+    const tilt = orientationFromAxisAngle(1, 1, 1, 0.6);
+    const o = createObb(1, 2, 3, 2, 3, 4, ...tilt);
+    const m = createMatrix4(0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 5, 6, 7, 1);
+    const out = createObb(0, 0, 0, 0, 0, 0, 0, 0, 0, 1);
+    transformObbByMatrix4(out, o, m);
+    expect(out.centerX).toBeCloseTo(5, 6);
+    expect(out.centerY).toBeCloseTo(6, 6);
+    expect(out.centerZ).toBeCloseTo(7, 6);
+    expect(out.halfExtentX).toBe(0);
+    expect(out.halfExtentY).toBe(0);
+    expect(out.halfExtentZ).toBe(0);
+    expect(out.orientationX).toBeCloseTo(tilt[0], 6);
+    expect(out.orientationY).toBeCloseTo(tilt[1], 6);
+    expect(out.orientationZ).toBeCloseTo(tilt[2], 6);
+    expect(out.orientationW).toBeCloseTo(tilt[3], 6);
+  });
+});

package/src/plane.test.ts

@@ -13,6 +13,8 @@ import {
   setPlaneFromPoints,
 } from '@flighthq/geometry/contract';
 
+import { initializePlane } from './plane';
+
 describe('clonePlane', () => {
   it('creates an independent copy', () => {
     const p = createPlane(1, 0, 0, -5);
@@ -124,6 +126,12 @@ describe('getPlaneSignedDistanceToPoint', () => {
   });
 });
 
+describe('initializePlane', () => {
+  it('is the construction initializer of createPlane', () => {
+    expect(typeof initializePlane).toBe('function');
+  });
+});
+
 describe('normalizePlane', () => {
   it('normalizes a plane with non-unit normal', () => {
     // Normal (0,2,0), d=4 → normalized (0,1,0), d=2.
@@ -201,7 +209,6 @@ describe('setPlaneFromNormalAndPoint', () => {
     expect(p.d).toBeCloseTo(-3, 6);
   });
 });
-
 describe('setPlaneFromPoints', () => {
   it('builds the XZ plane from three points in XZ', () => {
     const p = createPlane();

package/src/quaternion.test.ts

@@ -28,6 +28,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Quaternion } from '@flighthq/types/contract';
 
+import { initializeQuaternion } from './quaternion';
+
 function expectQuaternionClose(q: Readonly<Quaternion>, x: number, y: number, z: number, w: number): void {
   expect(q.x).toBeCloseTo(x, 6);
   expect(q.y).toBeCloseTo(y, 6);
@@ -277,6 +279,12 @@ describe('getQuaternionEuler', () => {
   });
 });
 
+describe('initializeQuaternion', () => {
+  it('is the construction initializer of createQuaternion', () => {
+    expect(typeof initializeQuaternion).toBe('function');
+  });
+});
+
 describe('inverseQuaternion', () => {
   it('q * inverse(q) ≈ identity for a unit quaternion', () => {
     const q = createQuaternion();
@@ -679,7 +687,6 @@ describe('setQuaternionLookRotation', () => {
     expectQuaternionClose(aliasedUp, distinct.x, distinct.y, distinct.z, distinct.w);
   });
 });
-
 describe('slerpQuaternion', () => {
   it('returns a at t=0 and b at t=1', () => {
     const a = createQuaternion();

package/src/ray3d.test.ts

@@ -9,6 +9,7 @@ import {
   getClosestPointBetweenRay3Ds,
   getClosestPointOnRay3D,
   getRay3DPointAt,
+  initializeRay3D,
   intersectRay3DAabb,
   intersectRay3DPlane,
   intersectRay3DSphere,
@@ -204,6 +205,12 @@ describe('getRay3DPointAt', () => {
   });
 });
 
+describe('initializeRay3D', () => {
+  it('is the construction initializer of createRay3D', () => {
+    expect(typeof initializeRay3D).toBe('function');
+  });
+});
+
 describe('intersectRay3DAabb', () => {
   it('returns t=0 when the ray origin is inside the box', () => {
     const ray = createRay3D(0, 0, 0, 0, 0, 1);
@@ -489,6 +496,18 @@ describe('setRay3D', () => {
   });
 });
 
+// Six rays, one per side of the box spanning (-1, -4, -6) to (2, 3, 5), each starting 10 units out
+// along an axis and aimed at it. The last number is where that ray meets the box: all six differ,
+// so a slab reading another axis' bound gives a different answer rather than the same one.
+// Reversing the direction of each turns the set into six rays aimed away from the box.
+const RAYS_INTO_LOPSIDED_BOX: ReadonlyArray<readonly [number, number, number, number, number, number, number]> = [
+  [-10, 0, 0, 1, 0, 0, 9],
+  [10, 0, 0, -1, 0, 0, 8],
+  [0, -10, 0, 0, 1, 0, 6],
+  [0, 10, 0, 0, -1, 0, 7],
+  [0, 0, -10, 0, 0, 1, 4],
+  [0, 0, 10, 0, 0, -1, 5],
+];
 describe('transformRay3DByMatrix4', () => {
   it('leaves the ray unchanged under the identity matrix', () => {
     const ray = createRay3D(1, 2, 3, 0, 0, 1);
@@ -530,16 +549,3 @@ describe('transformRay3DByMatrix4', () => {
     expect(ray.direction.z).toBeCloseTo(1);
   });
 });
-
-// Six rays, one per side of the box spanning (-1, -4, -6) to (2, 3, 5), each starting 10 units out
-// along an axis and aimed at it. The last number is where that ray meets the box: all six differ,
-// so a slab reading another axis' bound gives a different answer rather than the same one.
-// Reversing the direction of each turns the set into six rays aimed away from the box.
-const RAYS_INTO_LOPSIDED_BOX: ReadonlyArray<readonly [number, number, number, number, number, number, number]> = [
-  [-10, 0, 0, 1, 0, 0, 9],
-  [10, 0, 0, -1, 0, 0, 8],
-  [0, -10, 0, 0, 1, 0, 6],
-  [0, 10, 0, 0, -1, 0, 7],
-  [0, 0, -10, 0, 0, 1, 4],
-  [0, 0, 10, 0, 0, -1, 5],
-];

package/src/rectangle.test.ts

@@ -43,6 +43,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Rectangle } from '@flighthq/types/contract';
 
+import { initializeRectangle } from './rectangle';
+
 let r: Rectangle;
 let r2: Rectangle;
 
@@ -601,6 +603,12 @@ describe('inflateRectangle', () => {
   });
 });
 
+describe('initializeRectangle', () => {
+  it('is the construction initializer of createRectangle', () => {
+    expect(typeof initializeRectangle).toBe('function');
+  });
+});
+
 describe('intersectsRectangle', () => {
   it('returns true if rectangles overlap', () => {
     const r3 = createRectangle(5, 10, 10, 10);
@@ -1004,7 +1012,6 @@ describe('setRectangleTopLeft', () => {
     expect(r.y).toBe(10);
   });
 });
-
 describe('setTo', () => {
   it('assigns all four properties', () => {
     const rect = createRectangle();

package/src/transform2d.test.ts

@@ -1,7 +1,7 @@
 import { describe, expect, it } from 'vitest';
 
 import { createMatrix } from './matrix';
-import { createTransform2D, decomposeMatrixToTransform2D } from './transform2d';
+import { createTransform2D, decomposeMatrixToTransform2D, initializeTransform2D } from './transform2d';
 
 const DEG_TO_RAD = Math.PI / 180;
 
@@ -120,3 +120,8 @@ function expectRoundTrip(source: Readonly<ReturnType<typeof createMatrix>>) {
   expect(rebuilt.c).toBeCloseTo(source.c);
   expect(rebuilt.d).toBeCloseTo(source.d);
 }
+describe('initializeTransform2D', () => {
+  it('is the construction initializer of createTransform2D', () => {
+    expect(typeof initializeTransform2D).toBe('function');
+  });
+});

package/src/transform3d.test.ts

@@ -2,7 +2,12 @@ import { describe, expect, it } from 'vitest';
 
 import { createMatrix4 } from './matrix4';
 import { createQuaternion } from './quaternion';
-import { composeMatrix4FromTransform3D, createTransform3D, decomposeMatrix4ToTransform3D } from './transform3d';
+import {
+  composeMatrix4FromTransform3D,
+  createTransform3D,
+  decomposeMatrix4ToTransform3D,
+  initializeTransform3D,
+} from './transform3d';
 import { createVector3 } from './vector3';
 
 describe('composeMatrix4FromTransform3D', () => {
@@ -48,3 +53,8 @@ describe('decomposeMatrix4ToTransform3D', () => {
     expect(out.scale.z).toBeCloseTo(4);
   });
 });
+describe('initializeTransform3D', () => {
+  it('is the construction initializer of createTransform3D', () => {
+    expect(typeof initializeTransform3D).toBe('function');
+  });
+});

package/src/vector2.test.ts

@@ -36,6 +36,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Vector2 } from '@flighthq/types/contract';
 
+import { initializeVector2 } from './vector2';
+
 let pt: Vector2;
 let pt2: Vector2;
 
@@ -511,6 +513,12 @@ describe('getVector2LengthSquared', () => {
   });
 });
 
+describe('initializeVector2', () => {
+  it('is the construction initializer of createVector2', () => {
+    expect(typeof initializeVector2).toBe('function');
+  });
+});
+
 describe('interpolateVector2', () => {
   const cases = [
     { t: 0, expected: (a: number, _: number) => a },
@@ -1181,7 +1189,6 @@ describe('subtractVector2', () => {
     expect(result).not.toBe(pt2);
   });
 });
-
 describe('writeVector2ToFloat32Array', () => {
   it('writes 2 values at the offset', () => {
     const array = new Float32Array(6);

package/src/vector3.test.ts

@@ -39,6 +39,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Vector3 } from '@flighthq/types/contract';
 
+import { initializeVector3 } from './vector3';
+
 describe('addVector3', () => {
   it('returns a new vector when no target is passed', () => {
     const a = createVector3(1, 2, 3);
@@ -461,6 +463,12 @@ describe('getVector3Spherical', () => {
   });
 });
 
+describe('initializeVector3', () => {
+  it('is the construction initializer of createVector3', () => {
+    expect(typeof initializeVector3).toBe('function');
+  });
+});
+
 describe('interpolateVector3', () => {
   it('returns a at t=0 and b at t=1', () => {
     const out = createVector3();
@@ -947,7 +955,6 @@ describe('Y_AXIS', () => {
     expect(yAxis.z).toBe(0);
   });
 });
-
 describe('Z_AXIS', () => {
   it('returns the unit vector along the Z-axis', () => {
     const zAxis: Vector3 = VECTOR3_Z_AXIS;

package/src/vector4.test.ts

@@ -36,6 +36,8 @@ import {
 } from '@flighthq/geometry/contract';
 import type { Vector4 } from '@flighthq/types/contract';
 
+import { initializeVector4 } from './vector4';
+
 describe('addVector4', () => {
   it('returns a new vector when no target is passed', () => {
     const a = createVector4(1, 2, 3, 10);
@@ -398,6 +400,12 @@ describe('getVector4LengthSquared', () => {
   });
 });
 
+describe('initializeVector4', () => {
+  it('is the construction initializer of createVector4', () => {
+    expect(typeof initializeVector4).toBe('function');
+  });
+});
+
 describe('interpolateVector4', () => {
   it('returns a at t=0 and b at t=1', () => {
     const out = createVector4();
@@ -583,6 +591,8 @@ describe('nearEqualsVector4', () => {
   });
 });
 
+// Properties
+
 describe('negateVector4', () => {
   it('inverts the values of the vector components', () => {
     const v = createVector4(1, -2, 3, -4);
@@ -604,8 +614,6 @@ describe('negateVector4', () => {
   });
 });
 
-// Properties
-
 describe('normalizeVector4', () => {
   it('normalizes the vector', () => {
     const v = createVector4(3, 4, 0, 0);
@@ -878,7 +886,6 @@ describe('Y_AXIS', () => {
     expect(yAxis.w).toBe(0);
   });
 });
-
 describe('Z_AXIS', () => {
   it('returns the unit vector along the Z-axis', () => {
     const zAxis = VECTOR4_Z_AXIS;