@pawells/math-extended 2.0.0 → 3.0.0

package/build/vectors/core.d.ts

@@ -2,7 +2,7 @@
  * Core vector mathematics operations for linear algebra and geometric calculations.
  * Provides a comprehensive set of vector operations with type safety and error checking.
  */
-import { TAnyVector, TVectorResult, TVector, TVector2, TVector3 } from './types.js';
+import type { TAnyVector, TVectorResult, TVector, TVector2, TVector3 } from './types.js';
 /**
  * Creates a deep copy of a vector.
  * Essential for avoiding mutations when performing operations that should preserve the original vector.
@@ -12,11 +12,11 @@ import { TAnyVector, TVectorResult, TVector, TVector2, TVector3 } from './types.
  * @returns A new vector with identical components
  *
  * @example
-     * ```typescript
-     * const original = [1, 2, 3];
-     * const copy = VectorClone(original);
-     * copy[0] = 10; // original remains unchanged
-     * ```
+ * ```typescript
+ * const original = [1, 2, 3];
+ * const copy = VectorClone(original);
+ * copy[0] = 10; // original remains unchanged
+ * ```
  */
 export declare function VectorClone<T extends TAnyVector>(vector: T): TVectorResult<T>;
 /**
@@ -30,12 +30,12 @@ export declare function VectorClone<T extends TAnyVector>(vector: T): TVectorRes
  * @returns True if vectors are equal within tolerance, false otherwise
  *
  * @example
-     * ```typescript
-     * const a = [1.0001, 2.0001];
-     * const b = [1.0002, 2.0002];
-     * const exactlyEqual = VectorEquals(a, b); // false
-     * const approximatelyEqual = VectorEquals(a, b, 0.001); // true
-     * ```
+ * ```typescript
+ * const a = [1.0001, 2.0001];
+ * const b = [1.0002, 2.0002];
+ * const exactlyEqual = VectorEquals(a, b); // false
+ * const approximatelyEqual = VectorEquals(a, b, 0.001); // true
+ * ```
  */
 export declare function VectorEquals<T extends TAnyVector>(a: T, b: T, tolerance?: number): boolean;
 /**
@@ -47,11 +47,11 @@ export declare function VectorEquals<T extends TAnyVector>(a: T, b: T, tolerance
  * @returns String representation of the vector
  *
  * @example
-     * ```typescript
-     * const vec = [1, 2, 3];
-     * const parens = VectorToString(vec, 'parens'); // "(1, 2, 3)"
-     * const brackets = VectorToString(vec, 'brackets'); // "[1, 2, 3]"
-     * ```
+ * ```typescript
+ * const vec = [1, 2, 3];
+ * const parens = VectorToString(vec, 'parens'); // "(1, 2, 3)"
+ * const brackets = VectorToString(vec, 'brackets'); // "[1, 2, 3]"
+ * ```
  */
 export declare function VectorToString(vector: TVector, style?: 'parens' | 'brackets'): string;
 /**
@@ -64,11 +64,11 @@ export declare function VectorToString(vector: TVector, style?: 'parens' | 'brac
  * @returns New vector where each component is the sum of corresponding components
  *
  * @example
-     * ```typescript
-     * const position = [10, 20, 30];
-     * const velocity = [1, -2, 0.5];
-     * const newPosition = VectorAdd(position, velocity); // [11, 18, 30.5]
-     * ```
+ * ```typescript
+ * const position = [10, 20, 30];
+ * const velocity = [1, -2, 0.5];
+ * const newPosition = VectorAdd(position, velocity); // [11, 18, 30.5]
+ * ```
  */
 export declare function VectorAdd<T extends TAnyVector>(a: T, b: T): TVectorResult<T>;
 /**
@@ -81,11 +81,11 @@ export declare function VectorAdd<T extends TAnyVector>(a: T, b: T): TVectorResu
  * @returns New vector where each component is the difference of corresponding components
  *
  * @example
-     * ```typescript
-     * const target = [100, 50, 0];
-     * const current = [80, 30, 0];
-     * const direction = VectorSubtract(target, current); // [20, 20, 0]
-     * ```
+ * ```typescript
+ * const target = [100, 50, 0];
+ * const current = [80, 30, 0];
+ * const direction = VectorSubtract(target, current); // [20, 20, 0]
+ * ```
  */
 export declare function VectorSubtract<T extends TAnyVector>(a: T, b: T): TVectorResult<T>;
 /**
@@ -98,12 +98,12 @@ export declare function VectorSubtract<T extends TAnyVector>(a: T, b: T): TVecto
  * @returns New vector with multiplied components
  *
  * @example
-     * ```typescript
-     * const velocity = [10, 5, 0];
-     * const scaled = VectorMultiply(velocity, 2); // [20, 10, 0] - scalar multiplication
-     * const factors = [1, -1, 0.5];
-     * const componentWise = VectorMultiply(velocity, factors); // [10, -5, 0] - component-wise
-     * ```
+ * ```typescript
+ * const velocity = [10, 5, 0];
+ * const scaled = VectorMultiply(velocity, 2); // [20, 10, 0] - scalar multiplication
+ * const factors = [1, -1, 0.5];
+ * const componentWise = VectorMultiply(velocity, factors); // [10, -5, 0] - component-wise
+ * ```
  */
 export declare function VectorMultiply<T extends TAnyVector>(a: T, b: T | number): TVectorResult<T>;
 /**
@@ -115,11 +115,11 @@ export declare function VectorMultiply<T extends TAnyVector>(a: T, b: T | number
  * @returns The straight-line distance between the two points represented by the vectors
  *
  * @example
-     * ```typescript
-     * const pointA = [0, 0, 0];
-     * const pointB = [3, 4, 0];
-     * const distance = VectorDistance(pointA, pointB); // 5.0 (3-4-5 triangle)
-     * ```
+ * ```typescript
+ * const pointA = [0, 0, 0];
+ * const pointB = [3, 4, 0];
+ * const distance = VectorDistance(pointA, pointB); // 5.0 (3-4-5 triangle)
+ * ```
  */
 export declare function VectorDistance(a: TVector, b: TVector): number;
 /**
@@ -132,11 +132,11 @@ export declare function VectorDistance(a: TVector, b: TVector): number;
  * @returns The squared distance between vectors
  *
  * @example
-     * ```typescript
-     * const pointA = [1, 1];
-     * const pointB = [4, 5];
-     * const distSq = VectorDistanceSquared(pointA, pointB); // 25 (faster than distance comparison)
-     * ```
+ * ```typescript
+ * const pointA = [1, 1];
+ * const pointB = [4, 5];
+ * const distSq = VectorDistanceSquared(pointA, pointB); // 25 (faster than distance comparison)
+ * ```
  */
 export declare function VectorDistanceSquared(a: TVector, b: TVector): number;
 /**
@@ -149,13 +149,13 @@ export declare function VectorDistanceSquared(a: TVector, b: TVector): number;
  * @returns The dot product (scalar value)
  *
  * @example
-     * ```typescript
-     * const forward = [0, 0, 1];
-     * const direction = [0, 0, 2];
-     * const dot = VectorDot(forward, direction); // 2 (same direction)
-     * const perpendicular = [1, 0, 0];
-     * const dotPerp = VectorDot(forward, perpendicular); // 0 (perpendicular)
-     * ```
+ * ```typescript
+ * const forward = [0, 0, 1];
+ * const direction = [0, 0, 2];
+ * const dot = VectorDot(forward, direction); // 2 (same direction)
+ * const perpendicular = [1, 0, 0];
+ * const dotPerp = VectorDot(forward, perpendicular); // 0 (perpendicular)
+ * ```
  */
 export declare function VectorDot(a: TVector, b: TVector): number;
 /**
@@ -169,12 +169,12 @@ export declare function VectorDot(a: TVector, b: TVector): number;
  * @throws {VectorError} If the vector is zero or has infinite magnitude
  *
  * @example
-     * ```typescript
-     * const vector = [3, 4, 0];
-     * const normalized = VectorNormalize(vector); // [0.6, 0.8, 0] (magnitude = 1)
-     * const direction = [10, 0, 0];
-     * const unitDirection = VectorNormalize(direction); // [1, 0, 0]
-     * ```
+ * ```typescript
+ * const vector = [3, 4, 0];
+ * const normalized = VectorNormalize(vector); // [0.6, 0.8, 0] (magnitude = 1)
+ * const direction = [10, 0, 0];
+ * const unitDirection = VectorNormalize(direction); // [1, 0, 0]
+ * ```
  */
 export declare function VectorNormalize<T extends TAnyVector>(a: T): TVectorResult<T>;
 /**
@@ -185,12 +185,12 @@ export declare function VectorNormalize<T extends TAnyVector>(a: T): TVectorResu
  * @returns The magnitude (length) of the vector
  *
  * @example
-     * ```typescript
-     * const velocity = [3, 4, 0];
-     * const speed = VectorMagnitude(velocity); // 5.0
-     * const unitVector = [1, 0, 0];
-     * const unitLength = VectorMagnitude(unitVector); // 1.0
-     * ```
+ * ```typescript
+ * const velocity = [3, 4, 0];
+ * const speed = VectorMagnitude(velocity); // 5.0
+ * const unitVector = [1, 0, 0];
+ * const unitLength = VectorMagnitude(unitVector); // 1.0
+ * ```
  */
 export declare function VectorMagnitude(a: TVector): number;
 /**
@@ -202,12 +202,12 @@ export declare function VectorMagnitude(a: TVector): number;
  * @returns New vector with absolute values of all components
  *
  * @example
-     * ```typescript
-     * const vector = [-3, 4, -2];
-     * const absolute = VectorAbs(vector); // [3, 4, 2]
-     * const mixed = [1.5, -2.7, 0];
-     * const absValues = VectorAbs(mixed); // [1.5, 2.7, 0]
-     * ```
+ * ```typescript
+ * const vector = [-3, 4, -2];
+ * const absolute = VectorAbs(vector); // [3, 4, 2]
+ * const mixed = [1.5, -2.7, 0];
+ * const absValues = VectorAbs(mixed); // [1.5, 2.7, 0]
+ * ```
  */
 export declare function VectorAbs<T extends TAnyVector>(a: T): TVectorResult<T>;
 /**
@@ -218,12 +218,12 @@ export declare function VectorAbs<T extends TAnyVector>(a: T): TVectorResult<T>;
  * @returns True if all components are zero, false otherwise
  *
  * @example
-     * ```typescript
-     * const zero = [0, 0, 0];
-     * const isZero = VectorIsZero(zero); // true
-     * const notZero = [0, 0.001, 0];
-     * const isNotZero = VectorIsZero(notZero); // false
-     * ```
+ * ```typescript
+ * const zero = [0, 0, 0];
+ * const isZero = VectorIsZero(zero); // true
+ * const notZero = [0, 0.001, 0];
+ * const isNotZero = VectorIsZero(notZero); // false
+ * ```
  */
 export declare function VectorIsZero(vector: TVector): boolean;
 /**
@@ -237,16 +237,33 @@ export declare function VectorIsZero(vector: TVector): boolean;
  * @throws {VectorError} If either vector is zero
  *
  * @example
-     * ```typescript
-     * const right = [1, 0, 0];
-     * const up = [0, 1, 0];
-     * const angle = VectorAngle(right, up); // π/2 (90 degrees)
-     * const forward = [0, 0, 1];
-     * const backward = [0, 0, -1];
-     * const oppositeAngle = VectorAngle(forward, backward); // π (180 degrees)
-     * ```
+ * ```typescript
+ * const right = [1, 0, 0];
+ * const up = [0, 1, 0];
+ * const angle = VectorAngle(right, up); // π/2 (90 degrees)
+ * const forward = [0, 0, 1];
+ * const backward = [0, 0, -1];
+ * const oppositeAngle = VectorAngle(forward, backward); // π (180 degrees)
+ * ```
  */
 export declare function VectorAngle(a: TVector, b: TVector): number;
+/**
+ * Limits (clamps) the magnitude of a vector to a maximum value.
+ * If the vector's magnitude exceeds the maximum, scales it down proportionally while preserving direction.
+ * If the magnitude is below the maximum, the vector is returned unchanged.
+ *
+ * @param vector - The input vector to limit
+ * @param maxMagnitude - The maximum allowed magnitude (must be >= 0)
+ * @returns A new vector with magnitude limited to maxMagnitude
+ * @throws {VectorError} If vector is invalid or maxMagnitude is invalid
+ *
+ * @example
+ * ```typescript
+ *   const v = [3, 4]; // magnitude 5
+ *   VectorLimit(v, 3); // [1.8, 2.4] — magnitude now 3
+ * ```
+ */
+export declare function VectorLimit<T extends TAnyVector>(vector: T, maxMagnitude: number): TVectorResult<T>;
 /**
  * Rotates a 2D vector by the specified angle in radians.
  * Essential for 2D transformations, sprite rotations, and directional calculations.
@@ -256,11 +273,11 @@ export declare function VectorAngle(a: TVector, b: TVector): number;
  * @returns New rotated 2D vector
  *
  * @example
-     * ```typescript
-     * const right = [1, 0];
-     * const rotated90 = Vector2Rotate(right, Math.PI / 2); // [0, 1] (up)
-     * const rotated180 = Vector2Rotate(right, Math.PI); // [-1, 0] (left)
-     * ```
+ * ```typescript
+ * const right = [1, 0];
+ * const rotated90 = Vector2Rotate(right, Math.PI / 2); // [0, 1] (up)
+ * const rotated180 = Vector2Rotate(right, Math.PI); // [-1, 0] (left)
+ * ```
  */
 export declare function Vector2Rotate(vector: TVector2, radians: number): TVector2;
 /**
@@ -271,11 +288,11 @@ export declare function Vector2Rotate(vector: TVector2, radians: number): TVecto
  * @returns Unit vector pointing in the specified direction
  *
  * @example
-     * ```typescript
-     * const right = Vector2FromAngle(0); // [1, 0]
-     * const up = Vector2FromAngle(Math.PI / 2); // [0, 1]
-     * const diagonal = Vector2FromAngle(Math.PI / 4); // [0.707, 0.707]
-     * ```
+ * ```typescript
+ * const right = Vector2FromAngle(0); // [1, 0]
+ * const up = Vector2FromAngle(Math.PI / 2); // [0, 1]
+ * const diagonal = Vector2FromAngle(Math.PI / 4); // [0.707, 0.707]
+ * ```
  */
 export declare function Vector2FromAngle(radians: number): TVector2;
 /**
@@ -288,12 +305,12 @@ export declare function Vector2FromAngle(radians: number): TVector2;
  * @returns Scalar cross product (positive = counterclockwise, negative = clockwise)
  *
  * @example
-     * ```typescript
-     * const right = [1, 0];
-     * const up = [0, 1];
-     * const cross = Vector2Cross(right, up); // 1 (counterclockwise)
-     * const crossReverse = Vector2Cross(up, right); // -1 (clockwise)
-     * ```
+ * ```typescript
+ * const right = [1, 0];
+ * const up = [0, 1];
+ * const cross = Vector2Cross(right, up); // 1 (counterclockwise)
+ * const crossReverse = Vector2Cross(up, right); // -1 (clockwise)
+ * ```
  */
 export declare function Vector2Cross(a: TVector2, b: TVector2): number;
 /**
@@ -307,11 +324,11 @@ export declare function Vector2Cross(a: TVector2, b: TVector2): number;
  * @throws {VectorError} If vector b is zero
  *
  * @example
-     * ```typescript
-     * const force = [5, 3, 0];
-     * const surface = [1, 0, 0];
-     * const perpendicular = Vector3Reject(force, surface); // [0, 3, 0]
-     * ```
+ * ```typescript
+ * const force = [5, 3, 0];
+ * const surface = [1, 0, 0];
+ * const perpendicular = Vector3Reject(force, surface); // [0, 3, 0]
+ * ```
  */
 export declare function Vector3Reject(a: TVector3, b: TVector3): TVector3;
 /**
@@ -326,11 +343,11 @@ export declare function Vector3Reject(a: TVector3, b: TVector3): TVector3;
  * @throws {VectorError} If vector b is zero
  *
  * @example
-     * ```typescript
-     * const force = [5, 3, 0];
-     * const surface = [1, 0, 0];
-     * const parallel = VectorProject(force, surface); // [5, 0, 0]
-     * ```
+ * ```typescript
+ * const force = [5, 3, 0];
+ * const surface = [1, 0, 0];
+ * const parallel = VectorProject(force, surface); // [5, 0, 0]
+ * ```
  */
 export declare function VectorProject<T extends TAnyVector>(a: T, b: T): TVectorResult<T>;
 /**
@@ -344,11 +361,11 @@ export declare function VectorProject<T extends TAnyVector>(a: T, b: T): TVector
  * @throws {VectorError} If the normal is a zero vector
  *
  * @example
-     * ```typescript
-     * const incoming = [1, -1, 0];
-     * const normal = [0, 1, 0]; // surface normal (upward)
-     * const reflected = Vector3Reflect(incoming, normal); // [1, 1, 0]
-     * ```
+ * ```typescript
+ * const incoming = [1, -1, 0];
+ * const normal = [0, 1, 0]; // surface normal (upward)
+ * const reflected = Vector3Reflect(incoming, normal); // [1, 1, 0]
+ * ```
  */
 export declare function Vector3Reflect(incident: TVector3, normal: TVector3): TVector3;
 /**
@@ -361,12 +378,12 @@ export declare function Vector3Reflect(incident: TVector3, normal: TVector3): TV
  * @returns Vector perpendicular to both a and b (following right-hand rule)
  *
  * @example
-     * ```typescript
-     * const right = [1, 0, 0];
-     * const forward = [0, 0, 1];
-     * const up = Vector3Cross(right, forward); // [0, 1, 0]
-     * const normal = Vector3Cross([1, 0, 0], [0, 1, 0]); // [0, 0, 1]
-     * ```
+ * ```typescript
+ * const right = [1, 0, 0];
+ * const forward = [0, 0, 1];
+ * const up = Vector3Cross(right, forward); // [0, 1, 0]
+ * const normal = Vector3Cross([1, 0, 0], [0, 1, 0]); // [0, 0, 1]
+ * ```
  */
 export declare function Vector3Cross(a: TVector3, b: TVector3): TVector3;
 /**
@@ -379,11 +396,11 @@ export declare function Vector3Cross(a: TVector3, b: TVector3): TVector3;
  * @returns Magnitude of the cross product
  *
  * @example
-     * ```typescript
-     * const side1 = [3, 0, 0];
-     * const side2 = [0, 4, 0];
-     * const area = VectorCrossMagnitude(side1, side2); // 12 (area of rectangle)
-     * ```
+ * ```typescript
+ * const side1 = [3, 0, 0];
+ * const side2 = [0, 4, 0];
+ * const area = VectorCrossMagnitude(side1, side2); // 12 (area of rectangle)
+ * ```
  */
 export declare function VectorCrossMagnitude(a: TVector3, b: TVector3): number;
 /**
@@ -397,12 +414,12 @@ export declare function VectorCrossMagnitude(a: TVector3, b: TVector3): number;
  * @returns Signed volume (positive = right-handed orientation)
  *
  * @example
-     * ```typescript
-     * const x = [1, 0, 0];
-     * const y = [0, 1, 0];
-     * const z = [0, 0, 1];
-     * const volume = Vector3ScalarTripleProduct(x, y, z); // 1 (unit cube)
-     * ```
+ * ```typescript
+ * const x = [1, 0, 0];
+ * const y = [0, 1, 0];
+ * const z = [0, 0, 1];
+ * const volume = Vector3ScalarTripleProduct(x, y, z); // 1 (unit cube)
+ * ```
  */
 export declare function Vector3ScalarTripleProduct(a: TVector3, b: TVector3, c: TVector3): number;
 /**
@@ -416,12 +433,12 @@ export declare function Vector3ScalarTripleProduct(a: TVector3, b: TVector3, c:
  * @returns Vector result of a × (b × c)
  *
  * @example
-     * ```typescript
-     * const a = [1, 0, 0];
-     * const b = [0, 1, 0];
-     * const c = [0, 0, 1];
-     * const result = Vector3TripleProduct(a, b, c); // [0, 0, 0]
-     * ```
+ * ```typescript
+ * const a = [1, 0, 0];
+ * const b = [0, 1, 0];
+ * const c = [0, 0, 1];
+ * const result = Vector3TripleProduct(a, b, c); // [0, 0, 0]
+ * ```
  */
 export declare function Vector3TripleProduct(a: TVector3, b: TVector3, c: TVector3): TVector3;
 /**
@@ -435,11 +452,11 @@ export declare function Vector3TripleProduct(a: TVector3, b: TVector3, c: TVecto
  * @returns The reflected vector
  *
  * @example
-     * ```typescript
-     * const incoming = [1, -1, 0];
-     * const wall = [0, 1, 0]; // vertical wall normal
-     * const bounced = VectorReflect(incoming, wall); // [1, 1, 0]
-     * ```
+ * ```typescript
+ * const incoming = [1, -1, 0];
+ * const wall = [0, 1, 0]; // vertical wall normal
+ * const bounced = VectorReflect(incoming, wall); // [1, 1, 0]
+ * ```
  */
 export declare function VectorReflect<T extends TAnyVector>(incident: T, normal: T): TVectorResult<T>;
 /**
@@ -452,12 +469,12 @@ export declare function VectorReflect<T extends TAnyVector>(incident: T, normal:
  * @returns Vector with all components negated
  *
  * @example
-     * ```typescript
-     * const forward = [0, 0, 1];
-     * const backward = VectorNegate(forward); // [0, 0, -1]
-     * const velocity = [5, -3, 2];
-     * const opposite = VectorNegate(velocity); // [-5, 3, -2]
-     * ```
+ * ```typescript
+ * const forward = [0, 0, 1];
+ * const backward = VectorNegate(forward); // [0, 0, -1]
+ * const velocity = [5, -3, 2];
+ * const opposite = VectorNegate(velocity); // [-5, 3, -2]
+ * ```
  */
 export declare function VectorNegate<T extends TAnyVector>(a: T): TVectorResult<T>;
 /**
@@ -471,12 +488,12 @@ export declare function VectorNegate<T extends TAnyVector>(a: T): TVectorResult<
  * @throws {VectorError} If any divisor component is zero
  *
  * @example
-     * ```typescript
-     * const velocity = [20, 10, 0];
-     * const halved = VectorDivide(velocity, 2); // [10, 5, 0] - scalar division
-     * const factors = [2, 5, 1];
-     * const componentWise = VectorDivide(velocity, factors); // [10, 2, 0] - component-wise
-     * ```
+ * ```typescript
+ * const velocity = [20, 10, 0];
+ * const halved = VectorDivide(velocity, 2); // [10, 5, 0] - scalar division
+ * const factors = [2, 5, 1];
+ * const componentWise = VectorDivide(velocity, factors); // [10, 2, 0] - component-wise
+ * ```
  */
 export declare function VectorDivide<T extends TAnyVector>(a: T, b: T | number): TVectorResult<T>;
 /**
@@ -491,79 +508,15 @@ export declare function VectorDivide<T extends TAnyVector>(a: T, b: T | number):
  * @returns New vector with each component clamped between min and max
  *
  * @example
-     * ```typescript
-     * const v = [5, -3, 12, 0];
-     * VectorClamp(v, 0, 10); // [5, 0, 10, 0] - scalar bounds
-     * const mins = [0, -5, 0, -1];
-     * const maxs = [10, 5, 8, 1];
-     * VectorClamp(v, mins, maxs); // [5, -3, 8, 0] - per-component bounds
-     * ```
+ * ```typescript
+ * const v = [5, -3, 12, 0];
+ * VectorClamp(v, 0, 10); // [5, 0, 10, 0] - scalar bounds
+ * const mins = [0, -5, 0, -1];
+ * const maxs = [10, 5, 8, 1];
+ * VectorClamp(v, mins, maxs); // [5, -3, 8, 0] - per-component bounds
+ * ```
  */
 export declare function VectorClamp<T extends TAnyVector>(a: T, min: T | number, max: T | number): TVectorResult<T>;
-/**
- * Limits the magnitude of a vector to a maximum value.
- * If the vector's magnitude exceeds the limit, scales it down proportionally.
- * Preserves direction while constraining magnitude.
- *
- * @template T - The vector type extending TVector
- * @param a - Vector to limit
- * @param max - Maximum allowed magnitude
- * @returns Vector with magnitude limited to max
- * @throws {VectorError} If max is negative
- *
- * @example
-     * ```typescript
-     * const velocity = [15, 20, 0]; // magnitude ≈ 25
-     * const limited = VectorLimit(velocity, 10); // magnitude = 10, same direction
-     * const small = [1, 1, 0]; // magnitude ≈ 1.414
-     * const unchanged = VectorLimit(small, 5); // unchanged since already under limit
-     * ```
-export function VectorGramSchmidt<T extends TAnyVector>(vectors: T[], normalize: boolean = false): TVectorResult<T>[] {
-    if (vectors.length === 0) throw new VectorError('GramSchmidt: Empty Vector Set');
-    const [firstVector] = vectors;
-    if (!firstVector) throw new VectorError('GramSchmidt: Undefined First Vector');
-    const dimension = firstVector.length;
-    for (const [i, vector] of vectors.entries()) {
-        AssertVector(vector);
-        if (vector.length !== dimension) throw new VectorError(`GramSchmidt: Vector at index ${i} has different dimension than first vector. Expected ${dimension}, got ${vector.length}`);
-        if (VectorIsZero(vector)) throw new VectorError(`GramSchmidt: Vector at index ${i} is a zero vector. Cannot orthogonalize zero vectors.`);
-    }
-
-    const result: TVectorResult<T>[] = [];
-
-    for (const [i, currentVector] of vectors.entries()) {
-        AssertVector(currentVector);
-
-        // Start with a clone of the current vector
-        let orthogonalVector = VectorClone(currentVector);
-
-        // Subtract projections of all previous vectors
-        for (let j = 0; j < i; j++) {
-            const previousVector = result[j];
-            AssertVector(previousVector);
-
-            // Both vectors have the same dimension, projection result is compatible
-            const vectorT = currentVector as unknown as T;
-            const prevT = previousVector as unknown as T;
-            const projection = VectorProject(vectorT, prevT);
-            orthogonalVector = VectorSubtract(orthogonalVector as unknown as T, projection as unknown as T);
-        }
-        
-        if (VectorIsZero(orthogonalVector as unknown as T)) {
-            throw new VectorError(`GramSchmidt: Vector at index ${i} is linearly dependent on previous vectors. Cannot orthogonalize linearly dependent vectors.`);
-        }
-        
-        if (normalize) {
-            orthogonalVector = VectorNormalize(orthogonalVector as unknown as T);
-        }
-        
-        result.push(orthogonalVector as TVectorResult<T>);
-    }
-
-    return result;
-}}
-}
-
 /**
  * Performs Gram-Schmidt orthogonalization on a set of vectors.
  * Converts a set of linearly independent vectors into an orthogonal (or orthonormal) set.
@@ -576,11 +529,11 @@ export function VectorGramSchmidt<T extends TAnyVector>(vectors: T[], normalize:
  * @throws {VectorError} If vectors are linearly dependent or invalid
  *
  * @example
-     * ```typescript
-     * const vectors = [[1, 1, 0], [1, 0, 1], [0, 1, 1]];
-     * const orthogonal = VectorGramSchmidt(vectors); // Orthogonal set
-     * const orthonormal = VectorGramSchmidt(vectors, true); // Orthonormal set
-     * ```
+ * ```typescript
+ * const vectors = [[1, 1, 0], [1, 0, 1], [0, 1, 1]];
+ * const orthogonal = VectorGramSchmidt(vectors); // Orthogonal set
+ * const orthonormal = VectorGramSchmidt(vectors, true); // Orthonormal set
+ * ```
  */
 export declare function VectorGramSchmidt<T extends TAnyVector>(vectors: T[], normalize?: boolean): TVectorResult<T>[];
 //# sourceMappingURL=core.d.ts.map

package/build/vectors/core.d.ts.map

@@ -1 +1 @@
-{"version":3,"file":"core.d.ts","sourceRoot":"","sources":["../../src/vectors/core.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAIH,OAAO,EAAE,UAAU,EAAE,aAAa,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAEpF;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,WAAW,CAAC,CAAC,SAAS,UAAU,EAAE,MAAM,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAG7E;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,YAAY,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,GAAE,MAAU,GAAG,OAAO,CAgB7F;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,cAAc,CAAC,MAAM,EAAE,OAAO,EAAE,KAAK,GAAE,QAAQ,GAAG,UAAqB,GAAG,MAAM,CAQ/F;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,SAAS,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAe5E;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,cAAc,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAejF;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,cAAc,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CA0B1F;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAE7D;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,qBAAqB,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAgBpE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAexD;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,eAAe,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAgB5E;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,OAAO,GAAG,MAAM,CAWlD;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,SAAS,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAWtE;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,OAAO,GAAG,OAAO,CAErD;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAO1D;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,aAAa,CAAC,MAAM,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,GAAG,QAAQ,CAUzE;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,QAAQ,CAE1D;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,MAAM,CAI7D;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAMhE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,aAAa,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAehF;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,cAAc,CAAC,QAAQ,EAAE,QAAQ,EAAE,MAAM,EAAE,QAAQ,GAAG,QAAQ,CAM7E;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAQ/D;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,MAAM,CAMrE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,0BAA0B,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,MAAM,CAOxF;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAOpF;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,aAAa,CAAC,CAAC,SAAS,UAAU,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAkB5F;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,YAAY,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAczE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,YAAY,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CA4BxF;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,WAAW,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CAe1G;AAED;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;GAiFG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,SAAS,UAAU,EAAE,OAAO,EAAE,CAAC,EAAE,EAAE,SAAS,GAAE,OAAe,GAAG,aAAa,CAAC,CAAC,CAAC,EAAE,CA+BpH"}
+{"version":3,"file":"core.d.ts","sourceRoot":"","sources":["../../src/vectors/core.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAIH,OAAO,KAAK,EAAE,UAAU,EAAE,aAAa,EAAE,OAAO,EAAE,QAAQ,EAAE,QAAQ,EAAE,MAAM,YAAY,CAAC;AAEzF;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,WAAW,CAAC,CAAC,SAAS,UAAU,EAAE,MAAM,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAG7E;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,YAAY,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,SAAS,GAAE,MAAU,GAAG,OAAO,CAe7F;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,cAAc,CAAC,MAAM,EAAE,OAAO,EAAE,KAAK,GAAE,QAAQ,GAAG,UAAqB,GAAG,MAAM,CAM/F;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,SAAS,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAc5E;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,cAAc,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAcjF;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,cAAc,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CAuB1F;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,cAAc,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAE7D;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,qBAAqB,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAepE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,SAAS,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAcxD;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,eAAe,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAe5E;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,eAAe,CAAC,CAAC,EAAE,OAAO,GAAG,MAAM,CAQlD;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,SAAS,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAQtE;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,YAAY,CAAC,MAAM,EAAE,OAAO,GAAG,OAAO,CAErD;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,WAAW,CAAC,CAAC,EAAE,OAAO,EAAE,CAAC,EAAE,OAAO,GAAG,MAAM,CAS1D;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,WAAW,CAAC,CAAC,SAAS,UAAU,EAAE,MAAM,EAAE,CAAC,EAAE,YAAY,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CAUnG;AAED;;;;;;;;;;;;;;GAcG;AACH,wBAAgB,aAAa,CAAC,MAAM,EAAE,QAAQ,EAAE,OAAO,EAAE,MAAM,GAAG,QAAQ,CAUzE;AAED;;;;;;;;;;;;;GAaG;AACH,wBAAgB,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,QAAQ,CAE1D;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,MAAM,CAI7D;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,aAAa,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAMhE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,aAAa,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAchF;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,cAAc,CAAC,QAAQ,EAAE,QAAQ,EAAE,MAAM,EAAE,QAAQ,GAAG,QAAQ,CAM7E;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,YAAY,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAQ/D;AAED;;;;;;;;;;;;;;;GAeG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,MAAM,CAMrE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,0BAA0B,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,MAAM,CAOxF;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,oBAAoB,CAAC,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,EAAE,CAAC,EAAE,QAAQ,GAAG,QAAQ,CAOpF;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,aAAa,CAAC,CAAC,SAAS,UAAU,EAAE,QAAQ,EAAE,CAAC,EAAE,MAAM,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAiB5F;AAED;;;;;;;;;;;;;;;;GAgBG;AACH,wBAAgB,YAAY,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,GAAG,aAAa,CAAC,CAAC,CAAC,CAQzE;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,YAAY,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,EAAE,CAAC,GAAG,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CA0BxF;AAED;;;;;;;;;;;;;;;;;;;GAmBG;AACH,wBAAgB,WAAW,CAAC,CAAC,SAAS,UAAU,EAAE,CAAC,EAAE,CAAC,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,EAAE,GAAG,EAAE,CAAC,GAAG,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC,CAa1G;AAED;;;;;;;;;;;;;;;;;GAiBG;AACH,wBAAgB,iBAAiB,CAAC,CAAC,SAAS,UAAU,EAAE,OAAO,EAAE,CAAC,EAAE,EAAE,SAAS,GAAE,OAAe,GAAG,aAAa,CAAC,CAAC,CAAC,EAAE,CAkCpH"}