public override ECPoint Twice()
    {
        if (base.IsInfinity)
        {
            return(this);
        }
        ECCurve curve = Curve;
        SecP160R2FieldElement secP160R2FieldElement = (SecP160R2FieldElement)base.RawYCoord;

        if (secP160R2FieldElement.IsZero)
        {
            return(curve.Infinity);
        }
        SecP160R2FieldElement secP160R2FieldElement2 = (SecP160R2FieldElement)base.RawXCoord;
        SecP160R2FieldElement secP160R2FieldElement3 = (SecP160R2FieldElement)base.RawZCoords[0];

        uint[] array = Nat160.Create();
        SecP160R2Field.Square(secP160R2FieldElement.x, array);
        uint[] array2 = Nat160.Create();
        SecP160R2Field.Square(array, array2);
        uint[] array3 = Nat160.Create();
        SecP160R2Field.Square(secP160R2FieldElement2.x, array3);
        uint x = Nat160.AddBothTo(array3, array3, array3);

        SecP160R2Field.Reduce32(x, array3);
        uint[] array4 = array;
        SecP160R2Field.Multiply(array, secP160R2FieldElement2.x, array4);
        x = Nat.ShiftUpBits(5, array4, 2, 0u);
        SecP160R2Field.Reduce32(x, array4);
        uint[] array5 = Nat160.Create();
        x = Nat.ShiftUpBits(5, array2, 3, 0u, array5);
        SecP160R2Field.Reduce32(x, array5);
        SecP160R2FieldElement secP160R2FieldElement4 = new SecP160R2FieldElement(array2);

        SecP160R2Field.Square(array3, secP160R2FieldElement4.x);
        SecP160R2Field.Subtract(secP160R2FieldElement4.x, array4, secP160R2FieldElement4.x);
        SecP160R2Field.Subtract(secP160R2FieldElement4.x, array4, secP160R2FieldElement4.x);
        SecP160R2FieldElement secP160R2FieldElement5 = new SecP160R2FieldElement(array4);

        SecP160R2Field.Subtract(array4, secP160R2FieldElement4.x, secP160R2FieldElement5.x);
        SecP160R2Field.Multiply(secP160R2FieldElement5.x, array3, secP160R2FieldElement5.x);
        SecP160R2Field.Subtract(secP160R2FieldElement5.x, array5, secP160R2FieldElement5.x);
        SecP160R2FieldElement secP160R2FieldElement6 = new SecP160R2FieldElement(array3);

        SecP160R2Field.Twice(secP160R2FieldElement.x, secP160R2FieldElement6.x);
        if (!secP160R2FieldElement3.IsOne)
        {
            SecP160R2Field.Multiply(secP160R2FieldElement6.x, secP160R2FieldElement3.x, secP160R2FieldElement6.x);
        }
        return(new SecP160K1Point(curve, secP160R2FieldElement4, secP160R2FieldElement5, new ECFieldElement[1]
        {
            secP160R2FieldElement6
        }, base.IsCompressed));
    }
Beispiel #2
0
 public override ECFieldElement Subtract(ECFieldElement b)
 {
     uint[] z = Nat160.Create();
     SecP160R2Field.Subtract(x, ((SecP160R2FieldElement)b).x, z);
     return(new SecP160R2FieldElement(z));
 }
    public override ECPoint Add(ECPoint b)
    {
        if (base.IsInfinity)
        {
            return(b);
        }
        if (b.IsInfinity)
        {
            return(this);
        }
        if (this == b)
        {
            return(Twice());
        }
        ECCurve curve = Curve;
        SecP160R2FieldElement secP160R2FieldElement  = (SecP160R2FieldElement)base.RawXCoord;
        SecP160R2FieldElement secP160R2FieldElement2 = (SecP160R2FieldElement)base.RawYCoord;
        SecP160R2FieldElement secP160R2FieldElement3 = (SecP160R2FieldElement)b.RawXCoord;
        SecP160R2FieldElement secP160R2FieldElement4 = (SecP160R2FieldElement)b.RawYCoord;
        SecP160R2FieldElement secP160R2FieldElement5 = (SecP160R2FieldElement)base.RawZCoords[0];
        SecP160R2FieldElement secP160R2FieldElement6 = (SecP160R2FieldElement)b.RawZCoords[0];

        uint[] array  = Nat160.CreateExt();
        uint[] array2 = Nat160.Create();
        uint[] array3 = Nat160.Create();
        uint[] array4 = Nat160.Create();
        bool   isOne  = secP160R2FieldElement5.IsOne;

        uint[] array5;
        uint[] array6;
        if (isOne)
        {
            array5 = secP160R2FieldElement3.x;
            array6 = secP160R2FieldElement4.x;
        }
        else
        {
            array6 = array3;
            SecP160R2Field.Square(secP160R2FieldElement5.x, array6);
            array5 = array2;
            SecP160R2Field.Multiply(array6, secP160R2FieldElement3.x, array5);
            SecP160R2Field.Multiply(array6, secP160R2FieldElement5.x, array6);
            SecP160R2Field.Multiply(array6, secP160R2FieldElement4.x, array6);
        }
        bool isOne2 = secP160R2FieldElement6.IsOne;

        uint[] array7;
        uint[] array8;
        if (isOne2)
        {
            array7 = secP160R2FieldElement.x;
            array8 = secP160R2FieldElement2.x;
        }
        else
        {
            array8 = array4;
            SecP160R2Field.Square(secP160R2FieldElement6.x, array8);
            array7 = array;
            SecP160R2Field.Multiply(array8, secP160R2FieldElement.x, array7);
            SecP160R2Field.Multiply(array8, secP160R2FieldElement6.x, array8);
            SecP160R2Field.Multiply(array8, secP160R2FieldElement2.x, array8);
        }
        uint[] array9 = Nat160.Create();
        SecP160R2Field.Subtract(array7, array5, array9);
        uint[] array10 = array2;
        SecP160R2Field.Subtract(array8, array6, array10);
        if (Nat160.IsZero(array9))
        {
            if (Nat160.IsZero(array10))
            {
                return(Twice());
            }
            return(curve.Infinity);
        }
        uint[] array11 = array3;
        SecP160R2Field.Square(array9, array11);
        uint[] array12 = Nat160.Create();
        SecP160R2Field.Multiply(array11, array9, array12);
        uint[] array13 = array3;
        SecP160R2Field.Multiply(array11, array7, array13);
        SecP160R2Field.Negate(array12, array12);
        Nat160.Mul(array8, array12, array);
        uint x = Nat160.AddBothTo(array13, array13, array12);

        SecP160R2Field.Reduce32(x, array12);
        SecP160R2FieldElement secP160R2FieldElement7 = new SecP160R2FieldElement(array4);

        SecP160R2Field.Square(array10, secP160R2FieldElement7.x);
        SecP160R2Field.Subtract(secP160R2FieldElement7.x, array12, secP160R2FieldElement7.x);
        SecP160R2FieldElement secP160R2FieldElement8 = new SecP160R2FieldElement(array12);

        SecP160R2Field.Subtract(array13, secP160R2FieldElement7.x, secP160R2FieldElement8.x);
        SecP160R2Field.MultiplyAddToExt(secP160R2FieldElement8.x, array10, array);
        SecP160R2Field.Reduce(array, secP160R2FieldElement8.x);
        SecP160R2FieldElement secP160R2FieldElement9 = new SecP160R2FieldElement(array9);

        if (!isOne)
        {
            SecP160R2Field.Multiply(secP160R2FieldElement9.x, secP160R2FieldElement5.x, secP160R2FieldElement9.x);
        }
        if (!isOne2)
        {
            SecP160R2Field.Multiply(secP160R2FieldElement9.x, secP160R2FieldElement6.x, secP160R2FieldElement9.x);
        }
        ECFieldElement[] zs = new ECFieldElement[1]
        {
            secP160R2FieldElement9
        };
        return(new SecP160K1Point(curve, secP160R2FieldElement7, secP160R2FieldElement8, zs, base.IsCompressed));
    }