pub(super) trait Curve : Copy {
type Int: Int;
- // With const generics, both IntModP and IntModN can be replaced with a single IntMod.
- type IntModP: IntMod<I = Self::Int>;
- type IntModN: IntMod<I = Self::Int>;
+ // With const generics, both CurveField and ScalarField can be replaced with a single IntMod.
+ type CurveField: IntMod<I = Self::Int>;
+ type ScalarField: IntMod<I = Self::Int>;
- type P: PrimeModulus<Self::Int>;
- type N: PrimeModulus<Self::Int>;
+ type CurveModulus: PrimeModulus<Self::Int>;
+ type ScalarModulus: PrimeModulus<Self::Int>;
// Curve parameters y^2 = x^3 + ax + b
- const A: Self::IntModP;
- const B: Self::IntModP;
+ const A: Self::CurveField;
+ const B: Self::CurveField;
const G: Point<Self>;
}
#[derive(Clone, PartialEq, Eq)]
+/// A Point, stored in Jacobian coordinates
pub(super) struct Point<C: Curve + ?Sized> {
- x: C::IntModP,
- y: C::IntModP,
- z: C::IntModP,
+ x: C::CurveField,
+ y: C::CurveField,
+ z: C::CurveField,
}
impl<C: Curve + ?Sized> Point<C> {
- fn on_curve(x: &C::IntModP, y: &C::IntModP) -> Result<(), ()> {
+ fn check_curve_conditions() {
+ debug_assert!(C::ScalarModulus::PRIME < C::CurveModulus::PRIME, "N is < P");
+ }
+
+ fn on_curve(x: &C::CurveField, y: &C::CurveField) -> Result<(), ()> {
let x_2 = x.square();
let x_3 = x_2.mul(&x);
let v = x_3.add(&C::A.mul(&x)).add(&C::B);
}
#[cfg(debug_assertions)]
- fn on_curve_z(x: &C::IntModP, y: &C::IntModP, z: &C::IntModP) -> Result<(), ()> {
- let m = C::IntModP::from_modinv_of(z.clone().into_i())?;
+ fn on_curve_z(x: &C::CurveField, y: &C::CurveField, z: &C::CurveField) -> Result<(), ()> {
+ // m = 1 / z
+ // x_norm = x * m^2
+ // y_norm = y * m^3
+
+ let m = C::CurveField::from_modinv_of(z.clone().into_i())?;
let m_2 = m.square();
let m_3 = m_2.mul(&m);
let x_norm = x.mul(&m_2);
}
#[cfg(test)]
- fn normalize_x(&self) -> Result<C::IntModP, ()> {
- let m = C::IntModP::from_modinv_of(self.z.clone().into_i())?;
+ fn normalize_x(&self) -> Result<C::CurveField, ()> {
+ let m = C::CurveField::from_modinv_of(self.z.clone().into_i())?;
Ok(self.x.mul(&m.square()))
}
fn from_xy(x: C::Int, y: C::Int) -> Result<Self, ()> {
- let x = C::IntModP::from_i(x);
- let y = C::IntModP::from_i(y);
+ Self::check_curve_conditions();
+
+ let x = C::CurveField::from_i(x);
+ let y = C::CurveField::from_i(y);
Self::on_curve(&x, &y)?;
- Ok(Point { x, y, z: C::IntModP::ONE })
+ Ok(Point { x, y, z: C::CurveField::ONE })
}
- pub(super) const fn from_xy_assuming_on_curve(x: C::IntModP, y: C::IntModP) -> Self {
- Point { x, y, z: C::IntModP::ONE }
+ pub(super) const fn from_xy_assuming_on_curve(x: C::CurveField, y: C::CurveField) -> Self {
+ Point { x, y, z: C::CurveField::ONE }
}
/// Checks that `expected_x` is equal to our X affine coordinate (without modular inversion).
- fn eq_x(&self, expected_x: &C::IntModN) -> Result<(), ()> {
- debug_assert!(expected_x.clone().into_i() < C::P::PRIME, "N is < P");
-
+ fn eq_x(&self, expected_x: &C::ScalarField) -> Result<(), ()> {
// If x is between N and P the below calculations will fail and we'll spuriously reject a
// signature and the wycheproof tests will fail. We should in theory accept such
// signatures, but the probability of this happening at random is roughly 1/2^128, i.e. we
// really don't need to handle it in practice. Thus, we only bother to do this in tests.
+ debug_assert!(expected_x.clone().into_i() < C::CurveModulus::PRIME, "N is < P");
+ debug_assert!(C::ScalarModulus::PRIME < C::CurveModulus::PRIME, "N is < P");
+ #[cfg(debug_assertions)] {
+ // Check the above assertion - ensure the difference between the modulus of the scalar
+ // and curve fields is less than half the bit length of our integers, which are at
+ // least 256 bit long.
+ let scalar_mod_on_curve = C::CurveField::from_i(C::ScalarModulus::PRIME);
+ let diff = C::CurveField::ZERO.sub(&scalar_mod_on_curve);
+ assert!(C::Int::BYTES * 8 / 2 >= 128, "We assume 256-bit ints and longer");
+ assert!(C::CurveModulus::PRIME.limbs()[0] > (1 << 63), "PRIME should have the top bit set");
+ assert!(C::ScalarModulus::PRIME.limbs()[0] > (1 << 63), "PRIME should have the top bit set");
+ let mut half_bitlen = C::CurveField::ONE;
+ for _ in 0..C::Int::BYTES * 8 / 2 {
+ half_bitlen = half_bitlen.double();
+ }
+ assert!(diff.into_i() < half_bitlen.into_i());
+ }
+
#[allow(unused_mut, unused_assignments)]
let mut slow_check = None;
#[cfg(test)] {
- slow_check = Some(C::IntModN::from_i(self.normalize_x()?.into_i()) == *expected_x);
+ slow_check = Some(C::ScalarField::from_i(self.normalize_x()?.into_i()) == *expected_x);
}
- let e: C::IntModP = C::IntModP::from_i(expected_x.clone().into_i());
- if self.z == C::IntModP::ZERO { return Err(()); }
+ let e: C::CurveField = C::CurveField::from_i(expected_x.clone().into_i());
+ if self.z == C::CurveField::ZERO { return Err(()); }
let ezz = e.mul(&self.z).mul(&self.z);
- if self.x == ezz { Ok(()) } else {
- if slow_check == Some(true) { Ok(()) } else { Err(()) }
- }
+ if self.x == ezz || slow_check == Some(true) { Ok(()) } else { Err(()) }
}
fn double(&self) -> Result<Self, ()> {
- if self.y == C::IntModP::ZERO { return Err(()); }
- if self.z == C::IntModP::ZERO { return Err(()); }
-
- let s = self.x.times_four().mul(&self.y.square());
- let z_2 = self.z.square();
- let z_4 = z_2.square();
- let y_2 = self.y.square();
- let y_4 = y_2.square();
- let x_2 = self.x.square();
- let m = x_2.times_three().add(&C::A.mul(&z_4));
- let x = m.square().sub(&s.double());
- let y = m.mul(&s.sub(&x)).sub(&y_4.times_eight());
- let z = self.y.double().mul(&self.z);
+ if self.y == C::CurveField::ZERO { return Err(()); }
+ if self.z == C::CurveField::ZERO { return Err(()); }
+
+ // https://hyperelliptic.org/EFD/g1p/auto-shortw-jacobian-3.html#doubling-dbl-2001-b
+ // delta = Z1^2
+ // gamma = Y1^2
+ // beta = X1*gamma
+ // alpha = 3*(X1-delta)*(X1+delta)
+ // X3 = alpha^2-8*beta
+ // Z3 = (Y1+Z1)^2-gamma-delta
+ // Y3 = alpha*(4*beta-X3)-8*gamma^2
+
+ let delta = self.z.square();
+ let gamma = self.y.square();
+ let beta = self.x.mul(&gamma);
+ let alpha = self.x.sub(&delta).times_three().mul(&self.x.add(&delta));
+ let x = alpha.square().sub(&beta.times_eight());
+ let y = alpha.mul(&beta.times_four().sub(&x)).sub(&gamma.square().times_eight());
+ let z = self.y.add(&self.z).square().sub(&gamma).sub(&delta);
#[cfg(debug_assertions)] { assert!(Self::on_curve_z(&x, &y, &z).is_ok()); }
Ok(Point { x, y, z })
}
fn add(&self, o: &Self) -> Result<Self, ()> {
+ // https://hyperelliptic.org/EFD/g1p/auto-shortw-jacobian-3.html#addition-add-2007-bl
+ // Z1Z1 = Z1^2
+ // Z2Z2 = Z2^2
+ // U1 = X1*Z2Z2
+ // U2 = X2*Z1Z1
+ // S1 = Y1*Z2*Z2Z2
+ // S2 = Y2*Z1*Z1Z1
+ // H = U2-U1
+ // I = (2*H)^2
+ // J = H*I
+ // r = 2*(S2-S1)
+ // V = U1*I
+ // X3 = r^2-J-2*V
+ // Y3 = r*(V-X3)-2*S1*J
+ // Z3 = ((Z1+Z2)^2-Z1Z1-Z2Z2)*H
+
let o_z_2 = o.z.square();
let self_z_2 = self.z.square();
let s1 = self.y.mul(&o.z.mul(&o_z_2));
let s2 = o.y.mul(&self.z.mul(&self_z_2));
if u1 == u2 {
- if s1 != s2 { /* PAI */ return Err(()); }
+ if s1 != s2 { /* Point at Infinity */ return Err(()); }
return self.double();
}
let h = u2.sub(&u1);
- let h_2 = h.square();
- let h_3 = h.mul(&h_2);
- let r = s2.sub(&s1);
- let x = r.square().sub(&h_3).sub(&u1.double().mul(&h_2));
- let y = r.mul(&u1.mul(&h_2).sub(&x)).sub(&s1.mul(&h_3));
- let z = h.mul(&self.z).mul(&o.z);
+ let i = h.double().square();
+ let j = h.mul(&i);
+ let r = s2.sub(&s1).double();
+ let v = u1.mul(&i);
+ let x = r.square().sub(&j).sub(&v.double());
+ let y = r.mul(&v.sub(&x)).sub(&s1.double().mul(&j));
+ let z = self.z.add(&o.z).square().sub(&self_z_2).sub(&o_z_2).mul(&h);
#[cfg(debug_assertions)] { assert!(Self::on_curve_z(&x, &y, &z).is_ok()); }
Ok(Point { x, y, z})
/// Calculates i * I + j * J
#[allow(non_snake_case)]
-fn add_two_mul<C: Curve>(i: C::IntModN, I: &Point<C>, j: C::IntModN, J: &Point<C>) -> Result<Point<C>, ()> {
+fn add_two_mul<C: Curve>(i: C::ScalarField, I: &Point<C>, j: C::ScalarField, J: &Point<C>) -> Result<Point<C>, ()> {
let i = i.into_i();
let j = j.into_i();
}
if i_bit {
if let Ok(res) = res_opt.as_mut() {
- // The wycheproof tests expect to see signatures pass even if we hit PAI on an
- // intermediate result. While that's fine, I'm too lazy to go figure out if all
- // our PAI definitions are right and the probability of this happening at
- // random is, basically, the probability of guessing a private key anyway, so
- // its not really worth actually handling outside of tests.
+ // The wycheproof tests expect to see signatures pass even if we hit Point at
+ // Infinity (PAI) on an intermediate result. While that's fine, I'm too lazy to
+ // go figure out if all our PAI definitions are right and the probability of
+ // this happening at random is, basically, the probability of guessing a private
+ // key anyway, so its not really worth actually handling outside of tests.
#[cfg(test)] {
res_opt = res.add(I);
}
}
if j_bit {
if let Ok(res) = res_opt.as_mut() {
- // The wycheproof tests expect to see signatures pass even if we hit PAI on an
- // intermediate result. While that's fine, I'm too lazy to go figure out if all
- // our PAI definitions are right and the probability of this happening at
- // random is, basically, the probability of guessing a private key anyway, so
- // its not really worth actually handling outside of tests.
+ // The wycheproof tests expect to see signatures pass even if we hit Point at
+ // Infinity (PAI) on an intermediate result. While that's fine, I'm too lazy to
+ // go figure out if all our PAI definitions are right and the probability of
+ // this happening at random is, basically, the probability of guessing a private
+ // key anyway, so its not really worth actually handling outside of tests.
#[cfg(test)] {
res_opt = res.add(J);
}
// perfectly safe to do so, the wycheproof tests expect such signatures to be rejected, so we
// do so here.
let r_u256 = C::Int::from_be_bytes(r_bytes)?;
- if r_u256 > C::N::PRIME { return Err(()); }
+ if r_u256 > C::ScalarModulus::PRIME { return Err(()); }
let s_u256 = C::Int::from_be_bytes(s_bytes)?;
- if s_u256 > C::N::PRIME { return Err(()); }
+ if s_u256 > C::ScalarModulus::PRIME { return Err(()); }
- let r = C::IntModN::from_i(r_u256);
- let s_inv = C::IntModN::from_modinv_of(s_u256)?;
+ let r = C::ScalarField::from_i(r_u256);
+ let s_inv = C::ScalarField::from_modinv_of(s_u256)?;
- let z = C::IntModN::from_i(C::Int::from_be_bytes(hash_input)?);
+ let z = C::ScalarField::from_i(C::Int::from_be_bytes(hash_input)?);
let u_a = z.mul(&s_inv);
let u_b = r.mul(&s_inv);