curvy 1.0.9 → 2.0.0-alpha.2

package/dist/shared/vector4-5KbjoHVj.mjs

@@ -0,0 +1,109 @@
+import { t as __exportAll } from "./chunk-25Hzd3z8.mjs";
+import { a as hadamard$1, c as make$1, d as softmax$1, f as subtract$1, i as equals$1, l as normalize$1, n as components$1, o as isVector4$1, p as transpose$1, r as dot$1, s as magnitude$1, t as add$1, u as scale$1 } from "./vector4.internal-5XcaX12X.mjs";
+//#region src/vector/vector4.ts
+var vector4_exports = /* @__PURE__ */ __exportAll({
+	add: () => add,
+	components: () => components,
+	dot: () => dot,
+	equals: () => equals,
+	hadamard: () => hadamard,
+	isVector4: () => isVector4,
+	magnitude: () => magnitude,
+	make: () => make,
+	normalize: () => normalize,
+	scale: () => scale,
+	softmax: () => softmax,
+	subtract: () => subtract,
+	transpose: () => transpose,
+	unit: () => unit,
+	zero: () => zero
+});
+/**
+* Checks if a value is a `Vector4`.
+*
+* @param v - The value to check.
+* @returns `true` if the value is a `Vector4`, `false` otherwise.
+* @since 1.0.0
+*/
+const isVector4 = isVector4$1;
+const equals = equals$1;
+const make = make$1;
+/**
+* Transposes a 4-tuple of items into one or more `Vector4`s, one per channel
+* extracted by the projection function.
+*
+* For each output channel, takes the corresponding number from each of the
+* four projected items and packs them into a `Vector4` — the data-layout
+* transformation NumPy calls a transpose, FP calls `unzip`, GPUs call a
+* gather, and DSP people inexplicably call deinterleaving.
+*
+* The projection's return-tuple arity determines the number of output
+* vectors, and that arity is preserved in the return type.
+*
+* @example
+* ```ts
+* // Pack the x and y coordinates of four control points into two Vector4s.
+* const [xs, ys] = Vector4.transpose(
+*   [p0, p1, p2, p3],
+*   (p) => [p.x, p.y],
+* )
+* ```
+*
+* @param inputs - Exactly four items of any type.
+* @param project - A function returning a fixed-arity tuple of numbers — one per output channel.
+* @returns A tuple of `Vector4`s with the same arity as the projection's return tuple.
+* @since 2.1.0
+*/
+const transpose = transpose$1;
+/**
+* Gets the magnitude of a `Vector4`.
+*
+* @param vector - The vector to calculate the magnitude of.
+* @returns The magnitude of the vector.
+* @since 1.0.0
+*/
+const magnitude = magnitude$1;
+/**
+* Normalizes a `Vector4`.
+*
+* @param vector - The vector to normalize.
+* @returns A new `Vector4` instance with the same direction as the input vector but with a magnitude of 1.
+* @since 1.0.0
+*/
+const normalize = normalize$1;
+const dot = dot$1;
+/**
+* Returns the components of a `Vector4` as an array.
+*
+* @param v - The vector to get the components of.
+* @returns An array containing the x, y, z, and w components of the vector.
+* @since 1.0.0
+*/
+const components = components$1;
+/**
+* Applies the softmax function to a `Vector4`.
+*
+* @param v - The vector to calculate the softmax of.
+* @returns A new `Vector4` instance representing the softmax of the input vector.
+*/
+const softmax = softmax$1;
+/**
+* Creates a zero vector.
+*
+* @since 1.0.0
+*/
+const zero = make(0);
+/**
+* Creates a unit vector.
+*
+* @since 1.0.0
+*/
+const unit = make(1);
+const add = add$1;
+const subtract = subtract$1;
+const hadamard = hadamard$1;
+const scale = scale$1;
+//#endregion
+export { transpose as n, vector4_exports as r, make as t };
+
+//# sourceMappingURL=vector4-5KbjoHVj.mjs.map

package/dist/shared/vector4-5KbjoHVj.mjs.map

@@ -0,0 +1 @@
+{"version":3,"file":"vector4-5KbjoHVj.mjs","names":["internal.isVector4","internal.equals","internal.make","internal.transpose","internal.magnitude","internal.normalize","internal.dot","internal.components","internal.softmax","internal.add","internal.subtract","internal.hadamard","internal.scale"],"sources":["../../src/vector/vector4.ts"],"sourcesContent":["import type { FourDimensional } from '../dimensions'\nimport type { Pipeable } from '../utils'\nimport type { Vector4TypeId } from './vector4.internal'\nimport * as internal from './vector4.internal'\n\n/**\n * A 4D vector.\n *\n * All fields are readonly and immutable, and all operations create new instances.\n *\n * @since 1.0.0\n */\nexport interface Vector4 extends Pipeable, FourDimensional<number> {\n  readonly [Vector4TypeId]: Vector4TypeId\n}\n\n/**\n * Checks if a value is a `Vector4`.\n *\n * @param v - The value to check.\n * @returns `true` if the value is a `Vector4`, `false` otherwise.\n * @since 1.0.0\n */\nexport const isVector4: (v: unknown) => v is Vector4 = internal.isVector4\n\nexport const equals: {\n  /**\n   * Checks if two `Vector4` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param a - The first vector.\n   * @param b - The second vector.\n   * @returns `true` when each pair of components is within tolerance.\n   * @since 1.1.0\n   */\n  (a: Vector4, b: Vector4): boolean\n  /**\n   * Checks if two `Vector4` instances are approximately equal within the\n   * default absolute tolerance ({@link EPSILON}).\n   *\n   * @param b - The second vector.\n   * @returns A function that takes the first vector and returns the comparison result.\n   * @since 1.1.0\n   */\n  (b: Vector4): (a: Vector4) => boolean\n} = internal.equals\n\nexport const make: {\n  /**\n   * Creates a new `Vector4` instance.\n   *\n   * @param xyzw - The x, y, z, and w components of the vector.\n   * @returns A new `Vector4` instance.\n   * @since 1.0.0\n   */\n  (xyzw: number): Vector4\n  /**\n   * Creates a new `Vector4` instance.\n   *\n   * @param x - The x component of the vector.\n   * @param yzw - The y, z, and w components of the vector.\n   * @returns A new `Vector4` instance.\n   * @since 1.0.0\n   */\n  (x: number, yzw: number): Vector4\n  /**\n   * Creates a new `Vector4` instance.\n   *\n   * @param x - The x component of the vector.\n   * @param y - The y component of the vector.\n   * @param zw - The z and w components of the vector.\n   * @returns A new `Vector4` instance.\n   * @since 1.0.0\n   */\n  (x: number, y: number, zw: number): Vector4\n  /**\n   * Creates a new `Vector4` instance.\n   *\n   * @param x - The x component of the vector.\n   * @param y - The y component of the vector.\n   * @param z - The z component of the vector.\n   * @param w - The w component of the vector.\n   * @returns A new `Vector4` instance.\n   * @since 1.0.0\n   */\n  (x: number, y: number, z: number, w: number): Vector4\n} = internal.make\n\n/**\n * Transposes a 4-tuple of items into one or more `Vector4`s, one per channel\n * extracted by the projection function.\n *\n * For each output channel, takes the corresponding number from each of the\n * four projected items and packs them into a `Vector4` — the data-layout\n * transformation NumPy calls a transpose, FP calls `unzip`, GPUs call a\n * gather, and DSP people inexplicably call deinterleaving.\n *\n * The projection's return-tuple arity determines the number of output\n * vectors, and that arity is preserved in the return type.\n *\n * @example\n * ```ts\n * // Pack the x and y coordinates of four control points into two Vector4s.\n * const [xs, ys] = Vector4.transpose(\n *   [p0, p1, p2, p3],\n *   (p) => [p.x, p.y],\n * )\n * ```\n *\n * @param inputs - Exactly four items of any type.\n * @param project - A function returning a fixed-arity tuple of numbers — one per output channel.\n * @returns A tuple of `Vector4`s with the same arity as the projection's return tuple.\n * @since 2.1.0\n */\nexport const transpose: <T, const Channels extends ReadonlyArray<number>>(\n  inputs: readonly [T, T, T, T],\n  project: (item: T) => Channels,\n) => { readonly [K in keyof Channels]: Vector4 } = internal.transpose\n\n/**\n * Gets the magnitude of a `Vector4`.\n *\n * @param vector - The vector to calculate the magnitude of.\n * @returns The magnitude of the vector.\n * @since 1.0.0\n */\nexport const magnitude: (vector: Vector4) => number = internal.magnitude\n\n/**\n * Normalizes a `Vector4`.\n *\n * @param vector - The vector to normalize.\n * @returns A new `Vector4` instance with the same direction as the input vector but with a magnitude of 1.\n * @since 1.0.0\n */\nexport const normalize: (vector: Vector4) => Vector4 = internal.normalize\n\nexport const dot: {\n  /**\n   * Computes the dot product of two `Vector4` instances.\n   *\n   * @param a - The first vector.\n   * @param b - The second vector.\n   * @returns The dot product of the two vectors.\n   * @since 1.0.0\n   */\n  (a: Vector4, b: Vector4): number\n  /**\n   * Computes the dot product of a `Vector4` instance and another vector.\n   *\n   * @param b - The second vector.\n   * @returns A function that takes the first vector and returns the dot product.\n   * @since 1.0.0\n   */\n  (b: Vector4): (a: Vector4) => number\n} = internal.dot\n\n/**\n * Returns the components of a `Vector4` as an array.\n *\n * @param v - The vector to get the components of.\n * @returns An array containing the x, y, z, and w components of the vector.\n * @since 1.0.0\n */\nexport const components: (v: Vector4) => [number, number, number, number] = internal.components\n\n/**\n * Applies the softmax function to a `Vector4`.\n *\n * @param v - The vector to calculate the softmax of.\n * @returns A new `Vector4` instance representing the softmax of the input vector.\n */\nexport const softmax: (v: Vector4) => Vector4 = internal.softmax\n\n/**\n * Creates a zero vector.\n *\n * @since 1.0.0\n */\nexport const zero = make(0)\n\n/**\n * Creates a unit vector.\n *\n * @since 1.0.0\n */\nexport const unit = make(1)\n\nexport const add: {\n  /**\n   * Adds two `Vector4` instances.\n   *\n   * @param a - The first vector.\n   * @param b - The second vector.\n   * @returns A new `Vector4` instance representing the sum of the two vectors.\n   * @since 1.0.0\n   */\n  (a: Vector4, b: Vector4): Vector4\n\n  /**\n   * Adds a `Vector4` instance to another vector.\n   *\n   * @param b - The second vector.\n   * @returns A function that takes the first vector and returns the sum.\n   * @since 1.0.0\n   */\n  (b: Vector4): (a: Vector4) => Vector4\n} = internal.add\n\nexport const subtract: {\n  /**\n   * Subtracts one `Vector4` instance from another.\n   *\n   * @param a - The first vector.\n   * @param b - The second vector.\n   * @returns A new `Vector4` instance representing the difference of the two vectors.\n   * @since 1.0.0\n   */\n  (a: Vector4, b: Vector4): Vector4\n  /**\n   * Subtracts a `Vector4` instance from another vector.\n   *\n   * @param b - The second vector.\n   * @returns A function that takes the first vector and returns the difference.\n   * @since 1.0.0\n   */\n  (b: Vector4): (a: Vector4) => Vector4\n} = internal.subtract\n\nexport const hadamard: {\n  /**\n   * Computes the Hadamard product of two `Vector4` instances.\n   *\n   * @param a - The first vector.\n   * @param b - The second vector.\n   * @returns A new `Vector4` instance representing the Hadamard product of the two vectors.\n   * @since 1.0.0\n   */\n  (a: Vector4, b: Vector4): Vector4\n  /**\n   * Computes the Hadamard product of a `Vector4` instance and another vector.\n   *\n   * @param b - The second vector.\n   * @returns A function that takes the first vector and returns the Hadamard product.\n   * @since 1.0.0\n   */\n  (b: Vector4): (a: Vector4) => Vector4\n} = internal.hadamard\n\nexport const scale: {\n  /**\n   * Scales a `Vector4` by a scalar.\n   *\n   * @param s - The scalar to scale the vector by.\n   * @returns A function that takes the vector and returns the scaled vector.\n   * @since 1.0.0\n   */\n  (s: number): (v: Vector4) => Vector4\n  /**\n   * Scales a `Vector4` by a scalar.\n   *\n   * @param s - The scalar to scale the vector by.\n   * @param v - The vector to scale.\n   * @returns A new `Vector4` instance representing the scaled vector.\n   * @since 1.0.0\n   */\n  (v: Vector4, s: number): Vector4\n} = internal.scale\n"],"mappings":";;;;;;;;;;;;;;;;;;;;;;;;;;;AAuBA,MAAa,YAA0CA;AAEvD,MAAa,SAoBTC;AAEJ,MAAa,OAuCTC;;;;;;;;;;;;;;;;;;;;;;;;;;;AA4BJ,MAAa,YAGsCC;;;;;;;;AASnD,MAAa,YAAyCC;;;;;;;;AAStD,MAAa,YAA0CC;AAEvD,MAAa,MAkBTC;;;;;;;;AASJ,MAAa,aAA+DC;;;;;;;AAQ5E,MAAa,UAAmCC;;;;;;AAOhD,MAAa,OAAO,KAAK,EAAE;;;;;;AAO3B,MAAa,OAAO,KAAK,EAAE;AAE3B,MAAa,MAmBTC;AAEJ,MAAa,WAkBTC;AAEJ,MAAa,WAkBTC;AAEJ,MAAa,QAkBTC"}

package/dist/shared/vector4.internal-5XcaX12X.mjs

@@ -0,0 +1,163 @@
+import { Pipeable, dual } from "../utils.mjs";
+import { epsEquals } from "../number.mjs";
+//#region src/vector/vector3.internal.ts
+const Vector3TypeId = Symbol.for("curvy/vector3");
+var Vector3Impl = class extends Pipeable {
+	[Vector3TypeId] = Vector3TypeId;
+	x;
+	y;
+	z;
+	constructor(v0 = 0, v1 = 0, v2 = 0) {
+		super();
+		this.x = v0;
+		this.y = v1;
+		this.z = v2;
+	}
+	get [0]() {
+		return this.x;
+	}
+	get [1]() {
+		return this.y;
+	}
+	get [2]() {
+		return this.z;
+	}
+	get [Symbol.toStringTag]() {
+		return `Vector3(${this.x}, ${this.y}, ${this.z})`;
+	}
+	get [Symbol.for("nodejs.util.inspect.custom")]() {
+		return this[Symbol.toStringTag];
+	}
+};
+const isVector3 = (v) => typeof v === "object" && v !== null && Vector3TypeId in v;
+const equals$1 = dual(2, (a, b) => epsEquals(a.x, b.x) && epsEquals(a.y, b.y) && epsEquals(a.z, b.z));
+const make$1 = (v0, v1 = v0, v2 = v1) => new Vector3Impl(v0, v1, v2);
+const transpose$1 = (inputs, project) => {
+	const a = project(inputs[0]);
+	const b = project(inputs[1]);
+	const c = project(inputs[2]);
+	const result = [];
+	for (let i = 0; i < a.length; i++) result.push(make$1(a[i], b[i], c[i]));
+	return result;
+};
+const fromSpherical = (r, theta, phi) => new Vector3Impl(r * Math.sin(theta) * Math.cos(phi), r * Math.sin(theta) * Math.sin(phi), r * Math.cos(theta));
+const magnitude$1 = (vector) => Math.hypot(vector.x, vector.y, vector.z);
+const normalize$1 = (vector) => {
+	const m = magnitude$1(vector);
+	return make$1(vector.x / m, vector.y / m, vector.z / m);
+};
+const dot$1 = dual(2, (a, b) => a.x * b.x + a.y * b.y + a.z * b.z);
+const cross = dual(2, (a, b) => {
+	return make$1(a.y * b.z - a.z * b.y, a.z * b.x - a.x * b.z, a.x * b.y - a.y * b.x);
+});
+const components$1 = (v) => [
+	v.x,
+	v.y,
+	v.z
+];
+const softmax$1 = (v) => {
+	const max = Math.max(v.x, v.y, v.z);
+	const v0 = Math.exp(v.x - max);
+	const v1 = Math.exp(v.y - max);
+	const v2 = Math.exp(v.z - max);
+	const sum = v0 + v1 + v2;
+	return make$1(v0 / sum, v1 / sum, v2 / sum);
+};
+const add$1 = dual(2, (a, b) => make$1(a.x + b.x, a.y + b.y, a.z + b.z));
+const subtract$1 = dual(2, (a, b) => make$1(a.x - b.x, a.y - b.y, a.z - b.z));
+const hadamard$1 = dual(2, (a, b) => make$1(a.x * b.x, a.y * b.y, a.z * b.z));
+const scale$1 = dual(2, (v, s) => s === 1 ? v : make$1(v.x * s, v.y * s, v.z * s));
+const getX = (v) => v.x;
+const setX = dual(2, (v, x) => make$1(x, v.y, v.z));
+const getY = (v) => v.y;
+const setY = dual(2, (v, y) => make$1(v.x, y, v.z));
+const getZ = (v) => v.z;
+const setZ = dual(2, (v, z) => make$1(v.x, v.y, z));
+const setR = dual(2, (v, r) => scale$1(r / magnitude$1(v))(v));
+const getTheta = (v) => Math.acos(v.z / magnitude$1(v));
+const setTheta = dual(2, (v, phi) => {
+	const m = magnitude$1(v);
+	const sinTheta = Math.sin(phi);
+	return make$1(m * sinTheta * Math.cos(v.y), m * sinTheta * Math.sin(v.y), m * Math.cos(phi));
+});
+const getPhi = (v) => Math.atan2(v.y, v.x);
+const setPhi = dual(2, (v, theta) => {
+	const m = magnitude$1(v);
+	return make$1(m * Math.cos(theta), m * Math.sin(theta), v.z);
+});
+//#endregion
+//#region src/vector/vector4.internal.ts
+const Vector4TypeId = Symbol.for("curvy/vector4");
+var Vector4Impl = class extends Pipeable {
+	[Vector4TypeId] = Vector4TypeId;
+	x;
+	y;
+	z;
+	w;
+	constructor(v0 = 0, v1 = 0, v2 = 0, v3 = 0) {
+		super();
+		this.x = v0;
+		this.y = v1;
+		this.z = v2;
+		this.w = v3;
+	}
+	get [0]() {
+		return this.x;
+	}
+	get [1]() {
+		return this.y;
+	}
+	get [2]() {
+		return this.z;
+	}
+	get [3]() {
+		return this.w;
+	}
+	get [Symbol.toStringTag]() {
+		return `Vector4(${this.x}, ${this.y}, ${this.z}, ${this.w})`;
+	}
+	get [Symbol.for("nodejs.util.inspect.custom")]() {
+		return this[Symbol.toStringTag];
+	}
+};
+const isVector4 = (v) => typeof v === "object" && v !== null && Vector4TypeId in v;
+const equals = dual(2, (a, b) => epsEquals(a.x, b.x) && epsEquals(a.y, b.y) && epsEquals(a.z, b.z) && epsEquals(a.w, b.w));
+const make = (v0, v1 = v0, v2 = v1, v3 = v2) => new Vector4Impl(v0, v1, v2, v3);
+const transpose = (inputs, project) => {
+	const a = project(inputs[0]);
+	const b = project(inputs[1]);
+	const c = project(inputs[2]);
+	const d = project(inputs[3]);
+	const result = [];
+	for (let i = 0; i < a.length; i++) result.push(make(a[i], b[i], c[i], d[i]));
+	return result;
+};
+const magnitude = (vector) => Math.hypot(vector.x, vector.y, vector.z, vector.w);
+const normalize = (vector) => {
+	const m = magnitude(vector);
+	return make(vector.x / m, vector.y / m, vector.z / m, vector.w / m);
+};
+const dot = dual(2, (a, b) => a.x * b.x + a.y * b.y + a.z * b.z + a.w * b.w);
+const components = (v) => [
+	v.x,
+	v.y,
+	v.z,
+	v.w
+];
+const softmax = (v) => {
+	const max = Math.max(v.x, v.y, v.z, v.w);
+	const v0 = Math.exp(v.x - max);
+	const v1 = Math.exp(v.y - max);
+	const v2 = Math.exp(v.z - max);
+	const v3 = Math.exp(v.w - max);
+	const sum = v0 + v1 + v2 + v3;
+	return make(v0 / sum, v1 / sum, v2 / sum, v3 / sum);
+};
+const add = dual(2, (a, b) => make(a.x + b.x, a.y + b.y, a.z + b.z, a.w + b.w));
+const subtract = dual(2, (a, b) => make(a.x - b.x, a.y - b.y, a.z - b.z, a.w - b.w));
+const hadamard = dual(2, (a, b) => make(a.x * b.x, a.y * b.y, a.z * b.z, a.w * b.w));
+const scale = dual(2, (v, s) => s === 1 ? v : make(v.x * s, v.y * s, v.z * s, v.w * s));
+//#endregion
+export { scale$1 as A, getY as C, magnitude$1 as D, isVector3 as E, setY as F, setZ as I, softmax$1 as L, setR as M, setTheta as N, make$1 as O, setX as P, subtract$1 as R, getX as S, hadamard$1 as T, dot$1 as _, hadamard as a, getPhi as b, make as c, softmax as d, subtract as f, cross as g, components$1 as h, equals as i, setPhi as j, normalize$1 as k, normalize as l, add$1 as m, components as n, isVector4 as o, transpose as p, dot as r, magnitude as s, add as t, scale as u, equals$1 as v, getZ as w, getTheta as x, fromSpherical as y, transpose$1 as z };
+
+//# sourceMappingURL=vector4.internal-5XcaX12X.mjs.map

package/dist/shared/vector4.internal-5XcaX12X.mjs.map

@@ -0,0 +1 @@
+{"version":3,"file":"vector4.internal-5XcaX12X.mjs","names":["equals","make","transpose","magnitude","normalize","dot","components","softmax","add","subtract","hadamard","scale"],"sources":["../../src/vector/vector3.internal.ts","../../src/vector/vector4.internal.ts"],"sourcesContent":["import { dual, Pipeable } from '../utils'\nimport { epsEquals } from '../number'\nimport type { Vector3 } from './vector3'\n\nexport const Vector3TypeId: unique symbol = Symbol.for('curvy/vector3')\nexport type Vector3TypeId = typeof Vector3TypeId\n\nclass Vector3Impl extends Pipeable implements Vector3 {\n  readonly [Vector3TypeId]: Vector3TypeId = Vector3TypeId\n\n  readonly x: number\n  readonly y: number\n  readonly z: number\n\n  constructor(v0 = 0, v1 = 0, v2 = 0) {\n    super()\n    this.x = v0\n    this.y = v1\n    this.z = v2\n  }\n\n  get [0]() {\n    return this.x\n  }\n\n  get [1]() {\n    return this.y\n  }\n\n  get [2]() {\n    return this.z\n  }\n\n  get [Symbol.toStringTag]() {\n    return `Vector3(${this.x}, ${this.y}, ${this.z})`\n  }\n\n  get [Symbol.for('nodejs.util.inspect.custom')]() {\n    return this[Symbol.toStringTag]\n  }\n}\n\nexport const isVector3 = (v: unknown): v is Vector3 =>\n  typeof v === 'object' && v !== null && Vector3TypeId in v\n\nexport const equals = dual<\n  (b: Vector3) => (a: Vector3) => boolean,\n  (a: Vector3, b: Vector3) => boolean\n>(2, (a: Vector3, b: Vector3) => epsEquals(a.x, b.x) && epsEquals(a.y, b.y) && epsEquals(a.z, b.z))\n\nexport const make = (v0: number, v1 = v0, v2 = v1): Vector3 => new Vector3Impl(v0, v1, v2)\n\nexport const transpose = <T, const Channels extends ReadonlyArray<number>>(\n  inputs: readonly [T, T, T],\n  project: (item: T) => Channels,\n): { readonly [K in keyof Channels]: Vector3 } => {\n  const a = project(inputs[0])\n  const b = project(inputs[1])\n  const c = project(inputs[2])\n  const result: Array<Vector3> = []\n  for (let i = 0; i < a.length; i++) {\n    result.push(make(a[i] as number, b[i] as number, c[i] as number))\n  }\n  return result as never\n}\n\nexport const fromSpherical = (r: number, theta: number, phi: number) =>\n  new Vector3Impl(\n    r * Math.sin(theta) * Math.cos(phi),\n    r * Math.sin(theta) * Math.sin(phi),\n    r * Math.cos(theta),\n  )\n\nexport const magnitude = (vector: Vector3) => Math.hypot(vector.x, vector.y, vector.z)\n\nexport const normalize = (vector: Vector3) => {\n  const m = magnitude(vector)\n\n  return make(vector.x / m, vector.y / m, vector.z / m)\n}\n\nexport const dot = dual<(b: Vector3) => (a: Vector3) => number, (a: Vector3, b: Vector3) => number>(\n  2,\n  (a: Vector3, b: Vector3) => a.x * b.x + a.y * b.y + a.z * b.z,\n)\n\nexport const cross = dual<\n  (b: Vector3) => (a: Vector3) => Vector3,\n  (a: Vector3, b: Vector3) => Vector3\n>(2, (a: Vector3, b: Vector3) => {\n  const x = a.y * b.z - a.z * b.y\n  const y = a.z * b.x - a.x * b.z\n  const z = a.x * b.y - a.y * b.x\n\n  return make(x, y, z)\n})\n\nexport const components: (v: Vector3) => [number, number, number] = (v: Vector3) => [v.x, v.y, v.z]\n\nexport const softmax = (v: Vector3) => {\n  const max = Math.max(v.x, v.y, v.z)\n\n  const v0 = Math.exp(v.x - max)\n  const v1 = Math.exp(v.y - max)\n  const v2 = Math.exp(v.z - max)\n\n  const sum = v0 + v1 + v2\n\n  return make(v0 / sum, v1 / sum, v2 / sum)\n}\n\nexport const add = dual<\n  (b: Vector3) => (a: Vector3) => Vector3,\n  (a: Vector3, b: Vector3) => Vector3\n>(2, (a: Vector3, b: Vector3) => make(a.x + b.x, a.y + b.y, a.z + b.z))\n\nexport const subtract = dual<\n  (b: Vector3) => (a: Vector3) => Vector3,\n  (a: Vector3, b: Vector3) => Vector3\n>(2, (a: Vector3, b: Vector3) => make(a.x - b.x, a.y - b.y, a.z - b.z))\n\nexport const hadamard = dual<\n  (b: Vector3) => (a: Vector3) => Vector3,\n  (a: Vector3, b: Vector3) => Vector3\n>(2, (a: Vector3, b: Vector3) => make(a.x * b.x, a.y * b.y, a.z * b.z))\n\nexport const scale = dual<\n  (s: number) => (v: Vector3) => Vector3,\n  (v: Vector3, s: number) => Vector3\n>(2, (v: Vector3, s: number) => (s === 1 ? v : make(v.x * s, v.y * s, v.z * s)))\n\nexport const getX: (v: Vector3) => number = (v: Vector3) => v.x\nexport const setX = dual<\n  (x: number) => (v: Vector3) => Vector3,\n  (v: Vector3, x: number) => Vector3\n>(2, (v: Vector3, x: number) => make(x, v.y, v.z))\n\nexport const getY: (v: Vector3) => number = (v: Vector3) => v.y\nexport const setY = dual<\n  (y: number) => (v: Vector3) => Vector3,\n  (v: Vector3, y: number) => Vector3\n>(2, (v: Vector3, y: number) => make(v.x, y, v.z))\n\nexport const getZ: (v: Vector3) => number = (v: Vector3) => v.z\nexport const setZ = dual<\n  (z: number) => (v: Vector3) => Vector3,\n  (v: Vector3, z: number) => Vector3\n>(2, (v: Vector3, z: number) => make(v.x, v.y, z))\n\nexport const setR = dual<\n  (r: number) => (v: Vector3) => Vector3,\n  (v: Vector3, r: number) => Vector3\n>(2, (v: Vector3, r: number) => scale(r / magnitude(v))(v))\n\nexport const getTheta = (v: Vector3) => Math.acos(v.z / magnitude(v))\nexport const setTheta = dual<\n  (phi: number) => (v: Vector3) => Vector3,\n  (v: Vector3, phi: number) => Vector3\n>(2, (v: Vector3, phi: number) => {\n  const m = magnitude(v)\n  const sinTheta = Math.sin(phi)\n  return make(m * sinTheta * Math.cos(v.y), m * sinTheta * Math.sin(v.y), m * Math.cos(phi))\n})\n\nexport const getPhi = (v: Vector3) => Math.atan2(v.y, v.x)\nexport const setPhi = dual<\n  (theta: number) => (v: Vector3) => Vector3,\n  (v: Vector3, theta: number) => Vector3\n>(2, (v: Vector3, theta: number) => {\n  const m = magnitude(v)\n  return make(m * Math.cos(theta), m * Math.sin(theta), v.z)\n})\n","import { dual, Pipeable } from '../utils'\nimport { epsEquals } from '../number'\nimport type { Vector4 } from './vector4'\n\nexport const Vector4TypeId: unique symbol = Symbol.for('curvy/vector4')\nexport type Vector4TypeId = typeof Vector4TypeId\n\nclass Vector4Impl extends Pipeable implements Vector4 {\n  readonly [Vector4TypeId]: Vector4TypeId = Vector4TypeId\n\n  readonly x: number\n  readonly y: number\n  readonly z: number\n  readonly w: number\n\n  constructor(v0 = 0, v1 = 0, v2 = 0, v3 = 0) {\n    super()\n\n    this.x = v0\n    this.y = v1\n    this.z = v2\n    this.w = v3\n  }\n\n  get [0]() {\n    return this.x\n  }\n\n  get [1]() {\n    return this.y\n  }\n\n  get [2]() {\n    return this.z\n  }\n\n  get [3]() {\n    return this.w\n  }\n\n  get [Symbol.toStringTag]() {\n    return `Vector4(${this.x}, ${this.y}, ${this.z}, ${this.w})`\n  }\n\n  get [Symbol.for('nodejs.util.inspect.custom')]() {\n    return this[Symbol.toStringTag]\n  }\n}\n\nexport const isVector4 = (v: unknown): v is Vector4 =>\n  typeof v === 'object' && v !== null && Vector4TypeId in v\n\nexport const equals = dual<\n  (b: Vector4) => (a: Vector4) => boolean,\n  (a: Vector4, b: Vector4) => boolean\n>(\n  2,\n  (a: Vector4, b: Vector4) =>\n    epsEquals(a.x, b.x) && epsEquals(a.y, b.y) && epsEquals(a.z, b.z) && epsEquals(a.w, b.w),\n)\n\nexport const make = (v0: number, v1 = v0, v2 = v1, v3 = v2): Vector4 =>\n  new Vector4Impl(v0, v1, v2, v3)\n\nexport const transpose = <T, const Channels extends ReadonlyArray<number>>(\n  inputs: readonly [T, T, T, T],\n  project: (item: T) => Channels,\n): { readonly [K in keyof Channels]: Vector4 } => {\n  const a = project(inputs[0])\n  const b = project(inputs[1])\n  const c = project(inputs[2])\n  const d = project(inputs[3])\n  const result: Array<Vector4> = []\n  for (let i = 0; i < a.length; i++) {\n    result.push(make(a[i] as number, b[i] as number, c[i] as number, d[i] as number))\n  }\n  return result as never\n}\n\nexport const magnitude: (vector: Vector4) => number = (vector: Vector4) =>\n  Math.hypot(vector.x, vector.y, vector.z, vector.w)\n\nexport const normalize: (vector: Vector4) => Vector4 = (vector: Vector4) => {\n  const m = magnitude(vector)\n\n  return make(vector.x / m, vector.y / m, vector.z / m, vector.w / m)\n}\n\nexport const dot = dual<(b: Vector4) => (a: Vector4) => number, (a: Vector4, b: Vector4) => number>(\n  2,\n  (a: Vector4, b: Vector4) => a.x * b.x + a.y * b.y + a.z * b.z + a.w * b.w,\n)\n\nexport const components: (v: Vector4) => [number, number, number, number] = (v: Vector4) => [\n  v.x,\n  v.y,\n  v.z,\n  v.w,\n]\n\nexport const softmax: (v: Vector4) => Vector4 = (v: Vector4) => {\n  const max = Math.max(v.x, v.y, v.z, v.w)\n\n  const v0 = Math.exp(v.x - max)\n  const v1 = Math.exp(v.y - max)\n  const v2 = Math.exp(v.z - max)\n  const v3 = Math.exp(v.w - max)\n\n  const sum = v0 + v1 + v2 + v3\n\n  return make(v0 / sum, v1 / sum, v2 / sum, v3 / sum)\n}\n\nexport const add = dual<\n  (b: Vector4) => (a: Vector4) => Vector4,\n  (a: Vector4, b: Vector4) => Vector4\n>(2, (a: Vector4, b: Vector4) => make(a.x + b.x, a.y + b.y, a.z + b.z, a.w + b.w))\n\nexport const subtract = dual<\n  (b: Vector4) => (a: Vector4) => Vector4,\n  (a: Vector4, b: Vector4) => Vector4\n>(2, (a: Vector4, b: Vector4) => make(a.x - b.x, a.y - b.y, a.z - b.z, a.w - b.w))\n\nexport const hadamard = dual<\n  (b: Vector4) => (a: Vector4) => Vector4,\n  (a: Vector4, b: Vector4) => Vector4\n>(2, (a: Vector4, b: Vector4) => make(a.x * b.x, a.y * b.y, a.z * b.z, a.w * b.w))\nexport const scale = dual<\n  (s: number) => (v: Vector4) => Vector4,\n  (v: Vector4, s: number) => Vector4\n>(2, (v: Vector4, s: number) => (s === 1 ? v : make(v.x * s, v.y * s, v.z * s, v.w * s)))\n"],"mappings":";;;AAIA,MAAa,gBAA+B,OAAO,IAAI,gBAAgB;AAGvE,IAAM,cAAN,cAA0B,SAA4B;CACpD,CAAU,iBAAgC;CAE1C;CACA;CACA;CAEA,YAAY,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG;EAClC,OAAO;EACP,KAAK,IAAI;EACT,KAAK,IAAI;EACT,KAAK,IAAI;;CAGX,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,OAAO,eAAe;EACzB,OAAO,WAAW,KAAK,EAAE,IAAI,KAAK,EAAE,IAAI,KAAK,EAAE;;CAGjD,KAAK,OAAO,IAAI,6BAA6B,IAAI;EAC/C,OAAO,KAAK,OAAO;;;AAIvB,MAAa,aAAa,MACxB,OAAO,MAAM,YAAY,MAAM,QAAQ,iBAAiB;AAE1D,MAAaA,WAAS,KAGpB,IAAI,GAAY,MAAe,UAAU,EAAE,GAAG,EAAE,EAAE,IAAI,UAAU,EAAE,GAAG,EAAE,EAAE,IAAI,UAAU,EAAE,GAAG,EAAE,EAAE,CAAC;AAEnG,MAAaC,UAAQ,IAAY,KAAK,IAAI,KAAK,OAAgB,IAAI,YAAY,IAAI,IAAI,GAAG;AAE1F,MAAaC,eACX,QACA,YACgD;CAChD,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,SAAyB,EAAE;CACjC,KAAK,IAAI,IAAI,GAAG,IAAI,EAAE,QAAQ,KAC5B,OAAO,KAAKD,OAAK,EAAE,IAAc,EAAE,IAAc,EAAE,GAAa,CAAC;CAEnE,OAAO;;AAGT,MAAa,iBAAiB,GAAW,OAAe,QACtD,IAAI,YACF,IAAI,KAAK,IAAI,MAAM,GAAG,KAAK,IAAI,IAAI,EACnC,IAAI,KAAK,IAAI,MAAM,GAAG,KAAK,IAAI,IAAI,EACnC,IAAI,KAAK,IAAI,MAAM,CACpB;AAEH,MAAaE,eAAa,WAAoB,KAAK,MAAM,OAAO,GAAG,OAAO,GAAG,OAAO,EAAE;AAEtF,MAAaC,eAAa,WAAoB;CAC5C,MAAM,IAAID,YAAU,OAAO;CAE3B,OAAOF,OAAK,OAAO,IAAI,GAAG,OAAO,IAAI,GAAG,OAAO,IAAI,EAAE;;AAGvD,MAAaI,QAAM,KACjB,IACC,GAAY,MAAe,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,EAC7D;AAED,MAAa,QAAQ,KAGnB,IAAI,GAAY,MAAe;CAK/B,OAAOJ,OAJG,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,GACpB,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,GACpB,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,EAEV;EACpB;AAEF,MAAaK,gBAAwD,MAAe;CAAC,EAAE;CAAG,EAAE;CAAG,EAAE;CAAE;AAEnG,MAAaC,aAAW,MAAe;CACrC,MAAM,MAAM,KAAK,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE;CAEnC,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAC9B,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAC9B,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAE9B,MAAM,MAAM,KAAK,KAAK;CAEtB,OAAON,OAAK,KAAK,KAAK,KAAK,KAAK,KAAK,IAAI;;AAG3C,MAAaO,QAAM,KAGjB,IAAI,GAAY,MAAeP,OAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;AAEvE,MAAaQ,aAAW,KAGtB,IAAI,GAAY,MAAeR,OAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;AAEvE,MAAaS,aAAW,KAGtB,IAAI,GAAY,MAAeT,OAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;AAEvE,MAAaU,UAAQ,KAGnB,IAAI,GAAY,MAAe,MAAM,IAAI,IAAIV,OAAK,EAAE,IAAI,GAAG,EAAE,IAAI,GAAG,EAAE,IAAI,EAAE,CAAE;AAEhF,MAAa,QAAgC,MAAe,EAAE;AAC9D,MAAa,OAAO,KAGlB,IAAI,GAAY,MAAcA,OAAK,GAAG,EAAE,GAAG,EAAE,EAAE,CAAC;AAElD,MAAa,QAAgC,MAAe,EAAE;AAC9D,MAAa,OAAO,KAGlB,IAAI,GAAY,MAAcA,OAAK,EAAE,GAAG,GAAG,EAAE,EAAE,CAAC;AAElD,MAAa,QAAgC,MAAe,EAAE;AAC9D,MAAa,OAAO,KAGlB,IAAI,GAAY,MAAcA,OAAK,EAAE,GAAG,EAAE,GAAG,EAAE,CAAC;AAElD,MAAa,OAAO,KAGlB,IAAI,GAAY,MAAcU,QAAM,IAAIR,YAAU,EAAE,CAAC,CAAC,EAAE,CAAC;AAE3D,MAAa,YAAY,MAAe,KAAK,KAAK,EAAE,IAAIA,YAAU,EAAE,CAAC;AACrE,MAAa,WAAW,KAGtB,IAAI,GAAY,QAAgB;CAChC,MAAM,IAAIA,YAAU,EAAE;CACtB,MAAM,WAAW,KAAK,IAAI,IAAI;CAC9B,OAAOF,OAAK,IAAI,WAAW,KAAK,IAAI,EAAE,EAAE,EAAE,IAAI,WAAW,KAAK,IAAI,EAAE,EAAE,EAAE,IAAI,KAAK,IAAI,IAAI,CAAC;EAC1F;AAEF,MAAa,UAAU,MAAe,KAAK,MAAM,EAAE,GAAG,EAAE,EAAE;AAC1D,MAAa,SAAS,KAGpB,IAAI,GAAY,UAAkB;CAClC,MAAM,IAAIE,YAAU,EAAE;CACtB,OAAOF,OAAK,IAAI,KAAK,IAAI,MAAM,EAAE,IAAI,KAAK,IAAI,MAAM,EAAE,EAAE,EAAE;EAC1D;;;ACvKF,MAAa,gBAA+B,OAAO,IAAI,gBAAgB;AAGvE,IAAM,cAAN,cAA0B,SAA4B;CACpD,CAAU,iBAAgC;CAE1C;CACA;CACA;CACA;CAEA,YAAY,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG,KAAK,GAAG;EAC1C,OAAO;EAEP,KAAK,IAAI;EACT,KAAK,IAAI;EACT,KAAK,IAAI;EACT,KAAK,IAAI;;CAGX,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,KAAK;EACR,OAAO,KAAK;;CAGd,KAAK,OAAO,eAAe;EACzB,OAAO,WAAW,KAAK,EAAE,IAAI,KAAK,EAAE,IAAI,KAAK,EAAE,IAAI,KAAK,EAAE;;CAG5D,KAAK,OAAO,IAAI,6BAA6B,IAAI;EAC/C,OAAO,KAAK,OAAO;;;AAIvB,MAAa,aAAa,MACxB,OAAO,MAAM,YAAY,MAAM,QAAQ,iBAAiB;AAE1D,MAAa,SAAS,KAIpB,IACC,GAAY,MACX,UAAU,EAAE,GAAG,EAAE,EAAE,IAAI,UAAU,EAAE,GAAG,EAAE,EAAE,IAAI,UAAU,EAAE,GAAG,EAAE,EAAE,IAAI,UAAU,EAAE,GAAG,EAAE,EAAE,CAC3F;AAED,MAAa,QAAQ,IAAY,KAAK,IAAI,KAAK,IAAI,KAAK,OACtD,IAAI,YAAY,IAAI,IAAI,IAAI,GAAG;AAEjC,MAAa,aACX,QACA,YACgD;CAChD,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,IAAI,QAAQ,OAAO,GAAG;CAC5B,MAAM,SAAyB,EAAE;CACjC,KAAK,IAAI,IAAI,GAAG,IAAI,EAAE,QAAQ,KAC5B,OAAO,KAAK,KAAK,EAAE,IAAc,EAAE,IAAc,EAAE,IAAc,EAAE,GAAa,CAAC;CAEnF,OAAO;;AAGT,MAAa,aAA0C,WACrD,KAAK,MAAM,OAAO,GAAG,OAAO,GAAG,OAAO,GAAG,OAAO,EAAE;AAEpD,MAAa,aAA2C,WAAoB;CAC1E,MAAM,IAAI,UAAU,OAAO;CAE3B,OAAO,KAAK,OAAO,IAAI,GAAG,OAAO,IAAI,GAAG,OAAO,IAAI,GAAG,OAAO,IAAI,EAAE;;AAGrE,MAAa,MAAM,KACjB,IACC,GAAY,MAAe,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,EACzE;AAED,MAAa,cAAgE,MAAe;CAC1F,EAAE;CACF,EAAE;CACF,EAAE;CACF,EAAE;CACH;AAED,MAAa,WAAoC,MAAe;CAC9D,MAAM,MAAM,KAAK,IAAI,EAAE,GAAG,EAAE,GAAG,EAAE,GAAG,EAAE,EAAE;CAExC,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAC9B,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAC9B,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAC9B,MAAM,KAAK,KAAK,IAAI,EAAE,IAAI,IAAI;CAE9B,MAAM,MAAM,KAAK,KAAK,KAAK;CAE3B,OAAO,KAAK,KAAK,KAAK,KAAK,KAAK,KAAK,KAAK,KAAK,IAAI;;AAGrD,MAAa,MAAM,KAGjB,IAAI,GAAY,MAAe,KAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;AAElF,MAAa,WAAW,KAGtB,IAAI,GAAY,MAAe,KAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;AAElF,MAAa,WAAW,KAGtB,IAAI,GAAY,MAAe,KAAK,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,GAAG,EAAE,IAAI,EAAE,EAAE,CAAC;AAClF,MAAa,QAAQ,KAGnB,IAAI,GAAY,MAAe,MAAM,IAAI,IAAI,KAAK,EAAE,IAAI,GAAG,EAAE,IAAI,GAAG,EAAE,IAAI,GAAG,EAAE,IAAI,EAAE,CAAE"}

package/dist/solution/index.d.mts

@@ -0,0 +1,2 @@
+import { a as solution_d_exports } from "../shared/solution-sNQ-ZVhN.mjs";
+export { solution_d_exports as Solution };

package/dist/solution/index.mjs

@@ -0,0 +1,2 @@
+import { c as solution_exports } from "../shared/solution-CrtLfqEq.mjs";
+export { solution_exports as Solution };