@noble/curves 0.7.2 → 0.8.0
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- package/README.md +32 -13
- package/_shortw_utils.js +8 -13
- package/_shortw_utils.js.map +1 -1
- package/abstract/bls.js +12 -16
- package/abstract/bls.js.map +1 -1
- package/abstract/curve.js +7 -12
- package/abstract/curve.js.map +1 -1
- package/abstract/edwards.js +16 -20
- package/abstract/edwards.js.map +1 -1
- package/abstract/hash-to-curve.d.ts +16 -5
- package/abstract/hash-to-curve.d.ts.map +1 -1
- package/abstract/hash-to-curve.js +51 -56
- package/abstract/hash-to-curve.js.map +1 -1
- package/abstract/modular.d.ts.map +1 -1
- package/abstract/modular.js +25 -44
- package/abstract/modular.js.map +1 -1
- package/abstract/montgomery.js +11 -15
- package/abstract/montgomery.js.map +1 -1
- package/abstract/poseidon.js +6 -12
- package/abstract/poseidon.js.map +1 -1
- package/abstract/utils.js +19 -41
- package/abstract/utils.js.map +1 -1
- package/abstract/weierstrass.d.ts.map +1 -1
- package/abstract/weierstrass.js +25 -37
- package/abstract/weierstrass.js.map +1 -1
- package/bls12-381.js +63 -66
- package/bls12-381.js.map +1 -1
- package/bn.js +7 -10
- package/bn.js.map +1 -1
- package/ed25519.d.ts +3 -0
- package/ed25519.d.ts.map +1 -1
- package/ed25519.js +81 -74
- package/ed25519.js.map +1 -1
- package/ed448.js +37 -41
- package/ed448.js.map +1 -1
- package/jubjub.js +17 -22
- package/jubjub.js.map +1 -1
- package/p256.js +13 -17
- package/p256.js.map +1 -1
- package/p384.js +13 -17
- package/p384.js.map +1 -1
- package/p521.js +13 -17
- package/p521.js.map +1 -1
- package/package.json +4 -26
- package/pasta.js +16 -19
- package/pasta.js.map +1 -1
- package/secp256k1.d.ts +0 -6
- package/secp256k1.d.ts.map +1 -1
- package/secp256k1.js +58 -63
- package/secp256k1.js.map +1 -1
- package/src/abstract/hash-to-curve.ts +49 -47
- package/src/abstract/modular.ts +1 -0
- package/src/abstract/weierstrass.ts +4 -10
- package/src/ed25519.ts +17 -1
- package/src/secp256k1.ts +7 -7
- package/src/stark.ts +50 -27
- package/stark.d.ts +17 -19
- package/stark.d.ts.map +1 -1
- package/stark.js +76 -72
- package/stark.js.map +1 -1
- package/esm/_shortw_utils.js +0 -17
- package/esm/_shortw_utils.js.map +0 -1
- package/esm/abstract/bls.js +0 -226
- package/esm/abstract/bls.js.map +0 -1
- package/esm/abstract/curve.js +0 -152
- package/esm/abstract/curve.js.map +0 -1
- package/esm/abstract/edwards.js +0 -409
- package/esm/abstract/edwards.js.map +0 -1
- package/esm/abstract/hash-to-curve.js +0 -163
- package/esm/abstract/hash-to-curve.js.map +0 -1
- package/esm/abstract/modular.js +0 -345
- package/esm/abstract/modular.js.map +0 -1
- package/esm/abstract/montgomery.js +0 -157
- package/esm/abstract/montgomery.js.map +0 -1
- package/esm/abstract/poseidon.js +0 -110
- package/esm/abstract/poseidon.js.map +0 -1
- package/esm/abstract/utils.js +0 -222
- package/esm/abstract/utils.js.map +0 -1
- package/esm/abstract/weierstrass.js +0 -1016
- package/esm/abstract/weierstrass.js.map +0 -1
- package/esm/bls12-381.js +0 -1173
- package/esm/bls12-381.js.map +0 -1
- package/esm/bn.js +0 -22
- package/esm/bn.js.map +0 -1
- package/esm/ed25519.js +0 -385
- package/esm/ed25519.js.map +0 -1
- package/esm/ed448.js +0 -213
- package/esm/ed448.js.map +0 -1
- package/esm/index.js +0 -3
- package/esm/index.js.map +0 -1
- package/esm/jubjub.js +0 -54
- package/esm/jubjub.js.map +0 -1
- package/esm/p256.js +0 -42
- package/esm/p256.js.map +0 -1
- package/esm/p384.js +0 -47
- package/esm/p384.js.map +0 -1
- package/esm/p521.js +0 -48
- package/esm/p521.js.map +0 -1
- package/esm/package.json +0 -7
- package/esm/pasta.js +0 -30
- package/esm/pasta.js.map +0 -1
- package/esm/secp256k1.js +0 -253
- package/esm/secp256k1.js.map +0 -1
- package/esm/stark.js +0 -251
- package/esm/stark.js.map +0 -1
package/esm/bls12-381.js
DELETED
|
@@ -1,1173 +0,0 @@
|
|
|
1
|
-
/*! noble-curves - MIT License (c) 2022 Paul Miller (paulmillr.com) */
|
|
2
|
-
// bls12-381 pairing-friendly Barreto-Lynn-Scott elliptic curve construction allows to:
|
|
3
|
-
// - Construct zk-SNARKs at the 128-bit security
|
|
4
|
-
// - Use threshold signatures, which allows a user to sign lots of messages with one signature and
|
|
5
|
-
// verify them swiftly in a batch, using Boneh-Lynn-Shacham signature scheme.
|
|
6
|
-
//
|
|
7
|
-
// The library uses G1 for public keys and G2 for signatures. Support for G1 signatures is planned.
|
|
8
|
-
// Compatible with Algorand, Chia, Dfinity, Ethereum, FIL, Zcash. Matches specs
|
|
9
|
-
// [pairing-curves-10](https://tools.ietf.org/html/draft-irtf-cfrg-pairing-friendly-curves-10),
|
|
10
|
-
// [bls-sigs-04](https://tools.ietf.org/html/draft-irtf-cfrg-bls-signature-04),
|
|
11
|
-
// [hash-to-curve-12](https://tools.ietf.org/html/draft-irtf-cfrg-hash-to-curve-12).
|
|
12
|
-
//
|
|
13
|
-
// ### Summary
|
|
14
|
-
// 1. BLS Relies on Bilinear Pairing (expensive)
|
|
15
|
-
// 2. Private Keys: 32 bytes
|
|
16
|
-
// 3. Public Keys: 48 bytes: 381 bit affine x coordinate, encoded into 48 big-endian bytes.
|
|
17
|
-
// 4. Signatures: 96 bytes: two 381 bit integers (affine x coordinate), encoded into two 48 big-endian byte arrays.
|
|
18
|
-
// - The signature is a point on the G2 subgroup, which is defined over a finite field
|
|
19
|
-
// with elements twice as big as the G1 curve (G2 is over Fp2 rather than Fp. Fp2 is analogous to the complex numbers).
|
|
20
|
-
// 5. The 12 stands for the Embedding degree.
|
|
21
|
-
//
|
|
22
|
-
// ### Formulas
|
|
23
|
-
// - `P = pk x G` - public keys
|
|
24
|
-
// - `S = pk x H(m)` - signing
|
|
25
|
-
// - `e(P, H(m)) == e(G, S)` - verification using pairings
|
|
26
|
-
// - `e(G, S) = e(G, SUM(n)(Si)) = MUL(n)(e(G, Si))` - signature aggregation
|
|
27
|
-
// Filecoin uses little endian byte arrays for private keys -
|
|
28
|
-
// so ensure to reverse byte order if you'll use it with FIL.
|
|
29
|
-
//
|
|
30
|
-
// ### Resources
|
|
31
|
-
// - [BLS12-381 for the rest of us](https://hackmd.io/@benjaminion/bls12-381)
|
|
32
|
-
// - [Key concepts of pairings](https://medium.com/@alonmuroch_65570/bls-signatures-part-2-key-concepts-of-pairings-27a8a9533d0c)
|
|
33
|
-
// - Pairing over bls12-381:
|
|
34
|
-
// [part 1](https://research.nccgroup.com/2020/07/06/pairing-over-bls12-381-part-1-fields/),
|
|
35
|
-
// [part 2](https://research.nccgroup.com/2020/07/13/pairing-over-bls12-381-part-2-curves/),
|
|
36
|
-
// [part 3](https://research.nccgroup.com/2020/08/13/pairing-over-bls12-381-part-3-pairing/)
|
|
37
|
-
// - [Estimating the bit security of pairing-friendly curves](https://research.nccgroup.com/2022/02/03/estimating-the-bit-security-of-pairing-friendly-curves/)
|
|
38
|
-
//
|
|
39
|
-
// ### Differences from @noble/bls12-381 1.4
|
|
40
|
-
// - PointG1 -> G1.Point
|
|
41
|
-
// - PointG2 -> G2.Point
|
|
42
|
-
// - PointG2.fromSignature -> Signature.decode
|
|
43
|
-
// - PointG2.toSignature -> Signature.encode
|
|
44
|
-
// - Fixed Fp2 ORDER
|
|
45
|
-
// - Points now have only two coordinates
|
|
46
|
-
import { sha256 } from '@noble/hashes/sha256';
|
|
47
|
-
import { randomBytes } from '@noble/hashes/utils';
|
|
48
|
-
import { bls } from './abstract/bls.js';
|
|
49
|
-
import * as mod from './abstract/modular.js';
|
|
50
|
-
import { concatBytes as concatB, ensureBytes, numberToBytesBE, bytesToNumberBE, bitLen, bitSet, bitGet, bitMask, } from './abstract/utils.js';
|
|
51
|
-
// Types
|
|
52
|
-
import { mapToCurveSimpleSWU, } from './abstract/weierstrass.js';
|
|
53
|
-
import { isogenyMap } from './abstract/hash-to-curve.js';
|
|
54
|
-
// CURVE FIELDS
|
|
55
|
-
// Finite field over p.
|
|
56
|
-
const Fp = mod.Fp(0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaabn);
|
|
57
|
-
// Finite field over r.
|
|
58
|
-
// This particular field is not used anywhere in bls12-381, but it is still useful.
|
|
59
|
-
const Fr = mod.Fp(0x73eda753299d7d483339d80809a1d80553bda402fffe5bfeffffffff00000001n);
|
|
60
|
-
const Fp2Add = ({ c0, c1 }, { c0: r0, c1: r1 }) => ({
|
|
61
|
-
c0: Fp.add(c0, r0),
|
|
62
|
-
c1: Fp.add(c1, r1),
|
|
63
|
-
});
|
|
64
|
-
const Fp2Subtract = ({ c0, c1 }, { c0: r0, c1: r1 }) => ({
|
|
65
|
-
c0: Fp.sub(c0, r0),
|
|
66
|
-
c1: Fp.sub(c1, r1),
|
|
67
|
-
});
|
|
68
|
-
const Fp2Multiply = ({ c0, c1 }, rhs) => {
|
|
69
|
-
if (typeof rhs === 'bigint')
|
|
70
|
-
return { c0: Fp.mul(c0, rhs), c1: Fp.mul(c1, rhs) };
|
|
71
|
-
// (a+bi)(c+di) = (ac−bd) + (ad+bc)i
|
|
72
|
-
const { c0: r0, c1: r1 } = rhs;
|
|
73
|
-
let t1 = Fp.mul(c0, r0); // c0 * o0
|
|
74
|
-
let t2 = Fp.mul(c1, r1); // c1 * o1
|
|
75
|
-
// (T1 - T2) + ((c0 + c1) * (r0 + r1) - (T1 + T2))*i
|
|
76
|
-
const o0 = Fp.sub(t1, t2);
|
|
77
|
-
const o1 = Fp.sub(Fp.mul(Fp.add(c0, c1), Fp.add(r0, r1)), Fp.add(t1, t2));
|
|
78
|
-
return { c0: o0, c1: o1 };
|
|
79
|
-
};
|
|
80
|
-
const Fp2Square = ({ c0, c1 }) => {
|
|
81
|
-
const a = Fp.add(c0, c1);
|
|
82
|
-
const b = Fp.sub(c0, c1);
|
|
83
|
-
const c = Fp.add(c0, c0);
|
|
84
|
-
return { c0: Fp.mul(a, b), c1: Fp.mul(c, c1) };
|
|
85
|
-
};
|
|
86
|
-
// G2 is the order-q subgroup of E2(Fp²) : y² = x³+4(1+√−1),
|
|
87
|
-
// where Fp2 is Fp[√−1]/(x2+1). #E2(Fp2 ) = h2q, where
|
|
88
|
-
// G² - 1
|
|
89
|
-
// h2q
|
|
90
|
-
// NOTE: ORDER was wrong!
|
|
91
|
-
const FP2_ORDER = 0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaabn **
|
|
92
|
-
2n;
|
|
93
|
-
const Fp2 = {
|
|
94
|
-
ORDER: FP2_ORDER,
|
|
95
|
-
BITS: bitLen(FP2_ORDER),
|
|
96
|
-
BYTES: Math.ceil(bitLen(FP2_ORDER) / 8),
|
|
97
|
-
MASK: bitMask(bitLen(FP2_ORDER)),
|
|
98
|
-
ZERO: { c0: Fp.ZERO, c1: Fp.ZERO },
|
|
99
|
-
ONE: { c0: Fp.ONE, c1: Fp.ZERO },
|
|
100
|
-
create: (num) => num,
|
|
101
|
-
isValid: ({ c0, c1 }) => typeof c0 === 'bigint' && typeof c1 === 'bigint',
|
|
102
|
-
is0: ({ c0, c1 }) => Fp.is0(c0) && Fp.is0(c1),
|
|
103
|
-
eql: ({ c0, c1 }, { c0: r0, c1: r1 }) => Fp.eql(c0, r0) && Fp.eql(c1, r1),
|
|
104
|
-
neg: ({ c0, c1 }) => ({ c0: Fp.neg(c0), c1: Fp.neg(c1) }),
|
|
105
|
-
pow: (num, power) => mod.FpPow(Fp2, num, power),
|
|
106
|
-
invertBatch: (nums) => mod.FpInvertBatch(Fp2, nums),
|
|
107
|
-
// Normalized
|
|
108
|
-
add: Fp2Add,
|
|
109
|
-
sub: Fp2Subtract,
|
|
110
|
-
mul: Fp2Multiply,
|
|
111
|
-
sqr: Fp2Square,
|
|
112
|
-
// NonNormalized stuff
|
|
113
|
-
addN: Fp2Add,
|
|
114
|
-
subN: Fp2Subtract,
|
|
115
|
-
mulN: Fp2Multiply,
|
|
116
|
-
sqrN: Fp2Square,
|
|
117
|
-
// Why inversion for bigint inside Fp instead of Fp2? it is even used in that context?
|
|
118
|
-
div: (lhs, rhs) => Fp2.mul(lhs, typeof rhs === 'bigint' ? Fp.inv(Fp.create(rhs)) : Fp2.inv(rhs)),
|
|
119
|
-
inv: ({ c0: a, c1: b }) => {
|
|
120
|
-
// We wish to find the multiplicative inverse of a nonzero
|
|
121
|
-
// element a + bu in Fp2. We leverage an identity
|
|
122
|
-
//
|
|
123
|
-
// (a + bu)(a - bu) = a² + b²
|
|
124
|
-
//
|
|
125
|
-
// which holds because u² = -1. This can be rewritten as
|
|
126
|
-
//
|
|
127
|
-
// (a + bu)(a - bu)/(a² + b²) = 1
|
|
128
|
-
//
|
|
129
|
-
// because a² + b² = 0 has no nonzero solutions for (a, b).
|
|
130
|
-
// This gives that (a - bu)/(a² + b²) is the inverse
|
|
131
|
-
// of (a + bu). Importantly, this can be computing using
|
|
132
|
-
// only a single inversion in Fp.
|
|
133
|
-
const factor = Fp.inv(Fp.create(a * a + b * b));
|
|
134
|
-
return { c0: Fp.mul(factor, Fp.create(a)), c1: Fp.mul(factor, Fp.create(-b)) };
|
|
135
|
-
},
|
|
136
|
-
sqrt: (num) => {
|
|
137
|
-
if (Fp2.eql(num, Fp2.ZERO))
|
|
138
|
-
return Fp2.ZERO; // Algo doesn't handles this case
|
|
139
|
-
// TODO: Optimize this line. It's extremely slow.
|
|
140
|
-
// Speeding this up would boost aggregateSignatures.
|
|
141
|
-
// https://eprint.iacr.org/2012/685.pdf applicable?
|
|
142
|
-
// https://github.com/zkcrypto/bls12_381/blob/080eaa74ec0e394377caa1ba302c8c121df08b07/src/fp2.rs#L250
|
|
143
|
-
// https://github.com/supranational/blst/blob/aae0c7d70b799ac269ff5edf29d8191dbd357876/src/exp2.c#L1
|
|
144
|
-
// Inspired by https://github.com/dalek-cryptography/curve25519-dalek/blob/17698df9d4c834204f83a3574143abacb4fc81a5/src/field.rs#L99
|
|
145
|
-
const candidateSqrt = Fp2.pow(num, (Fp2.ORDER + 8n) / 16n);
|
|
146
|
-
const check = Fp2.div(Fp2.sqr(candidateSqrt), num); // candidateSqrt.square().div(this);
|
|
147
|
-
const R = FP2_ROOTS_OF_UNITY;
|
|
148
|
-
const divisor = [R[0], R[2], R[4], R[6]].find((r) => Fp2.eql(r, check));
|
|
149
|
-
if (!divisor)
|
|
150
|
-
throw new Error('No root');
|
|
151
|
-
const index = R.indexOf(divisor);
|
|
152
|
-
const root = R[index / 2];
|
|
153
|
-
if (!root)
|
|
154
|
-
throw new Error('Invalid root');
|
|
155
|
-
const x1 = Fp2.div(candidateSqrt, root);
|
|
156
|
-
const x2 = Fp2.neg(x1);
|
|
157
|
-
const { re: re1, im: im1 } = Fp2.reim(x1);
|
|
158
|
-
const { re: re2, im: im2 } = Fp2.reim(x2);
|
|
159
|
-
if (im1 > im2 || (im1 === im2 && re1 > re2))
|
|
160
|
-
return x1;
|
|
161
|
-
return x2;
|
|
162
|
-
},
|
|
163
|
-
// Same as sgn0_fp2 in draft-irtf-cfrg-hash-to-curve-16
|
|
164
|
-
isOdd: (x) => {
|
|
165
|
-
const { re: x0, im: x1 } = Fp2.reim(x);
|
|
166
|
-
const sign_0 = x0 % 2n;
|
|
167
|
-
const zero_0 = x0 === 0n;
|
|
168
|
-
const sign_1 = x1 % 2n;
|
|
169
|
-
return BigInt(sign_0 || (zero_0 && sign_1)) == 1n;
|
|
170
|
-
},
|
|
171
|
-
// Bytes util
|
|
172
|
-
fromBytes(b) {
|
|
173
|
-
if (b.length !== Fp2.BYTES)
|
|
174
|
-
throw new Error(`fromBytes wrong length=${b.length}`);
|
|
175
|
-
return { c0: Fp.fromBytes(b.subarray(0, Fp.BYTES)), c1: Fp.fromBytes(b.subarray(Fp.BYTES)) };
|
|
176
|
-
},
|
|
177
|
-
toBytes: ({ c0, c1 }) => concatB(Fp.toBytes(c0), Fp.toBytes(c1)),
|
|
178
|
-
cmov: ({ c0, c1 }, { c0: r0, c1: r1 }, c) => ({
|
|
179
|
-
c0: Fp.cmov(c0, r0, c),
|
|
180
|
-
c1: Fp.cmov(c1, r1, c),
|
|
181
|
-
}),
|
|
182
|
-
// Specific utils
|
|
183
|
-
// toString() {
|
|
184
|
-
// return `Fp2(${this.c0} + ${this.c1}×i)`;
|
|
185
|
-
// }
|
|
186
|
-
reim: ({ c0, c1 }) => ({ re: c0, im: c1 }),
|
|
187
|
-
// multiply by u + 1
|
|
188
|
-
mulByNonresidue: ({ c0, c1 }) => ({ c0: Fp.sub(c0, c1), c1: Fp.add(c0, c1) }),
|
|
189
|
-
multiplyByB: ({ c0, c1 }) => {
|
|
190
|
-
let t0 = Fp.mul(c0, 4n); // 4 * c0
|
|
191
|
-
let t1 = Fp.mul(c1, 4n); // 4 * c1
|
|
192
|
-
// (T0-T1) + (T0+T1)*i
|
|
193
|
-
return { c0: Fp.sub(t0, t1), c1: Fp.add(t0, t1) };
|
|
194
|
-
},
|
|
195
|
-
fromBigTuple: (tuple) => {
|
|
196
|
-
if (tuple.length !== 2)
|
|
197
|
-
throw new Error('Invalid tuple');
|
|
198
|
-
const fps = tuple.map((n) => Fp.create(n));
|
|
199
|
-
return { c0: fps[0], c1: fps[1] };
|
|
200
|
-
},
|
|
201
|
-
frobeniusMap: ({ c0, c1 }, power) => ({
|
|
202
|
-
c0,
|
|
203
|
-
c1: Fp.mul(c1, FP2_FROBENIUS_COEFFICIENTS[power % 2]),
|
|
204
|
-
}),
|
|
205
|
-
};
|
|
206
|
-
// Finite extension field over irreducible polynominal.
|
|
207
|
-
// Fp(u) / (u² - β) where β = -1
|
|
208
|
-
const FP2_FROBENIUS_COEFFICIENTS = [
|
|
209
|
-
0x1n,
|
|
210
|
-
0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaaan,
|
|
211
|
-
].map((item) => Fp.create(item));
|
|
212
|
-
// For Fp2 roots of unity.
|
|
213
|
-
const rv1 = 0x6af0e0437ff400b6831e36d6bd17ffe48395dabc2d3435e77f76e17009241c5ee67992f72ec05f4c81084fbede3cc09n;
|
|
214
|
-
// const ev1 =
|
|
215
|
-
// 0x699be3b8c6870965e5bf892ad5d2cc7b0e85a117402dfd83b7f4a947e02d978498255a2aaec0ac627b5afbdf1bf1c90n;
|
|
216
|
-
// const ev2 =
|
|
217
|
-
// 0x8157cd83046453f5dd0972b6e3949e4288020b5b8a9cc99ca07e27089a2ce2436d965026adad3ef7baba37f2183e9b5n;
|
|
218
|
-
// const ev3 =
|
|
219
|
-
// 0xab1c2ffdd6c253ca155231eb3e71ba044fd562f6f72bc5bad5ec46a0b7a3b0247cf08ce6c6317f40edbc653a72dee17n;
|
|
220
|
-
// const ev4 =
|
|
221
|
-
// 0xaa404866706722864480885d68ad0ccac1967c7544b447873cc37e0181271e006df72162a3d3e0287bf597fbf7f8fc1n;
|
|
222
|
-
// Eighth roots of unity, used for computing square roots in Fp2.
|
|
223
|
-
// To verify or re-calculate:
|
|
224
|
-
// Array(8).fill(new Fp2([1n, 1n])).map((fp2, k) => fp2.pow(Fp2.ORDER * BigInt(k) / 8n))
|
|
225
|
-
const FP2_ROOTS_OF_UNITY = [
|
|
226
|
-
[1n, 0n],
|
|
227
|
-
[rv1, -rv1],
|
|
228
|
-
[0n, 1n],
|
|
229
|
-
[rv1, rv1],
|
|
230
|
-
[-1n, 0n],
|
|
231
|
-
[-rv1, rv1],
|
|
232
|
-
[0n, -1n],
|
|
233
|
-
[-rv1, -rv1],
|
|
234
|
-
].map((pair) => Fp2.fromBigTuple(pair));
|
|
235
|
-
const Fp6Add = ({ c0, c1, c2 }, { c0: r0, c1: r1, c2: r2 }) => ({
|
|
236
|
-
c0: Fp2.add(c0, r0),
|
|
237
|
-
c1: Fp2.add(c1, r1),
|
|
238
|
-
c2: Fp2.add(c2, r2),
|
|
239
|
-
});
|
|
240
|
-
const Fp6Subtract = ({ c0, c1, c2 }, { c0: r0, c1: r1, c2: r2 }) => ({
|
|
241
|
-
c0: Fp2.sub(c0, r0),
|
|
242
|
-
c1: Fp2.sub(c1, r1),
|
|
243
|
-
c2: Fp2.sub(c2, r2),
|
|
244
|
-
});
|
|
245
|
-
const Fp6Multiply = ({ c0, c1, c2 }, rhs) => {
|
|
246
|
-
if (typeof rhs === 'bigint') {
|
|
247
|
-
return {
|
|
248
|
-
c0: Fp2.mul(c0, rhs),
|
|
249
|
-
c1: Fp2.mul(c1, rhs),
|
|
250
|
-
c2: Fp2.mul(c2, rhs),
|
|
251
|
-
};
|
|
252
|
-
}
|
|
253
|
-
const { c0: r0, c1: r1, c2: r2 } = rhs;
|
|
254
|
-
const t0 = Fp2.mul(c0, r0); // c0 * o0
|
|
255
|
-
const t1 = Fp2.mul(c1, r1); // c1 * o1
|
|
256
|
-
const t2 = Fp2.mul(c2, r2); // c2 * o2
|
|
257
|
-
return {
|
|
258
|
-
// t0 + (c1 + c2) * (r1 * r2) - (T1 + T2) * (u + 1)
|
|
259
|
-
c0: Fp2.add(t0, Fp2.mulByNonresidue(Fp2.sub(Fp2.mul(Fp2.add(c1, c2), Fp2.add(r1, r2)), Fp2.add(t1, t2)))),
|
|
260
|
-
// (c0 + c1) * (r0 + r1) - (T0 + T1) + T2 * (u + 1)
|
|
261
|
-
c1: Fp2.add(Fp2.sub(Fp2.mul(Fp2.add(c0, c1), Fp2.add(r0, r1)), Fp2.add(t0, t1)), Fp2.mulByNonresidue(t2)),
|
|
262
|
-
// T1 + (c0 + c2) * (r0 + r2) - T0 + T2
|
|
263
|
-
c2: Fp2.sub(Fp2.add(t1, Fp2.mul(Fp2.add(c0, c2), Fp2.add(r0, r2))), Fp2.add(t0, t2)),
|
|
264
|
-
};
|
|
265
|
-
};
|
|
266
|
-
const Fp6Square = ({ c0, c1, c2 }) => {
|
|
267
|
-
let t0 = Fp2.sqr(c0); // c0²
|
|
268
|
-
let t1 = Fp2.mul(Fp2.mul(c0, c1), 2n); // 2 * c0 * c1
|
|
269
|
-
let t3 = Fp2.mul(Fp2.mul(c1, c2), 2n); // 2 * c1 * c2
|
|
270
|
-
let t4 = Fp2.sqr(c2); // c2²
|
|
271
|
-
return {
|
|
272
|
-
c0: Fp2.add(Fp2.mulByNonresidue(t3), t0),
|
|
273
|
-
c1: Fp2.add(Fp2.mulByNonresidue(t4), t1),
|
|
274
|
-
// T1 + (c0 - c1 + c2)² + T3 - T0 - T4
|
|
275
|
-
c2: Fp2.sub(Fp2.sub(Fp2.add(Fp2.add(t1, Fp2.sqr(Fp2.add(Fp2.sub(c0, c1), c2))), t3), t0), t4),
|
|
276
|
-
};
|
|
277
|
-
};
|
|
278
|
-
const Fp6 = {
|
|
279
|
-
ORDER: Fp2.ORDER,
|
|
280
|
-
BITS: 3 * Fp2.BITS,
|
|
281
|
-
BYTES: 3 * Fp2.BYTES,
|
|
282
|
-
MASK: bitMask(3 * Fp2.BITS),
|
|
283
|
-
ZERO: { c0: Fp2.ZERO, c1: Fp2.ZERO, c2: Fp2.ZERO },
|
|
284
|
-
ONE: { c0: Fp2.ONE, c1: Fp2.ZERO, c2: Fp2.ZERO },
|
|
285
|
-
create: (num) => num,
|
|
286
|
-
isValid: ({ c0, c1, c2 }) => Fp2.isValid(c0) && Fp2.isValid(c1) && Fp2.isValid(c2),
|
|
287
|
-
is0: ({ c0, c1, c2 }) => Fp2.is0(c0) && Fp2.is0(c1) && Fp2.is0(c2),
|
|
288
|
-
neg: ({ c0, c1, c2 }) => ({ c0: Fp2.neg(c0), c1: Fp2.neg(c1), c2: Fp2.neg(c2) }),
|
|
289
|
-
eql: ({ c0, c1, c2 }, { c0: r0, c1: r1, c2: r2 }) => Fp2.eql(c0, r0) && Fp2.eql(c1, r1) && Fp2.eql(c2, r2),
|
|
290
|
-
sqrt: () => {
|
|
291
|
-
throw new Error('Not implemented');
|
|
292
|
-
},
|
|
293
|
-
// Do we need division by bigint at all? Should be done via order:
|
|
294
|
-
div: (lhs, rhs) => Fp6.mul(lhs, typeof rhs === 'bigint' ? Fp.inv(Fp.create(rhs)) : Fp6.inv(rhs)),
|
|
295
|
-
pow: (num, power) => mod.FpPow(Fp6, num, power),
|
|
296
|
-
invertBatch: (nums) => mod.FpInvertBatch(Fp6, nums),
|
|
297
|
-
// Normalized
|
|
298
|
-
add: Fp6Add,
|
|
299
|
-
sub: Fp6Subtract,
|
|
300
|
-
mul: Fp6Multiply,
|
|
301
|
-
sqr: Fp6Square,
|
|
302
|
-
// NonNormalized stuff
|
|
303
|
-
addN: Fp6Add,
|
|
304
|
-
subN: Fp6Subtract,
|
|
305
|
-
mulN: Fp6Multiply,
|
|
306
|
-
sqrN: Fp6Square,
|
|
307
|
-
inv: ({ c0, c1, c2 }) => {
|
|
308
|
-
let t0 = Fp2.sub(Fp2.sqr(c0), Fp2.mulByNonresidue(Fp2.mul(c2, c1))); // c0² - c2 * c1 * (u + 1)
|
|
309
|
-
let t1 = Fp2.sub(Fp2.mulByNonresidue(Fp2.sqr(c2)), Fp2.mul(c0, c1)); // c2² * (u + 1) - c0 * c1
|
|
310
|
-
let t2 = Fp2.sub(Fp2.sqr(c1), Fp2.mul(c0, c2)); // c1² - c0 * c2
|
|
311
|
-
// 1/(((c2 * T1 + c1 * T2) * v) + c0 * T0)
|
|
312
|
-
let t4 = Fp2.inv(Fp2.add(Fp2.mulByNonresidue(Fp2.add(Fp2.mul(c2, t1), Fp2.mul(c1, t2))), Fp2.mul(c0, t0)));
|
|
313
|
-
return { c0: Fp2.mul(t4, t0), c1: Fp2.mul(t4, t1), c2: Fp2.mul(t4, t2) };
|
|
314
|
-
},
|
|
315
|
-
// Bytes utils
|
|
316
|
-
fromBytes: (b) => {
|
|
317
|
-
if (b.length !== Fp6.BYTES)
|
|
318
|
-
throw new Error(`fromBytes wrong length=${b.length}`);
|
|
319
|
-
return {
|
|
320
|
-
c0: Fp2.fromBytes(b.subarray(0, Fp2.BYTES)),
|
|
321
|
-
c1: Fp2.fromBytes(b.subarray(Fp2.BYTES, 2 * Fp2.BYTES)),
|
|
322
|
-
c2: Fp2.fromBytes(b.subarray(2 * Fp2.BYTES)),
|
|
323
|
-
};
|
|
324
|
-
},
|
|
325
|
-
toBytes: ({ c0, c1, c2 }) => concatB(Fp2.toBytes(c0), Fp2.toBytes(c1), Fp2.toBytes(c2)),
|
|
326
|
-
cmov: ({ c0, c1, c2 }, { c0: r0, c1: r1, c2: r2 }, c) => ({
|
|
327
|
-
c0: Fp2.cmov(c0, r0, c),
|
|
328
|
-
c1: Fp2.cmov(c1, r1, c),
|
|
329
|
-
c2: Fp2.cmov(c2, r2, c),
|
|
330
|
-
}),
|
|
331
|
-
// Utils
|
|
332
|
-
// fromTriple(triple: [Fp2, Fp2, Fp2]) {
|
|
333
|
-
// return new Fp6(...triple);
|
|
334
|
-
// }
|
|
335
|
-
// toString() {
|
|
336
|
-
// return `Fp6(${this.c0} + ${this.c1} * v, ${this.c2} * v^2)`;
|
|
337
|
-
// }
|
|
338
|
-
fromBigSix: (t) => {
|
|
339
|
-
if (!Array.isArray(t) || t.length !== 6)
|
|
340
|
-
throw new Error('Invalid Fp6 usage');
|
|
341
|
-
return {
|
|
342
|
-
c0: Fp2.fromBigTuple(t.slice(0, 2)),
|
|
343
|
-
c1: Fp2.fromBigTuple(t.slice(2, 4)),
|
|
344
|
-
c2: Fp2.fromBigTuple(t.slice(4, 6)),
|
|
345
|
-
};
|
|
346
|
-
},
|
|
347
|
-
frobeniusMap: ({ c0, c1, c2 }, power) => ({
|
|
348
|
-
c0: Fp2.frobeniusMap(c0, power),
|
|
349
|
-
c1: Fp2.mul(Fp2.frobeniusMap(c1, power), FP6_FROBENIUS_COEFFICIENTS_1[power % 6]),
|
|
350
|
-
c2: Fp2.mul(Fp2.frobeniusMap(c2, power), FP6_FROBENIUS_COEFFICIENTS_2[power % 6]),
|
|
351
|
-
}),
|
|
352
|
-
mulByNonresidue: ({ c0, c1, c2 }) => ({ c0: Fp2.mulByNonresidue(c2), c1: c0, c2: c1 }),
|
|
353
|
-
// Sparse multiplication
|
|
354
|
-
multiplyBy1: ({ c0, c1, c2 }, b1) => ({
|
|
355
|
-
c0: Fp2.mulByNonresidue(Fp2.mul(c2, b1)),
|
|
356
|
-
c1: Fp2.mul(c0, b1),
|
|
357
|
-
c2: Fp2.mul(c1, b1),
|
|
358
|
-
}),
|
|
359
|
-
// Sparse multiplication
|
|
360
|
-
multiplyBy01({ c0, c1, c2 }, b0, b1) {
|
|
361
|
-
let t0 = Fp2.mul(c0, b0); // c0 * b0
|
|
362
|
-
let t1 = Fp2.mul(c1, b1); // c1 * b1
|
|
363
|
-
return {
|
|
364
|
-
// ((c1 + c2) * b1 - T1) * (u + 1) + T0
|
|
365
|
-
c0: Fp2.add(Fp2.mulByNonresidue(Fp2.sub(Fp2.mul(Fp2.add(c1, c2), b1), t1)), t0),
|
|
366
|
-
// (b0 + b1) * (c0 + c1) - T0 - T1
|
|
367
|
-
c1: Fp2.sub(Fp2.sub(Fp2.mul(Fp2.add(b0, b1), Fp2.add(c0, c1)), t0), t1),
|
|
368
|
-
// (c0 + c2) * b0 - T0 + T1
|
|
369
|
-
c2: Fp2.add(Fp2.sub(Fp2.mul(Fp2.add(c0, c2), b0), t0), t1),
|
|
370
|
-
};
|
|
371
|
-
},
|
|
372
|
-
multiplyByFp2: ({ c0, c1, c2 }, rhs) => ({
|
|
373
|
-
c0: Fp2.mul(c0, rhs),
|
|
374
|
-
c1: Fp2.mul(c1, rhs),
|
|
375
|
-
c2: Fp2.mul(c2, rhs),
|
|
376
|
-
}),
|
|
377
|
-
};
|
|
378
|
-
const FP6_FROBENIUS_COEFFICIENTS_1 = [
|
|
379
|
-
[0x1n, 0x0n],
|
|
380
|
-
[
|
|
381
|
-
0x0n,
|
|
382
|
-
0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaacn,
|
|
383
|
-
],
|
|
384
|
-
[
|
|
385
|
-
0x00000000000000005f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffefffen,
|
|
386
|
-
0x0n,
|
|
387
|
-
],
|
|
388
|
-
[0x0n, 0x1n],
|
|
389
|
-
[
|
|
390
|
-
0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaacn,
|
|
391
|
-
0x0n,
|
|
392
|
-
],
|
|
393
|
-
[
|
|
394
|
-
0x0n,
|
|
395
|
-
0x00000000000000005f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffefffen,
|
|
396
|
-
],
|
|
397
|
-
].map((pair) => Fp2.fromBigTuple(pair));
|
|
398
|
-
const FP6_FROBENIUS_COEFFICIENTS_2 = [
|
|
399
|
-
[0x1n, 0x0n],
|
|
400
|
-
[
|
|
401
|
-
0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaadn,
|
|
402
|
-
0x0n,
|
|
403
|
-
],
|
|
404
|
-
[
|
|
405
|
-
0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaacn,
|
|
406
|
-
0x0n,
|
|
407
|
-
],
|
|
408
|
-
[
|
|
409
|
-
0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaaan,
|
|
410
|
-
0x0n,
|
|
411
|
-
],
|
|
412
|
-
[
|
|
413
|
-
0x00000000000000005f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffefffen,
|
|
414
|
-
0x0n,
|
|
415
|
-
],
|
|
416
|
-
[
|
|
417
|
-
0x00000000000000005f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffeffffn,
|
|
418
|
-
0x0n,
|
|
419
|
-
],
|
|
420
|
-
].map((pair) => Fp2.fromBigTuple(pair));
|
|
421
|
-
// The BLS parameter x for BLS12-381
|
|
422
|
-
const BLS_X = 0xd201000000010000n;
|
|
423
|
-
const BLS_X_LEN = bitLen(BLS_X);
|
|
424
|
-
const Fp12Add = ({ c0, c1 }, { c0: r0, c1: r1 }) => ({
|
|
425
|
-
c0: Fp6.add(c0, r0),
|
|
426
|
-
c1: Fp6.add(c1, r1),
|
|
427
|
-
});
|
|
428
|
-
const Fp12Subtract = ({ c0, c1 }, { c0: r0, c1: r1 }) => ({
|
|
429
|
-
c0: Fp6.sub(c0, r0),
|
|
430
|
-
c1: Fp6.sub(c1, r1),
|
|
431
|
-
});
|
|
432
|
-
const Fp12Multiply = ({ c0, c1 }, rhs) => {
|
|
433
|
-
if (typeof rhs === 'bigint')
|
|
434
|
-
return { c0: Fp6.mul(c0, rhs), c1: Fp6.mul(c1, rhs) };
|
|
435
|
-
let { c0: r0, c1: r1 } = rhs;
|
|
436
|
-
let t1 = Fp6.mul(c0, r0); // c0 * r0
|
|
437
|
-
let t2 = Fp6.mul(c1, r1); // c1 * r1
|
|
438
|
-
return {
|
|
439
|
-
c0: Fp6.add(t1, Fp6.mulByNonresidue(t2)),
|
|
440
|
-
// (c0 + c1) * (r0 + r1) - (T1 + T2)
|
|
441
|
-
c1: Fp6.sub(Fp6.mul(Fp6.add(c0, c1), Fp6.add(r0, r1)), Fp6.add(t1, t2)),
|
|
442
|
-
};
|
|
443
|
-
};
|
|
444
|
-
const Fp12Square = ({ c0, c1 }) => {
|
|
445
|
-
let ab = Fp6.mul(c0, c1); // c0 * c1
|
|
446
|
-
return {
|
|
447
|
-
// (c1 * v + c0) * (c0 + c1) - AB - AB * v
|
|
448
|
-
c0: Fp6.sub(Fp6.sub(Fp6.mul(Fp6.add(Fp6.mulByNonresidue(c1), c0), Fp6.add(c0, c1)), ab), Fp6.mulByNonresidue(ab)),
|
|
449
|
-
c1: Fp6.add(ab, ab),
|
|
450
|
-
}; // AB + AB
|
|
451
|
-
};
|
|
452
|
-
function Fp4Square(a, b) {
|
|
453
|
-
const a2 = Fp2.sqr(a);
|
|
454
|
-
const b2 = Fp2.sqr(b);
|
|
455
|
-
return {
|
|
456
|
-
first: Fp2.add(Fp2.mulByNonresidue(b2), a2),
|
|
457
|
-
second: Fp2.sub(Fp2.sub(Fp2.sqr(Fp2.add(a, b)), a2), b2), // (a + b)² - a² - b²
|
|
458
|
-
};
|
|
459
|
-
}
|
|
460
|
-
const Fp12 = {
|
|
461
|
-
ORDER: Fp2.ORDER,
|
|
462
|
-
BITS: 2 * Fp2.BITS,
|
|
463
|
-
BYTES: 2 * Fp2.BYTES,
|
|
464
|
-
MASK: bitMask(2 * Fp2.BITS),
|
|
465
|
-
ZERO: { c0: Fp6.ZERO, c1: Fp6.ZERO },
|
|
466
|
-
ONE: { c0: Fp6.ONE, c1: Fp6.ZERO },
|
|
467
|
-
create: (num) => num,
|
|
468
|
-
isValid: ({ c0, c1 }) => Fp6.isValid(c0) && Fp6.isValid(c1),
|
|
469
|
-
is0: ({ c0, c1 }) => Fp6.is0(c0) && Fp6.is0(c1),
|
|
470
|
-
neg: ({ c0, c1 }) => ({ c0: Fp6.neg(c0), c1: Fp6.neg(c1) }),
|
|
471
|
-
eql: ({ c0, c1 }, { c0: r0, c1: r1 }) => Fp6.eql(c0, r0) && Fp6.eql(c1, r1),
|
|
472
|
-
sqrt: () => {
|
|
473
|
-
throw new Error('Not implemented');
|
|
474
|
-
},
|
|
475
|
-
inv: ({ c0, c1 }) => {
|
|
476
|
-
let t = Fp6.inv(Fp6.sub(Fp6.sqr(c0), Fp6.mulByNonresidue(Fp6.sqr(c1)))); // 1 / (c0² - c1² * v)
|
|
477
|
-
return { c0: Fp6.mul(c0, t), c1: Fp6.neg(Fp6.mul(c1, t)) }; // ((C0 * T) * T) + (-C1 * T) * w
|
|
478
|
-
},
|
|
479
|
-
div: (lhs, rhs) => Fp12.mul(lhs, typeof rhs === 'bigint' ? Fp.inv(Fp.create(rhs)) : Fp12.inv(rhs)),
|
|
480
|
-
pow: (num, power) => mod.FpPow(Fp12, num, power),
|
|
481
|
-
invertBatch: (nums) => mod.FpInvertBatch(Fp12, nums),
|
|
482
|
-
// Normalized
|
|
483
|
-
add: Fp12Add,
|
|
484
|
-
sub: Fp12Subtract,
|
|
485
|
-
mul: Fp12Multiply,
|
|
486
|
-
sqr: Fp12Square,
|
|
487
|
-
// NonNormalized stuff
|
|
488
|
-
addN: Fp12Add,
|
|
489
|
-
subN: Fp12Subtract,
|
|
490
|
-
mulN: Fp12Multiply,
|
|
491
|
-
sqrN: Fp12Square,
|
|
492
|
-
// Bytes utils
|
|
493
|
-
fromBytes: (b) => {
|
|
494
|
-
if (b.length !== Fp12.BYTES)
|
|
495
|
-
throw new Error(`fromBytes wrong length=${b.length}`);
|
|
496
|
-
return {
|
|
497
|
-
c0: Fp6.fromBytes(b.subarray(0, Fp6.BYTES)),
|
|
498
|
-
c1: Fp6.fromBytes(b.subarray(Fp6.BYTES)),
|
|
499
|
-
};
|
|
500
|
-
},
|
|
501
|
-
toBytes: ({ c0, c1 }) => concatB(Fp6.toBytes(c0), Fp6.toBytes(c1)),
|
|
502
|
-
cmov: ({ c0, c1 }, { c0: r0, c1: r1 }, c) => ({
|
|
503
|
-
c0: Fp6.cmov(c0, r0, c),
|
|
504
|
-
c1: Fp6.cmov(c1, r1, c),
|
|
505
|
-
}),
|
|
506
|
-
// Utils
|
|
507
|
-
// toString() {
|
|
508
|
-
// return `Fp12(${this.c0} + ${this.c1} * w)`;
|
|
509
|
-
// },
|
|
510
|
-
// fromTuple(c: [Fp6, Fp6]) {
|
|
511
|
-
// return new Fp12(...c);
|
|
512
|
-
// }
|
|
513
|
-
fromBigTwelve: (t) => ({
|
|
514
|
-
c0: Fp6.fromBigSix(t.slice(0, 6)),
|
|
515
|
-
c1: Fp6.fromBigSix(t.slice(6, 12)),
|
|
516
|
-
}),
|
|
517
|
-
// Raises to q**i -th power
|
|
518
|
-
frobeniusMap(lhs, power) {
|
|
519
|
-
const r0 = Fp6.frobeniusMap(lhs.c0, power);
|
|
520
|
-
const { c0, c1, c2 } = Fp6.frobeniusMap(lhs.c1, power);
|
|
521
|
-
const coeff = FP12_FROBENIUS_COEFFICIENTS[power % 12];
|
|
522
|
-
return {
|
|
523
|
-
c0: r0,
|
|
524
|
-
c1: Fp6.create({
|
|
525
|
-
c0: Fp2.mul(c0, coeff),
|
|
526
|
-
c1: Fp2.mul(c1, coeff),
|
|
527
|
-
c2: Fp2.mul(c2, coeff),
|
|
528
|
-
}),
|
|
529
|
-
};
|
|
530
|
-
},
|
|
531
|
-
// Sparse multiplication
|
|
532
|
-
multiplyBy014: ({ c0, c1 }, o0, o1, o4) => {
|
|
533
|
-
let t0 = Fp6.multiplyBy01(c0, o0, o1);
|
|
534
|
-
let t1 = Fp6.multiplyBy1(c1, o4);
|
|
535
|
-
return {
|
|
536
|
-
c0: Fp6.add(Fp6.mulByNonresidue(t1), t0),
|
|
537
|
-
// (c1 + c0) * [o0, o1+o4] - T0 - T1
|
|
538
|
-
c1: Fp6.sub(Fp6.sub(Fp6.multiplyBy01(Fp6.add(c1, c0), o0, Fp2.add(o1, o4)), t0), t1),
|
|
539
|
-
};
|
|
540
|
-
},
|
|
541
|
-
multiplyByFp2: ({ c0, c1 }, rhs) => ({
|
|
542
|
-
c0: Fp6.multiplyByFp2(c0, rhs),
|
|
543
|
-
c1: Fp6.multiplyByFp2(c1, rhs),
|
|
544
|
-
}),
|
|
545
|
-
conjugate: ({ c0, c1 }) => ({ c0, c1: Fp6.neg(c1) }),
|
|
546
|
-
// A cyclotomic group is a subgroup of Fp^n defined by
|
|
547
|
-
// GΦₙ(p) = {α ∈ Fpⁿ : α^Φₙ(p) = 1}
|
|
548
|
-
// The result of any pairing is in a cyclotomic subgroup
|
|
549
|
-
// https://eprint.iacr.org/2009/565.pdf
|
|
550
|
-
_cyclotomicSquare: ({ c0, c1 }) => {
|
|
551
|
-
const { c0: c0c0, c1: c0c1, c2: c0c2 } = c0;
|
|
552
|
-
const { c0: c1c0, c1: c1c1, c2: c1c2 } = c1;
|
|
553
|
-
const { first: t3, second: t4 } = Fp4Square(c0c0, c1c1);
|
|
554
|
-
const { first: t5, second: t6 } = Fp4Square(c1c0, c0c2);
|
|
555
|
-
const { first: t7, second: t8 } = Fp4Square(c0c1, c1c2);
|
|
556
|
-
let t9 = Fp2.mulByNonresidue(t8); // T8 * (u + 1)
|
|
557
|
-
return {
|
|
558
|
-
c0: Fp6.create({
|
|
559
|
-
c0: Fp2.add(Fp2.mul(Fp2.sub(t3, c0c0), 2n), t3),
|
|
560
|
-
c1: Fp2.add(Fp2.mul(Fp2.sub(t5, c0c1), 2n), t5),
|
|
561
|
-
c2: Fp2.add(Fp2.mul(Fp2.sub(t7, c0c2), 2n), t7),
|
|
562
|
-
}),
|
|
563
|
-
c1: Fp6.create({
|
|
564
|
-
c0: Fp2.add(Fp2.mul(Fp2.add(t9, c1c0), 2n), t9),
|
|
565
|
-
c1: Fp2.add(Fp2.mul(Fp2.add(t4, c1c1), 2n), t4),
|
|
566
|
-
c2: Fp2.add(Fp2.mul(Fp2.add(t6, c1c2), 2n), t6),
|
|
567
|
-
}),
|
|
568
|
-
}; // 2 * (T6 + c1c2) + T6
|
|
569
|
-
},
|
|
570
|
-
_cyclotomicExp(num, n) {
|
|
571
|
-
let z = Fp12.ONE;
|
|
572
|
-
for (let i = BLS_X_LEN - 1; i >= 0; i--) {
|
|
573
|
-
z = Fp12._cyclotomicSquare(z);
|
|
574
|
-
if (bitGet(n, i))
|
|
575
|
-
z = Fp12.mul(z, num);
|
|
576
|
-
}
|
|
577
|
-
return z;
|
|
578
|
-
},
|
|
579
|
-
// https://eprint.iacr.org/2010/354.pdf
|
|
580
|
-
// https://eprint.iacr.org/2009/565.pdf
|
|
581
|
-
finalExponentiate: (num) => {
|
|
582
|
-
const x = BLS_X;
|
|
583
|
-
// this^(q⁶) / this
|
|
584
|
-
const t0 = Fp12.div(Fp12.frobeniusMap(num, 6), num);
|
|
585
|
-
// t0^(q²) * t0
|
|
586
|
-
const t1 = Fp12.mul(Fp12.frobeniusMap(t0, 2), t0);
|
|
587
|
-
const t2 = Fp12.conjugate(Fp12._cyclotomicExp(t1, x));
|
|
588
|
-
const t3 = Fp12.mul(Fp12.conjugate(Fp12._cyclotomicSquare(t1)), t2);
|
|
589
|
-
const t4 = Fp12.conjugate(Fp12._cyclotomicExp(t3, x));
|
|
590
|
-
const t5 = Fp12.conjugate(Fp12._cyclotomicExp(t4, x));
|
|
591
|
-
const t6 = Fp12.mul(Fp12.conjugate(Fp12._cyclotomicExp(t5, x)), Fp12._cyclotomicSquare(t2));
|
|
592
|
-
const t7 = Fp12.conjugate(Fp12._cyclotomicExp(t6, x));
|
|
593
|
-
const t2_t5_pow_q2 = Fp12.frobeniusMap(Fp12.mul(t2, t5), 2);
|
|
594
|
-
const t4_t1_pow_q3 = Fp12.frobeniusMap(Fp12.mul(t4, t1), 3);
|
|
595
|
-
const t6_t1c_pow_q1 = Fp12.frobeniusMap(Fp12.mul(t6, Fp12.conjugate(t1)), 1);
|
|
596
|
-
const t7_t3c_t1 = Fp12.mul(Fp12.mul(t7, Fp12.conjugate(t3)), t1);
|
|
597
|
-
// (t2 * t5)^(q²) * (t4 * t1)^(q³) * (t6 * t1.conj)^(q^1) * t7 * t3.conj * t1
|
|
598
|
-
return Fp12.mul(Fp12.mul(Fp12.mul(t2_t5_pow_q2, t4_t1_pow_q3), t6_t1c_pow_q1), t7_t3c_t1);
|
|
599
|
-
},
|
|
600
|
-
};
|
|
601
|
-
const FP12_FROBENIUS_COEFFICIENTS = [
|
|
602
|
-
[0x1n, 0x0n],
|
|
603
|
-
[
|
|
604
|
-
0x1904d3bf02bb0667c231beb4202c0d1f0fd603fd3cbd5f4f7b2443d784bab9c4f67ea53d63e7813d8d0775ed92235fb8n,
|
|
605
|
-
0x00fc3e2b36c4e03288e9e902231f9fb854a14787b6c7b36fec0c8ec971f63c5f282d5ac14d6c7ec22cf78a126ddc4af3n,
|
|
606
|
-
],
|
|
607
|
-
[
|
|
608
|
-
0x00000000000000005f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffeffffn,
|
|
609
|
-
0x0n,
|
|
610
|
-
],
|
|
611
|
-
[
|
|
612
|
-
0x135203e60180a68ee2e9c448d77a2cd91c3dedd930b1cf60ef396489f61eb45e304466cf3e67fa0af1ee7b04121bdea2n,
|
|
613
|
-
0x06af0e0437ff400b6831e36d6bd17ffe48395dabc2d3435e77f76e17009241c5ee67992f72ec05f4c81084fbede3cc09n,
|
|
614
|
-
],
|
|
615
|
-
[
|
|
616
|
-
0x00000000000000005f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffefffen,
|
|
617
|
-
0x0n,
|
|
618
|
-
],
|
|
619
|
-
[
|
|
620
|
-
0x144e4211384586c16bd3ad4afa99cc9170df3560e77982d0db45f3536814f0bd5871c1908bd478cd1ee605167ff82995n,
|
|
621
|
-
0x05b2cfd9013a5fd8df47fa6b48b1e045f39816240c0b8fee8beadf4d8e9c0566c63a3e6e257f87329b18fae980078116n,
|
|
622
|
-
],
|
|
623
|
-
[
|
|
624
|
-
0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaaan,
|
|
625
|
-
0x0n,
|
|
626
|
-
],
|
|
627
|
-
[
|
|
628
|
-
0x00fc3e2b36c4e03288e9e902231f9fb854a14787b6c7b36fec0c8ec971f63c5f282d5ac14d6c7ec22cf78a126ddc4af3n,
|
|
629
|
-
0x1904d3bf02bb0667c231beb4202c0d1f0fd603fd3cbd5f4f7b2443d784bab9c4f67ea53d63e7813d8d0775ed92235fb8n,
|
|
630
|
-
],
|
|
631
|
-
[
|
|
632
|
-
0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaacn,
|
|
633
|
-
0x0n,
|
|
634
|
-
],
|
|
635
|
-
[
|
|
636
|
-
0x06af0e0437ff400b6831e36d6bd17ffe48395dabc2d3435e77f76e17009241c5ee67992f72ec05f4c81084fbede3cc09n,
|
|
637
|
-
0x135203e60180a68ee2e9c448d77a2cd91c3dedd930b1cf60ef396489f61eb45e304466cf3e67fa0af1ee7b04121bdea2n,
|
|
638
|
-
],
|
|
639
|
-
[
|
|
640
|
-
0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaadn,
|
|
641
|
-
0x0n,
|
|
642
|
-
],
|
|
643
|
-
[
|
|
644
|
-
0x05b2cfd9013a5fd8df47fa6b48b1e045f39816240c0b8fee8beadf4d8e9c0566c63a3e6e257f87329b18fae980078116n,
|
|
645
|
-
0x144e4211384586c16bd3ad4afa99cc9170df3560e77982d0db45f3536814f0bd5871c1908bd478cd1ee605167ff82995n,
|
|
646
|
-
],
|
|
647
|
-
].map((n) => Fp2.fromBigTuple(n));
|
|
648
|
-
// END OF CURVE FIELDS
|
|
649
|
-
// HashToCurve
|
|
650
|
-
// 3-isogeny map from E' to E
|
|
651
|
-
// https://datatracker.ietf.org/doc/html/draft-irtf-cfrg-hash-to-curve-11#appendix-E.3
|
|
652
|
-
const isogenyMapG2 = isogenyMap(Fp2, [
|
|
653
|
-
// xNum
|
|
654
|
-
[
|
|
655
|
-
[
|
|
656
|
-
'0x5c759507e8e333ebb5b7a9a47d7ed8532c52d39fd3a042a88b58423c50ae15d5c2638e343d9c71c6238aaaaaaaa97d6',
|
|
657
|
-
'0x5c759507e8e333ebb5b7a9a47d7ed8532c52d39fd3a042a88b58423c50ae15d5c2638e343d9c71c6238aaaaaaaa97d6',
|
|
658
|
-
],
|
|
659
|
-
[
|
|
660
|
-
'0x0',
|
|
661
|
-
'0x11560bf17baa99bc32126fced787c88f984f87adf7ae0c7f9a208c6b4f20a4181472aaa9cb8d555526a9ffffffffc71a',
|
|
662
|
-
],
|
|
663
|
-
[
|
|
664
|
-
'0x11560bf17baa99bc32126fced787c88f984f87adf7ae0c7f9a208c6b4f20a4181472aaa9cb8d555526a9ffffffffc71e',
|
|
665
|
-
'0x8ab05f8bdd54cde190937e76bc3e447cc27c3d6fbd7063fcd104635a790520c0a395554e5c6aaaa9354ffffffffe38d',
|
|
666
|
-
],
|
|
667
|
-
[
|
|
668
|
-
'0x171d6541fa38ccfaed6dea691f5fb614cb14b4e7f4e810aa22d6108f142b85757098e38d0f671c7188e2aaaaaaaa5ed1',
|
|
669
|
-
'0x0',
|
|
670
|
-
],
|
|
671
|
-
],
|
|
672
|
-
// xDen
|
|
673
|
-
[
|
|
674
|
-
[
|
|
675
|
-
'0x0',
|
|
676
|
-
'0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaa63',
|
|
677
|
-
],
|
|
678
|
-
[
|
|
679
|
-
'0xc',
|
|
680
|
-
'0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaa9f',
|
|
681
|
-
],
|
|
682
|
-
['0x1', '0x0'], // LAST 1
|
|
683
|
-
],
|
|
684
|
-
// yNum
|
|
685
|
-
[
|
|
686
|
-
[
|
|
687
|
-
'0x1530477c7ab4113b59a4c18b076d11930f7da5d4a07f649bf54439d87d27e500fc8c25ebf8c92f6812cfc71c71c6d706',
|
|
688
|
-
'0x1530477c7ab4113b59a4c18b076d11930f7da5d4a07f649bf54439d87d27e500fc8c25ebf8c92f6812cfc71c71c6d706',
|
|
689
|
-
],
|
|
690
|
-
[
|
|
691
|
-
'0x0',
|
|
692
|
-
'0x5c759507e8e333ebb5b7a9a47d7ed8532c52d39fd3a042a88b58423c50ae15d5c2638e343d9c71c6238aaaaaaaa97be',
|
|
693
|
-
],
|
|
694
|
-
[
|
|
695
|
-
'0x11560bf17baa99bc32126fced787c88f984f87adf7ae0c7f9a208c6b4f20a4181472aaa9cb8d555526a9ffffffffc71c',
|
|
696
|
-
'0x8ab05f8bdd54cde190937e76bc3e447cc27c3d6fbd7063fcd104635a790520c0a395554e5c6aaaa9354ffffffffe38f',
|
|
697
|
-
],
|
|
698
|
-
[
|
|
699
|
-
'0x124c9ad43b6cf79bfbf7043de3811ad0761b0f37a1e26286b0e977c69aa274524e79097a56dc4bd9e1b371c71c718b10',
|
|
700
|
-
'0x0',
|
|
701
|
-
],
|
|
702
|
-
],
|
|
703
|
-
// yDen
|
|
704
|
-
[
|
|
705
|
-
[
|
|
706
|
-
'0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffa8fb',
|
|
707
|
-
'0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffa8fb',
|
|
708
|
-
],
|
|
709
|
-
[
|
|
710
|
-
'0x0',
|
|
711
|
-
'0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffa9d3',
|
|
712
|
-
],
|
|
713
|
-
[
|
|
714
|
-
'0x12',
|
|
715
|
-
'0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaa99',
|
|
716
|
-
],
|
|
717
|
-
['0x1', '0x0'], // LAST 1
|
|
718
|
-
],
|
|
719
|
-
].map((i) => i.map((pair) => Fp2.fromBigTuple(pair.map(BigInt)))));
|
|
720
|
-
// 11-isogeny map from E' to E
|
|
721
|
-
const isogenyMapG1 = isogenyMap(Fp, [
|
|
722
|
-
// xNum
|
|
723
|
-
[
|
|
724
|
-
'0x11a05f2b1e833340b809101dd99815856b303e88a2d7005ff2627b56cdb4e2c85610c2d5f2e62d6eaeac1662734649b7',
|
|
725
|
-
'0x17294ed3e943ab2f0588bab22147a81c7c17e75b2f6a8417f565e33c70d1e86b4838f2a6f318c356e834eef1b3cb83bb',
|
|
726
|
-
'0xd54005db97678ec1d1048c5d10a9a1bce032473295983e56878e501ec68e25c958c3e3d2a09729fe0179f9dac9edcb0',
|
|
727
|
-
'0x1778e7166fcc6db74e0609d307e55412d7f5e4656a8dbf25f1b33289f1b330835336e25ce3107193c5b388641d9b6861',
|
|
728
|
-
'0xe99726a3199f4436642b4b3e4118e5499db995a1257fb3f086eeb65982fac18985a286f301e77c451154ce9ac8895d9',
|
|
729
|
-
'0x1630c3250d7313ff01d1201bf7a74ab5db3cb17dd952799b9ed3ab9097e68f90a0870d2dcae73d19cd13c1c66f652983',
|
|
730
|
-
'0xd6ed6553fe44d296a3726c38ae652bfb11586264f0f8ce19008e218f9c86b2a8da25128c1052ecaddd7f225a139ed84',
|
|
731
|
-
'0x17b81e7701abdbe2e8743884d1117e53356de5ab275b4db1a682c62ef0f2753339b7c8f8c8f475af9ccb5618e3f0c88e',
|
|
732
|
-
'0x80d3cf1f9a78fc47b90b33563be990dc43b756ce79f5574a2c596c928c5d1de4fa295f296b74e956d71986a8497e317',
|
|
733
|
-
'0x169b1f8e1bcfa7c42e0c37515d138f22dd2ecb803a0c5c99676314baf4bb1b7fa3190b2edc0327797f241067be390c9e',
|
|
734
|
-
'0x10321da079ce07e272d8ec09d2565b0dfa7dccdde6787f96d50af36003b14866f69b771f8c285decca67df3f1605fb7b',
|
|
735
|
-
'0x6e08c248e260e70bd1e962381edee3d31d79d7e22c837bc23c0bf1bc24c6b68c24b1b80b64d391fa9c8ba2e8ba2d229',
|
|
736
|
-
],
|
|
737
|
-
// xDen
|
|
738
|
-
[
|
|
739
|
-
'0x8ca8d548cff19ae18b2e62f4bd3fa6f01d5ef4ba35b48ba9c9588617fc8ac62b558d681be343df8993cf9fa40d21b1c',
|
|
740
|
-
'0x12561a5deb559c4348b4711298e536367041e8ca0cf0800c0126c2588c48bf5713daa8846cb026e9e5c8276ec82b3bff',
|
|
741
|
-
'0xb2962fe57a3225e8137e629bff2991f6f89416f5a718cd1fca64e00b11aceacd6a3d0967c94fedcfcc239ba5cb83e19',
|
|
742
|
-
'0x3425581a58ae2fec83aafef7c40eb545b08243f16b1655154cca8abc28d6fd04976d5243eecf5c4130de8938dc62cd8',
|
|
743
|
-
'0x13a8e162022914a80a6f1d5f43e7a07dffdfc759a12062bb8d6b44e833b306da9bd29ba81f35781d539d395b3532a21e',
|
|
744
|
-
'0xe7355f8e4e667b955390f7f0506c6e9395735e9ce9cad4d0a43bcef24b8982f7400d24bc4228f11c02df9a29f6304a5',
|
|
745
|
-
'0x772caacf16936190f3e0c63e0596721570f5799af53a1894e2e073062aede9cea73b3538f0de06cec2574496ee84a3a',
|
|
746
|
-
'0x14a7ac2a9d64a8b230b3f5b074cf01996e7f63c21bca68a81996e1cdf9822c580fa5b9489d11e2d311f7d99bbdcc5a5e',
|
|
747
|
-
'0xa10ecf6ada54f825e920b3dafc7a3cce07f8d1d7161366b74100da67f39883503826692abba43704776ec3a79a1d641',
|
|
748
|
-
'0x95fc13ab9e92ad4476d6e3eb3a56680f682b4ee96f7d03776df533978f31c1593174e4b4b7865002d6384d168ecdd0a',
|
|
749
|
-
'0x000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001', // LAST 1
|
|
750
|
-
],
|
|
751
|
-
// yNum
|
|
752
|
-
[
|
|
753
|
-
'0x90d97c81ba24ee0259d1f094980dcfa11ad138e48a869522b52af6c956543d3cd0c7aee9b3ba3c2be9845719707bb33',
|
|
754
|
-
'0x134996a104ee5811d51036d776fb46831223e96c254f383d0f906343eb67ad34d6c56711962fa8bfe097e75a2e41c696',
|
|
755
|
-
'0xcc786baa966e66f4a384c86a3b49942552e2d658a31ce2c344be4b91400da7d26d521628b00523b8dfe240c72de1f6',
|
|
756
|
-
'0x1f86376e8981c217898751ad8746757d42aa7b90eeb791c09e4a3ec03251cf9de405aba9ec61deca6355c77b0e5f4cb',
|
|
757
|
-
'0x8cc03fdefe0ff135caf4fe2a21529c4195536fbe3ce50b879833fd221351adc2ee7f8dc099040a841b6daecf2e8fedb',
|
|
758
|
-
'0x16603fca40634b6a2211e11db8f0a6a074a7d0d4afadb7bd76505c3d3ad5544e203f6326c95a807299b23ab13633a5f0',
|
|
759
|
-
'0x4ab0b9bcfac1bbcb2c977d027796b3ce75bb8ca2be184cb5231413c4d634f3747a87ac2460f415ec961f8855fe9d6f2',
|
|
760
|
-
'0x987c8d5333ab86fde9926bd2ca6c674170a05bfe3bdd81ffd038da6c26c842642f64550fedfe935a15e4ca31870fb29',
|
|
761
|
-
'0x9fc4018bd96684be88c9e221e4da1bb8f3abd16679dc26c1e8b6e6a1f20cabe69d65201c78607a360370e577bdba587',
|
|
762
|
-
'0xe1bba7a1186bdb5223abde7ada14a23c42a0ca7915af6fe06985e7ed1e4d43b9b3f7055dd4eba6f2bafaaebca731c30',
|
|
763
|
-
'0x19713e47937cd1be0dfd0b8f1d43fb93cd2fcbcb6caf493fd1183e416389e61031bf3a5cce3fbafce813711ad011c132',
|
|
764
|
-
'0x18b46a908f36f6deb918c143fed2edcc523559b8aaf0c2462e6bfe7f911f643249d9cdf41b44d606ce07c8a4d0074d8e',
|
|
765
|
-
'0xb182cac101b9399d155096004f53f447aa7b12a3426b08ec02710e807b4633f06c851c1919211f20d4c04f00b971ef8',
|
|
766
|
-
'0x245a394ad1eca9b72fc00ae7be315dc757b3b080d4c158013e6632d3c40659cc6cf90ad1c232a6442d9d3f5db980133',
|
|
767
|
-
'0x5c129645e44cf1102a159f748c4a3fc5e673d81d7e86568d9ab0f5d396a7ce46ba1049b6579afb7866b1e715475224b',
|
|
768
|
-
'0x15e6be4e990f03ce4ea50b3b42df2eb5cb181d8f84965a3957add4fa95af01b2b665027efec01c7704b456be69c8b604',
|
|
769
|
-
],
|
|
770
|
-
// yDen
|
|
771
|
-
[
|
|
772
|
-
'0x16112c4c3a9c98b252181140fad0eae9601a6de578980be6eec3232b5be72e7a07f3688ef60c206d01479253b03663c1',
|
|
773
|
-
'0x1962d75c2381201e1a0cbd6c43c348b885c84ff731c4d59ca4a10356f453e01f78a4260763529e3532f6102c2e49a03d',
|
|
774
|
-
'0x58df3306640da276faaae7d6e8eb15778c4855551ae7f310c35a5dd279cd2eca6757cd636f96f891e2538b53dbf67f2',
|
|
775
|
-
'0x16b7d288798e5395f20d23bf89edb4d1d115c5dbddbcd30e123da489e726af41727364f2c28297ada8d26d98445f5416',
|
|
776
|
-
'0xbe0e079545f43e4b00cc912f8228ddcc6d19c9f0f69bbb0542eda0fc9dec916a20b15dc0fd2ededda39142311a5001d',
|
|
777
|
-
'0x8d9e5297186db2d9fb266eaac783182b70152c65550d881c5ecd87b6f0f5a6449f38db9dfa9cce202c6477faaf9b7ac',
|
|
778
|
-
'0x166007c08a99db2fc3ba8734ace9824b5eecfdfa8d0cf8ef5dd365bc400a0051d5fa9c01a58b1fb93d1a1399126a775c',
|
|
779
|
-
'0x16a3ef08be3ea7ea03bcddfabba6ff6ee5a4375efa1f4fd7feb34fd206357132b920f5b00801dee460ee415a15812ed9',
|
|
780
|
-
'0x1866c8ed336c61231a1be54fd1d74cc4f9fb0ce4c6af5920abc5750c4bf39b4852cfe2f7bb9248836b233d9d55535d4a',
|
|
781
|
-
'0x167a55cda70a6e1cea820597d94a84903216f763e13d87bb5308592e7ea7d4fbc7385ea3d529b35e346ef48bb8913f55',
|
|
782
|
-
'0x4d2f259eea405bd48f010a01ad2911d9c6dd039bb61a6290e591b36e636a5c871a5c29f4f83060400f8b49cba8f6aa8',
|
|
783
|
-
'0xaccbb67481d033ff5852c1e48c50c477f94ff8aefce42d28c0f9a88cea7913516f968986f7ebbea9684b529e2561092',
|
|
784
|
-
'0xad6b9514c767fe3c3613144b45f1496543346d98adf02267d5ceef9a00d9b8693000763e3b90ac11e99b138573345cc',
|
|
785
|
-
'0x2660400eb2e4f3b628bdd0d53cd76f2bf565b94e72927c1cb748df27942480e420517bd8714cc80d1fadc1326ed06f7',
|
|
786
|
-
'0xe0fa1d816ddc03e6b24255e0d7819c171c40f65e273b853324efcd6356caa205ca2f570f13497804415473a1d634b8f',
|
|
787
|
-
'0x000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000000001', // LAST 1
|
|
788
|
-
],
|
|
789
|
-
].map((i) => i.map((j) => BigInt(j))));
|
|
790
|
-
// SWU Map - Fp2 to G2': y² = x³ + 240i * x + 1012 + 1012i
|
|
791
|
-
const G2_SWU = mapToCurveSimpleSWU(Fp2, {
|
|
792
|
-
A: Fp2.create({ c0: Fp.create(0n), c1: Fp.create(240n) }),
|
|
793
|
-
B: Fp2.create({ c0: Fp.create(1012n), c1: Fp.create(1012n) }),
|
|
794
|
-
Z: Fp2.create({ c0: Fp.create(-2n), c1: Fp.create(-1n) }), // Z: -(2 + I)
|
|
795
|
-
});
|
|
796
|
-
// Optimized SWU Map - Fp to G1
|
|
797
|
-
const G1_SWU = mapToCurveSimpleSWU(Fp, {
|
|
798
|
-
A: Fp.create(0x144698a3b8e9433d693a02c96d4982b0ea985383ee66a8d8e8981aefd881ac98936f8da0e0f97f5cf428082d584c1dn),
|
|
799
|
-
B: Fp.create(0x12e2908d11688030018b12e8753eee3b2016c1f0f24f4070a0b9c14fcef35ef55a23215a316ceaa5d1cc48e98e172be0n),
|
|
800
|
-
Z: Fp.create(11n),
|
|
801
|
-
});
|
|
802
|
-
// Endomorphisms (for fast cofactor clearing)
|
|
803
|
-
// Ψ(P) endomorphism
|
|
804
|
-
const ut_root = Fp6.create({ c0: Fp2.ZERO, c1: Fp2.ONE, c2: Fp2.ZERO });
|
|
805
|
-
const wsq = Fp12.create({ c0: ut_root, c1: Fp6.ZERO });
|
|
806
|
-
const wcu = Fp12.create({ c0: Fp6.ZERO, c1: ut_root });
|
|
807
|
-
const [wsq_inv, wcu_inv] = Fp12.invertBatch([wsq, wcu]);
|
|
808
|
-
function psi(x, y) {
|
|
809
|
-
// Untwist Fp2->Fp12 && frobenius(1) && twist back
|
|
810
|
-
const x2 = Fp12.mul(Fp12.frobeniusMap(Fp12.multiplyByFp2(wsq_inv, x), 1), wsq).c0.c0;
|
|
811
|
-
const y2 = Fp12.mul(Fp12.frobeniusMap(Fp12.multiplyByFp2(wcu_inv, y), 1), wcu).c0.c0;
|
|
812
|
-
return [x2, y2];
|
|
813
|
-
}
|
|
814
|
-
// Ψ endomorphism
|
|
815
|
-
function G2psi(c, P) {
|
|
816
|
-
const affine = P.toAffine();
|
|
817
|
-
const p = psi(affine.x, affine.y);
|
|
818
|
-
return new c(p[0], p[1], Fp2.ONE);
|
|
819
|
-
}
|
|
820
|
-
// Ψ²(P) endomorphism
|
|
821
|
-
// 1 / F2(2)^((p-1)/3) in GF(p²)
|
|
822
|
-
const PSI2_C1 = 0x1a0111ea397fe699ec02408663d4de85aa0d857d89759ad4897d29650fb85f9b409427eb4f49fffd8bfd00000000aaacn;
|
|
823
|
-
function psi2(x, y) {
|
|
824
|
-
return [Fp2.mul(x, PSI2_C1), Fp2.neg(y)];
|
|
825
|
-
}
|
|
826
|
-
function G2psi2(c, P) {
|
|
827
|
-
const affine = P.toAffine();
|
|
828
|
-
const p = psi2(affine.x, affine.y);
|
|
829
|
-
return new c(p[0], p[1], Fp2.ONE);
|
|
830
|
-
}
|
|
831
|
-
// Default hash_to_field options are for hash to G2.
|
|
832
|
-
//
|
|
833
|
-
// Parameter definitions are in section 5.3 of the spec unless otherwise noted.
|
|
834
|
-
// Parameter values come from section 8.8.2 of the spec.
|
|
835
|
-
// https://datatracker.ietf.org/doc/html/draft-irtf-cfrg-hash-to-curve-11#section-8.8.2
|
|
836
|
-
//
|
|
837
|
-
// Base field F is GF(p^m)
|
|
838
|
-
// p = 0x1a0111ea397fe69a4b1ba7b6434bacd764774b84f38512bf6730d2a0f6b0f6241eabfffeb153ffffb9feffffffffaaab
|
|
839
|
-
// m = 2 (or 1 for G1 see section 8.8.1)
|
|
840
|
-
// k = 128
|
|
841
|
-
const htfDefaults = Object.freeze({
|
|
842
|
-
// DST: a domain separation tag
|
|
843
|
-
// defined in section 2.2.5
|
|
844
|
-
// Use utils.getDSTLabel(), utils.setDSTLabel(value)
|
|
845
|
-
DST: 'BLS_SIG_BLS12381G2_XMD:SHA-256_SSWU_RO_NUL_',
|
|
846
|
-
encodeDST: 'BLS_SIG_BLS12381G2_XMD:SHA-256_SSWU_RO_NUL_',
|
|
847
|
-
// p: the characteristic of F
|
|
848
|
-
// where F is a finite field of characteristic p and order q = p^m
|
|
849
|
-
p: Fp.ORDER,
|
|
850
|
-
// m: the extension degree of F, m >= 1
|
|
851
|
-
// where F is a finite field of characteristic p and order q = p^m
|
|
852
|
-
m: 2,
|
|
853
|
-
// k: the target security level for the suite in bits
|
|
854
|
-
// defined in section 5.1
|
|
855
|
-
k: 128,
|
|
856
|
-
// option to use a message that has already been processed by
|
|
857
|
-
// expand_message_xmd
|
|
858
|
-
expand: 'xmd',
|
|
859
|
-
// Hash functions for: expand_message_xmd is appropriate for use with a
|
|
860
|
-
// wide range of hash functions, including SHA-2, SHA-3, BLAKE2, and others.
|
|
861
|
-
// BBS+ uses blake2: https://github.com/hyperledger/aries-framework-go/issues/2247
|
|
862
|
-
hash: sha256,
|
|
863
|
-
});
|
|
864
|
-
// Encoding utils
|
|
865
|
-
// Point on G1 curve: (x, y)
|
|
866
|
-
const C_BIT_POS = Fp.BITS; // C_bit, compression bit for serialization flag
|
|
867
|
-
const I_BIT_POS = Fp.BITS + 1; // I_bit, point-at-infinity bit for serialization flag
|
|
868
|
-
const S_BIT_POS = Fp.BITS + 2; // S_bit, sign bit for serialization flag
|
|
869
|
-
// Compressed point of infinity
|
|
870
|
-
const COMPRESSED_ZERO = Fp.toBytes(bitSet(bitSet(0n, I_BIT_POS, true), S_BIT_POS, true)); // set compressed & point-at-infinity bits
|
|
871
|
-
// To verify curve parameters, see pairing-friendly-curves spec:
|
|
872
|
-
// https://datatracker.ietf.org/doc/html/draft-irtf-cfrg-pairing-friendly-curves-09
|
|
873
|
-
// Basic math is done over finite fields over p.
|
|
874
|
-
// More complicated math is done over polynominal extension fields.
|
|
875
|
-
// To simplify calculations in Fp12, we construct extension tower:
|
|
876
|
-
// Fp₁₂ = Fp₆² => Fp₂³
|
|
877
|
-
// Fp(u) / (u² - β) where β = -1
|
|
878
|
-
// Fp₂(v) / (v³ - ξ) where ξ = u + 1
|
|
879
|
-
// Fp₆(w) / (w² - γ) where γ = v
|
|
880
|
-
// Here goes constants && point encoding format
|
|
881
|
-
export const bls12_381 = bls({
|
|
882
|
-
// Fields
|
|
883
|
-
Fr,
|
|
884
|
-
Fp,
|
|
885
|
-
Fp2,
|
|
886
|
-
Fp6,
|
|
887
|
-
Fp12,
|
|
888
|
-
// order; z⁴ − z² + 1
|
|
889
|
-
r: Fr.ORDER,
|
|
890
|
-
// G1 is the order-q subgroup of E1(Fp) : y² = x³ + 4, #E1(Fp) = h1q, where
|
|
891
|
-
// characteristic; z + (z⁴ - z² + 1)(z - 1)²/3
|
|
892
|
-
G1: {
|
|
893
|
-
Fp,
|
|
894
|
-
// cofactor; (z - 1)²/3
|
|
895
|
-
h: 0x396c8c005555e1568c00aaab0000aaabn,
|
|
896
|
-
// generator's coordinates
|
|
897
|
-
// x = 3685416753713387016781088315183077757961620795782546409894578378688607592378376318836054947676345821548104185464507
|
|
898
|
-
// y = 1339506544944476473020471379941921221584933875938349620426543736416511423956333506472724655353366534992391756441569
|
|
899
|
-
Gx: 0x17f1d3a73197d7942695638c4fa9ac0fc3688c4f9774b905a14e3a3f171bac586c55e83ff97a1aeffb3af00adb22c6bbn,
|
|
900
|
-
Gy: 0x08b3f481e3aaa0f1a09e30ed741d8ae4fcf5e095d5d00af600db18cb2c04b3edd03cc744a2888ae40caa232946c5e7e1n,
|
|
901
|
-
a: Fp.ZERO,
|
|
902
|
-
b: 4n,
|
|
903
|
-
htfDefaults: { ...htfDefaults, m: 1 },
|
|
904
|
-
wrapPrivateKey: true,
|
|
905
|
-
allowInfinityPoint: true,
|
|
906
|
-
// Checks is the point resides in prime-order subgroup.
|
|
907
|
-
// point.isTorsionFree() should return true for valid points
|
|
908
|
-
// It returns false for shitty points.
|
|
909
|
-
// https://eprint.iacr.org/2021/1130.pdf
|
|
910
|
-
isTorsionFree: (c, point) => {
|
|
911
|
-
// φ endomorphism
|
|
912
|
-
const cubicRootOfUnityModP = 0x5f19672fdf76ce51ba69c6076a0f77eaddb3a93be6f89688de17d813620a00022e01fffffffefffen;
|
|
913
|
-
const phi = new c(Fp.mul(point.px, cubicRootOfUnityModP), point.py, point.pz);
|
|
914
|
-
// todo: unroll
|
|
915
|
-
const xP = point.multiplyUnsafe(bls12_381.CURVE.x).negate(); // [x]P
|
|
916
|
-
const u2P = xP.multiplyUnsafe(bls12_381.CURVE.x); // [u2]P
|
|
917
|
-
return u2P.equals(phi);
|
|
918
|
-
// https://eprint.iacr.org/2019/814.pdf
|
|
919
|
-
// (z² − 1)/3
|
|
920
|
-
// const c1 = 0x396c8c005555e1560000000055555555n;
|
|
921
|
-
// const P = this;
|
|
922
|
-
// const S = P.sigma();
|
|
923
|
-
// const Q = S.double();
|
|
924
|
-
// const S2 = S.sigma();
|
|
925
|
-
// // [(z² − 1)/3](2σ(P) − P − σ²(P)) − σ²(P) = O
|
|
926
|
-
// const left = Q.subtract(P).subtract(S2).multiplyUnsafe(c1);
|
|
927
|
-
// const C = left.subtract(S2);
|
|
928
|
-
// return C.isZero();
|
|
929
|
-
},
|
|
930
|
-
// Clear cofactor of G1
|
|
931
|
-
// https://eprint.iacr.org/2019/403
|
|
932
|
-
clearCofactor: (c, point) => {
|
|
933
|
-
// return this.multiplyUnsafe(CURVE.h);
|
|
934
|
-
return point.multiplyUnsafe(bls12_381.CURVE.x).add(point); // x*P + P
|
|
935
|
-
},
|
|
936
|
-
mapToCurve: (scalars) => {
|
|
937
|
-
const { x, y } = G1_SWU(Fp.create(scalars[0]));
|
|
938
|
-
return isogenyMapG1(x, y);
|
|
939
|
-
},
|
|
940
|
-
fromBytes: (bytes) => {
|
|
941
|
-
if (bytes.length === 48) {
|
|
942
|
-
const P = Fp.ORDER;
|
|
943
|
-
const compressedValue = bytesToNumberBE(bytes);
|
|
944
|
-
const bflag = bitGet(compressedValue, I_BIT_POS);
|
|
945
|
-
// Zero
|
|
946
|
-
if (bflag === 1n)
|
|
947
|
-
return { x: 0n, y: 0n };
|
|
948
|
-
const x = Fp.create(compressedValue & Fp.MASK);
|
|
949
|
-
const right = Fp.add(Fp.pow(x, 3n), Fp.create(bls12_381.CURVE.G1.b)); // y² = x³ + b
|
|
950
|
-
let y = Fp.sqrt(right);
|
|
951
|
-
if (!y)
|
|
952
|
-
throw new Error('Invalid compressed G1 point');
|
|
953
|
-
const aflag = bitGet(compressedValue, C_BIT_POS);
|
|
954
|
-
if ((y * 2n) / P !== aflag)
|
|
955
|
-
y = Fp.neg(y);
|
|
956
|
-
return { x: Fp.create(x), y: Fp.create(y) };
|
|
957
|
-
}
|
|
958
|
-
else if (bytes.length === 96) {
|
|
959
|
-
// Check if the infinity flag is set
|
|
960
|
-
if ((bytes[0] & (1 << 6)) !== 0)
|
|
961
|
-
return bls12_381.G1.ProjectivePoint.ZERO.toAffine();
|
|
962
|
-
const x = bytesToNumberBE(bytes.slice(0, Fp.BYTES));
|
|
963
|
-
const y = bytesToNumberBE(bytes.slice(Fp.BYTES));
|
|
964
|
-
return { x: Fp.create(x), y: Fp.create(y) };
|
|
965
|
-
}
|
|
966
|
-
else {
|
|
967
|
-
throw new Error('Invalid point G1, expected 48/96 bytes');
|
|
968
|
-
}
|
|
969
|
-
},
|
|
970
|
-
toBytes: (c, point, isCompressed) => {
|
|
971
|
-
const isZero = point.equals(c.ZERO);
|
|
972
|
-
const { x, y } = point.toAffine();
|
|
973
|
-
if (isCompressed) {
|
|
974
|
-
if (isZero)
|
|
975
|
-
return COMPRESSED_ZERO.slice();
|
|
976
|
-
const P = Fp.ORDER;
|
|
977
|
-
let num;
|
|
978
|
-
num = bitSet(x, C_BIT_POS, Boolean((y * 2n) / P)); // set aflag
|
|
979
|
-
num = bitSet(num, S_BIT_POS, true);
|
|
980
|
-
return numberToBytesBE(num, Fp.BYTES);
|
|
981
|
-
}
|
|
982
|
-
else {
|
|
983
|
-
if (isZero) {
|
|
984
|
-
// 2x PUBLIC_KEY_LENGTH
|
|
985
|
-
const x = concatB(new Uint8Array([0x40]), new Uint8Array(2 * Fp.BYTES - 1));
|
|
986
|
-
return x;
|
|
987
|
-
}
|
|
988
|
-
else {
|
|
989
|
-
return concatB(numberToBytesBE(x, Fp.BYTES), numberToBytesBE(y, Fp.BYTES));
|
|
990
|
-
}
|
|
991
|
-
}
|
|
992
|
-
},
|
|
993
|
-
},
|
|
994
|
-
// G2 is the order-q subgroup of E2(Fp²) : y² = x³+4(1+√−1),
|
|
995
|
-
// where Fp2 is Fp[√−1]/(x2+1). #E2(Fp2 ) = h2q, where
|
|
996
|
-
// G² - 1
|
|
997
|
-
// h2q
|
|
998
|
-
G2: {
|
|
999
|
-
Fp: Fp2,
|
|
1000
|
-
// cofactor
|
|
1001
|
-
h: 0x5d543a95414e7f1091d50792876a202cd91de4547085abaa68a205b2e5a7ddfa628f1cb4d9e82ef21537e293a6691ae1616ec6e786f0c70cf1c38e31c7238e5n,
|
|
1002
|
-
Gx: Fp2.fromBigTuple([
|
|
1003
|
-
0x024aa2b2f08f0a91260805272dc51051c6e47ad4fa403b02b4510b647ae3d1770bac0326a805bbefd48056c8c121bdb8n,
|
|
1004
|
-
0x13e02b6052719f607dacd3a088274f65596bd0d09920b61ab5da61bbdc7f5049334cf11213945d57e5ac7d055d042b7en,
|
|
1005
|
-
]),
|
|
1006
|
-
// y =
|
|
1007
|
-
// 927553665492332455747201965776037880757740193453592970025027978793976877002675564980949289727957565575433344219582,
|
|
1008
|
-
// 1985150602287291935568054521177171638300868978215655730859378665066344726373823718423869104263333984641494340347905
|
|
1009
|
-
Gy: Fp2.fromBigTuple([
|
|
1010
|
-
0x0ce5d527727d6e118cc9cdc6da2e351aadfd9baa8cbdd3a76d429a695160d12c923ac9cc3baca289e193548608b82801n,
|
|
1011
|
-
0x0606c4a02ea734cc32acd2b02bc28b99cb3e287e85a763af267492ab572e99ab3f370d275cec1da1aaa9075ff05f79ben,
|
|
1012
|
-
]),
|
|
1013
|
-
a: Fp2.ZERO,
|
|
1014
|
-
b: Fp2.fromBigTuple([4n, 4n]),
|
|
1015
|
-
hEff: 0xbc69f08f2ee75b3584c6a0ea91b352888e2a8e9145ad7689986ff031508ffe1329c2f178731db956d82bf015d1212b02ec0ec69d7477c1ae954cbc06689f6a359894c0adebbf6b4e8020005aaa95551n,
|
|
1016
|
-
htfDefaults: { ...htfDefaults },
|
|
1017
|
-
wrapPrivateKey: true,
|
|
1018
|
-
allowInfinityPoint: true,
|
|
1019
|
-
mapToCurve: (scalars) => {
|
|
1020
|
-
const { x, y } = G2_SWU(Fp2.fromBigTuple(scalars));
|
|
1021
|
-
return isogenyMapG2(x, y);
|
|
1022
|
-
},
|
|
1023
|
-
// Checks is the point resides in prime-order subgroup.
|
|
1024
|
-
// point.isTorsionFree() should return true for valid points
|
|
1025
|
-
// It returns false for shitty points.
|
|
1026
|
-
// https://eprint.iacr.org/2021/1130.pdf
|
|
1027
|
-
isTorsionFree: (c, P) => {
|
|
1028
|
-
return P.multiplyUnsafe(bls12_381.CURVE.x).negate().equals(G2psi(c, P)); // ψ(P) == [u](P)
|
|
1029
|
-
// Older version: https://eprint.iacr.org/2019/814.pdf
|
|
1030
|
-
// Ψ²(P) => Ψ³(P) => [z]Ψ³(P) where z = -x => [z]Ψ³(P) - Ψ²(P) + P == O
|
|
1031
|
-
// return P.psi2().psi().mulNegX().subtract(psi2).add(P).isZero();
|
|
1032
|
-
},
|
|
1033
|
-
// Maps the point into the prime-order subgroup G2.
|
|
1034
|
-
// clear_cofactor_bls12381_g2 from cfrg-hash-to-curve-11
|
|
1035
|
-
// https://eprint.iacr.org/2017/419.pdf
|
|
1036
|
-
// prettier-ignore
|
|
1037
|
-
clearCofactor: (c, P) => {
|
|
1038
|
-
const { x } = bls12_381.CURVE;
|
|
1039
|
-
let t1 = P.multiplyUnsafe(x).negate(); // [-x]P
|
|
1040
|
-
let t2 = G2psi(c, P); // Ψ(P)
|
|
1041
|
-
let t3 = P.double(); // 2P
|
|
1042
|
-
t3 = G2psi2(c, t3); // Ψ²(2P)
|
|
1043
|
-
t3 = t3.subtract(t2); // Ψ²(2P) - Ψ(P)
|
|
1044
|
-
t2 = t1.add(t2); // [-x]P + Ψ(P)
|
|
1045
|
-
t2 = t2.multiplyUnsafe(x).negate(); // [x²]P - [x]Ψ(P)
|
|
1046
|
-
t3 = t3.add(t2); // Ψ²(2P) - Ψ(P) + [x²]P - [x]Ψ(P)
|
|
1047
|
-
t3 = t3.subtract(t1); // Ψ²(2P) - Ψ(P) + [x²]P - [x]Ψ(P) + [x]P
|
|
1048
|
-
const Q = t3.subtract(P); // Ψ²(2P) - Ψ(P) + [x²]P - [x]Ψ(P) + [x]P - 1P
|
|
1049
|
-
return Q; // [x²-x-1]P + [x-1]Ψ(P) + Ψ²(2P)
|
|
1050
|
-
},
|
|
1051
|
-
fromBytes: (bytes) => {
|
|
1052
|
-
const m_byte = bytes[0] & 0xe0;
|
|
1053
|
-
if (m_byte === 0x20 || m_byte === 0x60 || m_byte === 0xe0) {
|
|
1054
|
-
throw new Error('Invalid encoding flag: ' + m_byte);
|
|
1055
|
-
}
|
|
1056
|
-
const bitC = m_byte & 0x80; // compression bit
|
|
1057
|
-
const bitI = m_byte & 0x40; // point at infinity bit
|
|
1058
|
-
const bitS = m_byte & 0x20; // sign bit
|
|
1059
|
-
const L = Fp.BYTES;
|
|
1060
|
-
const slc = (b, from, to) => bytesToNumberBE(b.slice(from, to));
|
|
1061
|
-
if (bytes.length === 96 && bitC) {
|
|
1062
|
-
const { b } = bls12_381.CURVE.G2;
|
|
1063
|
-
const P = Fp.ORDER;
|
|
1064
|
-
bytes[0] = bytes[0] & 0x1f; // clear flags
|
|
1065
|
-
if (bitI) {
|
|
1066
|
-
// check that all bytes are 0
|
|
1067
|
-
if (bytes.reduce((p, c) => (p !== 0 ? c + 1 : c), 0) > 0) {
|
|
1068
|
-
throw new Error('Invalid compressed G2 point');
|
|
1069
|
-
}
|
|
1070
|
-
return { x: Fp2.ZERO, y: Fp2.ZERO };
|
|
1071
|
-
}
|
|
1072
|
-
const x_1 = slc(bytes, 0, L);
|
|
1073
|
-
const x_0 = slc(bytes, L, 2 * L);
|
|
1074
|
-
const x = Fp2.create({ c0: Fp.create(x_0), c1: Fp.create(x_1) });
|
|
1075
|
-
const right = Fp2.add(Fp2.pow(x, 3n), b); // y² = x³ + 4 * (u+1) = x³ + b
|
|
1076
|
-
let y = Fp2.sqrt(right);
|
|
1077
|
-
const Y_bit = y.c1 === 0n ? (y.c0 * 2n) / P : (y.c1 * 2n) / P ? 1n : 0n;
|
|
1078
|
-
y = bitS > 0 && Y_bit > 0 ? y : Fp2.neg(y);
|
|
1079
|
-
return { x, y };
|
|
1080
|
-
}
|
|
1081
|
-
else if (bytes.length === 192 && !bitC) {
|
|
1082
|
-
// Check if the infinity flag is set
|
|
1083
|
-
if ((bytes[0] & (1 << 6)) !== 0) {
|
|
1084
|
-
return { x: Fp2.ZERO, y: Fp2.ZERO };
|
|
1085
|
-
}
|
|
1086
|
-
const x1 = slc(bytes, 0, L);
|
|
1087
|
-
const x0 = slc(bytes, L, 2 * L);
|
|
1088
|
-
const y1 = slc(bytes, 2 * L, 3 * L);
|
|
1089
|
-
const y0 = slc(bytes, 3 * L, 4 * L);
|
|
1090
|
-
return { x: Fp2.fromBigTuple([x0, x1]), y: Fp2.fromBigTuple([y0, y1]) };
|
|
1091
|
-
}
|
|
1092
|
-
else {
|
|
1093
|
-
throw new Error('Invalid point G2, expected 96/192 bytes');
|
|
1094
|
-
}
|
|
1095
|
-
},
|
|
1096
|
-
toBytes: (c, point, isCompressed) => {
|
|
1097
|
-
const isZero = point.equals(c.ZERO);
|
|
1098
|
-
const { x, y } = point.toAffine();
|
|
1099
|
-
if (isCompressed) {
|
|
1100
|
-
const P = Fp.ORDER;
|
|
1101
|
-
if (isZero)
|
|
1102
|
-
return concatB(COMPRESSED_ZERO, numberToBytesBE(0n, Fp.BYTES));
|
|
1103
|
-
const flag = Boolean(y.c1 === 0n ? (y.c0 * 2n) / P : (y.c1 * 2n) / P);
|
|
1104
|
-
// set compressed & sign bits (looks like different offsets than for G1/Fp?)
|
|
1105
|
-
let x_1 = bitSet(x.c1, C_BIT_POS, flag);
|
|
1106
|
-
x_1 = bitSet(x_1, S_BIT_POS, true);
|
|
1107
|
-
return concatB(numberToBytesBE(x_1, Fp.BYTES), numberToBytesBE(x.c0, Fp.BYTES));
|
|
1108
|
-
}
|
|
1109
|
-
else {
|
|
1110
|
-
if (isZero)
|
|
1111
|
-
return concatB(new Uint8Array([0x40]), new Uint8Array(4 * Fp.BYTES - 1)); // bytes[0] |= 1 << 6;
|
|
1112
|
-
const { re: x0, im: x1 } = Fp2.reim(x);
|
|
1113
|
-
const { re: y0, im: y1 } = Fp2.reim(y);
|
|
1114
|
-
return concatB(numberToBytesBE(x1, Fp.BYTES), numberToBytesBE(x0, Fp.BYTES), numberToBytesBE(y1, Fp.BYTES), numberToBytesBE(y0, Fp.BYTES));
|
|
1115
|
-
}
|
|
1116
|
-
},
|
|
1117
|
-
Signature: {
|
|
1118
|
-
// TODO: Optimize, it's very slow because of sqrt.
|
|
1119
|
-
decode(hex) {
|
|
1120
|
-
hex = ensureBytes('signatureHex', hex);
|
|
1121
|
-
const P = Fp.ORDER;
|
|
1122
|
-
const half = hex.length / 2;
|
|
1123
|
-
if (half !== 48 && half !== 96)
|
|
1124
|
-
throw new Error('Invalid compressed signature length, must be 96 or 192');
|
|
1125
|
-
const z1 = bytesToNumberBE(hex.slice(0, half));
|
|
1126
|
-
const z2 = bytesToNumberBE(hex.slice(half));
|
|
1127
|
-
// Indicates the infinity point
|
|
1128
|
-
const bflag1 = bitGet(z1, I_BIT_POS);
|
|
1129
|
-
if (bflag1 === 1n)
|
|
1130
|
-
return bls12_381.G2.ProjectivePoint.ZERO;
|
|
1131
|
-
const x1 = Fp.create(z1 & Fp.MASK);
|
|
1132
|
-
const x2 = Fp.create(z2);
|
|
1133
|
-
const x = Fp2.create({ c0: x2, c1: x1 });
|
|
1134
|
-
const y2 = Fp2.add(Fp2.pow(x, 3n), bls12_381.CURVE.G2.b); // y² = x³ + 4
|
|
1135
|
-
// The slow part
|
|
1136
|
-
let y = Fp2.sqrt(y2);
|
|
1137
|
-
if (!y)
|
|
1138
|
-
throw new Error('Failed to find a square root');
|
|
1139
|
-
// Choose the y whose leftmost bit of the imaginary part is equal to the a_flag1
|
|
1140
|
-
// If y1 happens to be zero, then use the bit of y0
|
|
1141
|
-
const { re: y0, im: y1 } = Fp2.reim(y);
|
|
1142
|
-
const aflag1 = bitGet(z1, 381);
|
|
1143
|
-
const isGreater = y1 > 0n && (y1 * 2n) / P !== aflag1;
|
|
1144
|
-
const isZero = y1 === 0n && (y0 * 2n) / P !== aflag1;
|
|
1145
|
-
if (isGreater || isZero)
|
|
1146
|
-
y = Fp2.neg(y);
|
|
1147
|
-
const point = bls12_381.G2.ProjectivePoint.fromAffine({ x, y });
|
|
1148
|
-
point.assertValidity();
|
|
1149
|
-
return point;
|
|
1150
|
-
},
|
|
1151
|
-
encode(point) {
|
|
1152
|
-
// NOTE: by some reasons it was missed in bls12-381, looks like bug
|
|
1153
|
-
point.assertValidity();
|
|
1154
|
-
if (point.equals(bls12_381.G2.ProjectivePoint.ZERO))
|
|
1155
|
-
return concatB(COMPRESSED_ZERO, numberToBytesBE(0n, Fp.BYTES));
|
|
1156
|
-
const a = point.toAffine();
|
|
1157
|
-
const { re: x0, im: x1 } = Fp2.reim(a.x);
|
|
1158
|
-
const { re: y0, im: y1 } = Fp2.reim(a.y);
|
|
1159
|
-
const tmp = y1 > 0n ? y1 * 2n : y0 * 2n;
|
|
1160
|
-
const aflag1 = Boolean((tmp / Fp.ORDER) & 1n);
|
|
1161
|
-
const z1 = bitSet(bitSet(x1, 381, aflag1), S_BIT_POS, true);
|
|
1162
|
-
const z2 = x0;
|
|
1163
|
-
return concatB(numberToBytesBE(z1, Fp.BYTES), numberToBytesBE(z2, Fp.BYTES));
|
|
1164
|
-
},
|
|
1165
|
-
},
|
|
1166
|
-
},
|
|
1167
|
-
// The BLS parameter x for BLS12-381
|
|
1168
|
-
x: BLS_X,
|
|
1169
|
-
htfDefaults,
|
|
1170
|
-
hash: sha256,
|
|
1171
|
-
randomBytes,
|
|
1172
|
-
});
|
|
1173
|
-
//# sourceMappingURL=bls12-381.js.map
|