@fluixi/utils 1.0.0-alpha.53

package/dist/lib/maths/fn/complex.js

@@ -0,0 +1,383 @@
+/**
+ * Complex number operations.
+ * @module ComplexFunctions
+ */
+// Constants
+const EPSILON = 1e-10;
+const TWO_PI = 2 * Math.PI;
+/**
+ * Checks if a value is approximately zero within epsilon.
+ * @param value - The value to check
+ * @returns True if the value is effectively zero
+ */
+function isApproximatelyZero(value) {
+    return Math.abs(value) < EPSILON;
+}
+/**
+ * Safely compares two numbers with epsilon precision.
+ * @param a - First number
+ * @param b - Second number
+ * @returns True if numbers are approximately equal
+ */
+function approximatelyEqual(a, b) {
+    return Math.abs(a - b) < EPSILON;
+}
+/**
+ * Normalizes a complex number by handling near-zero values.
+ * @param c - The complex number to normalize
+ * @returns Normalized complex number
+ */
+function normalizeComplex(c) {
+    return {
+        real: isApproximatelyZero(c.real) ? 0 : c.real,
+        imag: isApproximatelyZero(c.imag) ? 0 : c.imag,
+    };
+}
+/**
+ * Validates that a value is a finite number.
+ * @param value - The value to validate
+ * @param name - Name of the value for error messages
+ * @throws {Error} If value is not a finite number
+ */
+function validateNumber(value, name) {
+    if (typeof value !== 'number' || !Number.isFinite(value)) {
+        throw new Error(`${name} must be a finite number, got ${value}`);
+    }
+}
+/**
+ * Validates a complex number.
+ * @param c - The complex number to validate
+ * @param name - Name of the complex number for error messages
+ * @throws {Error} If complex number is invalid
+ */
+function validateComplexNumber(c, name = 'Complex number') {
+    if (!c || typeof c !== 'object') {
+        throw new Error(`${name} must be an object`);
+    }
+    validateNumber(c.real, `${name}.real`);
+    validateNumber(c.imag, `${name}.imag`);
+}
+/**
+ * Creates a complex number with validation and normalization.
+ * @param real - Real part
+ * @param imag - Imaginary part
+ * @returns Validated and normalized complex number
+ */
+export function createComplex(real, imag = 0) {
+    validateNumber(real, 'real');
+    validateNumber(imag, 'imag');
+    return normalizeComplex({ real, imag });
+}
+/**
+ * Checks if two complex numbers are approximately equal.
+ * @param a - First complex number
+ * @param b - Second complex number
+ * @returns True if numbers are approximately equal
+ */
+export function complexEquals(a, b) {
+    validateComplexNumber(a, 'First complex number');
+    validateComplexNumber(b, 'Second complex number');
+    return (approximatelyEqual(a.real, b.real) && approximatelyEqual(a.imag, b.imag));
+}
+/**
+ * Adds two complex numbers.
+ * @param a - First complex number
+ * @param b - Second complex number
+ * @returns The sum as a complex number
+ */
+export function addComplex(a, b) {
+    validateComplexNumber(a, 'First complex number');
+    validateComplexNumber(b, 'Second complex number');
+    return normalizeComplex({
+        real: a.real + b.real,
+        imag: a.imag + b.imag,
+    });
+}
+/**
+ * Subtracts two complex numbers.
+ * @param a - First complex number
+ * @param b - Second complex number
+ * @returns The difference as a complex number (a - b)
+ */
+export function subtractComplex(a, b) {
+    validateComplexNumber(a, 'First complex number');
+    validateComplexNumber(b, 'Second complex number');
+    return normalizeComplex({
+        real: a.real - b.real,
+        imag: a.imag - b.imag,
+    });
+}
+/**
+ * Multiplies two complex numbers.
+ * @param a - First complex number
+ * @param b - Second complex number
+ * @returns The product as a complex number
+ */
+export function multiplyComplex(a, b) {
+    validateComplexNumber(a, 'First complex number');
+    validateComplexNumber(b, 'Second complex number');
+    return normalizeComplex({
+        real: a.real * b.real - a.imag * b.imag,
+        imag: a.real * b.imag + a.imag * b.real,
+    });
+}
+/**
+ * Divides two complex numbers.
+ * @param a - Numerator complex number
+ * @param b - Denominator complex number
+ * @returns The quotient as a complex number
+ * @throws {Error} If dividing by zero or near-zero complex number
+ */
+export function divideComplex(a, b) {
+    validateComplexNumber(a, 'Numerator complex number');
+    validateComplexNumber(b, 'Denominator complex number');
+    const denominator = b.real * b.real + b.imag * b.imag;
+    if (isApproximatelyZero(denominator)) {
+        throw new Error('Cannot divide by zero or near-zero complex number');
+    }
+    return normalizeComplex({
+        real: (a.real * b.real + a.imag * b.imag) / denominator,
+        imag: (a.imag * b.real - a.real * b.imag) / denominator,
+    });
+}
+/**
+ * Calculates the magnitude (absolute value) of a complex number.
+ * @param c - The complex number
+ * @returns The magnitude (non-negative)
+ */
+export function magnitudeComplex(c) {
+    validateComplexNumber(c);
+    const magnitude = Math.sqrt(c.real * c.real + c.imag * c.imag);
+    return isApproximatelyZero(magnitude) ? 0 : magnitude;
+}
+/**
+ * Calculates the conjugate of a complex number.
+ * @param c - The complex number
+ * @returns The complex conjugate
+ */
+export function conjugateComplex(c) {
+    validateComplexNumber(c);
+    return normalizeComplex({
+        real: c.real,
+        imag: -c.imag,
+    });
+}
+/**
+ * Calculates the argument (angle) of a complex number in radians.
+ * @param c - The complex number
+ * @returns The principal argument in radians (-π, π]
+ */
+export function argumentComplex(c) {
+    validateComplexNumber(c);
+    if (isApproximatelyZero(c.real) && isApproximatelyZero(c.imag)) {
+        return 0; // Argument of zero is conventionally 0
+    }
+    return Math.atan2(c.imag, c.real);
+}
+/**
+ * Converts a complex number to polar form.
+ * @param c - The complex number
+ * @returns Polar representation {magnitude, angle}
+ */
+export function toPolarComplex(c) {
+    validateComplexNumber(c);
+    return {
+        magnitude: magnitudeComplex(c),
+        angle: argumentComplex(c),
+    };
+}
+/**
+ * Converts polar form to a complex number.
+ * @param magnitude - The magnitude (must be non-negative)
+ * @param angle - The angle in radians
+ * @returns The complex number
+ * @throws {Error} If magnitude is negative
+ */
+export function fromPolarComplex(magnitude, angle) {
+    validateNumber(magnitude, 'magnitude');
+    validateNumber(angle, 'angle');
+    if (magnitude < -EPSILON) {
+        throw new Error('Magnitude must be non-negative');
+    }
+    // Normalize angle to (-π, π]
+    let normalizedAngle = angle % TWO_PI;
+    if (normalizedAngle > Math.PI)
+        normalizedAngle -= TWO_PI;
+    if (normalizedAngle <= -Math.PI)
+        normalizedAngle += TWO_PI;
+    const mag = Math.max(0, magnitude); // Ensure non-negative
+    return normalizeComplex({
+        real: mag * Math.cos(normalizedAngle),
+        imag: mag * Math.sin(normalizedAngle),
+    });
+}
+/**
+ * Calculates the square root of a complex number.
+ * @param c - The complex number
+ * @returns Array of two complex square roots (principal and secondary)
+ */
+export function sqrtComplex(c) {
+    validateComplexNumber(c);
+    // Special case for zero
+    if (isApproximatelyZero(c.real) && isApproximatelyZero(c.imag)) {
+        const zero = createComplex(0, 0);
+        return [zero, zero];
+    }
+    const magnitude = magnitudeComplex(c);
+    const angle = argumentComplex(c);
+    const rootMagnitude = Math.sqrt(magnitude);
+    const rootAngle1 = angle / 2;
+    const rootAngle2 = rootAngle1 + Math.PI;
+    return [
+        fromPolarComplex(rootMagnitude, rootAngle1),
+        fromPolarComplex(rootMagnitude, rootAngle2),
+    ];
+}
+/**
+ * Calculates the power of a complex number with a real exponent.
+ * @param c - The complex number
+ * @param exponent - The exponent (real number)
+ * @returns The complex number raised to the power
+ */
+export function powerComplex(c, exponent) {
+    validateComplexNumber(c);
+    validateNumber(exponent, 'exponent');
+    // Special cases for optimization
+    if (isApproximatelyZero(exponent)) {
+        return createComplex(1, 0); // c^0 = 1
+    }
+    if (approximatelyEqual(exponent, 1)) {
+        return createComplex(c.real, c.imag); // c^1 = c
+    }
+    if (approximatelyEqual(exponent, 2)) {
+        // Square directly for better performance
+        return multiplyComplex(c, c);
+    }
+    if (isApproximatelyZero(c.real) && isApproximatelyZero(c.imag)) {
+        if (exponent > 0) {
+            return createComplex(0, 0); // 0^n = 0 for n > 0
+        }
+        else if (exponent < 0) {
+            throw new Error('Cannot raise zero to a negative power');
+        }
+    }
+    const polar = toPolarComplex(c);
+    const magnitude = Math.pow(polar.magnitude, exponent);
+    // Handle negative magnitude due to floating point errors
+    const safeMagnitude = magnitude < -EPSILON ? Math.abs(magnitude) : Math.max(0, magnitude);
+    const angle = polar.angle * exponent;
+    return fromPolarComplex(safeMagnitude, angle);
+}
+/**
+ * Calculates the natural logarithm of a complex number.
+ * @param c - The complex number
+ * @returns The principal value of natural logarithm as a complex number
+ * @throws {Error} If the complex number is zero or near-zero
+ */
+export function lnComplex(c) {
+    validateComplexNumber(c);
+    const magnitude = magnitudeComplex(c);
+    if (isApproximatelyZero(magnitude)) {
+        throw new Error('Cannot take logarithm of zero or near-zero complex number');
+    }
+    const angle = argumentComplex(c);
+    return normalizeComplex({
+        real: Math.log(magnitude),
+        imag: angle,
+    });
+}
+/**
+ * Calculates the exponential of a complex number.
+ * @param c - The complex number
+ * @returns The exponential as a complex number
+ */
+export function expComplex(c) {
+    validateComplexNumber(c);
+    const expReal = Math.exp(c.real);
+    return normalizeComplex({
+        real: expReal * Math.cos(c.imag),
+        imag: expReal * Math.sin(c.imag),
+    });
+}
+/**
+ * Calculates the sine of a complex number.
+ * @param c - The complex number
+ * @returns The sine as a complex number
+ */
+export function sinComplex(c) {
+    validateComplexNumber(c);
+    // sin(a + bi) = sin(a)cosh(b) + i cos(a)sinh(b)
+    return normalizeComplex({
+        real: Math.sin(c.real) * Math.cosh(c.imag),
+        imag: Math.cos(c.real) * Math.sinh(c.imag),
+    });
+}
+/**
+ * Calculates the cosine of a complex number.
+ * @param c - The complex number
+ * @returns The cosine as a complex number
+ */
+export function cosComplex(c) {
+    validateComplexNumber(c);
+    // cos(a + bi) = cos(a)cosh(b) - i sin(a)sinh(b)
+    return normalizeComplex({
+        real: Math.cos(c.real) * Math.cosh(c.imag),
+        imag: -Math.sin(c.real) * Math.sinh(c.imag),
+    });
+}
+/**
+ * Checks if a complex number is zero (or effectively zero).
+ * @param c - The complex number to check
+ * @returns True if the complex number is effectively zero
+ */
+export function isZeroComplex(c) {
+    validateComplexNumber(c);
+    return isApproximatelyZero(c.real) && isApproximatelyZero(c.imag);
+}
+/**
+ * Returns the real part of a complex number.
+ * @param c - The complex number
+ * @returns The real part
+ */
+export function getReal(c) {
+    validateComplexNumber(c);
+    return c.real;
+}
+/**
+ * Returns the imaginary part of a complex number.
+ * @param c - The complex number
+ * @returns The imaginary part
+ */
+export function getImag(c) {
+    validateComplexNumber(c);
+    return c.imag;
+}
+/**
+ * Negates a complex number.
+ * @param c - The complex number
+ * @returns The negated complex number
+ */
+export function negateComplex(c) {
+    validateComplexNumber(c);
+    return normalizeComplex({
+        real: -c.real,
+        imag: -c.imag,
+    });
+}
+/**
+ * Calculates the reciprocal of a complex number.
+ * @param c - The complex number
+ * @returns The reciprocal (1/c)
+ * @throws {Error} If the complex number is zero or near-zero
+ */
+export function reciprocalComplex(c) {
+    validateComplexNumber(c);
+    const denominator = c.real * c.real + c.imag * c.imag;
+    if (isApproximatelyZero(denominator)) {
+        throw new Error('Cannot calculate reciprocal of zero or near-zero complex number');
+    }
+    return normalizeComplex({
+        real: c.real / denominator,
+        imag: -c.imag / denominator,
+    });
+}

package/dist/lib/maths/fn/complex.mjs

@@ -0,0 +1 @@
+var g=2*Math.PI;function n(e){return Math.abs(e)<1e-10}function p(e,r){return Math.abs(e-r)<1e-10}function t(e){return{real:n(e.real)?0:e.real,imag:n(e.imag)?0:e.imag}}function l(e,r){if(typeof e!="number"||!Number.isFinite(e))throw new Error(`${r} must be a finite number, got ${e}`)}function m(e,r="Complex number"){if(!e||typeof e!="object")throw new Error(`${r} must be an object`);l(e.real,`${r}.real`),l(e.imag,`${r}.imag`)}function i(e,r=0){return l(e,"real"),l(r,"imag"),t({real:e,imag:r})}function M(e,r){return m(e,"First complex number"),m(r,"Second complex number"),p(e.real,r.real)&&p(e.imag,r.imag)}function c(e,r){return m(e,"First complex number"),m(r,"Second complex number"),t({real:e.real+r.real,imag:e.imag+r.imag})}function d(e,r){return m(e,"First complex number"),m(r,"Second complex number"),t({real:e.real-r.real,imag:e.imag-r.imag})}function f(e,r){return m(e,"First complex number"),m(r,"Second complex number"),t({real:e.real*r.real-e.imag*r.imag,imag:e.real*r.imag+e.imag*r.real})}function w(e,r){m(e,"Numerator complex number"),m(r,"Denominator complex number");let o=r.real*r.real+r.imag*r.imag;if(n(o))throw new Error("Cannot divide by zero or near-zero complex number");return t({real:(e.real*r.real+e.imag*r.imag)/o,imag:(e.imag*r.real-e.real*r.imag)/o})}function C(e){m(e);let r=Math.sqrt(e.real*e.real+e.imag*e.imag);return n(r)?0:r}function E(e){return m(e),t({real:e.real,imag:-e.imag})}function s(e){return m(e),n(e.real)&&n(e.imag)?0:Math.atan2(e.imag,e.real)}function h(e){return m(e),{magnitude:C(e),angle:s(e)}}function b(e,r){if(l(e,"magnitude"),l(r,"angle"),e<-1e-10)throw new Error("Magnitude must be non-negative");let o=r%g;o>Math.PI&&(o-=g),o<=-Math.PI&&(o+=g);let a=Math.max(0,e);return t({real:a*Math.cos(o),imag:a*Math.sin(o)})}function P(e){if(m(e),n(e.real)&&n(e.imag)){let N=i(0,0);return[N,N]}let r=C(e),o=s(e),a=Math.sqrt(r),u=o/2,x=u+Math.PI;return[b(a,u),b(a,x)]}function I(e,r){if(m(e),l(r,"exponent"),n(r))return i(1,0);if(p(r,1))return i(e.real,e.imag);if(p(r,2))return f(e,e);if(n(e.real)&&n(e.imag)){if(r>0)return i(0,0);if(r<0)throw new Error("Cannot raise zero to a negative power")}let o=h(e),a=Math.pow(o.magnitude,r),u=a<-1e-10?Math.abs(a):Math.max(0,a),x=o.angle*r;return b(u,x)}function z(e){m(e);let r=C(e);if(n(r))throw new Error("Cannot take logarithm of zero or near-zero complex number");let o=s(e);return t({real:Math.log(r),imag:o})}function S(e){m(e);let r=Math.exp(e.real);return t({real:r*Math.cos(e.imag),imag:r*Math.sin(e.imag)})}function y(e){return m(e),t({real:Math.sin(e.real)*Math.cosh(e.imag),imag:Math.cos(e.real)*Math.sinh(e.imag)})}function v(e){return m(e),t({real:Math.cos(e.real)*Math.cosh(e.imag),imag:-Math.sin(e.real)*Math.sinh(e.imag)})}function O(e){return m(e),n(e.real)&&n(e.imag)}function q(e){return m(e),e.real}function F(e){return m(e),e.imag}function L(e){return m(e),t({real:-e.real,imag:-e.imag})}function $(e){m(e);let r=e.real*e.real+e.imag*e.imag;if(n(r))throw new Error("Cannot calculate reciprocal of zero or near-zero complex number");return t({real:e.real/r,imag:-e.imag/r})}export{c as addComplex,s as argumentComplex,M as complexEquals,E as conjugateComplex,v as cosComplex,i as createComplex,w as divideComplex,S as expComplex,b as fromPolarComplex,F as getImag,q as getReal,O as isZeroComplex,z as lnComplex,C as magnitudeComplex,f as multiplyComplex,L as negateComplex,I as powerComplex,$ as reciprocalComplex,y as sinComplex,P as sqrtComplex,d as subtractComplex,h as toPolarComplex};

package/dist/lib/maths/fn/index.cjs

@@ -0,0 +1 @@
+"use strict";var S=Object.defineProperty;var ae=Object.getOwnPropertyDescriptor;var se=Object.getOwnPropertyNames;var xe=Object.prototype.hasOwnProperty;var pe=(r,e)=>{for(var n in e)S(r,n,{get:e[n],enumerable:!0})},le=(r,e,n,t)=>{if(e&&typeof e=="object"||typeof e=="function")for(let o of se(e))!xe.call(r,o)&&o!==n&&S(r,o,{get:()=>e[o],enumerable:!(t=ae(e,o))||t.enumerable});return r};var fe=r=>le(S({},"__esModule",{value:!0}),r);var Fe={};pe(Fe,{abs:()=>O,acos:()=>Jr,acosh:()=>oe,add:()=>_,addComplex:()=>Pr,addFractions:()=>xr,addVectors2D:()=>Cr,angleBetween2D:()=>Sr,argumentComplex:()=>w,asin:()=>Wr,asinh:()=>te,atan:()=>Kr,atan2:()=>vr,atanh:()=>ue,average:()=>K,cbrt:()=>cr,ceil:()=>k,clamp:()=>M,conjugateComplex:()=>Lr,cos:()=>Hr,cosh:()=>ee,cot:()=>me,crossProduct3D:()=>Mr,csc:()=>ce,degToGrad:()=>kr,degToRad:()=>Br,distance2D:()=>Rr,distance3D:()=>Fr,divide:()=>$,divideComplex:()=>_r,dotProduct2D:()=>T,dotProduct3D:()=>zr,evaluatePolynomial:()=>sr,exp:()=>ir,factorial:()=>Y,floor:()=>Q,fromPolarComplex:()=>g,gcd:()=>h,gradToDeg:()=>Qr,invLerp:()=>F,lcm:()=>H,lerp:()=>R,lerp2D:()=>Ir,ln:()=>nr,lnComplex:()=>Gr,log:()=>ur,log10:()=>tr,log2:()=>or,magnitude2D:()=>y,magnitude3D:()=>I,magnitude4D:()=>gr,magnitudeComplex:()=>V,mapRange:()=>U,max:()=>J,median:()=>rr,min:()=>W,modulo:()=>j,multiplyComplex:()=>P,multiplyFractions:()=>pr,normalize:()=>hr,normalize2D:()=>Vr,normalize3D:()=>wr,normalizeDegrees:()=>Xr,normalizeRadians:()=>Zr,nthRoot:()=>mr,percentage:()=>lr,power:()=>G,powerComplex:()=>jr,product:()=>er,radToDeg:()=>Or,range:()=>br,rotate2D:()=>Tr,round:()=>b,scaleVector2D:()=>Nr,sec:()=>ie,sign:()=>X,simplifyFraction:()=>C,sin:()=>Yr,sinh:()=>re,solveQuadratic:()=>ar,sqrt:()=>B,sqrtComplex:()=>$r,subtract:()=>L,subtractComplex:()=>qr,subtractVectors2D:()=>Er,sum:()=>v,tan:()=>Ur,tanh:()=>ne,toBase:()=>fr,toOrdinal:()=>Dr,toPolarComplex:()=>q,trunc:()=>Z,truncateDecimals:()=>dr,wrap:()=>yr});module.exports=fe(Fe);function _(r,e){return r+e}function L(r,e){return r-e}var $=(r,e,n)=>e===0?n!==void 0?n:NaN:r/e;function j(r,e){return r%e}function G(r,e){return Math.pow(r,e)}function B(r){if(r<0)throw new Error("Cannot calculate square root of negative number");return Math.sqrt(r)}function O(r){return Math.abs(r)}var b=(r,e=0,n="round")=>{if(n==="david")return De(r,e,n);if(!Number.isFinite(r))return r;if(e===0)return Math[n](r);if(e<0){let z=10**-e;return Math[n](r/z)*z}let t=10**e,o=Math.pow(2,-52),u=r*t,s=Math.sign(u)*o*Math.abs(u),l;switch(n){case"round":l=Math.round(u+s);break;case"ceil":l=Math.ceil(u-o);break;case"floor":l=Math.floor(u+o);break;case"trunc":l=u<0?Math.ceil(u):Math.floor(u);break;default:l=Math.round(u+s)}return l/t},De=(r,e=0,n="round")=>{let t=Math.abs(r),o=r<0?-1:1;if(e===0)return Math.round(t)*o;if(!Number.isFinite(t))return t*o;let u=10**e,s=t*u,l=Math.pow(2,-53),z=Math.sign(s)*l*Math.abs(s);return s+=z,Math.round(s)/u*o},k=(r,e=0)=>b(r,e,"ceil"),Q=(r,e=0)=>b(r,e,"floor"),Z=(r,e=0)=>b(r,e,"trunc"),M=(r,e,n)=>(e>n&&([e,n]=[n,e]),r<e?e:r>n?n:r);function X(r){return r>0?1:r<0?-1:0}function Y(r){if(!Number.isInteger(r)||r<0)throw new Error("Factorial is only defined for non-negative integers");if(r===0||r===1)return 1;let e=1;for(let n=2;n<=r;n++)e*=n;return e}function h(r,e){for(r=Math.abs(r),e=Math.abs(e);e!==0;){let n=e;e=r%e,r=n}return r}function H(r,e){return Math.abs(r*e)/h(r,e)}function R(r,e,n){return r+(e-r)*M(n,0,1)}var F=(r,e,n)=>M((n-r)/(e-r),0,1),U=(r,e,n,t,o)=>{let u=F(e,n,r);return R(t,o,u)};function W(...r){return Math.min(...r)}function J(...r){return Math.max(...r)}var K=(...r)=>{if(!Array.isArray(r)||r.length===0)return 0;let e=0,n=0;for(let t of r)Number.isFinite(t)&&(e+=t,n++);return n===0?0:e/n},v=(...r)=>Array.isArray(r)?r.reduce((e,n)=>e+(Number.isFinite(n)?n:0),0):0,rr=(...r)=>{if(!Array.isArray(r)||r.length===0)return 0;let e=r.filter(Number.isFinite);if(e.length===0)return 0;e.sort((t,o)=>t-o);let n=Math.floor(e.length/2);return e.length%2===0?(e[n-1]+e[n])/2:e[n]};function er(...r){return r.reduce((e,n)=>e*n,1)}function nr(r){if(r<=0)throw new Error("Natural logarithm is only defined for positive numbers");return Math.log(r)}function tr(r){if(r<=0)throw new Error("Logarithm is only defined for positive numbers");return Math.log10(r)}function or(r){if(r<=0)throw new Error("Logarithm is only defined for positive numbers");return Math.log2(r)}function ur(r,e){if(r<=0||e<=0)throw new Error("Logarithm is only defined for positive numbers and bases");return Math.log(r)/Math.log(e)}function ir(r){return Math.exp(r)}function cr(r){return Math.cbrt(r)}function mr(r,e){if(e===0)throw new Error("Cannot calculate 0th root");if(r<0&&e%2===0)throw new Error("Cannot calculate even root of negative number");let n=1/e;return Math.sign(r)*Math.pow(Math.abs(r),n)}function ar(r,e,n){if(r===0)return e===0?[]:[-n/e];let t=e*e-4*r*n;if(t<0)return[];if(t===0)return[-e/(2*r)];{let o=Math.sqrt(t);return[(-e-o)/(2*r),(-e+o)/(2*r)]}}function sr(r,e){let n=0;for(let t=0;t<r.length;t++)n=n*e+r[t];return n}function C(r){let e=h(r.numerator,r.denominator);return{numerator:r.numerator/e,denominator:r.denominator/e}}function xr(r,e){let n=r.numerator*e.denominator+e.numerator*r.denominator,t=r.denominator*e.denominator;return C({numerator:n,denominator:t})}function pr(r,e){let n=r.numerator*e.numerator,t=r.denominator*e.denominator;return C({numerator:n,denominator:t})}var lr=(r,e,n=2)=>{if(e===0)return 0;let t=r/e*100;return b(t,n)},fr=(r,e,n=1)=>{if(e<2||e>36)throw new Error("Base must be between 2 and 36");return r.toString(e).padStart(n,"0")},Dr=r=>{let e=["th","st","nd","rd"],n=r%100;return r+(e[(n-20)%10]||e[n]||e[0])},br=(r,e,n=1)=>{let t=[];if(n===0)return t;if(n>0)for(let o=r;o<e;o+=n)t.push(o);else for(let o=r;o>e;o+=n)t.push(o);return t},dr=(r,e)=>{if(!Number.isFinite(r))return r;let n=10**e;return Math.trunc(r*n)/n},hr=(r,e,n)=>e===n?0:(r-e)/(n-e),yr=(r,e,n)=>{let t=n-e;return((r-e)%t+t)%t+e};function m(r,e){if(typeof r!="number"||!Number.isFinite(r))throw new Error(`${e} must be a finite number, got ${r}`)}function p(r){return Math.abs(r)<1e-9}function i(r,e="Vector2D"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);m(r.x,`${e}.x`),m(r.y,`${e}.y`)}function f(r,e="Vector3D"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);m(r.x,`${e}.x`),m(r.y,`${e}.y`),m(r.z,`${e}.z`)}function be(r,e="Vector4D"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);m(r.x,`${e}.x`),m(r.y,`${e}.y`),m(r.z,`${e}.z`),m(r.w,`${e}.w`)}function y(r){i(r);let e=Math.sqrt(r.x*r.x+r.y*r.y);return p(e)?0:e}function I(r){f(r);let e=Math.sqrt(r.x*r.x+r.y*r.y+r.z*r.z);return p(e)?0:e}function gr(r){be(r);let e=Math.sqrt(r.x*r.x+r.y*r.y+r.z*r.z+r.w*r.w);return p(e)?0:e}function Vr(r){i(r);let e=y(r);if(p(e))throw new Error("Cannot normalize zero or near-zero vector");return{x:r.x/e,y:r.y/e}}function wr(r){f(r);let e=I(r);if(p(e))throw new Error("Cannot normalize zero or near-zero vector");return{x:r.x/e,y:r.y/e,z:r.z/e}}function T(r,e){return i(r,"First vector"),i(e,"Second vector"),r.x*e.x+r.y*e.y}function zr(r,e){return f(r,"First vector"),f(e,"Second vector"),r.x*e.x+r.y*e.y+r.z*e.z}function Mr(r,e){return f(r,"First vector"),f(e,"Second vector"),{x:r.y*e.z-r.z*e.y,y:r.z*e.x-r.x*e.z,z:r.x*e.y-r.y*e.x}}function Cr(r,e){return i(r,"First vector"),i(e,"Second vector"),{x:r.x+e.x,y:r.y+e.y}}function Er(r,e){return i(r,"First vector"),i(e,"Second vector"),{x:r.x-e.x,y:r.y-e.y}}function Nr(r,e){return i(r),m(e,"scalar"),{x:r.x*e,y:r.y*e}}function Sr(r,e){i(r,"First vector"),i(e,"Second vector");let n=T(r,e),t=y(r),o=y(e);if(p(t)||p(o))throw new Error("Cannot calculate angle with zero or near-zero vector");let u=Math.max(-1,Math.min(1,n/(t*o)));return Math.acos(u)}function Rr(r,e){i(r,"First point"),i(e,"Second point");let n=e.x-r.x,t=e.y-r.y,o=Math.sqrt(n*n+t*t);return p(o)?0:o}function Fr(r,e){f(r,"First point"),f(e,"Second point");let n=e.x-r.x,t=e.y-r.y,o=e.z-r.z,u=Math.sqrt(n*n+t*t+o*o);return p(u)?0:u}function Ir(r,e,n){return i(r,"Start vector"),i(e,"End vector"),m(n,"t"),{x:r.x+(e.x-r.x)*n,y:r.y+(e.y-r.y)*n}}function Tr(r,e){i(r),m(e,"angle");let n=Math.cos(e),t=Math.sin(e);return{x:r.x*n-r.y*t,y:r.x*t+r.y*n}}var A=2*Math.PI;function a(r){return Math.abs(r)<1e-10}function Ar(r,e){return Math.abs(r-e)<1e-10}function D(r){return{real:a(r.real)?0:r.real,imag:a(r.imag)?0:r.imag}}function d(r,e){if(typeof r!="number"||!Number.isFinite(r))throw new Error(`${e} must be a finite number, got ${r}`)}function c(r,e="Complex number"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);d(r.real,`${e}.real`),d(r.imag,`${e}.imag`)}function E(r,e=0){return d(r,"real"),d(e,"imag"),D({real:r,imag:e})}function Pr(r,e){return c(r,"First complex number"),c(e,"Second complex number"),D({real:r.real+e.real,imag:r.imag+e.imag})}function qr(r,e){return c(r,"First complex number"),c(e,"Second complex number"),D({real:r.real-e.real,imag:r.imag-e.imag})}function P(r,e){return c(r,"First complex number"),c(e,"Second complex number"),D({real:r.real*e.real-r.imag*e.imag,imag:r.real*e.imag+r.imag*e.real})}function _r(r,e){c(r,"Numerator complex number"),c(e,"Denominator complex number");let n=e.real*e.real+e.imag*e.imag;if(a(n))throw new Error("Cannot divide by zero or near-zero complex number");return D({real:(r.real*e.real+r.imag*e.imag)/n,imag:(r.imag*e.real-r.real*e.imag)/n})}function V(r){c(r);let e=Math.sqrt(r.real*r.real+r.imag*r.imag);return a(e)?0:e}function Lr(r){return c(r),D({real:r.real,imag:-r.imag})}function w(r){return c(r),a(r.real)&&a(r.imag)?0:Math.atan2(r.imag,r.real)}function q(r){return c(r),{magnitude:V(r),angle:w(r)}}function g(r,e){if(d(r,"magnitude"),d(e,"angle"),r<-1e-10)throw new Error("Magnitude must be non-negative");let n=e%A;n>Math.PI&&(n-=A),n<=-Math.PI&&(n+=A);let t=Math.max(0,r);return D({real:t*Math.cos(n),imag:t*Math.sin(n)})}function $r(r){if(c(r),a(r.real)&&a(r.imag)){let s=E(0,0);return[s,s]}let e=V(r),n=w(r),t=Math.sqrt(e),o=n/2,u=o+Math.PI;return[g(t,o),g(t,u)]}function jr(r,e){if(c(r),d(e,"exponent"),a(e))return E(1,0);if(Ar(e,1))return E(r.real,r.imag);if(Ar(e,2))return P(r,r);if(a(r.real)&&a(r.imag)){if(e>0)return E(0,0);if(e<0)throw new Error("Cannot raise zero to a negative power")}let n=q(r),t=Math.pow(n.magnitude,e),o=t<-1e-10?Math.abs(t):Math.max(0,t),u=n.angle*e;return g(o,u)}function Gr(r){c(r);let e=V(r);if(a(e))throw new Error("Cannot take logarithm of zero or near-zero complex number");let n=w(r);return D({real:Math.log(e),imag:n})}var N=3.141592653589793,de=2.718281828459045,he=1.618033988749894,ye=1.4142135623730951,ge=.7071067811865476,Ve=.6931471805599453,we=2.302585092994046,ze=1.4426950408889634,Me=.4342944819032518,Ce=360,Ee=2*N,Ne=400,Se=180/N,Re=N/180,x={PI:N,E:de,PHI:he,SQRT2:ye,SQRT1_2:ge,LN2:Ve,LN10:we,LOG2E:ze,LOG10E:Me,DEGREES_IN_CIRCLE:Ce,RADIANS_IN_CIRCLE:Ee,GRADIANS_IN_CIRCLE:Ne,DEGREES_PER_RADIAN:Se,RADIANS_PER_DEGREE:Re};function Br(r){return r*x.RADIANS_PER_DEGREE}function Or(r){return r*x.DEGREES_PER_RADIAN}function kr(r){return r*(x.GRADIANS_IN_CIRCLE/x.DEGREES_IN_CIRCLE)}function Qr(r){return r*(x.DEGREES_IN_CIRCLE/x.GRADIANS_IN_CIRCLE)}function Zr(r){let e=x.RADIANS_IN_CIRCLE,n=r%e;return n<0&&(n+=e),n}function Xr(r){let e=r%x.DEGREES_IN_CIRCLE;return e<0&&(e+=x.DEGREES_IN_CIRCLE),e}function Yr(r){return Math.sin(r)}function Hr(r){return Math.cos(r)}function Ur(r){let e=Math.cos(r);if(Math.abs(e)<1e-10)throw new Error("Tangent undefined for this angle (cosine is zero)");return Math.tan(r)}function Wr(r){if(r<-1||r>1)throw new Error("Arcsine is only defined for values between -1 and 1");return Math.asin(r)}function Jr(r){if(r<-1||r>1)throw new Error("Arccosine is only defined for values between -1 and 1");return Math.acos(r)}function Kr(r){return Math.atan(r)}function vr(r,e){return Math.atan2(r,e)}function re(r){return Math.sinh(r)}function ee(r){return Math.cosh(r)}function ne(r){return Math.tanh(r)}function te(r){return Math.asinh(r)}function oe(r){if(r<1)throw new Error("Inverse hyperbolic cosine is only defined for x >= 1");return Math.acosh(r)}function ue(r){if(Math.abs(r)>=1)throw new Error("Inverse hyperbolic tangent is only defined for |x| < 1");return Math.atanh(r)}function ie(r){let e=Math.cos(r);if(Math.abs(e)<1e-10)throw new Error("Secant undefined for this angle (cosine is zero)");return 1/e}function ce(r){let e=Math.sin(r);if(Math.abs(e)<1e-10)throw new Error("Cosecant undefined for this angle (sine is zero)");return 1/e}function me(r){let e=Math.tan(r);if(Math.abs(e)<1e-10)throw new Error("Cotangent undefined for this angle (tangent is zero)");return 1/e}

package/dist/lib/maths/fn/index.d.ts

@@ -0,0 +1,10 @@
+/**
+ * Re-export all mathematical functions.
+ * @module MathFunctions
+ */
+export { add, subtract, divide, modulo, power, sqrt, abs, round, floor, ceil, trunc, clamp, sign, factorial, gcd, lcm, lerp, min, max, average, sum, product, invLerp, mapRange, median } from './basic.js';
+export { ln, log10, log2, log, exp, cbrt, nthRoot, solveQuadratic, evaluatePolynomial, simplifyFraction, addFractions, multiplyFractions, percentage, normalize, range, toBase, toOrdinal, truncateDecimals, wrap } from './advanced.js';
+export { magnitude2D, magnitude3D, magnitude4D, normalize2D, normalize3D, dotProduct2D, dotProduct3D, crossProduct3D, addVectors2D, subtractVectors2D, scaleVector2D, angleBetween2D, distance2D, distance3D, lerp2D, rotate2D, } from './vector.js';
+export { addComplex, subtractComplex, multiplyComplex, divideComplex, magnitudeComplex, conjugateComplex, argumentComplex, toPolarComplex, fromPolarComplex, sqrtComplex, powerComplex, lnComplex, } from './complex.js';
+export { degToRad, radToDeg, degToGrad, gradToDeg, normalizeRadians, normalizeDegrees, sin, cos, tan, asin, acos, atan, atan2, sinh, cosh, tanh, asinh, acosh, atanh, sec, csc, cot, } from './trigo.js';
+//# sourceMappingURL=index.d.ts.map

package/dist/lib/maths/fn/index.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"index.d.ts","sourceRoot":"","sources":["../../../../src/lib/maths/fn/index.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAGH,OAAO,EACL,GAAG,EACH,QAAQ,EACR,MAAM,EACN,MAAM,EACN,KAAK,EACL,IAAI,EACJ,GAAG,EACH,KAAK,EACL,KAAK,EACL,IAAI,EACJ,KAAK,EACL,KAAK,EACL,IAAI,EACJ,SAAS,EACT,GAAG,EACH,GAAG,EACH,IAAI,EACJ,GAAG,EACH,GAAG,EACH,OAAO,EACP,GAAG,EACH,OAAO,EACP,OAAO,EACP,QAAQ,EACR,MAAM,EACP,MAAM,YAAY,CAAC;AAGpB,OAAO,EACL,EAAE,EACF,KAAK,EACL,IAAI,EACJ,GAAG,EACH,GAAG,EACH,IAAI,EACJ,OAAO,EACP,cAAc,EACd,kBAAkB,EAClB,gBAAgB,EAChB,YAAY,EACZ,iBAAiB,EACjB,UAAU,EACV,SAAS,EACT,KAAK,EACL,MAAM,EACN,SAAS,EACT,gBAAgB,EAChB,IAAI,EACL,MAAM,eAAe,CAAC;AAGvB,OAAO,EACL,WAAW,EACX,WAAW,EACX,WAAW,EACX,WAAW,EACX,WAAW,EACX,YAAY,EACZ,YAAY,EACZ,cAAc,EACd,YAAY,EACZ,iBAAiB,EACjB,aAAa,EACb,cAAc,EACd,UAAU,EACV,UAAU,EACV,MAAM,EACN,QAAQ,GACT,MAAM,aAAa,CAAC;AAGrB,OAAO,EACL,UAAU,EACV,eAAe,EACf,eAAe,EACf,aAAa,EACb,gBAAgB,EAChB,gBAAgB,EAChB,eAAe,EACf,cAAc,EACd,gBAAgB,EAChB,WAAW,EACX,YAAY,EACZ,SAAS,GACV,MAAM,cAAc,CAAC;AAGtB,OAAO,EACL,QAAQ,EACR,QAAQ,EACR,SAAS,EACT,SAAS,EACT,gBAAgB,EAChB,gBAAgB,EAChB,GAAG,EACH,GAAG,EACH,GAAG,EACH,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,KAAK,EACL,IAAI,EACJ,IAAI,EACJ,IAAI,EACJ,KAAK,EACL,KAAK,EACL,KAAK,EACL,GAAG,EACH,GAAG,EACH,GAAG,GACJ,MAAM,YAAY,CAAC"}

package/dist/lib/maths/fn/index.js

@@ -0,0 +1,9 @@
+/**
+ * Re-export all mathematical functions.
+ * @module MathFunctions
+ */
+export { add, subtract, divide, modulo, power, sqrt, abs, round, floor, ceil, trunc, clamp, sign, factorial, gcd, lcm, lerp, min, max, average, sum, product, invLerp, mapRange, median } from './basic.js';
+export { ln, log10, log2, log, exp, cbrt, nthRoot, solveQuadratic, evaluatePolynomial, simplifyFraction, addFractions, multiplyFractions, percentage, normalize, range, toBase, toOrdinal, truncateDecimals, wrap } from './advanced.js';
+export { magnitude2D, magnitude3D, magnitude4D, normalize2D, normalize3D, dotProduct2D, dotProduct3D, crossProduct3D, addVectors2D, subtractVectors2D, scaleVector2D, angleBetween2D, distance2D, distance3D, lerp2D, rotate2D, } from './vector.js';
+export { addComplex, subtractComplex, multiplyComplex, divideComplex, magnitudeComplex, conjugateComplex, argumentComplex, toPolarComplex, fromPolarComplex, sqrtComplex, powerComplex, lnComplex, } from './complex.js';
+export { degToRad, radToDeg, degToGrad, gradToDeg, normalizeRadians, normalizeDegrees, sin, cos, tan, asin, acos, atan, atan2, sinh, cosh, tanh, asinh, acosh, atanh, sec, csc, cot, } from './trigo.js';

package/dist/lib/maths/fn/index.mjs

@@ -0,0 +1 @@
+function _(r,e){return r+e}function L(r,e){return r-e}var $=(r,e,n)=>e===0?n!==void 0?n:NaN:r/e;function j(r,e){return r%e}function G(r,e){return Math.pow(r,e)}function B(r){if(r<0)throw new Error("Cannot calculate square root of negative number");return Math.sqrt(r)}function O(r){return Math.abs(r)}var d=(r,e=0,n="round")=>{if(n==="david")return k(r,e,n);if(!Number.isFinite(r))return r;if(e===0)return Math[n](r);if(e<0){let h=10**-e;return Math[n](r/h)*h}let t=10**e,o=Math.pow(2,-52),u=r*t,s=Math.sign(u)*o*Math.abs(u),l;switch(n){case"round":l=Math.round(u+s);break;case"ceil":l=Math.ceil(u-o);break;case"floor":l=Math.floor(u+o);break;case"trunc":l=u<0?Math.ceil(u):Math.floor(u);break;default:l=Math.round(u+s)}return l/t},k=(r,e=0,n="round")=>{let t=Math.abs(r),o=r<0?-1:1;if(e===0)return Math.round(t)*o;if(!Number.isFinite(t))return t*o;let u=10**e,s=t*u,l=Math.pow(2,-53),h=Math.sign(s)*l*Math.abs(s);return s+=h,Math.round(s)/u*o},Q=(r,e=0)=>d(r,e,"ceil"),Z=(r,e=0)=>d(r,e,"floor"),X=(r,e=0)=>d(r,e,"trunc"),E=(r,e,n)=>(e>n&&([e,n]=[n,e]),r<e?e:r>n?n:r);function Y(r){return r>0?1:r<0?-1:0}function H(r){if(!Number.isInteger(r)||r<0)throw new Error("Factorial is only defined for non-negative integers");if(r===0||r===1)return 1;let e=1;for(let n=2;n<=r;n++)e*=n;return e}function y(r,e){for(r=Math.abs(r),e=Math.abs(e);e!==0;){let n=e;e=r%e,r=n}return r}function U(r,e){return Math.abs(r*e)/y(r,e)}function R(r,e,n){return r+(e-r)*E(n,0,1)}var F=(r,e,n)=>E((n-r)/(e-r),0,1),W=(r,e,n,t,o)=>{let u=F(e,n,r);return R(t,o,u)};function J(...r){return Math.min(...r)}function K(...r){return Math.max(...r)}var v=(...r)=>{if(!Array.isArray(r)||r.length===0)return 0;let e=0,n=0;for(let t of r)Number.isFinite(t)&&(e+=t,n++);return n===0?0:e/n},rr=(...r)=>Array.isArray(r)?r.reduce((e,n)=>e+(Number.isFinite(n)?n:0),0):0,er=(...r)=>{if(!Array.isArray(r)||r.length===0)return 0;let e=r.filter(Number.isFinite);if(e.length===0)return 0;e.sort((t,o)=>t-o);let n=Math.floor(e.length/2);return e.length%2===0?(e[n-1]+e[n])/2:e[n]};function nr(...r){return r.reduce((e,n)=>e*n,1)}function tr(r){if(r<=0)throw new Error("Natural logarithm is only defined for positive numbers");return Math.log(r)}function or(r){if(r<=0)throw new Error("Logarithm is only defined for positive numbers");return Math.log10(r)}function ur(r){if(r<=0)throw new Error("Logarithm is only defined for positive numbers");return Math.log2(r)}function ir(r,e){if(r<=0||e<=0)throw new Error("Logarithm is only defined for positive numbers and bases");return Math.log(r)/Math.log(e)}function cr(r){return Math.exp(r)}function mr(r){return Math.cbrt(r)}function ar(r,e){if(e===0)throw new Error("Cannot calculate 0th root");if(r<0&&e%2===0)throw new Error("Cannot calculate even root of negative number");let n=1/e;return Math.sign(r)*Math.pow(Math.abs(r),n)}function sr(r,e,n){if(r===0)return e===0?[]:[-n/e];let t=e*e-4*r*n;if(t<0)return[];if(t===0)return[-e/(2*r)];{let o=Math.sqrt(t);return[(-e-o)/(2*r),(-e+o)/(2*r)]}}function xr(r,e){let n=0;for(let t=0;t<r.length;t++)n=n*e+r[t];return n}function N(r){let e=y(r.numerator,r.denominator);return{numerator:r.numerator/e,denominator:r.denominator/e}}function pr(r,e){let n=r.numerator*e.denominator+e.numerator*r.denominator,t=r.denominator*e.denominator;return N({numerator:n,denominator:t})}function lr(r,e){let n=r.numerator*e.numerator,t=r.denominator*e.denominator;return N({numerator:n,denominator:t})}var fr=(r,e,n=2)=>{if(e===0)return 0;let t=r/e*100;return d(t,n)},Dr=(r,e,n=1)=>{if(e<2||e>36)throw new Error("Base must be between 2 and 36");return r.toString(e).padStart(n,"0")},br=r=>{let e=["th","st","nd","rd"],n=r%100;return r+(e[(n-20)%10]||e[n]||e[0])},dr=(r,e,n=1)=>{let t=[];if(n===0)return t;if(n>0)for(let o=r;o<e;o+=n)t.push(o);else for(let o=r;o>e;o+=n)t.push(o);return t},hr=(r,e)=>{if(!Number.isFinite(r))return r;let n=10**e;return Math.trunc(r*n)/n},yr=(r,e,n)=>e===n?0:(r-e)/(n-e),gr=(r,e,n)=>{let t=n-e;return((r-e)%t+t)%t+e};function m(r,e){if(typeof r!="number"||!Number.isFinite(r))throw new Error(`${e} must be a finite number, got ${r}`)}function p(r){return Math.abs(r)<1e-9}function i(r,e="Vector2D"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);m(r.x,`${e}.x`),m(r.y,`${e}.y`)}function f(r,e="Vector3D"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);m(r.x,`${e}.x`),m(r.y,`${e}.y`),m(r.z,`${e}.z`)}function Vr(r,e="Vector4D"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);m(r.x,`${e}.x`),m(r.y,`${e}.y`),m(r.z,`${e}.z`),m(r.w,`${e}.w`)}function g(r){i(r);let e=Math.sqrt(r.x*r.x+r.y*r.y);return p(e)?0:e}function I(r){f(r);let e=Math.sqrt(r.x*r.x+r.y*r.y+r.z*r.z);return p(e)?0:e}function wr(r){Vr(r);let e=Math.sqrt(r.x*r.x+r.y*r.y+r.z*r.z+r.w*r.w);return p(e)?0:e}function zr(r){i(r);let e=g(r);if(p(e))throw new Error("Cannot normalize zero or near-zero vector");return{x:r.x/e,y:r.y/e}}function Mr(r){f(r);let e=I(r);if(p(e))throw new Error("Cannot normalize zero or near-zero vector");return{x:r.x/e,y:r.y/e,z:r.z/e}}function T(r,e){return i(r,"First vector"),i(e,"Second vector"),r.x*e.x+r.y*e.y}function Cr(r,e){return f(r,"First vector"),f(e,"Second vector"),r.x*e.x+r.y*e.y+r.z*e.z}function Er(r,e){return f(r,"First vector"),f(e,"Second vector"),{x:r.y*e.z-r.z*e.y,y:r.z*e.x-r.x*e.z,z:r.x*e.y-r.y*e.x}}function Nr(r,e){return i(r,"First vector"),i(e,"Second vector"),{x:r.x+e.x,y:r.y+e.y}}function Sr(r,e){return i(r,"First vector"),i(e,"Second vector"),{x:r.x-e.x,y:r.y-e.y}}function Rr(r,e){return i(r),m(e,"scalar"),{x:r.x*e,y:r.y*e}}function Fr(r,e){i(r,"First vector"),i(e,"Second vector");let n=T(r,e),t=g(r),o=g(e);if(p(t)||p(o))throw new Error("Cannot calculate angle with zero or near-zero vector");let u=Math.max(-1,Math.min(1,n/(t*o)));return Math.acos(u)}function Ir(r,e){i(r,"First point"),i(e,"Second point");let n=e.x-r.x,t=e.y-r.y,o=Math.sqrt(n*n+t*t);return p(o)?0:o}function Tr(r,e){f(r,"First point"),f(e,"Second point");let n=e.x-r.x,t=e.y-r.y,o=e.z-r.z,u=Math.sqrt(n*n+t*t+o*o);return p(u)?0:u}function Ar(r,e,n){return i(r,"Start vector"),i(e,"End vector"),m(n,"t"),{x:r.x+(e.x-r.x)*n,y:r.y+(e.y-r.y)*n}}function Pr(r,e){i(r),m(e,"angle");let n=Math.cos(e),t=Math.sin(e);return{x:r.x*n-r.y*t,y:r.x*t+r.y*n}}var S=2*Math.PI;function a(r){return Math.abs(r)<1e-10}function A(r,e){return Math.abs(r-e)<1e-10}function D(r){return{real:a(r.real)?0:r.real,imag:a(r.imag)?0:r.imag}}function b(r,e){if(typeof r!="number"||!Number.isFinite(r))throw new Error(`${e} must be a finite number, got ${r}`)}function c(r,e="Complex number"){if(!r||typeof r!="object")throw new Error(`${e} must be an object`);b(r.real,`${e}.real`),b(r.imag,`${e}.imag`)}function V(r,e=0){return b(r,"real"),b(e,"imag"),D({real:r,imag:e})}function qr(r,e){return c(r,"First complex number"),c(e,"Second complex number"),D({real:r.real+e.real,imag:r.imag+e.imag})}function _r(r,e){return c(r,"First complex number"),c(e,"Second complex number"),D({real:r.real-e.real,imag:r.imag-e.imag})}function P(r,e){return c(r,"First complex number"),c(e,"Second complex number"),D({real:r.real*e.real-r.imag*e.imag,imag:r.real*e.imag+r.imag*e.real})}function Lr(r,e){c(r,"Numerator complex number"),c(e,"Denominator complex number");let n=e.real*e.real+e.imag*e.imag;if(a(n))throw new Error("Cannot divide by zero or near-zero complex number");return D({real:(r.real*e.real+r.imag*e.imag)/n,imag:(r.imag*e.real-r.real*e.imag)/n})}function z(r){c(r);let e=Math.sqrt(r.real*r.real+r.imag*r.imag);return a(e)?0:e}function $r(r){return c(r),D({real:r.real,imag:-r.imag})}function M(r){return c(r),a(r.real)&&a(r.imag)?0:Math.atan2(r.imag,r.real)}function q(r){return c(r),{magnitude:z(r),angle:M(r)}}function w(r,e){if(b(r,"magnitude"),b(e,"angle"),r<-1e-10)throw new Error("Magnitude must be non-negative");let n=e%S;n>Math.PI&&(n-=S),n<=-Math.PI&&(n+=S);let t=Math.max(0,r);return D({real:t*Math.cos(n),imag:t*Math.sin(n)})}function jr(r){if(c(r),a(r.real)&&a(r.imag)){let s=V(0,0);return[s,s]}let e=z(r),n=M(r),t=Math.sqrt(e),o=n/2,u=o+Math.PI;return[w(t,o),w(t,u)]}function Gr(r,e){if(c(r),b(e,"exponent"),a(e))return V(1,0);if(A(e,1))return V(r.real,r.imag);if(A(e,2))return P(r,r);if(a(r.real)&&a(r.imag)){if(e>0)return V(0,0);if(e<0)throw new Error("Cannot raise zero to a negative power")}let n=q(r),t=Math.pow(n.magnitude,e),o=t<-1e-10?Math.abs(t):Math.max(0,t),u=n.angle*e;return w(o,u)}function Br(r){c(r);let e=z(r);if(a(e))throw new Error("Cannot take logarithm of zero or near-zero complex number");let n=M(r);return D({real:Math.log(e),imag:n})}var C=3.141592653589793,Or=2.718281828459045,kr=1.618033988749894,Qr=1.4142135623730951,Zr=.7071067811865476,Xr=.6931471805599453,Yr=2.302585092994046,Hr=1.4426950408889634,Ur=.4342944819032518,Wr=360,Jr=2*C,Kr=400,vr=180/C,re=C/180,x={PI:C,E:Or,PHI:kr,SQRT2:Qr,SQRT1_2:Zr,LN2:Xr,LN10:Yr,LOG2E:Hr,LOG10E:Ur,DEGREES_IN_CIRCLE:Wr,RADIANS_IN_CIRCLE:Jr,GRADIANS_IN_CIRCLE:Kr,DEGREES_PER_RADIAN:vr,RADIANS_PER_DEGREE:re};function ee(r){return r*x.RADIANS_PER_DEGREE}function ne(r){return r*x.DEGREES_PER_RADIAN}function te(r){return r*(x.GRADIANS_IN_CIRCLE/x.DEGREES_IN_CIRCLE)}function oe(r){return r*(x.DEGREES_IN_CIRCLE/x.GRADIANS_IN_CIRCLE)}function ue(r){let e=x.RADIANS_IN_CIRCLE,n=r%e;return n<0&&(n+=e),n}function ie(r){let e=r%x.DEGREES_IN_CIRCLE;return e<0&&(e+=x.DEGREES_IN_CIRCLE),e}function ce(r){return Math.sin(r)}function me(r){return Math.cos(r)}function ae(r){let e=Math.cos(r);if(Math.abs(e)<1e-10)throw new Error("Tangent undefined for this angle (cosine is zero)");return Math.tan(r)}function se(r){if(r<-1||r>1)throw new Error("Arcsine is only defined for values between -1 and 1");return Math.asin(r)}function xe(r){if(r<-1||r>1)throw new Error("Arccosine is only defined for values between -1 and 1");return Math.acos(r)}function pe(r){return Math.atan(r)}function le(r,e){return Math.atan2(r,e)}function fe(r){return Math.sinh(r)}function De(r){return Math.cosh(r)}function be(r){return Math.tanh(r)}function de(r){return Math.asinh(r)}function he(r){if(r<1)throw new Error("Inverse hyperbolic cosine is only defined for x >= 1");return Math.acosh(r)}function ye(r){if(Math.abs(r)>=1)throw new Error("Inverse hyperbolic tangent is only defined for |x| < 1");return Math.atanh(r)}function ge(r){let e=Math.cos(r);if(Math.abs(e)<1e-10)throw new Error("Secant undefined for this angle (cosine is zero)");return 1/e}function Ve(r){let e=Math.sin(r);if(Math.abs(e)<1e-10)throw new Error("Cosecant undefined for this angle (sine is zero)");return 1/e}function we(r){let e=Math.tan(r);if(Math.abs(e)<1e-10)throw new Error("Cotangent undefined for this angle (tangent is zero)");return 1/e}export{O as abs,xe as acos,he as acosh,_ as add,qr as addComplex,pr as addFractions,Nr as addVectors2D,Fr as angleBetween2D,M as argumentComplex,se as asin,de as asinh,pe as atan,le as atan2,ye as atanh,v as average,mr as cbrt,Q as ceil,E as clamp,$r as conjugateComplex,me as cos,De as cosh,we as cot,Er as crossProduct3D,Ve as csc,te as degToGrad,ee as degToRad,Ir as distance2D,Tr as distance3D,$ as divide,Lr as divideComplex,T as dotProduct2D,Cr as dotProduct3D,xr as evaluatePolynomial,cr as exp,H as factorial,Z as floor,w as fromPolarComplex,y as gcd,oe as gradToDeg,F as invLerp,U as lcm,R as lerp,Ar as lerp2D,tr as ln,Br as lnComplex,ir as log,or as log10,ur as log2,g as magnitude2D,I as magnitude3D,wr as magnitude4D,z as magnitudeComplex,W as mapRange,K as max,er as median,J as min,j as modulo,P as multiplyComplex,lr as multiplyFractions,yr as normalize,zr as normalize2D,Mr as normalize3D,ie as normalizeDegrees,ue as normalizeRadians,ar as nthRoot,fr as percentage,G as power,Gr as powerComplex,nr as product,ne as radToDeg,dr as range,Pr as rotate2D,d as round,Rr as scaleVector2D,ge as sec,Y as sign,N as simplifyFraction,ce as sin,fe as sinh,sr as solveQuadratic,B as sqrt,jr as sqrtComplex,L as subtract,_r as subtractComplex,Sr as subtractVectors2D,rr as sum,ae as tan,be as tanh,Dr as toBase,br as toOrdinal,q as toPolarComplex,X as trunc,hr as truncateDecimals,gr as wrap};

package/dist/lib/maths/fn/trigo.cjs

@@ -0,0 +1 @@
+"use strict";var a=Object.defineProperty;var u=Object.getOwnPropertyDescriptor;var c=Object.getOwnPropertyNames;var E=Object.prototype.hasOwnProperty;var f=(n,e)=>{for(var r in e)a(n,r,{get:e[r],enumerable:!0})},b=(n,e,r,i)=>{if(e&&typeof e=="object"||typeof e=="function")for(let o of c(e))!E.call(n,o)&&o!==r&&a(n,o,{get:()=>e[o],enumerable:!(i=u(e,o))||i.enumerable});return n};var h=n=>b(a({},"__esModule",{value:!0}),n);var W={};f(W,{acos:()=>Q,acosh:()=>F,asin:()=>O,asinh:()=>B,atan:()=>j,atan2:()=>H,atanh:()=>J,cos:()=>P,cosh:()=>q,cot:()=>V,csc:()=>U,degToGrad:()=>G,degToRad:()=>w,gradToDeg:()=>l,normalizeDegrees:()=>y,normalizeRadians:()=>L,radToDeg:()=>S,sec:()=>K,sin:()=>z,sinh:()=>k,tan:()=>T,tanh:()=>v});module.exports=h(W);var s=3.141592653589793,R=2.718281828459045,m=1.618033988749894,p=1.4142135623730951,d=.7071067811865476,x=.6931471805599453,I=2.302585092994046,_=1.4426950408889634,D=.4342944819032518,C=360,M=2*s,g=400,N=180/s,A=s/180,t={PI:s,E:R,PHI:m,SQRT2:p,SQRT1_2:d,LN2:x,LN10:I,LOG2E:_,LOG10E:D,DEGREES_IN_CIRCLE:C,RADIANS_IN_CIRCLE:M,GRADIANS_IN_CIRCLE:g,DEGREES_PER_RADIAN:N,RADIANS_PER_DEGREE:A};function w(n){return n*t.RADIANS_PER_DEGREE}function S(n){return n*t.DEGREES_PER_RADIAN}function G(n){return n*(t.GRADIANS_IN_CIRCLE/t.DEGREES_IN_CIRCLE)}function l(n){return n*(t.DEGREES_IN_CIRCLE/t.GRADIANS_IN_CIRCLE)}function L(n){let e=t.RADIANS_IN_CIRCLE,r=n%e;return r<0&&(r+=e),r}function y(n){let e=n%t.DEGREES_IN_CIRCLE;return e<0&&(e+=t.DEGREES_IN_CIRCLE),e}function z(n){return Math.sin(n)}function P(n){return Math.cos(n)}function T(n){let e=Math.cos(n);if(Math.abs(e)<1e-10)throw new Error("Tangent undefined for this angle (cosine is zero)");return Math.tan(n)}function O(n){if(n<-1||n>1)throw new Error("Arcsine is only defined for values between -1 and 1");return Math.asin(n)}function Q(n){if(n<-1||n>1)throw new Error("Arccosine is only defined for values between -1 and 1");return Math.acos(n)}function j(n){return Math.atan(n)}function H(n,e){return Math.atan2(n,e)}function k(n){return Math.sinh(n)}function q(n){return Math.cosh(n)}function v(n){return Math.tanh(n)}function B(n){return Math.asinh(n)}function F(n){if(n<1)throw new Error("Inverse hyperbolic cosine is only defined for x >= 1");return Math.acosh(n)}function J(n){if(Math.abs(n)>=1)throw new Error("Inverse hyperbolic tangent is only defined for |x| < 1");return Math.atanh(n)}function K(n){let e=Math.cos(n);if(Math.abs(e)<1e-10)throw new Error("Secant undefined for this angle (cosine is zero)");return 1/e}function U(n){let e=Math.sin(n);if(Math.abs(e)<1e-10)throw new Error("Cosecant undefined for this angle (sine is zero)");return 1/e}function V(n){let e=Math.tan(n);if(Math.abs(e)<1e-10)throw new Error("Cotangent undefined for this angle (tangent is zero)");return 1/e}

package/dist/lib/maths/fn/trigo.d.ts

@@ -0,0 +1,147 @@
+/**
+ * Trigonometric functions and angle conversions.
+ * @module TrigonometryFunctions
+ */
+import type { Degrees, Radians } from '../types/basic.js';
+/**
+ * Converts degrees to radians.
+ * @param degrees - Angle in degrees
+ * @returns Angle in radians
+ */
+export declare function degToRad(degrees: number): Radians;
+/**
+ * Converts radians to degrees.
+ * @param radians - Angle in radians
+ * @returns Angle in degrees
+ */
+export declare function radToDeg(radians: number): Degrees;
+/**
+ * Converts degrees to gradians.
+ * @param degrees - Angle in degrees
+ * @returns Angle in gradians
+ */
+export declare function degToGrad(degrees: number): number;
+/**
+ * Converts gradians to degrees.
+ * @param gradians - Angle in gradians
+ * @returns Angle in degrees
+ */
+export declare function gradToDeg(gradians: number): Degrees;
+/**
+ * Normalizes an angle to the range [0, 2π) radians.
+ * @param radians - Angle in radians
+ * @returns Normalized angle in radians
+ */
+export declare function normalizeRadians(radians: number): Radians;
+/**
+ * Normalizes an angle to the range [0, 360) degrees.
+ * @param degrees - Angle in degrees
+ * @returns Normalized angle in degrees
+ */
+export declare function normalizeDegrees(degrees: number): Degrees;
+/**
+ * Calculates the sine of an angle in radians.
+ * @param radians - Angle in radians
+ * @returns The sine value
+ */
+export declare function sin(radians: number): number;
+/**
+ * Calculates the cosine of an angle in radians.
+ * @param radians - Angle in radians
+ * @returns The cosine value
+ */
+export declare function cos(radians: number): number;
+/**
+ * Calculates the tangent of an angle in radians.
+ * @param radians - Angle in radians
+ * @returns The tangent value
+ * @throws {Error} If cosine is zero (angle is π/2 + nπ)
+ */
+export declare function tan(radians: number): number;
+/**
+ * Calculates the arcsine (inverse sine) in radians.
+ * @param value - Sine value between -1 and 1
+ * @returns Angle in radians
+ * @throws {Error} If value is outside [-1, 1]
+ */
+export declare function asin(value: number): Radians;
+/**
+ * Calculates the arccosine (inverse cosine) in radians.
+ * @param value - Cosine value between -1 and 1
+ * @returns Angle in radians
+ * @throws {Error} If value is outside [-1, 1]
+ */
+export declare function acos(value: number): Radians;
+/**
+ * Calculates the arctangent (inverse tangent) in radians.
+ * @param value - Tangent value
+ * @returns Angle in radians between -π/2 and π/2
+ */
+export declare function atan(value: number): Radians;
+/**
+ * Calculates the arctangent of y/x in radians, handling all quadrants.
+ * @param y - Numerator
+ * @param x - Denominator
+ * @returns Angle in radians between -π and π
+ */
+export declare function atan2(y: number, x: number): Radians;
+/**
+ * Calculates the hyperbolic sine.
+ * @param x - Input value
+ * @returns Hyperbolic sine
+ */
+export declare function sinh(x: number): number;
+/**
+ * Calculates the hyperbolic cosine.
+ * @param x - Input value
+ * @returns Hyperbolic cosine
+ */
+export declare function cosh(x: number): number;
+/**
+ * Calculates the hyperbolic tangent.
+ * @param x - Input value
+ * @returns Hyperbolic tangent
+ */
+export declare function tanh(x: number): number;
+/**
+ * Calculates the inverse hyperbolic sine.
+ * @param x - Input value
+ * @returns Inverse hyperbolic sine
+ */
+export declare function asinh(x: number): number;
+/**
+ * Calculates the inverse hyperbolic cosine.
+ * @param x - Input value >= 1
+ * @returns Inverse hyperbolic cosine
+ * @throws {Error} If x < 1
+ */
+export declare function acosh(x: number): number;
+/**
+ * Calculates the inverse hyperbolic tangent.
+ * @param x - Input value between -1 and 1 (exclusive)
+ * @returns Inverse hyperbolic tangent
+ * @throws {Error} If |x| >= 1
+ */
+export declare function atanh(x: number): number;
+/**
+ * Calculates the secant of an angle.
+ * @param radians - Angle in radians
+ * @returns The secant value
+ * @throws {Error} If cosine is zero
+ */
+export declare function sec(radians: number): number;
+/**
+ * Calculates the cosecant of an angle.
+ * @param radians - Angle in radians
+ * @returns The cosecant value
+ * @throws {Error} If sine is zero
+ */
+export declare function csc(radians: number): number;
+/**
+ * Calculates the cotangent of an angle.
+ * @param radians - Angle in radians
+ * @returns The cotangent value
+ * @throws {Error} If tangent is zero
+ */
+export declare function cot(radians: number): number;
+//# sourceMappingURL=trigo.d.ts.map

package/dist/lib/maths/fn/trigo.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"trigo.d.ts","sourceRoot":"","sources":["../../../../src/lib/maths/fn/trigo.ts"],"names":[],"mappings":"AAAA;;;GAGG;AAGH,OAAO,KAAK,EAAE,OAAO,EAAE,OAAO,EAAE,MAAM,mBAAmB,CAAC;AAE1D;;;;GAIG;AACH,wBAAgB,QAAQ,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO,CAGjD;AAED;;;;GAIG;AACH,wBAAgB,QAAQ,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO,CAGjD;AAED;;;;GAIG;AACH,wBAAgB,SAAS,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAEjD;AAED;;;;GAIG;AACH,wBAAgB,SAAS,CAAC,QAAQ,EAAE,MAAM,GAAG,OAAO,CAEnD;AAED;;;;GAIG;AACH,wBAAgB,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO,CAKzD;AAED;;;;GAIG;AACH,wBAAgB,gBAAgB,CAAC,OAAO,EAAE,MAAM,GAAG,OAAO,CAIzD;AAED;;;;GAIG;AACH,wBAAgB,GAAG,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAE3C;AAED;;;;GAIG;AACH,wBAAgB,GAAG,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAE3C;AAED;;;;;GAKG;AACH,wBAAgB,GAAG,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAM3C;AAED;;;;;GAKG;AACH,wBAAgB,IAAI,CAAC,KAAK,EAAE,MAAM,GAAG,OAAO,CAK3C;AAED;;;;;GAKG;AACH,wBAAgB,IAAI,CAAC,KAAK,EAAE,MAAM,GAAG,OAAO,CAK3C;AAED;;;;GAIG;AACH,wBAAgB,IAAI,CAAC,KAAK,EAAE,MAAM,GAAG,OAAO,CAE3C;AAED;;;;;GAKG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,EAAE,CAAC,EAAE,MAAM,GAAG,OAAO,CAEnD;AAED;;;;GAIG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtC;AAED;;;;GAIG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtC;AAED;;;;GAIG;AACH,wBAAgB,IAAI,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEtC;AAED;;;;GAIG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAEvC;AAED;;;;;GAKG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKvC;AAED;;;;;GAKG;AACH,wBAAgB,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,MAAM,CAKvC;AAED;;;;;GAKG;AACH,wBAAgB,GAAG,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAM3C;AAED;;;;;GAKG;AACH,wBAAgB,GAAG,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAM3C;AAED;;;;;GAKG;AACH,wBAAgB,GAAG,CAAC,OAAO,EAAE,MAAM,GAAG,MAAM,CAM3C"}