@bsv/sdk 1.6.18 → 1.6.20

package/dist/cjs/src/primitives/Point.js

@@ -3,11 +3,216 @@ var __importDefault = (this && this.__importDefault) || function (mod) {
     return (mod && mod.__esModule) ? mod : { "default": mod };
 };
 Object.defineProperty(exports, "__esModule", { value: true });
+exports.modInvN = exports.modMulN = exports.modN = exports.scalarMultiplyWNAF = exports.jpNeg = exports.jpAdd = exports.jpDouble = exports.GY_BIGINT = exports.GX_BIGINT = exports.biModSqrt = exports.P_PLUS1_DIV4 = exports.biModPow = exports.biModSqr = exports.biModInv = exports.biModAdd = exports.biModMul = exports.biModSub = exports.biMod = exports.red = exports.MASK_256 = exports.N_BIGINT = exports.P_BIGINT = exports.BI_EIGHT = exports.BI_FOUR = exports.BI_THREE = exports.BI_TWO = exports.BI_ONE = exports.BI_ZERO = void 0;
 const BasePoint_js_1 = __importDefault(require("./BasePoint.js"));
 const JacobianPoint_js_1 = __importDefault(require("./JacobianPoint.js"));
 const BigNumber_js_1 = __importDefault(require("./BigNumber.js"));
 const utils_js_1 = require("./utils.js");
-const ReductionContext_js_1 = __importDefault(require("./ReductionContext.js"));
+// -----------------------------------------------------------------------------
+// BigInt helpers & constants (secp256k1) – hoisted so we don't recreate them on
+// every Point.mul() call.
+// -----------------------------------------------------------------------------
+exports.BI_ZERO = 0n;
+exports.BI_ONE = 1n;
+exports.BI_TWO = 2n;
+exports.BI_THREE = 3n;
+exports.BI_FOUR = 4n;
+exports.BI_EIGHT = 8n;
+exports.P_BIGINT = 0xfffffffffffffffffffffffffffffffffffffffffffffffffffffffefffffc2fn;
+exports.N_BIGINT = 0xfffffffffffffffffffffffffffffffebaaedce6af48a03bbfd25e8cd0364141n;
+exports.MASK_256 = (1n << 256n) - 1n; // 0xffff…ffff (256 sones)
+function red(x) {
+    // first fold
+    let hi = x >> 256n;
+    x = (x & exports.MASK_256) + (hi << 32n) + hi * 977n;
+    // second fold  (hi ≤ 2³² + 977 here, so one more pass is enough)
+    hi = x >> 256n;
+    x = (x & exports.MASK_256) + (hi << 32n) + hi * 977n;
+    // final conditional subtraction
+    if (x >= exports.P_BIGINT)
+        x -= exports.P_BIGINT;
+    return x;
+}
+exports.red = red;
+const biMod = (a) => red((a % exports.P_BIGINT + exports.P_BIGINT) % exports.P_BIGINT);
+exports.biMod = biMod;
+const biModSub = (a, b) => (a >= b ? a - b : exports.P_BIGINT - (b - a));
+exports.biModSub = biModSub;
+const biModMul = (a, b) => red(a * b);
+exports.biModMul = biModMul;
+const biModAdd = (a, b) => red(a + b);
+exports.biModAdd = biModAdd;
+const biModInv = (a) => {
+    let lm = exports.BI_ONE;
+    let hm = exports.BI_ZERO;
+    let low = (0, exports.biMod)(a);
+    let high = exports.P_BIGINT;
+    while (low > exports.BI_ONE) {
+        const r = high / low;
+        [lm, hm] = [hm - lm * r, lm];
+        [low, high] = [high - low * r, low];
+    }
+    return (0, exports.biMod)(lm);
+};
+exports.biModInv = biModInv;
+const biModSqr = (a) => (0, exports.biModMul)(a, a);
+exports.biModSqr = biModSqr;
+const biModPow = (base, exp) => {
+    let result = exports.BI_ONE;
+    base = (0, exports.biMod)(base);
+    let e = exp;
+    while (e > exports.BI_ZERO) {
+        if ((e & exports.BI_ONE) === exports.BI_ONE)
+            result = (0, exports.biModMul)(result, base);
+        base = (0, exports.biModMul)(base, base);
+        e >>= exports.BI_ONE;
+    }
+    return result;
+};
+exports.biModPow = biModPow;
+exports.P_PLUS1_DIV4 = (exports.P_BIGINT + 1n) >> 2n;
+const biModSqrt = (a) => {
+    const r = (0, exports.biModPow)(a, exports.P_PLUS1_DIV4);
+    return (0, exports.biModMul)(r, r) === (0, exports.biMod)(a) ? r : null;
+};
+exports.biModSqrt = biModSqrt;
+const toBigInt = (x) => {
+    if (BigNumber_js_1.default.isBN(x))
+        return BigInt('0x' + x.toString(16));
+    if (typeof x === 'string')
+        return BigInt('0x' + x);
+    if (Array.isArray(x))
+        return BigInt('0x' + (0, utils_js_1.toHex)(x));
+    return BigInt(x);
+};
+// Generator point coordinates as bigint constants
+exports.GX_BIGINT = BigInt('0x79be667ef9dcbbac55a06295ce870b07029bfcdb2dce28d959f2815b16f81798');
+exports.GY_BIGINT = BigInt('0x483ada7726a3c4655da4fbfc0e1108a8fd17b448a68554199c47d08ffb10d4b8');
+// Cache for precomputed windowed tables keyed by 'window:x:y'
+const WNAF_TABLE_CACHE = new Map();
+const jpDouble = (P) => {
+    const { X: X1, Y: Y1, Z: Z1 } = P;
+    if (Y1 === exports.BI_ZERO)
+        return { X: exports.BI_ZERO, Y: exports.BI_ONE, Z: exports.BI_ZERO };
+    const Y1sq = (0, exports.biModMul)(Y1, Y1);
+    const S = (0, exports.biModMul)(exports.BI_FOUR, (0, exports.biModMul)(X1, Y1sq));
+    const M = (0, exports.biModMul)(exports.BI_THREE, (0, exports.biModMul)(X1, X1));
+    const X3 = (0, exports.biModSub)((0, exports.biModMul)(M, M), (0, exports.biModMul)(exports.BI_TWO, S));
+    const Y3 = (0, exports.biModSub)((0, exports.biModMul)(M, (0, exports.biModSub)(S, X3)), (0, exports.biModMul)(exports.BI_EIGHT, (0, exports.biModMul)(Y1sq, Y1sq)));
+    const Z3 = (0, exports.biModMul)(exports.BI_TWO, (0, exports.biModMul)(Y1, Z1));
+    return { X: X3, Y: Y3, Z: Z3 };
+};
+exports.jpDouble = jpDouble;
+const jpAdd = (P, Q) => {
+    if (P.Z === exports.BI_ZERO)
+        return Q;
+    if (Q.Z === exports.BI_ZERO)
+        return P;
+    const Z1Z1 = (0, exports.biModMul)(P.Z, P.Z);
+    const Z2Z2 = (0, exports.biModMul)(Q.Z, Q.Z);
+    const U1 = (0, exports.biModMul)(P.X, Z2Z2);
+    const U2 = (0, exports.biModMul)(Q.X, Z1Z1);
+    const S1 = (0, exports.biModMul)(P.Y, (0, exports.biModMul)(Z2Z2, Q.Z));
+    const S2 = (0, exports.biModMul)(Q.Y, (0, exports.biModMul)(Z1Z1, P.Z));
+    const H = (0, exports.biModSub)(U2, U1);
+    const r = (0, exports.biModSub)(S2, S1);
+    if (H === exports.BI_ZERO) {
+        if (r === exports.BI_ZERO)
+            return (0, exports.jpDouble)(P);
+        return { X: exports.BI_ZERO, Y: exports.BI_ONE, Z: exports.BI_ZERO }; // Infinity
+    }
+    const HH = (0, exports.biModMul)(H, H);
+    const HHH = (0, exports.biModMul)(H, HH);
+    const V = (0, exports.biModMul)(U1, HH);
+    const X3 = (0, exports.biModSub)((0, exports.biModSub)((0, exports.biModMul)(r, r), HHH), (0, exports.biModMul)(exports.BI_TWO, V));
+    const Y3 = (0, exports.biModSub)((0, exports.biModMul)(r, (0, exports.biModSub)(V, X3)), (0, exports.biModMul)(S1, HHH));
+    const Z3 = (0, exports.biModMul)(H, (0, exports.biModMul)(P.Z, Q.Z));
+    return { X: X3, Y: Y3, Z: Z3 };
+};
+exports.jpAdd = jpAdd;
+const jpNeg = (P) => {
+    if (P.Z === exports.BI_ZERO)
+        return P;
+    return { X: P.X, Y: exports.P_BIGINT - P.Y, Z: P.Z };
+};
+exports.jpNeg = jpNeg;
+// Fast windowed-NAF scalar multiplication (default window = 5) in Jacobian
+// coordinates.  Returns Q = k * P0 as a JacobianPoint.
+const scalarMultiplyWNAF = (k, P0, window = 5) => {
+    const key = `${window}:${P0.x.toString(16)}:${P0.y.toString(16)}`;
+    let tbl = WNAF_TABLE_CACHE.get(key);
+    let P;
+    if (tbl === undefined) {
+        // Convert affine to Jacobian and pre-compute odd multiples
+        const tblSize = 1 << (window - 1); // e.g. w=5 → 16 entries
+        tbl = new Array(tblSize);
+        P = { X: P0.x, Y: P0.y, Z: exports.BI_ONE };
+        tbl[0] = P;
+        const twoP = (0, exports.jpDouble)(P);
+        for (let i = 1; i < tblSize; i++) {
+            tbl[i] = (0, exports.jpAdd)(tbl[i - 1], twoP);
+        }
+        WNAF_TABLE_CACHE.set(key, tbl);
+    }
+    else {
+        P = tbl[0];
+    }
+    // Build wNAF representation of k
+    const wnaf = [];
+    const wBig = 1n << BigInt(window);
+    const wHalf = wBig >> 1n;
+    let kTmp = k;
+    while (kTmp > 0n) {
+        if ((kTmp & exports.BI_ONE) === exports.BI_ZERO) {
+            wnaf.push(0);
+            kTmp >>= exports.BI_ONE;
+        }
+        else {
+            let z = kTmp & (wBig - 1n); // kTmp mod 2^w
+            if (z > wHalf)
+                z -= wBig; // make it odd & within (-2^{w-1}, 2^{w-1})
+            wnaf.push(Number(z));
+            kTmp -= z;
+            kTmp >>= exports.BI_ONE;
+        }
+    }
+    // Accumulate from MSB to LSB
+    let Q = { X: exports.BI_ZERO, Y: exports.BI_ONE, Z: exports.BI_ZERO }; // infinity
+    for (let i = wnaf.length - 1; i >= 0; i--) {
+        Q = (0, exports.jpDouble)(Q);
+        const di = wnaf[i];
+        if (di !== 0) {
+            const idx = Math.abs(di) >> 1; // (|di|-1)/2  because di is odd
+            const addend = di > 0 ? tbl[idx] : (0, exports.jpNeg)(tbl[idx]);
+            Q = (0, exports.jpAdd)(Q, addend);
+        }
+    }
+    return Q;
+};
+exports.scalarMultiplyWNAF = scalarMultiplyWNAF;
+const modN = (a) => {
+    let r = a % exports.N_BIGINT;
+    if (r < 0n)
+        r += exports.N_BIGINT;
+    return r;
+};
+exports.modN = modN;
+const modMulN = (a, b) => (0, exports.modN)(a * b);
+exports.modMulN = modMulN;
+/** modular inverse modulo n with plain extended‑gcd (not constant‑time) */
+const modInvN = (a) => {
+    let lm = 1n;
+    let hm = 0n;
+    let low = (0, exports.modN)(a);
+    let high = exports.N_BIGINT;
+    while (low > 1n) {
+        const q = high / low;
+        [lm, hm] = [hm - lm * q, lm];
+        [low, high] = [high - low * q, low];
+    }
+    return (0, exports.modN)(lm);
+};
+exports.modInvN = modInvN;
 /**
  * `Point` class is a representation of an elliptic curve point with affine coordinates.
  * It extends the functionality of BasePoint and carries x, y coordinates of point on the curve.
@@ -82,12 +287,6 @@ class Point extends BasePoint_js_1.default {
         const bytes = (0, utils_js_1.toArray)(str, 'hex');
         return Point.fromDER(bytes);
     }
-    static redSqrtOptimized(y2) {
-        const red = Point.red;
-        const p = red.m; // The modulus
-        const exponent = p.addn(1).iushrn(2); // (p + 1) / 4
-        return y2.redPow(exponent);
-    }
     /**
      * Generates a point from an x coordinate and a boolean indicating whether the corresponding
      * y coordinate is odd.
@@ -104,62 +303,18 @@ class Point extends BasePoint_js_1.default {
      * const point = Point.fromX(xCoordinate, true);
      */
     static fromX(x, odd) {
-        function mod(a, n) {
-            return ((a % n) + n) % n;
-        }
-        function modPow(base, exponent, modulus) {
-            let result = BigInt(1);
-            base = mod(base, modulus);
-            while (exponent > BigInt(0)) {
-                if ((exponent & BigInt(1)) === BigInt(1)) {
-                    result = mod(result * base, modulus);
-                }
-                exponent >>= BigInt(1);
-                base = mod(base * base, modulus);
-            }
-            return result;
-        }
-        function sqrtMod(a, p) {
-            const exponent = (p + BigInt(1)) >> BigInt(2);
-            const sqrtCandidate = modPow(a, exponent, p);
-            if (mod(sqrtCandidate * sqrtCandidate, p) === mod(a, p)) {
-                return sqrtCandidate;
-            }
-            else {
-                return null;
-            }
-        }
-        // Curve parameters for secp256k1
-        const p = BigInt('0xFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFC2F');
-        const b = BigInt(7);
-        let xBigInt;
-        if (x instanceof BigNumber_js_1.default) {
-            xBigInt = BigInt('0x' + x.toString(16));
-        }
-        else if (typeof x === 'string') {
-            xBigInt = BigInt('0x' + x);
-        }
-        else if (Array.isArray(x)) {
-            xBigInt = BigInt('0x' + (0, utils_js_1.toHex)(x).padStart(64, '0'));
-        }
-        else if (typeof x === 'number') {
-            xBigInt = BigInt(x);
-        }
-        else {
-            throw new Error('Invalid x-coordinate type');
-        }
-        xBigInt = mod(xBigInt, p);
-        const y2 = mod(modPow(xBigInt, BigInt(3), p) + b, p);
-        let y = sqrtMod(y2, p);
-        if (y === null) {
+        let xBigInt = toBigInt(x);
+        xBigInt = (0, exports.biMod)(xBigInt);
+        const y2 = (0, exports.biModAdd)((0, exports.biModMul)((0, exports.biModSqr)(xBigInt), xBigInt), 7n);
+        const y = (0, exports.biModSqrt)(y2);
+        if (y === null)
             throw new Error('Invalid point');
-        }
-        const isYOdd = y % BigInt(2) === BigInt(1);
-        if ((odd && !isYOdd) || (!odd && isYOdd)) {
-            y = p - y;
+        let yBig = y;
+        if ((yBig & exports.BI_ONE) !== (odd ? exports.BI_ONE : exports.BI_ZERO)) {
+            yBig = (0, exports.biModSub)(exports.P_BIGINT, yBig);
         }
         const xBN = new BigNumber_js_1.default(xBigInt.toString(16), 16);
-        const yBN = new BigNumber_js_1.default(y.toString(16), 16);
+        const yBN = new BigNumber_js_1.default(yBig.toString(16), 16);
         return new Point(xBN, yBN);
     }
     /**
@@ -401,21 +556,30 @@ class Point extends BasePoint_js_1.default {
         }
         // P + (-P) = O
         if (this.neg().eq(p)) {
-            return new Point(new BigNumber_js_1.default(0), new BigNumber_js_1.default(0));
+            return new Point(null, null);
         }
         // P + Q = O
         if (this.x?.cmp(p.x ?? new BigNumber_js_1.default(0)) === 0) {
-            return new Point(new BigNumber_js_1.default(0), new BigNumber_js_1.default(0));
-        }
-        let c = this.y?.redSub(p.y ?? new BigNumber_js_1.default(0)) ?? new BigNumber_js_1.default(0);
-        if (c.cmpn(0) !== 0) {
-            c = c.redMul(this.x?.redSub(p.x ?? new BigNumber_js_1.default(0)).redInvm() ?? new BigNumber_js_1.default(1));
+            return new Point(null, null);
         }
-        const nx = c?.redSqr().redISub(this.x ?? new BigNumber_js_1.default(0)).redISub(p.x ?? new BigNumber_js_1.default(0));
-        const ny = (c ?? new BigNumber_js_1.default(1))
-            .redMul((this.x ?? new BigNumber_js_1.default(0)).redSub(nx ?? new BigNumber_js_1.default(0)))
-            .redISub(this.y ?? new BigNumber_js_1.default(0));
-        return new Point(nx ?? new BigNumber_js_1.default(0), ny ?? new BigNumber_js_1.default(0));
+        const P1 = {
+            X: BigInt('0x' + this.x.fromRed().toString(16)),
+            Y: BigInt('0x' + this.y.fromRed().toString(16)),
+            Z: exports.BI_ONE
+        };
+        const Q1 = {
+            X: BigInt('0x' + p.x.fromRed().toString(16)),
+            Y: BigInt('0x' + p.y.fromRed().toString(16)),
+            Z: exports.BI_ONE
+        };
+        const R = (0, exports.jpAdd)(P1, Q1);
+        if (R.Z === exports.BI_ZERO)
+            return new Point(null, null);
+        const zInv = (0, exports.biModInv)(R.Z);
+        const zInv2 = (0, exports.biModMul)(zInv, zInv);
+        const xRes = (0, exports.biModMul)(R.X, zInv2);
+        const yRes = (0, exports.biModMul)(R.Y, (0, exports.biModMul)(zInv2, zInv));
+        return new Point(xRes.toString(16), yRes.toString(16));
     }
     /**
      * Doubles the current point.
@@ -427,21 +591,21 @@ class Point extends BasePoint_js_1.default {
      * const result = P.dbl();
      * */
     dbl() {
-        if (this.inf) {
+        if (this.inf)
             return this;
+        if (this.x === null || this.y === null) {
+            throw new Error('Point coordinates cannot be null');
         }
-        // 2P = O
-        const ys1 = (this.y ?? new BigNumber_js_1.default(0)).redAdd(this.y ?? new BigNumber_js_1.default(0));
-        if (ys1.cmpn(0) === 0) {
-            return new Point(new BigNumber_js_1.default(0), new BigNumber_js_1.default(0));
-        }
-        const a = this.curve.a;
-        const x2 = (this.x ?? new BigNumber_js_1.default(0)).redSqr();
-        const dyinv = ys1.redInvm();
-        const c = x2.redAdd(x2).redIAdd(x2).redIAdd(a).redMul(dyinv);
-        const nx = c.redSqr().redISub((this.x ?? new BigNumber_js_1.default(0)).redAdd(this.x ?? new BigNumber_js_1.default(0)));
-        const ny = c.redMul((this.x ?? new BigNumber_js_1.default(0)).redSub(nx)).redISub(this.y ?? new BigNumber_js_1.default(0));
-        return new Point(nx, ny);
+        const X = BigInt('0x' + this.x.fromRed().toString(16));
+        const Y = BigInt('0x' + this.y.fromRed().toString(16));
+        if (Y === exports.BI_ZERO)
+            return new Point(null, null);
+        const R = (0, exports.jpDouble)({ X, Y, Z: exports.BI_ONE });
+        const zInv = (0, exports.biModInv)(R.Z);
+        const zInv2 = (0, exports.biModMul)(zInv, zInv);
+        const xRes = (0, exports.biModMul)(R.X, zInv2);
+        const yRes = (0, exports.biModMul)(R.Y, (0, exports.biModMul)(zInv2, zInv));
+        return new Point(xRes.toString(16), yRes.toString(16));
     }
     /**
      * Returns X coordinate of point
@@ -479,86 +643,45 @@ class Point extends BasePoint_js_1.default {
             k = new BigNumber_js_1.default(k, 16);
         }
         k = k;
-        if (typeof BigInt === 'function') {
-            if (this.inf) {
-                return this;
-            }
-            const zero = 0n;
-            const one = 1n;
-            const p = BigInt('0xFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFC2F');
-            const n = BigInt('0xFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEBAAEDCE6AF48A03BBFD25E8CD0364141');
-            let kBig = BigInt('0x' + k.toString(16));
-            const isNeg = kBig < zero;
-            if (isNeg)
-                kBig = -kBig;
-            kBig = ((kBig % n) + n) % n;
-            if (kBig === zero) {
-                return new Point(null, null);
-            }
-            if (this.x === null || this.y === null) {
-                throw new Error('Point coordinates cannot be null');
-            }
-            let Px;
-            let Py;
-            if (this === this.curve.g) {
-                Px = GX_BIGINT;
-                Py = GY_BIGINT;
-            }
-            else {
-                Px = BigInt('0x' + this.x.fromRed().toString(16));
-                Py = BigInt('0x' + this.y.fromRed().toString(16));
-            }
-            const mod = (a, m) => ((a % m) + m) % m;
-            const modMul = (a, b, m) => mod(a * b, m);
-            const modInv = (a, m) => {
-                let lm = one;
-                let hm = zero;
-                let low = mod(a, m);
-                let high = m;
-                while (low > one) {
-                    const r = high / low;
-                    const nm = hm - lm * r;
-                    const neww = high - low * r;
-                    hm = lm;
-                    lm = nm;
-                    high = low;
-                    low = neww;
-                }
-                return mod(lm, m);
-            };
-            const scalarMultiply = (kVal, P0) => {
-                // Delegate to the hoisted windowed-NAF implementation above.  We
-                // keep the wrapper so that the rest of the mul() code remains
-                // untouched while providing a massive speed-up (≈4-6×).
-                return scalarMultiplyWNAF(kVal, P0);
-            };
-            const R = scalarMultiply(kBig, { x: Px, y: Py });
-            if (R.Z === zero) {
-                return new Point(null, null);
-            }
-            const zInv = modInv(R.Z, p);
-            const zInv2 = modMul(zInv, zInv, p);
-            const xRes = modMul(R.X, zInv2, p);
-            const yRes = modMul(R.Y, modMul(zInv2, zInv, p), p);
-            const xBN = new BigNumber_js_1.default(xRes.toString(16), 16);
-            const yBN = new BigNumber_js_1.default(yRes.toString(16), 16);
-            const result = new Point(xBN, yBN);
-            if (isNeg) {
-                return result.neg();
-            }
-            return result;
+        if (this.inf) {
+            return this;
+        }
+        let kBig = BigInt('0x' + k.toString(16));
+        const isNeg = kBig < exports.BI_ZERO;
+        if (isNeg)
+            kBig = -kBig;
+        kBig = (0, exports.biMod)(kBig);
+        if (kBig === exports.BI_ZERO) {
+            return new Point(null, null);
+        }
+        if (this.x === null || this.y === null) {
+            throw new Error('Point coordinates cannot be null');
+        }
+        let Px;
+        let Py;
+        if (this === this.curve.g) {
+            Px = exports.GX_BIGINT;
+            Py = exports.GY_BIGINT;
         }
         else {
-            if (this.isInfinity()) {
-                return this;
-            }
-            else if (this._hasDoubles(k)) {
-                return this._fixedNafMul(k);
-            }
-            else {
-                return this._endoWnafMulAdd([this], [k]);
-            }
+            Px = BigInt('0x' + this.x.fromRed().toString(16));
+            Py = BigInt('0x' + this.y.fromRed().toString(16));
+        }
+        const R = (0, exports.scalarMultiplyWNAF)(kBig, { x: Px, y: Py });
+        if (R.Z === exports.BI_ZERO) {
+            return new Point(null, null);
         }
+        const zInv = (0, exports.biModInv)(R.Z);
+        const zInv2 = (0, exports.biModMul)(zInv, zInv);
+        const xRes = (0, exports.biModMul)(R.X, zInv2);
+        const yRes = (0, exports.biModMul)(R.Y, (0, exports.biModMul)(zInv2, zInv));
+        const xBN = new BigNumber_js_1.default(xRes.toString(16), 16);
+        const yBN = new BigNumber_js_1.default(yRes.toString(16), 16);
+        const result = new Point(xBN, yBN);
+        if (isNeg) {
+            return result.neg();
+        }
+        return result;
     }
     /**
      * Performs a multiplication and addition operation in a single step.
@@ -892,7 +1015,7 @@ class Point extends BasePoint_js_1.default {
         for (i = 0; i < points.length; i++) {
             const split = this.curve._endoSplit(coeffs[i]);
             let p = points[i];
-            let beta = p._getBeta() ?? new Point(new BigNumber_js_1.default(0), new BigNumber_js_1.default(0));
+            let beta = p._getBeta() ?? new Point(null, null);
             if (split.k1.negative !== 0) {
                 split.k1.ineg();
                 p = p.neg(true);
@@ -966,136 +1089,5 @@ class Point extends BasePoint_js_1.default {
         };
     }
 }
-Point.red = new ReductionContext_js_1.default('k256');
-Point.a = new BigNumber_js_1.default(0).toRed(Point.red);
-Point.b = new BigNumber_js_1.default(7).toRed(Point.red);
-Point.zero = new BigNumber_js_1.default(0).toRed(Point.red);
 exports.default = Point;
-// -----------------------------------------------------------------------------
-// BigInt helpers & constants (secp256k1) – hoisted so we don't recreate them on
-// every Point.mul() call.
-// -----------------------------------------------------------------------------
-const BI_ZERO = 0n;
-const BI_ONE = 1n;
-const BI_TWO = 2n;
-const BI_THREE = 3n;
-const BI_FOUR = 4n;
-const BI_EIGHT = 8n;
-const P_BIGINT = 0xfffffffffffffffffffffffffffffffffffffffffffffffffffffffefffffc2fn;
-const MASK_256 = (1n << 256n) - 1n; // 0xffff…ffff (256 sones)
-function red(x) {
-    // first fold
-    let hi = x >> 256n;
-    x = (x & MASK_256) + (hi << 32n) + hi * 977n;
-    // second fold  (hi ≤ 2³² + 977 here, so one more pass is enough)
-    hi = x >> 256n;
-    x = (x & MASK_256) + (hi << 32n) + hi * 977n;
-    // final conditional subtraction
-    if (x >= P_BIGINT)
-        x -= P_BIGINT;
-    return x;
-}
-const biModSub = (a, b) => (a >= b ? a - b : P_BIGINT - (b - a));
-const biModMul = (a, b) => red(a * b);
-// Generator point coordinates as bigint constants
-const GX_BIGINT = BigInt('0x79be667ef9dcbbac55a06295ce870b07029bfcdb2dce28d959f2815b16f81798');
-const GY_BIGINT = BigInt('0x483ada7726a3c4655da4fbfc0e1108a8fd17b448a68554199c47d08ffb10d4b8');
-// Cache for precomputed windowed tables keyed by 'window:x:y'
-const WNAF_TABLE_CACHE = new Map();
-const jpDouble = (P) => {
-    const { X: X1, Y: Y1, Z: Z1 } = P;
-    if (Y1 === BI_ZERO)
-        return { X: BI_ZERO, Y: BI_ONE, Z: BI_ZERO };
-    const Y1sq = biModMul(Y1, Y1);
-    const S = biModMul(BI_FOUR, biModMul(X1, Y1sq));
-    const M = biModMul(BI_THREE, biModMul(X1, X1));
-    const X3 = biModSub(biModMul(M, M), biModMul(BI_TWO, S));
-    const Y3 = biModSub(biModMul(M, biModSub(S, X3)), biModMul(BI_EIGHT, biModMul(Y1sq, Y1sq)));
-    const Z3 = biModMul(BI_TWO, biModMul(Y1, Z1));
-    return { X: X3, Y: Y3, Z: Z3 };
-};
-const jpAdd = (P, Q) => {
-    if (P.Z === BI_ZERO)
-        return Q;
-    if (Q.Z === BI_ZERO)
-        return P;
-    const Z1Z1 = biModMul(P.Z, P.Z);
-    const Z2Z2 = biModMul(Q.Z, Q.Z);
-    const U1 = biModMul(P.X, Z2Z2);
-    const U2 = biModMul(Q.X, Z1Z1);
-    const S1 = biModMul(P.Y, biModMul(Z2Z2, Q.Z));
-    const S2 = biModMul(Q.Y, biModMul(Z1Z1, P.Z));
-    const H = biModSub(U2, U1);
-    const r = biModSub(S2, S1);
-    if (H === BI_ZERO) {
-        if (r === BI_ZERO)
-            return jpDouble(P);
-        return { X: BI_ZERO, Y: BI_ONE, Z: BI_ZERO }; // Infinity
-    }
-    const HH = biModMul(H, H);
-    const HHH = biModMul(H, HH);
-    const V = biModMul(U1, HH);
-    const X3 = biModSub(biModSub(biModMul(r, r), HHH), biModMul(BI_TWO, V));
-    const Y3 = biModSub(biModMul(r, biModSub(V, X3)), biModMul(S1, HHH));
-    const Z3 = biModMul(H, biModMul(P.Z, Q.Z));
-    return { X: X3, Y: Y3, Z: Z3 };
-};
-const jpNeg = (P) => {
-    if (P.Z === BI_ZERO)
-        return P;
-    return { X: P.X, Y: P_BIGINT - P.Y, Z: P.Z };
-};
-// Fast windowed-NAF scalar multiplication (default window = 5) in Jacobian
-// coordinates.  Returns Q = k * P0 as a JacobianPoint.
-const scalarMultiplyWNAF = (k, P0, window = 5) => {
-    const key = `${window}:${P0.x.toString(16)}:${P0.y.toString(16)}`;
-    let tbl = WNAF_TABLE_CACHE.get(key);
-    let P;
-    if (tbl === undefined) {
-        // Convert affine to Jacobian and pre-compute odd multiples
-        const tblSize = 1 << (window - 1); // e.g. w=5 → 16 entries
-        tbl = new Array(tblSize);
-        P = { X: P0.x, Y: P0.y, Z: BI_ONE };
-        tbl[0] = P;
-        const twoP = jpDouble(P);
-        for (let i = 1; i < tblSize; i++) {
-            tbl[i] = jpAdd(tbl[i - 1], twoP);
-        }
-        WNAF_TABLE_CACHE.set(key, tbl);
-    }
-    else {
-        P = tbl[0];
-    }
-    // Build wNAF representation of k
-    const wnaf = [];
-    const wBig = 1n << BigInt(window);
-    const wHalf = wBig >> 1n;
-    let kTmp = k;
-    while (kTmp > 0n) {
-        if ((kTmp & BI_ONE) === BI_ZERO) {
-            wnaf.push(0);
-            kTmp >>= BI_ONE;
-        }
-        else {
-            let z = kTmp & (wBig - 1n); // kTmp mod 2^w
-            if (z > wHalf)
-                z -= wBig; // make it odd & within (-2^{w-1}, 2^{w-1})
-            wnaf.push(Number(z));
-            kTmp -= z;
-            kTmp >>= BI_ONE;
-        }
-    }
-    // Accumulate from MSB to LSB
-    let Q = { X: BI_ZERO, Y: BI_ONE, Z: BI_ZERO }; // infinity
-    for (let i = wnaf.length - 1; i >= 0; i--) {
-        Q = jpDouble(Q);
-        const di = wnaf[i];
-        if (di !== 0) {
-            const idx = Math.abs(di) >> 1; // (|di|-1)/2  because di is odd
-            const addend = di > 0 ? tbl[idx] : jpNeg(tbl[idx]);
-            Q = jpAdd(Q, addend);
-        }
-    }
-    return Q;
-};
 //# sourceMappingURL=Point.js.map