private static Fp12 MillerLoop(Bn128Fp g1, Bn128Fp2 g2)
{
// convert to affine coordinates
g1 = g1.ToAffine();
g2 = g2.ToAffine();
// calculate Ell coefficients
List <EllCoeffs> coeffs = CalcEllCoeffs(g2);
Fp12 f = Fp12.One;
int idx = 0;
// for each bit except most significant one
for (int i = LoopCount.BitLength() - 2; i >= 0; i--)
{
EllCoeffs cInLoop = coeffs[idx++];
f = f.Squared();
f = f.MulBy024(cInLoop.Ell0, g1.Y.Mul(cInLoop.EllVw), g1.X.Mul(cInLoop.EllVv));
if (LoopCount.TestBit(i))
{
cInLoop = coeffs[idx++];
f = f.MulBy024(cInLoop.Ell0, g1.Y.Mul(cInLoop.EllVw), g1.X.Mul(cInLoop.EllVv));
}
}
EllCoeffs c = coeffs[idx++];
f = f.MulBy024(c.Ell0, g1.Y.Mul(c.EllVw), g1.X.Mul(c.EllVv));
c = coeffs[idx];
f = f.MulBy024(c.Ell0, g1.Y.Mul(c.EllVw), g1.X.Mul(c.EllVv));
return(f);
}
public void AddPair(Bn128Fp g1, Bn128Fp2 g2)
{
_pairs.Add(new Pair(g1, g2));
}
public Pair(Bn128Fp g1, Bn128Fp2 g2)
{
G1 = g1;
G2 = g2;
}