Skip to content

Commit

Permalink
Merge pull request #308 from cuviper/less-wide
Browse files Browse the repository at this point in the history
Avoid some widening operations
  • Loading branch information
cuviper committed May 11, 2024
2 parents f511841 + 98bea13 commit 71b78bc
Show file tree
Hide file tree
Showing 4 changed files with 17 additions and 24 deletions.
12 changes: 6 additions & 6 deletions src/biguint/addition.rs
Expand Up @@ -25,14 +25,14 @@ fn adc(carry: u8, a: u32, b: u32, out: &mut u32) -> u8 {
}

// fallback for environments where we don't have an addcarry intrinsic
// (copied from the standard library's `carrying_add`)
#[cfg(not(any(target_arch = "x86", target_arch = "x86_64")))]
#[inline]
fn adc(carry: u8, a: BigDigit, b: BigDigit, out: &mut BigDigit) -> u8 {
use crate::big_digit::DoubleBigDigit;

let sum = DoubleBigDigit::from(a) + DoubleBigDigit::from(b) + DoubleBigDigit::from(carry);
*out = sum as BigDigit;
(sum >> big_digit::BITS) as u8
fn adc(carry: u8, lhs: BigDigit, rhs: BigDigit, out: &mut BigDigit) -> u8 {
let (a, b) = lhs.overflowing_add(rhs);
let (c, d) = a.overflowing_add(carry as BigDigit);
*out = c;
u8::from(b || d)
}

/// Two argument addition of raw slices, `a += b`, returning the carry.
Expand Down
9 changes: 5 additions & 4 deletions src/biguint/monty.rs
Expand Up @@ -3,7 +3,7 @@ use core::mem;
use core::ops::Shl;
use num_traits::One;

use crate::big_digit::{self, BigDigit, DoubleBigDigit, SignedDoubleBigDigit};
use crate::big_digit::{self, BigDigit, DoubleBigDigit};
use crate::biguint::BigUint;

struct MontyReducer {
Expand All @@ -15,16 +15,17 @@ struct MontyReducer {
fn inv_mod_alt(b: BigDigit) -> BigDigit {
assert_ne!(b & 1, 0);

let mut k0 = 2 - b as SignedDoubleBigDigit;
let mut t = (b - 1) as SignedDoubleBigDigit;
let mut k0 = BigDigit::wrapping_sub(2, b);
let mut t = b - 1;
let mut i = 1;
while i < big_digit::BITS {
t = t.wrapping_mul(t);
k0 = k0.wrapping_mul(t + 1);

i <<= 1;
}
-k0 as BigDigit
debug_assert_eq!(k0.wrapping_mul(b), 1);
k0.wrapping_neg()
}

impl MontyReducer {
Expand Down
14 changes: 6 additions & 8 deletions src/biguint/subtraction.rs
Expand Up @@ -25,16 +25,14 @@ fn sbb(borrow: u8, a: u32, b: u32, out: &mut u32) -> u8 {
}

// fallback for environments where we don't have a subborrow intrinsic
// (copied from the standard library's `borrowing_sub`)
#[cfg(not(any(target_arch = "x86", target_arch = "x86_64")))]
#[inline]
fn sbb(borrow: u8, a: BigDigit, b: BigDigit, out: &mut BigDigit) -> u8 {
use crate::big_digit::SignedDoubleBigDigit;

let difference = SignedDoubleBigDigit::from(a)
- SignedDoubleBigDigit::from(b)
- SignedDoubleBigDigit::from(borrow);
*out = difference as BigDigit;
u8::from(difference < 0)
fn sbb(borrow: u8, lhs: BigDigit, rhs: BigDigit, out: &mut BigDigit) -> u8 {
let (a, b) = lhs.overflowing_sub(rhs);
let (c, d) = a.overflowing_sub(borrow as BigDigit);
*out = c;
u8::from(b || d)
}

pub(super) fn sub2(a: &mut [BigDigit], b: &[BigDigit]) {
Expand Down
6 changes: 0 additions & 6 deletions src/lib.rs
Expand Up @@ -238,12 +238,6 @@ mod big_digit {
pub(crate) type DoubleBigDigit = u128;
);

// A [`SignedDoubleBigDigit`] is the signed version of [`DoubleBigDigit`].
cfg_digit!(
pub(crate) type SignedDoubleBigDigit = i64;
pub(crate) type SignedDoubleBigDigit = i128;
);

pub(crate) const BITS: u8 = BigDigit::BITS as u8;
pub(crate) const HALF_BITS: u8 = BITS / 2;
pub(crate) const HALF: BigDigit = (1 << HALF_BITS) - 1;
Expand Down

0 comments on commit 71b78bc

Please sign in to comment.