using Org.BouncyCastle.Math.Raw; namespace Org.BouncyCastle.Math.EC.Custom.Sec { internal class SecP521R1Point : AbstractFpPoint { internal SecP521R1Point(ECCurve curve, ECFieldElement x, ECFieldElement y) : base(curve, x, y) { } internal SecP521R1Point(ECCurve curve, ECFieldElement x, ECFieldElement y, ECFieldElement[] zs) : base(curve, x, y, zs) { } protected override ECPoint Detach() { return new SecP521R1Point(null, AffineXCoord, AffineYCoord); } public override ECPoint Add(ECPoint b) { if (base.IsInfinity) return b; if (b.IsInfinity) return this; if (this == b) return Twice(); ECCurve curve = Curve; SecP521R1FieldElement secP521R1FieldElement = (SecP521R1FieldElement)base.RawXCoord; SecP521R1FieldElement secP521R1FieldElement2 = (SecP521R1FieldElement)base.RawYCoord; SecP521R1FieldElement secP521R1FieldElement3 = (SecP521R1FieldElement)b.RawXCoord; SecP521R1FieldElement secP521R1FieldElement4 = (SecP521R1FieldElement)b.RawYCoord; SecP521R1FieldElement secP521R1FieldElement5 = (SecP521R1FieldElement)base.RawZCoords[0]; SecP521R1FieldElement secP521R1FieldElement6 = (SecP521R1FieldElement)b.RawZCoords[0]; uint[] tt = Nat.Create(33); uint[] array = Nat.Create(17); uint[] array2 = Nat.Create(17); uint[] array3 = Nat.Create(17); uint[] array4 = Nat.Create(17); bool isOne = secP521R1FieldElement5.IsOne; uint[] array5; uint[] array6; if (isOne) { array5 = secP521R1FieldElement3.x; array6 = secP521R1FieldElement4.x; } else { array6 = array3; SecP521R1Field.Square(secP521R1FieldElement5.x, array6, tt); array5 = array2; SecP521R1Field.Multiply(array6, secP521R1FieldElement3.x, array5, tt); SecP521R1Field.Multiply(array6, secP521R1FieldElement5.x, array6, tt); SecP521R1Field.Multiply(array6, secP521R1FieldElement4.x, array6, tt); } bool isOne2 = secP521R1FieldElement6.IsOne; uint[] array7; uint[] array8; if (isOne2) { array7 = secP521R1FieldElement.x; array8 = secP521R1FieldElement2.x; } else { array8 = array4; SecP521R1Field.Square(secP521R1FieldElement6.x, array8, tt); array7 = array; SecP521R1Field.Multiply(array8, secP521R1FieldElement.x, array7, tt); SecP521R1Field.Multiply(array8, secP521R1FieldElement6.x, array8, tt); SecP521R1Field.Multiply(array8, secP521R1FieldElement2.x, array8, tt); } uint[] array9 = Nat.Create(17); SecP521R1Field.Subtract(array7, array5, array9); uint[] array10 = array2; SecP521R1Field.Subtract(array8, array6, array10); if (Nat.IsZero(17, array9)) { if (Nat.IsZero(17, array10)) return Twice(); return curve.Infinity; } uint[] array11 = array3; SecP521R1Field.Square(array9, array11, tt); uint[] array12 = Nat.Create(17); SecP521R1Field.Multiply(array11, array9, array12, tt); uint[] array13 = array3; SecP521R1Field.Multiply(array11, array7, array13, tt); SecP521R1Field.Multiply(array8, array12, array, tt); SecP521R1FieldElement secP521R1FieldElement7 = new SecP521R1FieldElement(array4); SecP521R1Field.Square(array10, secP521R1FieldElement7.x, tt); SecP521R1Field.Add(secP521R1FieldElement7.x, array12, secP521R1FieldElement7.x); SecP521R1Field.Subtract(secP521R1FieldElement7.x, array13, secP521R1FieldElement7.x); SecP521R1Field.Subtract(secP521R1FieldElement7.x, array13, secP521R1FieldElement7.x); SecP521R1FieldElement secP521R1FieldElement8 = new SecP521R1FieldElement(array12); SecP521R1Field.Subtract(array13, secP521R1FieldElement7.x, secP521R1FieldElement8.x); SecP521R1Field.Multiply(secP521R1FieldElement8.x, array10, array2, tt); SecP521R1Field.Subtract(array2, array, secP521R1FieldElement8.x); SecP521R1FieldElement secP521R1FieldElement9 = new SecP521R1FieldElement(array9); if (!isOne) SecP521R1Field.Multiply(secP521R1FieldElement9.x, secP521R1FieldElement5.x, secP521R1FieldElement9.x, tt); if (!isOne2) SecP521R1Field.Multiply(secP521R1FieldElement9.x, secP521R1FieldElement6.x, secP521R1FieldElement9.x, tt); ECFieldElement[] zs = new ECFieldElement[1] { secP521R1FieldElement9 }; return new SecP521R1Point(curve, secP521R1FieldElement7, secP521R1FieldElement8, zs); } public override ECPoint Twice() { if (base.IsInfinity) return this; ECCurve curve = Curve; SecP521R1FieldElement secP521R1FieldElement = (SecP521R1FieldElement)base.RawYCoord; if (secP521R1FieldElement.IsZero) return curve.Infinity; SecP521R1FieldElement secP521R1FieldElement2 = (SecP521R1FieldElement)base.RawXCoord; SecP521R1FieldElement secP521R1FieldElement3 = (SecP521R1FieldElement)base.RawZCoords[0]; uint[] tt = Nat.Create(33); uint[] array = Nat.Create(17); uint[] array2 = Nat.Create(17); uint[] array3 = Nat.Create(17); SecP521R1Field.Square(secP521R1FieldElement.x, array3, tt); uint[] array4 = Nat.Create(17); SecP521R1Field.Square(array3, array4, tt); bool isOne = secP521R1FieldElement3.IsOne; uint[] array5 = secP521R1FieldElement3.x; if (!isOne) { array5 = array2; SecP521R1Field.Square(secP521R1FieldElement3.x, array5, tt); } SecP521R1Field.Subtract(secP521R1FieldElement2.x, array5, array); uint[] array6 = array2; SecP521R1Field.Add(secP521R1FieldElement2.x, array5, array6); SecP521R1Field.Multiply(array6, array, array6, tt); Nat.AddBothTo(17, array6, array6, array6); SecP521R1Field.Reduce23(array6); uint[] array7 = array3; SecP521R1Field.Multiply(array3, secP521R1FieldElement2.x, array7, tt); Nat.ShiftUpBits(17, array7, 2, 0); SecP521R1Field.Reduce23(array7); Nat.ShiftUpBits(17, array4, 3, 0, array); SecP521R1Field.Reduce23(array); SecP521R1FieldElement secP521R1FieldElement4 = new SecP521R1FieldElement(array4); SecP521R1Field.Square(array6, secP521R1FieldElement4.x, tt); SecP521R1Field.Subtract(secP521R1FieldElement4.x, array7, secP521R1FieldElement4.x); SecP521R1Field.Subtract(secP521R1FieldElement4.x, array7, secP521R1FieldElement4.x); SecP521R1FieldElement secP521R1FieldElement5 = new SecP521R1FieldElement(array7); SecP521R1Field.Subtract(array7, secP521R1FieldElement4.x, secP521R1FieldElement5.x); SecP521R1Field.Multiply(secP521R1FieldElement5.x, array6, secP521R1FieldElement5.x, tt); SecP521R1Field.Subtract(secP521R1FieldElement5.x, array, secP521R1FieldElement5.x); SecP521R1FieldElement secP521R1FieldElement6 = new SecP521R1FieldElement(array6); SecP521R1Field.Twice(secP521R1FieldElement.x, secP521R1FieldElement6.x); if (!isOne) SecP521R1Field.Multiply(secP521R1FieldElement6.x, secP521R1FieldElement3.x, secP521R1FieldElement6.x, tt); return new SecP521R1Point(curve, secP521R1FieldElement4, secP521R1FieldElement5, new ECFieldElement[1] { secP521R1FieldElement6 }); } public override ECPoint TwicePlus(ECPoint b) { if (this == b) return ThreeTimes(); if (base.IsInfinity) return b; if (b.IsInfinity) return Twice(); if (base.RawYCoord.IsZero) return b; return Twice().Add(b); } public override ECPoint ThreeTimes() { if (base.IsInfinity || base.RawYCoord.IsZero) return this; return Twice().Add(this); } public override ECPoint Negate() { if (base.IsInfinity) return this; return new SecP521R1Point(Curve, base.RawXCoord, base.RawYCoord.Negate(), base.RawZCoords); } } }