public override ECFieldElement Add(ECFieldElement b)
 {
     uint[] z = Nat.Create(12);
     SecP384R1Field.Add(x, ((SecP384R1FieldElement)b).x, z);
     return(new SecP384R1FieldElement(z));
 }
        public override ECPoint Twice()
        {
            if (this.IsInfinity)
            {
                return(this);
            }

            ECCurve curve = this.Curve;

            SecP384R1FieldElement Y1 = (SecP384R1FieldElement)this.RawYCoord;

            if (Y1.IsZero)
            {
                return(curve.Infinity);
            }

            SecP384R1FieldElement X1 = (SecP384R1FieldElement)this.RawXCoord, Z1 = (SecP384R1FieldElement)this.RawZCoords[0];

            uint c;

            uint[] t1 = Nat.Create(12);
            uint[] t2 = Nat.Create(12);

            uint[] Y1Squared = Nat.Create(12);
            SecP384R1Field.Square(Y1.x, Y1Squared);

            uint[] T = Nat.Create(12);
            SecP384R1Field.Square(Y1Squared, T);

            bool Z1IsOne = Z1.IsOne;

            uint[] Z1Squared = Z1.x;
            if (!Z1IsOne)
            {
                Z1Squared = t2;
                SecP384R1Field.Square(Z1.x, Z1Squared);
            }

            SecP384R1Field.Subtract(X1.x, Z1Squared, t1);

            uint[] M = t2;
            SecP384R1Field.Add(X1.x, Z1Squared, M);
            SecP384R1Field.Multiply(M, t1, M);
            c = Nat.AddBothTo(12, M, M, M);
            SecP384R1Field.Reduce32(c, M);

            uint[] S = Y1Squared;
            SecP384R1Field.Multiply(Y1Squared, X1.x, S);
            c = Nat.ShiftUpBits(12, S, 2, 0);
            SecP384R1Field.Reduce32(c, S);

            c = Nat.ShiftUpBits(12, T, 3, 0, t1);
            SecP384R1Field.Reduce32(c, t1);

            SecP384R1FieldElement X3 = new SecP384R1FieldElement(T);

            SecP384R1Field.Square(M, X3.x);
            SecP384R1Field.Subtract(X3.x, S, X3.x);
            SecP384R1Field.Subtract(X3.x, S, X3.x);

            SecP384R1FieldElement Y3 = new SecP384R1FieldElement(S);

            SecP384R1Field.Subtract(S, X3.x, Y3.x);
            SecP384R1Field.Multiply(Y3.x, M, Y3.x);
            SecP384R1Field.Subtract(Y3.x, t1, Y3.x);

            SecP384R1FieldElement Z3 = new SecP384R1FieldElement(M);

            SecP384R1Field.Twice(Y1.x, Z3.x);
            if (!Z1IsOne)
            {
                SecP384R1Field.Multiply(Z3.x, Z1.x, Z3.x);
            }

            return(new SecP384R1Point(curve, X3, Y3, new ECFieldElement[] { Z3 }, IsCompressed));
        }