You can not select more than 25 topics
Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
355 lines
15 KiB
355 lines
15 KiB
"use strict";
|
|
Object.defineProperty(exports, "__esModule", { value: true });
|
|
exports.bls = bls;
|
|
/**
|
|
* BLS (Barreto-Lynn-Scott) family of pairing-friendly curves.
|
|
* BLS != BLS.
|
|
* The file implements BLS (Boneh-Lynn-Shacham) signatures.
|
|
* Used in both BLS (Barreto-Lynn-Scott) and BN (Barreto-Naehrig)
|
|
* families of pairing-friendly curves.
|
|
* Consists of two curves: G1 and G2:
|
|
* - G1 is a subgroup of (x, y) E(Fq) over y² = x³ + 4.
|
|
* - G2 is a subgroup of ((x₁, x₂+i), (y₁, y₂+i)) E(Fq²) over y² = x³ + 4(1 + i) where i is √-1
|
|
* - Gt, created by bilinear (ate) pairing e(G1, G2), consists of p-th roots of unity in
|
|
* Fq^k where k is embedding degree. Only degree 12 is currently supported, 24 is not.
|
|
* Pairing is used to aggregate and verify signatures.
|
|
* There are two main ways to use it:
|
|
* 1. Fp for short private keys, Fp₂ for signatures
|
|
* 2. Fp for short signatures, Fp₂ for private keys
|
|
* @module
|
|
**/
|
|
/*! noble-curves - MIT License (c) 2022 Paul Miller (paulmillr.com) */
|
|
// TODO: import { AffinePoint } from './curve.ts';
|
|
const hash_to_curve_ts_1 = require("./hash-to-curve.js");
|
|
const modular_ts_1 = require("./modular.js");
|
|
const utils_ts_1 = require("./utils.js");
|
|
const weierstrass_ts_1 = require("./weierstrass.js");
|
|
// prettier-ignore
|
|
const _0n = BigInt(0), _1n = BigInt(1), _2n = BigInt(2), _3n = BigInt(3);
|
|
// Not used with BLS12-381 (no sequential `11` in X). Useful for other curves.
|
|
function NAfDecomposition(a) {
|
|
const res = [];
|
|
// a>1 because of marker bit
|
|
for (; a > _1n; a >>= _1n) {
|
|
if ((a & _1n) === _0n)
|
|
res.unshift(0);
|
|
else if ((a & _3n) === _3n) {
|
|
res.unshift(-1);
|
|
a += _1n;
|
|
}
|
|
else
|
|
res.unshift(1);
|
|
}
|
|
return res;
|
|
}
|
|
function bls(CURVE) {
|
|
// Fields are specific for curve, so for now we'll need to pass them with opts
|
|
const { Fp, Fr, Fp2, Fp6, Fp12 } = CURVE.fields;
|
|
const BLS_X_IS_NEGATIVE = CURVE.params.xNegative;
|
|
const TWIST = CURVE.params.twistType;
|
|
// Point on G1 curve: (x, y)
|
|
const G1_ = (0, weierstrass_ts_1.weierstrassPoints)({ n: Fr.ORDER, ...CURVE.G1 });
|
|
const G1 = Object.assign(G1_, (0, hash_to_curve_ts_1.createHasher)(G1_.ProjectivePoint, CURVE.G1.mapToCurve, {
|
|
...CURVE.htfDefaults,
|
|
...CURVE.G1.htfDefaults,
|
|
}));
|
|
// Point on G2 curve (complex numbers): (x₁, x₂+i), (y₁, y₂+i)
|
|
const G2_ = (0, weierstrass_ts_1.weierstrassPoints)({ n: Fr.ORDER, ...CURVE.G2 });
|
|
const G2 = Object.assign(G2_, (0, hash_to_curve_ts_1.createHasher)(G2_.ProjectivePoint, CURVE.G2.mapToCurve, {
|
|
...CURVE.htfDefaults,
|
|
...CURVE.G2.htfDefaults,
|
|
}));
|
|
// Applies sparse multiplication as line function
|
|
let lineFunction;
|
|
if (TWIST === 'multiplicative') {
|
|
lineFunction = (c0, c1, c2, f, Px, Py) => Fp12.mul014(f, c0, Fp2.mul(c1, Px), Fp2.mul(c2, Py));
|
|
}
|
|
else if (TWIST === 'divisive') {
|
|
// NOTE: it should be [c0, c1, c2], but we use different order here to reduce complexity of
|
|
// precompute calculations.
|
|
lineFunction = (c0, c1, c2, f, Px, Py) => Fp12.mul034(f, Fp2.mul(c2, Py), Fp2.mul(c1, Px), c0);
|
|
}
|
|
else
|
|
throw new Error('bls: unknown twist type');
|
|
const Fp2div2 = Fp2.div(Fp2.ONE, Fp2.mul(Fp2.ONE, _2n));
|
|
function pointDouble(ell, Rx, Ry, Rz) {
|
|
const t0 = Fp2.sqr(Ry); // Ry²
|
|
const t1 = Fp2.sqr(Rz); // Rz²
|
|
const t2 = Fp2.mulByB(Fp2.mul(t1, _3n)); // 3 * T1 * B
|
|
const t3 = Fp2.mul(t2, _3n); // 3 * T2
|
|
const t4 = Fp2.sub(Fp2.sub(Fp2.sqr(Fp2.add(Ry, Rz)), t1), t0); // (Ry + Rz)² - T1 - T0
|
|
const c0 = Fp2.sub(t2, t0); // T2 - T0 (i)
|
|
const c1 = Fp2.mul(Fp2.sqr(Rx), _3n); // 3 * Rx²
|
|
const c2 = Fp2.neg(t4); // -T4 (-h)
|
|
ell.push([c0, c1, c2]);
|
|
Rx = Fp2.mul(Fp2.mul(Fp2.mul(Fp2.sub(t0, t3), Rx), Ry), Fp2div2); // ((T0 - T3) * Rx * Ry) / 2
|
|
Ry = Fp2.sub(Fp2.sqr(Fp2.mul(Fp2.add(t0, t3), Fp2div2)), Fp2.mul(Fp2.sqr(t2), _3n)); // ((T0 + T3) / 2)² - 3 * T2²
|
|
Rz = Fp2.mul(t0, t4); // T0 * T4
|
|
return { Rx, Ry, Rz };
|
|
}
|
|
function pointAdd(ell, Rx, Ry, Rz, Qx, Qy) {
|
|
// Addition
|
|
const t0 = Fp2.sub(Ry, Fp2.mul(Qy, Rz)); // Ry - Qy * Rz
|
|
const t1 = Fp2.sub(Rx, Fp2.mul(Qx, Rz)); // Rx - Qx * Rz
|
|
const c0 = Fp2.sub(Fp2.mul(t0, Qx), Fp2.mul(t1, Qy)); // T0 * Qx - T1 * Qy == Ry * Qx - Rx * Qy
|
|
const c1 = Fp2.neg(t0); // -T0 == Qy * Rz - Ry
|
|
const c2 = t1; // == Rx - Qx * Rz
|
|
ell.push([c0, c1, c2]);
|
|
const t2 = Fp2.sqr(t1); // T1²
|
|
const t3 = Fp2.mul(t2, t1); // T2 * T1
|
|
const t4 = Fp2.mul(t2, Rx); // T2 * Rx
|
|
const t5 = Fp2.add(Fp2.sub(t3, Fp2.mul(t4, _2n)), Fp2.mul(Fp2.sqr(t0), Rz)); // T3 - 2 * T4 + T0² * Rz
|
|
Rx = Fp2.mul(t1, t5); // T1 * T5
|
|
Ry = Fp2.sub(Fp2.mul(Fp2.sub(t4, t5), t0), Fp2.mul(t3, Ry)); // (T4 - T5) * T0 - T3 * Ry
|
|
Rz = Fp2.mul(Rz, t3); // Rz * T3
|
|
return { Rx, Ry, Rz };
|
|
}
|
|
// Pre-compute coefficients for sparse multiplication
|
|
// Point addition and point double calculations is reused for coefficients
|
|
// pointAdd happens only if bit set, so wNAF is reasonable. Unfortunately we cannot combine
|
|
// add + double in windowed precomputes here, otherwise it would be single op (since X is static)
|
|
const ATE_NAF = NAfDecomposition(CURVE.params.ateLoopSize);
|
|
const calcPairingPrecomputes = (0, utils_ts_1.memoized)((point) => {
|
|
const p = point;
|
|
const { x, y } = p.toAffine();
|
|
// prettier-ignore
|
|
const Qx = x, Qy = y, negQy = Fp2.neg(y);
|
|
// prettier-ignore
|
|
let Rx = Qx, Ry = Qy, Rz = Fp2.ONE;
|
|
const ell = [];
|
|
for (const bit of ATE_NAF) {
|
|
const cur = [];
|
|
({ Rx, Ry, Rz } = pointDouble(cur, Rx, Ry, Rz));
|
|
if (bit)
|
|
({ Rx, Ry, Rz } = pointAdd(cur, Rx, Ry, Rz, Qx, bit === -1 ? negQy : Qy));
|
|
ell.push(cur);
|
|
}
|
|
if (CURVE.postPrecompute) {
|
|
const last = ell[ell.length - 1];
|
|
CURVE.postPrecompute(Rx, Ry, Rz, Qx, Qy, pointAdd.bind(null, last));
|
|
}
|
|
return ell;
|
|
});
|
|
function millerLoopBatch(pairs, withFinalExponent = false) {
|
|
let f12 = Fp12.ONE;
|
|
if (pairs.length) {
|
|
const ellLen = pairs[0][0].length;
|
|
for (let i = 0; i < ellLen; i++) {
|
|
f12 = Fp12.sqr(f12); // This allows us to do sqr only one time for all pairings
|
|
// NOTE: we apply multiple pairings in parallel here
|
|
for (const [ell, Px, Py] of pairs) {
|
|
for (const [c0, c1, c2] of ell[i])
|
|
f12 = lineFunction(c0, c1, c2, f12, Px, Py);
|
|
}
|
|
}
|
|
}
|
|
if (BLS_X_IS_NEGATIVE)
|
|
f12 = Fp12.conjugate(f12);
|
|
return withFinalExponent ? Fp12.finalExponentiate(f12) : f12;
|
|
}
|
|
// Calculates product of multiple pairings
|
|
// This up to x2 faster than just `map(({g1, g2})=>pairing({g1,g2}))`
|
|
function pairingBatch(pairs, withFinalExponent = true) {
|
|
const res = [];
|
|
// Cache precomputed toAffine for all points
|
|
G1.ProjectivePoint.normalizeZ(pairs.map(({ g1 }) => g1));
|
|
G2.ProjectivePoint.normalizeZ(pairs.map(({ g2 }) => g2));
|
|
for (const { g1, g2 } of pairs) {
|
|
if (g1.equals(G1.ProjectivePoint.ZERO) || g2.equals(G2.ProjectivePoint.ZERO))
|
|
throw new Error('pairing is not available for ZERO point');
|
|
// This uses toAffine inside
|
|
g1.assertValidity();
|
|
g2.assertValidity();
|
|
const Qa = g1.toAffine();
|
|
res.push([calcPairingPrecomputes(g2), Qa.x, Qa.y]);
|
|
}
|
|
return millerLoopBatch(res, withFinalExponent);
|
|
}
|
|
// Calculates bilinear pairing
|
|
function pairing(Q, P, withFinalExponent = true) {
|
|
return pairingBatch([{ g1: Q, g2: P }], withFinalExponent);
|
|
}
|
|
const utils = {
|
|
randomPrivateKey: () => {
|
|
const length = (0, modular_ts_1.getMinHashLength)(Fr.ORDER);
|
|
return (0, modular_ts_1.mapHashToField)(CURVE.randomBytes(length), Fr.ORDER);
|
|
},
|
|
calcPairingPrecomputes,
|
|
};
|
|
const { ShortSignature } = CURVE.G1;
|
|
const { Signature } = CURVE.G2;
|
|
function normP1(point) {
|
|
return point instanceof G1.ProjectivePoint ? point : G1.ProjectivePoint.fromHex(point);
|
|
}
|
|
function normP1Hash(point, htfOpts) {
|
|
return point instanceof G1.ProjectivePoint
|
|
? point
|
|
: G1.hashToCurve((0, utils_ts_1.ensureBytes)('point', point), htfOpts);
|
|
}
|
|
function normP2(point) {
|
|
return point instanceof G2.ProjectivePoint ? point : Signature.fromHex(point);
|
|
}
|
|
function normP2Hash(point, htfOpts) {
|
|
return point instanceof G2.ProjectivePoint
|
|
? point
|
|
: G2.hashToCurve((0, utils_ts_1.ensureBytes)('point', point), htfOpts);
|
|
}
|
|
// Multiplies generator (G1) by private key.
|
|
// P = pk x G
|
|
function getPublicKey(privateKey) {
|
|
return G1.ProjectivePoint.fromPrivateKey(privateKey).toRawBytes(true);
|
|
}
|
|
// Multiplies generator (G2) by private key.
|
|
// P = pk x G
|
|
function getPublicKeyForShortSignatures(privateKey) {
|
|
return G2.ProjectivePoint.fromPrivateKey(privateKey).toRawBytes(true);
|
|
}
|
|
function sign(message, privateKey, htfOpts) {
|
|
const msgPoint = normP2Hash(message, htfOpts);
|
|
msgPoint.assertValidity();
|
|
const sigPoint = msgPoint.multiply(G1.normPrivateKeyToScalar(privateKey));
|
|
if (message instanceof G2.ProjectivePoint)
|
|
return sigPoint;
|
|
return Signature.toRawBytes(sigPoint);
|
|
}
|
|
function signShortSignature(message, privateKey, htfOpts) {
|
|
const msgPoint = normP1Hash(message, htfOpts);
|
|
msgPoint.assertValidity();
|
|
const sigPoint = msgPoint.multiply(G1.normPrivateKeyToScalar(privateKey));
|
|
if (message instanceof G1.ProjectivePoint)
|
|
return sigPoint;
|
|
return ShortSignature.toRawBytes(sigPoint);
|
|
}
|
|
// Checks if pairing of public key & hash is equal to pairing of generator & signature.
|
|
// e(P, H(m)) == e(G, S)
|
|
function verify(signature, message, publicKey, htfOpts) {
|
|
const P = normP1(publicKey);
|
|
const Hm = normP2Hash(message, htfOpts);
|
|
const G = G1.ProjectivePoint.BASE;
|
|
const S = normP2(signature);
|
|
const exp = pairingBatch([
|
|
{ g1: P.negate(), g2: Hm }, // ePHM = pairing(P.negate(), Hm, false);
|
|
{ g1: G, g2: S }, // eGS = pairing(G, S, false);
|
|
]);
|
|
return Fp12.eql(exp, Fp12.ONE);
|
|
}
|
|
// Checks if pairing of public key & hash is equal to pairing of generator & signature.
|
|
// e(S, G) == e(H(m), P)
|
|
function verifyShortSignature(signature, message, publicKey, htfOpts) {
|
|
const P = normP2(publicKey);
|
|
const Hm = normP1Hash(message, htfOpts);
|
|
const G = G2.ProjectivePoint.BASE;
|
|
const S = normP1(signature);
|
|
const exp = pairingBatch([
|
|
{ g1: Hm, g2: P }, // eHmP = pairing(Hm, P, false);
|
|
{ g1: S, g2: G.negate() }, // eSG = pairing(S, G.negate(), false);
|
|
]);
|
|
return Fp12.eql(exp, Fp12.ONE);
|
|
}
|
|
function aNonEmpty(arr) {
|
|
if (!Array.isArray(arr) || arr.length === 0)
|
|
throw new Error('expected non-empty array');
|
|
}
|
|
function aggregatePublicKeys(publicKeys) {
|
|
aNonEmpty(publicKeys);
|
|
const agg = publicKeys.map(normP1).reduce((sum, p) => sum.add(p), G1.ProjectivePoint.ZERO);
|
|
const aggAffine = agg; //.toAffine();
|
|
if (publicKeys[0] instanceof G1.ProjectivePoint) {
|
|
aggAffine.assertValidity();
|
|
return aggAffine;
|
|
}
|
|
// toRawBytes ensures point validity
|
|
return aggAffine.toRawBytes(true);
|
|
}
|
|
function aggregateSignatures(signatures) {
|
|
aNonEmpty(signatures);
|
|
const agg = signatures.map(normP2).reduce((sum, s) => sum.add(s), G2.ProjectivePoint.ZERO);
|
|
const aggAffine = agg; //.toAffine();
|
|
if (signatures[0] instanceof G2.ProjectivePoint) {
|
|
aggAffine.assertValidity();
|
|
return aggAffine;
|
|
}
|
|
return Signature.toRawBytes(aggAffine);
|
|
}
|
|
function aggregateShortSignatures(signatures) {
|
|
aNonEmpty(signatures);
|
|
const agg = signatures.map(normP1).reduce((sum, s) => sum.add(s), G1.ProjectivePoint.ZERO);
|
|
const aggAffine = agg; //.toAffine();
|
|
if (signatures[0] instanceof G1.ProjectivePoint) {
|
|
aggAffine.assertValidity();
|
|
return aggAffine;
|
|
}
|
|
return ShortSignature.toRawBytes(aggAffine);
|
|
}
|
|
// https://ethresear.ch/t/fast-verification-of-multiple-bls-signatures/5407
|
|
// e(G, S) = e(G, SUM(n)(Si)) = MUL(n)(e(G, Si))
|
|
function verifyBatch(signature,
|
|
// TODO: maybe `{message: G2Hex, publicKey: G1Hex}[]` instead?
|
|
messages, publicKeys, htfOpts) {
|
|
aNonEmpty(messages);
|
|
if (publicKeys.length !== messages.length)
|
|
throw new Error('amount of public keys and messages should be equal');
|
|
const sig = normP2(signature);
|
|
const nMessages = messages.map((i) => normP2Hash(i, htfOpts));
|
|
const nPublicKeys = publicKeys.map(normP1);
|
|
// NOTE: this works only for exact same object
|
|
const messagePubKeyMap = new Map();
|
|
for (let i = 0; i < nPublicKeys.length; i++) {
|
|
const pub = nPublicKeys[i];
|
|
const msg = nMessages[i];
|
|
let keys = messagePubKeyMap.get(msg);
|
|
if (keys === undefined) {
|
|
keys = [];
|
|
messagePubKeyMap.set(msg, keys);
|
|
}
|
|
keys.push(pub);
|
|
}
|
|
const paired = [];
|
|
try {
|
|
for (const [msg, keys] of messagePubKeyMap) {
|
|
const groupPublicKey = keys.reduce((acc, msg) => acc.add(msg));
|
|
paired.push({ g1: groupPublicKey, g2: msg });
|
|
}
|
|
paired.push({ g1: G1.ProjectivePoint.BASE.negate(), g2: sig });
|
|
return Fp12.eql(pairingBatch(paired), Fp12.ONE);
|
|
}
|
|
catch {
|
|
return false;
|
|
}
|
|
}
|
|
G1.ProjectivePoint.BASE._setWindowSize(4);
|
|
return {
|
|
getPublicKey,
|
|
getPublicKeyForShortSignatures,
|
|
sign,
|
|
signShortSignature,
|
|
verify,
|
|
verifyBatch,
|
|
verifyShortSignature,
|
|
aggregatePublicKeys,
|
|
aggregateSignatures,
|
|
aggregateShortSignatures,
|
|
millerLoopBatch,
|
|
pairing,
|
|
pairingBatch,
|
|
G1,
|
|
G2,
|
|
Signature,
|
|
ShortSignature,
|
|
fields: {
|
|
Fr,
|
|
Fp,
|
|
Fp2,
|
|
Fp6,
|
|
Fp12,
|
|
},
|
|
params: {
|
|
ateLoopSize: CURVE.params.ateLoopSize,
|
|
r: CURVE.params.r,
|
|
G1b: CURVE.G1.b,
|
|
G2b: CURVE.G2.b,
|
|
},
|
|
utils,
|
|
};
|
|
}
|
|
//# sourceMappingURL=bls.js.map
|