saymrwulf/curve25519-dalek-source
1
0
Fork 0
mirror of https://github.com/saymrwulf/curve25519-dalek-source.git synced 2026-09-06 20:41:14 +00:00
Code Issues Projects Releases Packages Wiki Activity Actions

Merge pull request #245 from dalek-cryptography/optimize-variable-base

Browse source 
Save 2.5% on variable-base scmul by squeezing some multiplications.
This commit is contained in:
Henry de Valence 2019-05-21 08:52:53 -07:00 committed by GitHub
commit e726147af8
No known key found for this signature in database
GPG key ID: 4AEE18F83AFDEB23
1 changed files with 20 additions and 6 deletions
Show stats Download patch file Download diff file Expand all files Collapse all files

26
src/backend/serial/scalar_mul/variable_base.rs
View file

@ -3,7 +3,7 @@
use traits::Identity; use traits::Identity;
use scalar::Scalar; use scalar::Scalar;
use edwards::EdwardsPoint; use edwards::EdwardsPoint;
use backend::serial::curve_models::ProjectiveNielsPoint; use backend::serial::curve_models::{ProjectiveNielsPoint, ProjectivePoint};
use window::LookupTable; use window::LookupTable;
/// Perform constant-time, variable-base scalar multiplication. /// Perform constant-time, variable-base scalar multiplication.
@ -23,10 +23,24 @@ pub(crate) fn mul(point: &EdwardsPoint, scalar: &Scalar) -> EdwardsPoint {
// s*P = P*s_0 + 16*(P*s_1 + 16*(P*s_2 + 16*( ... + P*s_63)...)) // s*P = P*s_0 + 16*(P*s_1 + 16*(P*s_2 + 16*( ... + P*s_63)...))
// //
// We sum right-to-left. // We sum right-to-left.
let mut Q = EdwardsPoint::identity();
for i in (0..64).rev() { // Unwrap first loop iteration to save computing 16*identity
Q = Q.mul_by_pow_2(4); let mut tmp2 = ProjectivePoint::identity();
Q = (&Q + &lookup_table.select(scalar_digits[i])).to_extended(); let mut tmp3 = EdwardsPoint::identity();
let mut tmp1 = &tmp3 + &lookup_table.select(scalar_digits[63]);
// Now tmp1 = s_63*P in P1xP1 coords
for i in (0..63).rev() {
tmp2 = tmp1.to_projective(); // tmp2 = (prev) in P2 coords
tmp1 = tmp2.double(); // tmp1 = 2*(prev) in P1xP1 coords
tmp2 = tmp1.to_projective(); // tmp2 = 2*(prev) in P2 coords
tmp1 = tmp2.double(); // tmp1 = 4*(prev) in P1xP1 coords
tmp2 = tmp1.to_projective(); // tmp2 = 4*(prev) in P2 coords
tmp1 = tmp2.double(); // tmp1 = 8*(prev) in P1xP1 coords
tmp2 = tmp1.to_projective(); // tmp2 = 8*(prev) in P2 coords
tmp1 = tmp2.double(); // tmp1 = 16*(prev) in P1xP1 coords
tmp3 = tmp1.to_extended(); // tmp3 = 16*(prev) in P3 coords
tmp1 = &tmp3 + &lookup_table.select(scalar_digits[i]);
// Now tmp1 = s_i*P + 16*(prev) in P1xP1 coords
} }
Q tmp1.to_extended()
} }