mirror of
https://github.com/saymrwulf/risc0-curve25519-dalek-source.git
synced 2026-09-04 20:03:40 +00:00
This doesn't (yet) give any speedup over the non-precomputed multiscalar multiplication, and it's not clear that it's a good idea to commit to supporting it in the future. Removing it means that it's not committed-to as part of the public API, but the source is still there in the tree if we want to revisit it later.
377 lines
13 KiB
Rust
377 lines
13 KiB
Rust
// -*- mode: rust; -*-
|
|
//
|
|
// This file is part of curve25519-dalek.
|
|
// Copyright (c) 2016-2018 Isis Lovecruft, Henry de Valence
|
|
// See LICENSE for licensing information.
|
|
//
|
|
// Authors:
|
|
// - Isis Agora Lovecruft <isis@patternsinthevoid.net>
|
|
// - Henry de Valence <hdevalence@hdevalence.ca>
|
|
|
|
//! Module for common traits.
|
|
|
|
use core::borrow::Borrow;
|
|
|
|
use subtle;
|
|
|
|
use scalar::Scalar;
|
|
|
|
// ------------------------------------------------------------------------
|
|
// Public Traits
|
|
// ------------------------------------------------------------------------
|
|
|
|
/// Trait for getting the identity element of a point type.
|
|
pub trait Identity {
|
|
/// Returns the identity element of the curve.
|
|
/// Can be used as a constructor.
|
|
fn identity() -> Self;
|
|
}
|
|
|
|
/// Trait for testing if a curve point is equivalent to the identity point.
|
|
pub trait IsIdentity {
|
|
/// Return true if this element is the identity element of the curve.
|
|
fn is_identity(&self) -> bool;
|
|
}
|
|
|
|
/// Implement generic identity equality testing for a point representations
|
|
/// which have constant-time equality testing and a defined identity
|
|
/// constructor.
|
|
impl<T> IsIdentity for T
|
|
where
|
|
T: subtle::ConstantTimeEq + Identity,
|
|
{
|
|
fn is_identity(&self) -> bool {
|
|
self.ct_eq(&T::identity()).unwrap_u8() == 1u8
|
|
}
|
|
}
|
|
|
|
/// A trait for constant-time multiscalar multiplication without precomputation.
|
|
pub trait MultiscalarMul {
|
|
/// The type of point being multiplied, e.g., `RistrettoPoint`.
|
|
type Point;
|
|
|
|
/// Given an iterator of (possibly secret) scalars and an iterator of
|
|
/// public points, compute
|
|
/// $$
|
|
/// Q = c\_1 P\_1 + \cdots + c\_n P\_n.
|
|
/// $$
|
|
///
|
|
/// It is an error to call this function with two iterators of different lengths.
|
|
///
|
|
/// # Examples
|
|
///
|
|
/// The trait bound aims for maximum flexibility: the inputs must be
|
|
/// convertable to iterators (`I: IntoIter`), and the iterator's items
|
|
/// must be `Borrow<Scalar>` (or `Borrow<Point>`), to allow
|
|
/// iterators returning either `Scalar`s or `&Scalar`s.
|
|
///
|
|
/// ```
|
|
/// use curve25519_dalek::constants;
|
|
/// use curve25519_dalek::traits::MultiscalarMul;
|
|
/// use curve25519_dalek::ristretto::RistrettoPoint;
|
|
/// use curve25519_dalek::scalar::Scalar;
|
|
///
|
|
/// // Some scalars
|
|
/// let a = Scalar::from(87329482u64);
|
|
/// let b = Scalar::from(37264829u64);
|
|
/// let c = Scalar::from(98098098u64);
|
|
///
|
|
/// // Some points
|
|
/// let P = constants::RISTRETTO_BASEPOINT_POINT;
|
|
/// let Q = P + P;
|
|
/// let R = P + Q;
|
|
///
|
|
/// // A1 = a*P + b*Q + c*R
|
|
/// let abc = [a,b,c];
|
|
/// let A1 = RistrettoPoint::multiscalar_mul(&abc, &[P,Q,R]);
|
|
/// // Note: (&abc).into_iter(): Iterator<Item=&Scalar>
|
|
///
|
|
/// // A2 = (-a)*P + (-b)*Q + (-c)*R
|
|
/// let minus_abc = abc.iter().map(|x| -x);
|
|
/// let A2 = RistrettoPoint::multiscalar_mul(minus_abc, &[P,Q,R]);
|
|
/// // Note: minus_abc.into_iter(): Iterator<Item=Scalar>
|
|
///
|
|
/// assert_eq!(A1.compress(), (-A2).compress());
|
|
/// ```
|
|
fn multiscalar_mul<I, J>(scalars: I, points: J) -> Self::Point
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Scalar>,
|
|
J: IntoIterator,
|
|
J::Item: Borrow<Self::Point>;
|
|
}
|
|
|
|
/// A trait for variable-time multiscalar multiplication without precomputation.
|
|
pub trait VartimeMultiscalarMul {
|
|
/// The type of point being multiplied, e.g., `RistrettoPoint`.
|
|
type Point;
|
|
|
|
/// Given an iterator of public scalars and an iterator of
|
|
/// `Option`s of points, compute either `Some(Q)`, where
|
|
/// $$
|
|
/// Q = c\_1 P\_1 + \cdots + c\_n P\_n,
|
|
/// $$
|
|
/// if all points were `Some(P_i)`, or else return `None`.
|
|
///
|
|
/// This function is particularly useful when verifying statements
|
|
/// involving compressed points. Accepting `Option<Point>` allows
|
|
/// inlining point decompression into the multiscalar call,
|
|
/// avoiding the need for temporary buffers.
|
|
/// ```
|
|
/// use curve25519_dalek::constants;
|
|
/// use curve25519_dalek::traits::VartimeMultiscalarMul;
|
|
/// use curve25519_dalek::ristretto::RistrettoPoint;
|
|
/// use curve25519_dalek::scalar::Scalar;
|
|
///
|
|
/// // Some scalars
|
|
/// let a = Scalar::from(87329482u64);
|
|
/// let b = Scalar::from(37264829u64);
|
|
/// let c = Scalar::from(98098098u64);
|
|
/// let abc = [a,b,c];
|
|
///
|
|
/// // Some points
|
|
/// let P = constants::RISTRETTO_BASEPOINT_POINT;
|
|
/// let Q = P + P;
|
|
/// let R = P + Q;
|
|
/// let PQR = [P, Q, R];
|
|
///
|
|
/// let compressed = [P.compress(), Q.compress(), R.compress()];
|
|
///
|
|
/// // Now we can compute A1 = a*P + b*Q + c*R using P, Q, R:
|
|
/// let A1 = RistrettoPoint::vartime_multiscalar_mul(&abc, &PQR);
|
|
///
|
|
/// // Or using the compressed points:
|
|
/// let A2 = RistrettoPoint::optional_multiscalar_mul(
|
|
/// &abc,
|
|
/// compressed.iter().map(|pt| pt.decompress()),
|
|
/// );
|
|
///
|
|
/// assert_eq!(A2, Some(A1));
|
|
///
|
|
/// // It's also possible to mix compressed and uncompressed points:
|
|
/// let A3 = RistrettoPoint::optional_multiscalar_mul(
|
|
/// abc.iter()
|
|
/// .chain(abc.iter()),
|
|
/// compressed.iter().map(|pt| pt.decompress())
|
|
/// .chain(PQR.iter().map(|&pt| Some(pt))),
|
|
/// );
|
|
///
|
|
/// assert_eq!(A3, Some(A1+A1));
|
|
/// ```
|
|
fn optional_multiscalar_mul<I, J>(scalars: I, points: J) -> Option<Self::Point>
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Scalar>,
|
|
J: IntoIterator<Item = Option<Self::Point>>;
|
|
|
|
/// Given an iterator of public scalars and an iterator of
|
|
/// public points, compute
|
|
/// $$
|
|
/// Q = c\_1 P\_1 + \cdots + c\_n P\_n,
|
|
/// $$
|
|
/// using variable-time operations.
|
|
///
|
|
/// It is an error to call this function with two iterators of different lengths.
|
|
///
|
|
/// # Examples
|
|
///
|
|
/// The trait bound aims for maximum flexibility: the inputs must be
|
|
/// convertable to iterators (`I: IntoIter`), and the iterator's items
|
|
/// must be `Borrow<Scalar>` (or `Borrow<Point>`), to allow
|
|
/// iterators returning either `Scalar`s or `&Scalar`s.
|
|
///
|
|
/// ```
|
|
/// use curve25519_dalek::constants;
|
|
/// use curve25519_dalek::traits::VartimeMultiscalarMul;
|
|
/// use curve25519_dalek::ristretto::RistrettoPoint;
|
|
/// use curve25519_dalek::scalar::Scalar;
|
|
///
|
|
/// // Some scalars
|
|
/// let a = Scalar::from(87329482u64);
|
|
/// let b = Scalar::from(37264829u64);
|
|
/// let c = Scalar::from(98098098u64);
|
|
///
|
|
/// // Some points
|
|
/// let P = constants::RISTRETTO_BASEPOINT_POINT;
|
|
/// let Q = P + P;
|
|
/// let R = P + Q;
|
|
///
|
|
/// // A1 = a*P + b*Q + c*R
|
|
/// let abc = [a,b,c];
|
|
/// let A1 = RistrettoPoint::vartime_multiscalar_mul(&abc, &[P,Q,R]);
|
|
/// // Note: (&abc).into_iter(): Iterator<Item=&Scalar>
|
|
///
|
|
/// // A2 = (-a)*P + (-b)*Q + (-c)*R
|
|
/// let minus_abc = abc.iter().map(|x| -x);
|
|
/// let A2 = RistrettoPoint::vartime_multiscalar_mul(minus_abc, &[P,Q,R]);
|
|
/// // Note: minus_abc.into_iter(): Iterator<Item=Scalar>
|
|
///
|
|
/// assert_eq!(A1.compress(), (-A2).compress());
|
|
/// ```
|
|
#[allow(non_snake_case)]
|
|
fn vartime_multiscalar_mul<I, J>(scalars: I, points: J) -> Self::Point
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Scalar>,
|
|
J: IntoIterator,
|
|
J::Item: Borrow<Self::Point>,
|
|
Self::Point: Clone,
|
|
{
|
|
Self::optional_multiscalar_mul(
|
|
scalars,
|
|
points.into_iter().map(|P| Some(P.borrow().clone())),
|
|
)
|
|
.unwrap()
|
|
}
|
|
}
|
|
|
|
/// A trait for variable-time multiscalar multiplication with precomputation.
|
|
///
|
|
/// A general multiscalar multiplication with precomputation can be written as
|
|
/// $$
|
|
/// Q = a_1 A_1 + \cdots + a_n A_n + b_1 B_1 + \cdots + b_m B_m,
|
|
/// $$
|
|
/// where the \\(B_i\\) are *static* points, for which precomputation
|
|
/// is possible, and the \\(A_j\\) are *dynamic* points, for which
|
|
/// precomputation is not possible.
|
|
///
|
|
/// This trait has three methods for performing this computation:
|
|
///
|
|
/// * [`vartime_multiscalar_mul`], which handles the special case
|
|
/// where \\(n = 0\\) and there are no dynamic points;
|
|
///
|
|
/// * [`vartime_mixed_multiscalar_mul`], which takes the dynamic
|
|
/// points as already-validated `Point`s and is infallible;
|
|
///
|
|
/// * [`optional_mixed_multiscalar_mul`], which takes the dynamic
|
|
/// points as `Option<Point>`s and returns an `Option<Point>`,
|
|
/// allowing decompression to be composed into the input iterators.
|
|
///
|
|
/// All methods require that the lengths of the input iterators be
|
|
/// known and matching, as if they were `ExactSizeIterator`s. (It
|
|
/// does not require `ExactSizeIterator` only because that trait is
|
|
/// broken).
|
|
pub trait VartimePrecomputedMultiscalarMul: Sized {
|
|
/// The type of point to be multiplied, e.g., `RistrettoPoint`.
|
|
type Point: Clone;
|
|
|
|
/// Given the static points \\( B_i \\), perform precomputation
|
|
/// and return the precomputation data.
|
|
fn new<I>(static_points: I) -> Self
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Self::Point>;
|
|
|
|
/// Given `static_scalars`, an iterator of public scalars
|
|
/// \\(b_i\\), compute
|
|
/// $$
|
|
/// Q = b_1 B_1 + \cdots + b_m B_m,
|
|
/// $$
|
|
/// where the \\(B_j\\) are the points that were supplied to `new`.
|
|
///
|
|
/// It is an error to call this function with iterators of
|
|
/// inconsistent lengths.
|
|
///
|
|
/// The trait bound aims for maximum flexibility: the input must
|
|
/// be convertable to iterators (`I: IntoIter`), and the
|
|
/// iterator's items must be `Borrow<Scalar>`, to allow iterators
|
|
/// returning either `Scalar`s or `&Scalar`s.
|
|
fn vartime_multiscalar_mul<I>(&self, static_scalars: I) -> Self::Point
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Scalar>,
|
|
{
|
|
use core::iter;
|
|
|
|
Self::vartime_mixed_multiscalar_mul(
|
|
self,
|
|
static_scalars,
|
|
iter::empty::<Scalar>(),
|
|
iter::empty::<Self::Point>(),
|
|
)
|
|
}
|
|
|
|
/// Given `static_scalars`, an iterator of public scalars
|
|
/// \\(b_i\\), `dynamic_scalars`, an iterator of public scalars
|
|
/// \\(a_i\\), and `dynamic_points`, an iterator of points
|
|
/// \\(A_i\\), compute
|
|
/// $$
|
|
/// Q = a_1 A_1 + \cdots + a_n A_n + b_1 B_1 + \cdots + b_m B_m,
|
|
/// $$
|
|
/// where the \\(B_j\\) are the points that were supplied to `new`.
|
|
///
|
|
/// It is an error to call this function with iterators of
|
|
/// inconsistent lengths.
|
|
///
|
|
/// The trait bound aims for maximum flexibility: the inputs must be
|
|
/// convertable to iterators (`I: IntoIter`), and the iterator's items
|
|
/// must be `Borrow<Scalar>` (or `Borrow<Point>`), to allow
|
|
/// iterators returning either `Scalar`s or `&Scalar`s.
|
|
fn vartime_mixed_multiscalar_mul<I, J, K>(
|
|
&self,
|
|
static_scalars: I,
|
|
dynamic_scalars: J,
|
|
dynamic_points: K,
|
|
) -> Self::Point
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Scalar>,
|
|
J: IntoIterator,
|
|
J::Item: Borrow<Scalar>,
|
|
K: IntoIterator,
|
|
K::Item: Borrow<Self::Point>,
|
|
{
|
|
Self::optional_mixed_multiscalar_mul(
|
|
self,
|
|
static_scalars,
|
|
dynamic_scalars,
|
|
dynamic_points.into_iter().map(|P| Some(P.borrow().clone())),
|
|
)
|
|
.unwrap()
|
|
}
|
|
|
|
/// Given `static_scalars`, an iterator of public scalars
|
|
/// \\(b_i\\), `dynamic_scalars`, an iterator of public scalars
|
|
/// \\(a_i\\), and `dynamic_points`, an iterator of points
|
|
/// \\(A_i\\), compute
|
|
/// $$
|
|
/// Q = a_1 A_1 + \cdots + a_n A_n + b_1 B_1 + \cdots + b_m B_m,
|
|
/// $$
|
|
/// where the \\(B_j\\) are the points that were supplied to `new`.
|
|
///
|
|
/// If any of the dynamic points were `None`, return `None`.
|
|
///
|
|
/// It is an error to call this function with iterators of
|
|
/// inconsistent lengths.
|
|
///
|
|
/// This function is particularly useful when verifying statements
|
|
/// involving compressed points. Accepting `Option<Point>` allows
|
|
/// inlining point decompression into the multiscalar call,
|
|
/// avoiding the need for temporary buffers.
|
|
fn optional_mixed_multiscalar_mul<I, J, K>(
|
|
&self,
|
|
static_scalars: I,
|
|
dynamic_scalars: J,
|
|
dynamic_points: K,
|
|
) -> Option<Self::Point>
|
|
where
|
|
I: IntoIterator,
|
|
I::Item: Borrow<Scalar>,
|
|
J: IntoIterator,
|
|
J::Item: Borrow<Scalar>,
|
|
K: IntoIterator<Item = Option<Self::Point>>;
|
|
}
|
|
|
|
// ------------------------------------------------------------------------
|
|
// Private Traits
|
|
// ------------------------------------------------------------------------
|
|
|
|
/// Trait for checking whether a point is on the curve.
|
|
///
|
|
/// This trait is only for debugging/testing, since it should be
|
|
/// impossible for a `curve25519-dalek` user to construct an invalid
|
|
/// point.
|
|
pub(crate) trait ValidityCheck {
|
|
/// Checks whether the point is on the curve. Not CT.
|
|
fn is_valid(&self) -> bool;
|
|
}
|