@johnhenry/math 0.0.0

package/dist/ComplexMath.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ComplexMath.d.ts","sourceRoot":"","sources":["../src/ComplexMath.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AAEnD,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAGrC;;;;;;;;;;;;;;;;;;;GAmBG;AAEH,KAAK,OAAO,GAAG,aAAa,GAAG,MAAM,CAAC;AACtC,KAAK,IAAI,GAAG,MAAM,CAAC,aAAa,CAAC,CAAC;AAClC,KAAK,MAAM,GAAG,CAAC,CAAC,EAAE,aAAa,KAAK,aAAa,CAAC;AAMlD,wBAAgB,SAAS,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,IAAI,CAEvD;AAED,wBAAgB,WAAW,CAAC,KAAK,EAAE,IAAI,EAAE,MAAM,EAAE,OAAO,GAAG,IAAI,CAG9D;AAED,wBAAgB,cAAc,CAAC,KAAK,EAAE,IAAI,GAAG,IAAI,CAEhD;AAED,wBAAgB,cAAc,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,IAAI,CAE5D;AAED,wBAAgB,UAAU,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,aAAa,CAGjE;AAED,2FAA2F;AAC3F,wBAAgB,YAAY,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,IAAI,CAW1D;AAED,wBAAgB,gBAAgB,CAAC,CAAC,EAAE,IAAI,EAAE,CAAC,EAAE,IAAI,GAAG,IAAI,CAIvD;AAED,wBAAgB,KAAK,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,GAAE,OAAW,GAAG,aAAa,CASnE;AAED,wBAAgB,cAAc,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAE,OAAW,GAAG,aAAa,CAExF;AAED,wBAAgB,eAAe,CAAC,KAAK,EAAE,IAAI,GAAG,aAAa,CAE1D;AAED,wBAAgB,YAAY,CAAC,KAAK,EAAE,IAAI,EAAE,IAAI,EAAE,IAAI,GAAG,aAAa,CAInE;AAID,wBAAgB,MAAM,CAAC,CAAC,EAAE,OAAO,EAAE,IAAI,GAAE,OAAW,EAAE,IAAI,GAAE,OAAW,GAAG,aAAa,CAEtF;AAED,wBAAgB,YAAY,CAAC,CAAC,EAAE,OAAO,EAAE,IAAI,GAAE,OAAW,EAAE,IAAI,GAAE,OAAW,GAAG,aAAa,CAE5F;AAED,iDAAiD;AACjD,wBAAgB,kBAAkB,CAAC,IAAI,EAAE,OAAO,EAAE,IAAI,EAAE,OAAO,GAAG,MAAM,CAMvE;AAED,wBAAgB,0BAA0B,CAAC,CAAC,EAAE,OAAO,GAAG,aAAa,CAEpE;AAkBD,wBAAgB,iBAAiB,CAC/B,CAAC,EAAE,OAAO,EACV,IAAI,EAAE,OAAO,EACb,IAAI,EAAE,OAAO,EACb,6BAA6B,SAAK,EAClC,mBAAmB,SAAS,GAC3B,aAAa,CAOf;AAED,wBAAgB,yBAAyB,CACvC,CAAC,EAAE,OAAO,EACV,6BAA6B,SAAK,EAClC,QAAQ,SAAS,GAChB,aAAa,CAEf"}

package/dist/ComplexMath.js

@@ -0,0 +1,98 @@
+import { ComplexNumber } from "./ComplexNumber.js";
+import { integrateN } from "./RealMath.js";
+import { Vector } from "./Vector.js";
+import { VectorUtils } from "./VectorUtils.js";
+const cn = (x) => (x instanceof ComplexNumber ? x : new ComplexNumber(x));
+// -- vectors ------------------------------------------------------------------
+export function addVector(alpha, beta) {
+    return VectorUtils.combine(alpha, beta, (a, b) => a.add(b));
+}
+export function scaleVector(alpha, scalar) {
+    const s = cn(scalar);
+    return VectorUtils.transform(alpha, (x) => s.multiply(x));
+}
+export function negativeVector(alpha) {
+    return scaleVector(alpha, -1);
+}
+export function subtractVector(alpha, beta) {
+    return addVector(alpha, negativeVector(beta));
+}
+export function dotProduct(alpha, beta) {
+    const collapsable = VectorUtils.combine(alpha, beta, (a, b) => a.multiply(b), ComplexNumber.Zero);
+    return VectorUtils.collapse(collapsable, (a, b) => a.add(b));
+}
+/** 3D cross product (bug fix: index 2 for z, presence not truthiness — as in RealMath). */
+export function crossProduct(alpha, beta) {
+    const a = [ComplexNumber.Zero, ComplexNumber.Zero, ComplexNumber.Zero];
+    const b = [ComplexNumber.Zero, ComplexNumber.Zero, ComplexNumber.Zero];
+    for (let idx = 0; idx < 3; idx++) {
+        if (idx < alpha.length && alpha[idx] != null)
+            a[idx] = alpha[idx];
+        if (idx < beta.length && beta[idx] != null)
+            b[idx] = beta[idx];
+    }
+    const i = a[1].multiply(b[2]).subtract(a[2].multiply(b[1]));
+    const j = a[2].multiply(b[0]).subtract(a[0].multiply(b[2]));
+    const k = a[0].multiply(b[1]).subtract(a[1].multiply(b[0]));
+    return Vector.fromArray([i, j, k]);
+}
+export function kroneckerProduct(A, B) {
+    let kronecker = VectorUtils.constantVector(A.length, null);
+    kronecker = VectorUtils.fillByIndex(kronecker, (i) => scaleVector(B, A[i]));
+    return VectorUtils.flattenSDLevels(kronecker, 2);
+}
+export function pNorm(alpha, norm = 2) {
+    const p = cn(norm);
+    if (p.value === 0) {
+        const magnitudes = VectorUtils.transform(alpha, (x) => new ComplexNumber(x.magnitude()));
+        return VectorUtils.collapse(magnitudes, (a, b) => (a.magCompare(b) >= 0 ? a : b));
+    }
+    const summable = VectorUtils.transform(alpha, (x) => new ComplexNumber(x.magnitude()).power(p.value));
+    const sum = VectorUtils.collapse(summable, (a, b) => a.add(b), ComplexNumber.Zero);
+    return sum.power(1 / p.value);
+}
+export function distanceVector(alpha, beta, norm = 2) {
+    return pNorm(subtractVector(alpha, beta), norm);
+}
+export function magnitudeVector(alpha) {
+    return pNorm(alpha, 2);
+}
+export function angleBetween(alpha, beta) {
+    return dotProduct(alpha, beta)
+        .divide(magnitudeVector(alpha).multiply(magnitudeVector(beta)))
+        .arcCosine();
+}
+// -- probability ----------------------------------------------------------------
+export function zScore(x, mean = 0, sDev = 1) {
+    return cn(x).subtract(mean).divide(sDev);
+}
+export function invertZScore(z, mean = 0, sDev = 1) {
+    return cn(z).multiply(sDev).add(mean);
+}
+/** Normal PDF closure over complex arguments. */
+export function normalDistribution(mean, sDev) {
+    return (x) => cn(sDev)
+        .multiply(Math.sqrt(2 * Math.PI))
+        .reciprocal()
+        .multiply(ComplexNumber.E.power(zScore(x, mean, sDev).square().multiply(-0.5)));
+}
+export function standardNormalDistribution(x) {
+    return normalDistribution(0, 1)(cn(x));
+}
+/**
+ * Numeric integral of a complex-valued function over a real interval, via the
+ * real {@link integrateN} applied component-wise (bug fix: same midpoint-rule
+ * fix as `RealMath.integrateN`).
+ */
+function integrateComplexN(unaryOperation, lowerLimit, upperLimit, interval = 0.0001) {
+    const re = integrateN((t) => unaryOperation(new ComplexNumber(t, 0)).value, lowerLimit, upperLimit, interval);
+    const im = integrateN((t) => unaryOperation(new ComplexNumber(t, 0)).iValue, lowerLimit, upperLimit, interval);
+    return new ComplexNumber(re, im);
+}
+export function normalProbability(x, mean, sDev, negativeInfinityApproximation = -5, integrationInterval = 0.0001) {
+    return integrateComplexN(normalDistribution(mean, sDev), negativeInfinityApproximation, cn(x).value, integrationInterval);
+}
+export function standardNormalProbability(z, negativeInfinityApproximation = -5, interval = 0.0001) {
+    return integrateComplexN(standardNormalDistribution, negativeInfinityApproximation, cn(z).value, interval);
+}
+//# sourceMappingURL=ComplexMath.js.map

package/dist/ComplexMath.js.map

@@ -0,0 +1 @@
+{"version":3,"file":"ComplexMath.js","sourceRoot":"","sources":["../src/ComplexMath.ts"],"names":[],"mappings":"AAAA,OAAO,EAAE,aAAa,EAAE,MAAM,oBAAoB,CAAC;AACnD,OAAO,EAAE,UAAU,EAAE,MAAM,eAAe,CAAC;AAC3C,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AACrC,OAAO,EAAE,WAAW,EAAE,MAAM,kBAAkB,CAAC;AA2B/C,MAAM,EAAE,GAAG,CAAC,CAAU,EAAiB,EAAE,CAAC,CAAC,CAAC,YAAY,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,IAAI,aAAa,CAAC,CAAC,CAAC,CAAC,CAAC;AAElG,gFAAgF;AAEhF,MAAM,UAAU,SAAS,CAAC,KAAW,EAAE,IAAU;IAC/C,OAAO,WAAW,CAAC,OAAO,CAAC,KAAK,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAS,CAAC;AACtE,CAAC;AAED,MAAM,UAAU,WAAW,CAAC,KAAW,EAAE,MAAe;IACtD,MAAM,CAAC,GAAG,EAAE,CAAC,MAAM,CAAC,CAAC;IACrB,OAAO,WAAW,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC;AAC5D,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,KAAW;IACxC,OAAO,WAAW,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,CAAC;AAChC,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,KAAW,EAAE,IAAU;IACpD,OAAO,SAAS,CAAC,KAAK,EAAE,cAAc,CAAC,IAAI,CAAC,CAAC,CAAC;AAChD,CAAC;AAED,MAAM,UAAU,UAAU,CAAC,KAAW,EAAE,IAAU;IAChD,MAAM,WAAW,GAAG,WAAW,CAAC,OAAO,CAAC,KAAK,EAAE,IAAI,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,aAAa,CAAC,IAAI,CAAS,CAAC;IAC1G,OAAO,WAAW,CAAC,QAAQ,CAAC,WAAW,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAkB,CAAC;AAChF,CAAC;AAED,2FAA2F;AAC3F,MAAM,UAAU,YAAY,CAAC,KAAW,EAAE,IAAU;IAClD,MAAM,CAAC,GAAG,CAAC,aAAa,CAAC,IAAI,EAAE,aAAa,CAAC,IAAI,EAAE,aAAa,CAAC,IAAI,CAAC,CAAC;IACvE,MAAM,CAAC,GAAG,CAAC,aAAa,CAAC,IAAI,EAAE,aAAa,CAAC,IAAI,EAAE,aAAa,CAAC,IAAI,CAAC,CAAC;IACvE,KAAK,IAAI,GAAG,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,GAAG,EAAE,EAAE,CAAC;QACjC,IAAI,GAAG,GAAG,KAAK,CAAC,MAAM,IAAI,KAAK,CAAC,GAAG,CAAC,IAAI,IAAI;YAAE,CAAC,CAAC,GAAG,CAAC,GAAG,KAAK,CAAC,GAAG,CAAkB,CAAC;QACnF,IAAI,GAAG,GAAG,IAAI,CAAC,MAAM,IAAI,IAAI,CAAC,GAAG,CAAC,IAAI,IAAI;YAAE,CAAC,CAAC,GAAG,CAAC,GAAG,IAAI,CAAC,GAAG,CAAkB,CAAC;IAClF,CAAC;IACD,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,MAAM,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC;IAC5D,OAAO,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC,EAAE,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC;AACrC,CAAC;AAED,MAAM,UAAU,gBAAgB,CAAC,CAAO,EAAE,CAAO;IAC/C,IAAI,SAAS,GAAG,WAAW,CAAC,cAAc,CAAC,CAAC,CAAC,MAAM,EAAE,IAAe,CAAoB,CAAC;IACzF,SAAS,GAAG,WAAW,CAAC,WAAW,CAAC,SAAS,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,WAAW,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAkB,CAAC,CAAC,CAAC;IAC7F,OAAO,WAAW,CAAC,eAAe,CAAC,SAAS,EAAE,CAAC,CAAS,CAAC;AAC3D,CAAC;AAED,MAAM,UAAU,KAAK,CAAC,KAAW,EAAE,OAAgB,CAAC;IAClD,MAAM,CAAC,GAAG,EAAE,CAAC,IAAI,CAAC,CAAC;IACnB,IAAI,CAAC,CAAC,KAAK,KAAK,CAAC,EAAE,CAAC;QAClB,MAAM,UAAU,GAAG,WAAW,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,aAAa,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC,CAAC;QACzF,OAAO,WAAW,CAAC,QAAQ,CAAC,UAAU,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,UAAU,CAAC,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAC,CAAkB,CAAC;IACrG,CAAC;IACD,MAAM,QAAQ,GAAG,WAAW,CAAC,SAAS,CAAC,KAAK,EAAE,CAAC,CAAC,EAAE,EAAE,CAAC,IAAI,aAAa,CAAC,CAAC,CAAC,SAAS,EAAE,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC,KAAK,CAAC,CAAC,CAAC;IACtG,MAAM,GAAG,GAAG,WAAW,CAAC,QAAQ,CAAC,QAAQ,EAAE,CAAC,CAAC,EAAE,CAAC,EAAE,EAAE,CAAC,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,EAAE,aAAa,CAAC,IAAI,CAAkB,CAAC;IACpG,OAAO,GAAG,CAAC,KAAK,CAAC,CAAC,GAAG,CAAC,CAAC,KAAK,CAAC,CAAC;AAChC,CAAC;AAED,MAAM,UAAU,cAAc,CAAC,KAAW,EAAE,IAAU,EAAE,OAAgB,CAAC;IACvE,OAAO,KAAK,CAAC,cAAc,CAAC,KAAK,EAAE,IAAI,CAAC,EAAE,IAAI,CAAC,CAAC;AAClD,CAAC;AAED,MAAM,UAAU,eAAe,CAAC,KAAW;IACzC,OAAO,KAAK,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC;AACzB,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,KAAW,EAAE,IAAU;IAClD,OAAO,UAAU,CAAC,KAAK,EAAE,IAAI,CAAC;SAC3B,MAAM,CAAC,eAAe,CAAC,KAAK,CAAC,CAAC,QAAQ,CAAC,eAAe,CAAC,IAAI,CAAC,CAAC,CAAC;SAC9D,SAAS,EAAE,CAAC;AACjB,CAAC;AAED,kFAAkF;AAElF,MAAM,UAAU,MAAM,CAAC,CAAU,EAAE,OAAgB,CAAC,EAAE,OAAgB,CAAC;IACrE,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,MAAM,CAAC,IAAI,CAAC,CAAC;AAC3C,CAAC;AAED,MAAM,UAAU,YAAY,CAAC,CAAU,EAAE,OAAgB,CAAC,EAAE,OAAgB,CAAC;IAC3E,OAAO,EAAE,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,IAAI,CAAC,CAAC;AACxC,CAAC;AAED,iDAAiD;AACjD,MAAM,UAAU,kBAAkB,CAAC,IAAa,EAAE,IAAa;IAC7D,OAAO,CAAC,CAAgB,EAAE,EAAE,CAC1B,EAAE,CAAC,IAAI,CAAC;SACL,QAAQ,CAAC,IAAI,CAAC,IAAI,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,CAAC,CAAC;SAChC,UAAU,EAAE;SACZ,QAAQ,CAAC,aAAa,CAAC,CAAC,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,IAAI,EAAE,IAAI,CAAC,CAAC,MAAM,EAAE,CAAC,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;AACtF,CAAC;AAED,MAAM,UAAU,0BAA0B,CAAC,CAAU;IACnD,OAAO,kBAAkB,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;AACzC,CAAC;AAED;;;;GAIG;AACH,SAAS,iBAAiB,CACxB,cAAsB,EACtB,UAAkB,EAClB,UAAkB,EAClB,QAAQ,GAAG,MAAM;IAEjB,MAAM,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,cAAc,CAAC,IAAI,aAAa,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;IAC9G,MAAM,EAAE,GAAG,UAAU,CAAC,CAAC,CAAC,EAAE,EAAE,CAAC,cAAc,CAAC,IAAI,aAAa,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,MAAM,EAAE,UAAU,EAAE,UAAU,EAAE,QAAQ,CAAC,CAAC;IAC/G,OAAO,IAAI,aAAa,CAAC,EAAE,EAAE,EAAE,CAAC,CAAC;AACnC,CAAC;AAED,MAAM,UAAU,iBAAiB,CAC/B,CAAU,EACV,IAAa,EACb,IAAa,EACb,6BAA6B,GAAG,CAAC,CAAC,EAClC,mBAAmB,GAAG,MAAM;IAE5B,OAAO,iBAAiB,CACtB,kBAAkB,CAAC,IAAI,EAAE,IAAI,CAAC,EAC9B,6BAA6B,EAC7B,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,EACX,mBAAmB,CACpB,CAAC;AACJ,CAAC;AAED,MAAM,UAAU,yBAAyB,CACvC,CAAU,EACV,6BAA6B,GAAG,CAAC,CAAC,EAClC,QAAQ,GAAG,MAAM;IAEjB,OAAO,iBAAiB,CAAC,0BAA0B,EAAE,6BAA6B,EAAE,EAAE,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,QAAQ,CAAC,CAAC;AAC7G,CAAC"}

package/dist/ComplexNumber.d.ts

@@ -0,0 +1,174 @@
+import { Vector } from "./Vector.ts";
+/** Anything that can stand in for a complex operand in a fluent method call. */
+export type CNInput = ComplexNumber | number;
+/**
+ * ComplexNumber — an immutable-ish complex value `value + iValue·i`, with a
+ * fluent arithmetic/trigonometric API matching the style of `Quaternion`/
+ * `Rational`/`DualNumber`/`Decimal`.
+ *
+ * Ported from Mallory's ActionScript `ComplexNumber`. The AS3 class extended
+ * `Array` and stashed its two coordinates in a private `Vector`; that indirection
+ * bought nothing, so here the real and imaginary parts are plain number fields.
+ * The historical accessor names `value` (real) and `iValue` (imaginary) are kept
+ * because the rest of the library reads them hundreds of times; `re`/`im` are
+ * provided as modern aliases.
+ *
+ * The arithmetic/trig/log methods and the `E`/`I`/`PI`/`PHI`/`Zero`/`One`
+ * constants and `random()` factory were previously the static-method "namespace
+ * class" `ComplexMath` — folded in here so a complex value carries its own
+ * operations rather than requiring a second class to operate on it.
+ *
+ * Bug fixes relative to the AS3 original:
+ *  - `parse` (was `fromString`) no longer throws when `String.match` returns
+ *    `null`, and it correctly handles negative imaginary parts and the bare
+ *    `i` / `-i` forms, so it round-trips with `toString`.
+ *  - `toString` renders a negative imaginary part as `a-b*i` rather than the
+ *    old `a+-b*i`.
+ *  - `divide`'s zero-divisor branches compared `alpha.ivalue` (lowercase v),
+ *    which is always `undefined`, so the eight directed infinities were never
+ *    returned. Fixed to `iValue`.
+ *  - `power`: a zero base produced `NaN` (via `ln 0`); now handled directly, so
+ *    `square(0)`, `squareRoot(0)` and e.g. `arcSine(0)` are correct.
+ */
+export declare class ComplexNumber {
+    value: number;
+    iValue: number;
+    /** "Not a Complex Number" — both parts NaN. */
+    static readonly NaCN: ComplexNumber;
+    static readonly PositiveInfinity: ComplexNumber;
+    static readonly InfinityQ1: ComplexNumber;
+    static readonly PositiveInfinityI: ComplexNumber;
+    static readonly InfinityQ2: ComplexNumber;
+    static readonly NegativeInfinity: ComplexNumber;
+    static readonly InfinityQ3: ComplexNumber;
+    static readonly NegativeInfinityI: ComplexNumber;
+    static readonly InfinityQ4: ComplexNumber;
+    static readonly Zero: ComplexNumber;
+    static readonly One: ComplexNumber;
+    static readonly I: ComplexNumber;
+    static readonly E: ComplexNumber;
+    static readonly PI: ComplexNumber;
+    static readonly PHI: ComplexNumber;
+    /**
+     * Construct a complex number. Accepts:
+     *  - `()` → 0
+     *  - `(re, im)` → re + im·i
+     *  - `(ComplexNumber)` → a copy
+     *  - `(number)` → a real number
+     *  - `(string)` → parsed (see {@link ComplexNumber.parse})
+     */
+    constructor(value?: number | ComplexNumber | string, iValue?: number);
+    /** Modern alias for {@link value} (the real part). */
+    get re(): number;
+    set re(v: number);
+    /** Modern alias for {@link iValue} (the imaginary part). */
+    get im(): number;
+    set im(v: number);
+    /**
+     * Coerce an arbitrary value into a ComplexNumber, returning {@link NaCN} for
+     * anything that cannot be interpreted numerically. Safe replacement for the
+     * AS3 pattern `new ComplexNumber(element)` used by {@link isComplex}.
+     */
+    static from(input: unknown): ComplexNumber;
+    /** True when `input` is NOT a finite/definite complex number. */
+    static isNotComplex(input: unknown): boolean;
+    /** True when `input` can be interpreted as a complex number. */
+    static isComplex(input: unknown): boolean;
+    /**
+     * Wrap arguments into a ComplexNumber. Mirrors AS3 `Wrap`: one argument is
+     * coerced, two are taken as (re, im). Returns {@link NaCN} for other arities
+     * (the AS3 version returned `false`; a ComplexNumber is friendlier and typed).
+     */
+    static wrap(...args: Array<number | ComplexNumber | string>): ComplexNumber;
+    /** Build a real complex number from a plain number. */
+    static fromNumber(num: number): ComplexNumber;
+    /** Build a complex number from polar form `r·e^(iθ)` (inverse of magnitude/angle). */
+    static fromPolar(magnitude: number, angle: number): ComplexNumber;
+    /** Build a complex number from a `[re, im]` vector-like (inverse of {@link toVector}). */
+    static fromVector(coordinates: ArrayLike<number>): ComplexNumber;
+    /** Parse the `<complexNumber>` form produced by {@link toXML} (inverse of it). */
+    static fromXML(xml: string): ComplexNumber;
+    /**
+     * Parse a string such as `"3"`, `"-2*i"`, `"i"`, `"3+2*i"`, `"3-2i"`.
+     * The `*` between coefficient and `i` is optional. Whitespace is ignored.
+     * Returns {@link NaCN} when the string is not a recognisable complex literal.
+     */
+    static parse(expression: string): ComplexNumber;
+    /** Uniform random complex number, real part in `[min(r1,r2), max(r1,r2)]`, imaginary in `[min(i1,i2), max(i1,i2)]`. */
+    static random(r1?: number, r2?: number, i1?: number, i2?: number, _inclusive?: boolean): ComplexNumber;
+    /** XML-ish serialisation kept for API compatibility. */
+    toXML(): string;
+    /** The two coordinates as a fresh {@link Vector}. */
+    toVector(): Vector<number>;
+    /**
+     * String representation. In `fullMode` the raw `a+bi` form is emitted (used
+     * for debugging); otherwise a canonical, sign-correct form is produced that
+     * round-trips through {@link parse}.
+     */
+    toString(fullMode?: boolean): string;
+    /** Additive inverse `-z`. */
+    neg(): ComplexNumber;
+    /** Multiplicative inverse `1/z` (returns {@link NaCN} for zero). */
+    reciprocal(): ComplexNumber;
+    /** Complex conjugate `a - b·i`. */
+    conjugate(): ComplexNumber;
+    /** Swap the real and imaginary parts (`a + b·i` → `b + a·i`). */
+    flip(): ComplexNumber;
+    /** Structural equality (NaN parts compare unequal, matching IEEE semantics). */
+    equals(other: ComplexNumber): boolean;
+    /** Lexicographic comparison: real part first, then imaginary. */
+    lexCompare(other: CNInput): number;
+    /** Comparison by magnitude. */
+    magCompare(other: CNInput): number;
+    add(other?: CNInput): ComplexNumber;
+    subtract(other?: CNInput): ComplexNumber;
+    multiply(other?: CNInput): ComplexNumber;
+    /** Division, returning one of the eight directed infinities on a zero divisor. */
+    divide(other?: CNInput): ComplexNumber;
+    magnitude(): number;
+    angle(): number;
+    private static power2;
+    private selectedNaturalLogarithm;
+    /** `this ** exponent`, with branch selection for the multivalued complex logarithm. */
+    power(exponent?: CNInput, selector?: number): ComplexNumber;
+    square(): ComplexNumber;
+    squareRoot(): ComplexNumber;
+    /** Logarithm in an arbitrary `base` (natural log by default), with branch selection. */
+    logarithm(base?: CNInput, selectA?: number, selectB?: number): ComplexNumber;
+    sine(): ComplexNumber;
+    cosine(): ComplexNumber;
+    tangent(): ComplexNumber;
+    hyperbolicSine(): ComplexNumber;
+    hyperbolicCosine(): ComplexNumber;
+    hyperbolicTangent(): ComplexNumber;
+    arcSine(): ComplexNumber;
+    arcCosine(): ComplexNumber;
+    arcTangent(): ComplexNumber;
+    /** `asinh(z) = log(z + sqrt(z^2 + 1))`. */
+    arcHyperbolicSine(): ComplexNumber;
+    /**
+     * `acosh(z) = log(z + sqrt(z+1)*sqrt(z-1))` — the numerically-careful
+     * factored form, not the naive `log(z + sqrt(z^2-1))`. The two are
+     * algebraically identical but the factored form keeps the two square
+     * roots' branch cuts separate (one for `z+1`, one for `z-1`), which
+     * matches the principal-branch convention every serious complex `acosh`
+     * implementation uses (e.g. the C99/IEEE 754-2008 `cacosh` reference
+     * algorithm) — the naive form can pick the wrong branch/sign near its
+     * own cut along the real axis below 1.
+     */
+    arcHyperbolicCosine(): ComplexNumber;
+    /** `atanh(z) = 0.5 * log((1+z)/(1-z))`. */
+    arcHyperbolicTangent(): ComplexNumber;
+    /**
+     * Coercion safety: a genuinely complex value (nonzero imaginary part)
+     * cannot be implicitly converted to a `Number` — without this override,
+     * `+z`/`z * 2`/etc. on a non-real `ComplexNumber` silently produce `NaN`
+     * (the default `Object`-to-primitive path falls through to parsing
+     * `toString()`'s `"a+bi"` output as a number), which is exactly the kind
+     * of silent-wrong-answer failure this library avoids elsewhere. A purely
+     * real `ComplexNumber` (iValue === 0) still converts safely, matching
+     * `ℝ ⊂ ℂ`.
+     */
+    valueOf(): number;
+}
+//# sourceMappingURL=ComplexNumber.d.ts.map

package/dist/ComplexNumber.d.ts.map

@@ -0,0 +1 @@
+{"version":3,"file":"ComplexNumber.d.ts","sourceRoot":"","sources":["../src/ComplexNumber.ts"],"names":[],"mappings":"AACA,OAAO,EAAE,MAAM,EAAE,MAAM,aAAa,CAAC;AAErC,gFAAgF;AAChF,MAAM,MAAM,OAAO,GAAG,aAAa,GAAG,MAAM,CAAC;AAW7C;;;;;;;;;;;;;;;;;;;;;;;;;;;;GA4BG;AACH,qBAAa,aAAa;IACxB,KAAK,EAAE,MAAM,CAAC;IACd,MAAM,EAAE,MAAM,CAAC;IAEf,+CAA+C;IAC/C,MAAM,CAAC,QAAQ,CAAC,IAAI,gBAA+B;IAGnD,MAAM,CAAC,QAAQ,CAAC,gBAAgB,gBAAkC;IAClE,MAAM,CAAC,QAAQ,CAAC,UAAU,gBAAyC;IACnE,MAAM,CAAC,QAAQ,CAAC,iBAAiB,gBAAkC;IACnE,MAAM,CAAC,QAAQ,CAAC,UAAU,gBAA0C;IACpE,MAAM,CAAC,QAAQ,CAAC,gBAAgB,gBAAmC;IACnE,MAAM,CAAC,QAAQ,CAAC,UAAU,gBAA2C;IACrE,MAAM,CAAC,QAAQ,CAAC,iBAAiB,gBAAmC;IACpE,MAAM,CAAC,QAAQ,CAAC,UAAU,gBAA0C;IAEpE,MAAM,CAAC,QAAQ,CAAC,IAAI,gBAA2B;IAC/C,MAAM,CAAC,QAAQ,CAAC,GAAG,gBAA2B;IAC9C,MAAM,CAAC,QAAQ,CAAC,CAAC,gBAA2B;IAC5C,MAAM,CAAC,QAAQ,CAAC,CAAC,gBAAgC;IACjD,MAAM,CAAC,QAAQ,CAAC,EAAE,gBAAiC;IACnD,MAAM,CAAC,QAAQ,CAAC,GAAG,gBAAgD;IAEnE;;;;;;;OAOG;gBACS,KAAK,GAAE,MAAM,GAAG,aAAa,GAAG,MAAU,EAAE,MAAM,SAAI;IAclE,sDAAsD;IACtD,IAAI,EAAE,IAAI,MAAM,CAEf;IACD,IAAI,EAAE,CAAC,CAAC,EAAE,MAAM,EAEf;IAED,4DAA4D;IAC5D,IAAI,EAAE,IAAI,MAAM,CAEf;IACD,IAAI,EAAE,CAAC,CAAC,EAAE,MAAM,EAEf;IAED;;;;OAIG;IACH,MAAM,CAAC,IAAI,CAAC,KAAK,EAAE,OAAO,GAAG,aAAa;IAO1C,iEAAiE;IACjE,MAAM,CAAC,YAAY,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAK5C,gEAAgE;IAChE,MAAM,CAAC,SAAS,CAAC,KAAK,EAAE,OAAO,GAAG,OAAO;IAKzC;;;;OAIG;IACH,MAAM,CAAC,IAAI,CAAC,GAAG,IAAI,EAAE,KAAK,CAAC,MAAM,GAAG,aAAa,GAAG,MAAM,CAAC,GAAG,aAAa;IAM3E,uDAAuD;IACvD,MAAM,CAAC,UAAU,CAAC,GAAG,EAAE,MAAM,GAAG,aAAa;IAI7C,sFAAsF;IACtF,MAAM,CAAC,SAAS,CAAC,SAAS,EAAE,MAAM,EAAE,KAAK,EAAE,MAAM,GAAG,aAAa;IAIjE,0FAA0F;IAC1F,MAAM,CAAC,UAAU,CAAC,WAAW,EAAE,SAAS,CAAC,MAAM,CAAC,GAAG,aAAa;IAIhE,kFAAkF;IAClF,MAAM,CAAC,OAAO,CAAC,GAAG,EAAE,MAAM,GAAG,aAAa;IAM1C;;;;OAIG;IACH,MAAM,CAAC,KAAK,CAAC,UAAU,EAAE,MAAM,GAAG,aAAa;IA4C/C,uHAAuH;IACvH,MAAM,CAAC,MAAM,CAAC,EAAE,SAAI,EAAE,EAAE,SAAI,EAAE,EAAE,SAAI,EAAE,EAAE,SAAI,EAAE,UAAU,UAAO,GAAG,aAAa;IAQ/E,wDAAwD;IACxD,KAAK,IAAI,MAAM;IAIf,qDAAqD;IACrD,QAAQ,IAAI,MAAM,CAAC,MAAM,CAAC;IAI1B;;;;OAIG;IACH,QAAQ,CAAC,QAAQ,UAAQ,GAAG,MAAM;IAalC,6BAA6B;IAC7B,GAAG,IAAI,aAAa;IAIpB,oEAAoE;IACpE,UAAU,IAAI,aAAa;IAO3B,mCAAmC;IACnC,SAAS,IAAI,aAAa;IAI1B,iEAAiE;IACjE,IAAI,IAAI,aAAa;IAIrB,gFAAgF;IAChF,MAAM,CAAC,KAAK,EAAE,aAAa,GAAG,OAAO;IAIrC,iEAAiE;IACjE,UAAU,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IASlC,+BAA+B;IAC/B,UAAU,CAAC,KAAK,EAAE,OAAO,GAAG,MAAM;IAQlC,GAAG,CAAC,KAAK,GAAE,OAAW,GAAG,aAAa;IAKtC,QAAQ,CAAC,KAAK,GAAE,OAAW,GAAG,aAAa;IAK3C,QAAQ,CAAC,KAAK,GAAE,OAAW,GAAG,aAAa;IAQ3C,kFAAkF;IAClF,MAAM,CAAC,KAAK,GAAE,OAAW,GAAG,aAAa;IAoBzC,SAAS,IAAI,MAAM;IAInB,KAAK,IAAI,MAAM;IAIf,OAAO,CAAC,MAAM,CAAC,MAAM;IAcrB,OAAO,CAAC,wBAAwB;IAMhC,uFAAuF;IACvF,KAAK,CAAC,QAAQ,GAAE,OAAW,EAAE,QAAQ,SAAI,GAAG,aAAa;IAYzD,MAAM,IAAI,aAAa;IAIvB,UAAU,IAAI,aAAa;IAI3B,wFAAwF;IACxF,SAAS,CAAC,IAAI,GAAE,OAAgB,EAAE,OAAO,SAAI,EAAE,OAAO,SAAI,GAAG,aAAa;IAK1E,IAAI,IAAI,aAAa;IAOrB,MAAM,IAAI,aAAa;IAOvB,OAAO,IAAI,aAAa;IAIxB,cAAc,IAAI,aAAa;IAO/B,gBAAgB,IAAI,aAAa;IAOjC,iBAAiB,IAAI,aAAa;IAIlC,OAAO,IAAI,aAAa;IAKxB,SAAS,IAAI,aAAa;IAK1B,UAAU,IAAI,aAAa;IAO3B,2CAA2C;IAC3C,iBAAiB,IAAI,aAAa;IAKlC;;;;;;;;;OASG;IACH,mBAAmB,IAAI,aAAa;IAKpC,2CAA2C;IAC3C,oBAAoB,IAAI,aAAa;IAMrC;;;;;;;;;OASG;IACH,OAAO,IAAI,MAAM;CASlB"}