public override ECPoint Add(ECPoint b)
    {
        if (base.IsInfinity)
        {
            return(b);
        }
        if (b.IsInfinity)
        {
            return(this);
        }
        if (this == b)
        {
            return(Twice());
        }
        ECCurve curve = Curve;
        SecP160R1FieldElement secP160R1FieldElement  = (SecP160R1FieldElement)base.RawXCoord;
        SecP160R1FieldElement secP160R1FieldElement2 = (SecP160R1FieldElement)base.RawYCoord;
        SecP160R1FieldElement secP160R1FieldElement3 = (SecP160R1FieldElement)b.RawXCoord;
        SecP160R1FieldElement secP160R1FieldElement4 = (SecP160R1FieldElement)b.RawYCoord;
        SecP160R1FieldElement secP160R1FieldElement5 = (SecP160R1FieldElement)base.RawZCoords[0];
        SecP160R1FieldElement secP160R1FieldElement6 = (SecP160R1FieldElement)b.RawZCoords[0];

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

        uint[] array5;
        uint[] array6;
        if (isOne)
        {
            array5 = secP160R1FieldElement3.x;
            array6 = secP160R1FieldElement4.x;
        }
        else
        {
            array6 = array3;
            SecP160R1Field.Square(secP160R1FieldElement5.x, array6);
            array5 = array2;
            SecP160R1Field.Multiply(array6, secP160R1FieldElement3.x, array5);
            SecP160R1Field.Multiply(array6, secP160R1FieldElement5.x, array6);
            SecP160R1Field.Multiply(array6, secP160R1FieldElement4.x, array6);
        }
        bool isOne2 = secP160R1FieldElement6.IsOne;

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

        SecP160R1Field.Reduce32(x, array12);
        SecP160R1FieldElement secP160R1FieldElement7 = new SecP160R1FieldElement(array4);

        SecP160R1Field.Square(array10, secP160R1FieldElement7.x);
        SecP160R1Field.Subtract(secP160R1FieldElement7.x, array12, secP160R1FieldElement7.x);
        SecP160R1FieldElement secP160R1FieldElement8 = new SecP160R1FieldElement(array12);

        SecP160R1Field.Subtract(array13, secP160R1FieldElement7.x, secP160R1FieldElement8.x);
        SecP160R1Field.MultiplyAddToExt(secP160R1FieldElement8.x, array10, array);
        SecP160R1Field.Reduce(array, secP160R1FieldElement8.x);
        SecP160R1FieldElement secP160R1FieldElement9 = new SecP160R1FieldElement(array9);

        if (!isOne)
        {
            SecP160R1Field.Multiply(secP160R1FieldElement9.x, secP160R1FieldElement5.x, secP160R1FieldElement9.x);
        }
        if (!isOne2)
        {
            SecP160R1Field.Multiply(secP160R1FieldElement9.x, secP160R1FieldElement6.x, secP160R1FieldElement9.x);
        }
        ECFieldElement[] zs = new ECFieldElement[1]
        {
            secP160R1FieldElement9
        };
        return(new SecP160R1Point(curve, secP160R1FieldElement7, secP160R1FieldElement8, zs, base.IsCompressed));
    }