/* return TRUE if this==a */
public bool Equals(FP a)
{
FP f = new FP(this);
FP s = new FP(a);
f.Reduce();
s.Reduce();
if (BIG.Comp(f.x, s.x) == 0)
{
return(true);
}
return(false);
}
/* r=ck^a.cl^n using XTR double exponentiation method on traces of FP12s. See Stam thesis. */
public FP4 Xtr_Pow2(FP4 ck, FP4 ckml, FP4 ckm2l, BIG a, BIG b)
{
a.Norm();
b.Norm();
BIG e = new BIG(a);
BIG d = new BIG(b);
BIG w = new BIG(0);
FP4 cu = new FP4(ck); // can probably be passed in w/o copying
FP4 cv = new FP4(this);
FP4 cumv = new FP4(ckml);
FP4 cum2v = new FP4(ckm2l);
FP4 r = new FP4(0);
FP4 t = new FP4(0);
int f2 = 0;
while (d.Parity() == 0 && e.Parity() == 0)
{
d.FShr(1);
e.FShr(1);
f2++;
}
while (BIG.Comp(d, e) != 0)
{
if (BIG.Comp(d, e) > 0)
{
w.Copy(e);
w.IMul(4);
w.Norm();
if (BIG.Comp(d, w) <= 0)
{
w.Copy(d);
d.Copy(e);
e.RSub(w);
e.Norm();
t.Copy(cv);
t.Xtr_A(cu, cumv, cum2v);
cum2v.Copy(cumv);
cum2v.Conj();
cumv.Copy(cv);
cv.Copy(cu);
cu.Copy(t);
}
else if (d.Parity() == 0)
{
d.FShr(1);
r.Copy(cum2v);
r.Conj();
t.Copy(cumv);
t.Xtr_A(cu, cv, r);
cum2v.Copy(cumv);
cum2v.Xtr_D();
cumv.Copy(t);
cu.Xtr_D();
}
else if (e.Parity() == 1)
{
d.Sub(e);
d.Norm();
d.FShr(1);
t.Copy(cv);
t.Xtr_A(cu, cumv, cum2v);
cu.Xtr_D();
cum2v.Copy(cv);
cum2v.Xtr_D();
cum2v.Conj();
cv.Copy(t);
}
else
{
w.Copy(d);
d.Copy(e);
d.FShr(1);
e.Copy(w);
t.Copy(cumv);
t.Xtr_D();
cumv.Copy(cum2v);
cumv.Conj();
cum2v.Copy(t);
cum2v.Conj();
t.Copy(cv);
t.Xtr_D();
cv.Copy(cu);
cu.Copy(t);
}
}
if (BIG.Comp(d, e) < 0)
{
w.Copy(d);
w.IMul(4);
w.Norm();
if (BIG.Comp(e, w) <= 0)
{
e.Sub(d);
e.Norm();
t.Copy(cv);
t.Xtr_A(cu, cumv, cum2v);
cum2v.Copy(cumv);
cumv.Copy(cu);
cu.Copy(t);
}
else if (e.Parity() == 0)
{
w.Copy(d);
d.Copy(e);
d.FShr(1);
e.Copy(w);
t.Copy(cumv);
t.Xtr_D();
cumv.Copy(cum2v);
cumv.Conj();
cum2v.Copy(t);
cum2v.Conj();
t.Copy(cv);
t.Xtr_D();
cv.Copy(cu);
cu.Copy(t);
}
else if (d.Parity() == 1)
{
w.Copy(e);
e.Copy(d);
w.Sub(d);
w.Norm();
d.Copy(w);
d.FShr(1);
t.Copy(cv);
t.Xtr_A(cu, cumv, cum2v);
cumv.Conj();
cum2v.Copy(cu);
cum2v.Xtr_D();
cum2v.Conj();
cu.Copy(cv);
cu.Xtr_D();
cv.Copy(t);
}
else
{
d.FShr(1);
r.Copy(cum2v);
r.Conj();
t.Copy(cumv);
t.Xtr_A(cu, cv, r);
cum2v.Copy(cumv);
cum2v.Xtr_D();
cumv.Copy(t);
cu.Xtr_D();
}
}
}
r.Copy(cv);
r.Xtr_A(cu, cumv, cum2v);
for (int i = 0; i < f2; i++)
{
r.Xtr_D();
}
r = r.Xtr_Pow(d);
return(r);
}
