rscrypto 0.1.1

Pure Rust cryptography, hardware-accelerated: BLAKE3, SHA-2/3, AES-GCM, ChaCha20-Poly1305, Ed25519, X25519, HMAC, HKDF, Argon2, CRC. no_std, WASM, ten CPU architectures.
Documentation
use core::arch::wasm32::*;

use super::{BLOCK_SIZE, KEY_SIZE, NONCE_SIZE, load_u32_le, xor_keystream_portable};

const BLOCKS_PER_BATCH: usize = 4;

#[inline]
pub(super) fn xor_keystream(key: &[u8; KEY_SIZE], initial_counter: u32, nonce: &[u8; NONCE_SIZE], buffer: &mut [u8]) {
  // SAFETY: Backend selection guarantees simd128 is available before this wrapper is chosen.
  unsafe { xor_keystream_impl(key, initial_counter, nonce, buffer) }
}

#[target_feature(enable = "simd128")]
unsafe fn xor_keystream_impl(key: &[u8; KEY_SIZE], initial_counter: u32, nonce: &[u8; NONCE_SIZE], buffer: &mut [u8]) {
  let mut counter = initial_counter;
  let mut batches = buffer.chunks_exact_mut(BLOCK_SIZE * BLOCKS_PER_BATCH);
  for chunk in &mut batches {
    debug_assert!(counter.checked_add((BLOCKS_PER_BATCH - 1) as u32).is_some());

    let mut x0 = u32x4_splat(0x6170_7865);
    let mut x1 = u32x4_splat(0x3320_646e);
    let mut x2 = u32x4_splat(0x7962_2d32);
    let mut x3 = u32x4_splat(0x6b20_6574);
    let mut x4 = u32x4_splat(load_u32_le(&key[0..4]));
    let mut x5 = u32x4_splat(load_u32_le(&key[4..8]));
    let mut x6 = u32x4_splat(load_u32_le(&key[8..12]));
    let mut x7 = u32x4_splat(load_u32_le(&key[12..16]));
    let mut x8 = u32x4_splat(load_u32_le(&key[16..20]));
    let mut x9 = u32x4_splat(load_u32_le(&key[20..24]));
    let mut x10 = u32x4_splat(load_u32_le(&key[24..28]));
    let mut x11 = u32x4_splat(load_u32_le(&key[28..32]));
    let mut x12 = u32x4(
      counter,
      counter.wrapping_add(1),
      counter.wrapping_add(2),
      counter.wrapping_add(3),
    );
    let mut x13 = u32x4_splat(load_u32_le(&nonce[0..4]));
    let mut x14 = u32x4_splat(load_u32_le(&nonce[4..8]));
    let mut x15 = u32x4_splat(load_u32_le(&nonce[8..12]));

    let o0 = x0;
    let o1 = x1;
    let o2 = x2;
    let o3 = x3;
    let o4 = x4;
    let o5 = x5;
    let o6 = x6;
    let o7 = x7;
    let o8 = x8;
    let o9 = x9;
    let o10 = x10;
    let o11 = x11;
    let o12 = x12;
    let o13 = x13;
    let o14 = x14;
    let o15 = x15;

    let mut round = 0usize;
    while round < 10 {
      quarter_round(&mut x0, &mut x4, &mut x8, &mut x12);
      quarter_round(&mut x1, &mut x5, &mut x9, &mut x13);
      quarter_round(&mut x2, &mut x6, &mut x10, &mut x14);
      quarter_round(&mut x3, &mut x7, &mut x11, &mut x15);

      quarter_round(&mut x0, &mut x5, &mut x10, &mut x15);
      quarter_round(&mut x1, &mut x6, &mut x11, &mut x12);
      quarter_round(&mut x2, &mut x7, &mut x8, &mut x13);
      quarter_round(&mut x3, &mut x4, &mut x9, &mut x14);

      round = round.strict_add(1);
    }

    x0 = u32x4_add(x0, o0);
    x1 = u32x4_add(x1, o1);
    x2 = u32x4_add(x2, o2);
    x3 = u32x4_add(x3, o3);
    x4 = u32x4_add(x4, o4);
    x5 = u32x4_add(x5, o5);
    x6 = u32x4_add(x6, o6);
    x7 = u32x4_add(x7, o7);
    x8 = u32x4_add(x8, o8);
    x9 = u32x4_add(x9, o9);
    x10 = u32x4_add(x10, o10);
    x11 = u32x4_add(x11, o11);
    x12 = u32x4_add(x12, o12);
    x13 = u32x4_add(x13, o13);
    x14 = u32x4_add(x14, o14);
    x15 = u32x4_add(x15, o15);

    let mut words = [[0u32; BLOCKS_PER_BATCH]; 16];
    // SAFETY: each destination is a valid four-lane `u32` array for one unaligned `v128` store.
    unsafe {
      v128_store(words[0].as_mut_ptr() as *mut v128, x0);
      v128_store(words[1].as_mut_ptr() as *mut v128, x1);
      v128_store(words[2].as_mut_ptr() as *mut v128, x2);
      v128_store(words[3].as_mut_ptr() as *mut v128, x3);
      v128_store(words[4].as_mut_ptr() as *mut v128, x4);
      v128_store(words[5].as_mut_ptr() as *mut v128, x5);
      v128_store(words[6].as_mut_ptr() as *mut v128, x6);
      v128_store(words[7].as_mut_ptr() as *mut v128, x7);
      v128_store(words[8].as_mut_ptr() as *mut v128, x8);
      v128_store(words[9].as_mut_ptr() as *mut v128, x9);
      v128_store(words[10].as_mut_ptr() as *mut v128, x10);
      v128_store(words[11].as_mut_ptr() as *mut v128, x11);
      v128_store(words[12].as_mut_ptr() as *mut v128, x12);
      v128_store(words[13].as_mut_ptr() as *mut v128, x13);
      v128_store(words[14].as_mut_ptr() as *mut v128, x14);
      v128_store(words[15].as_mut_ptr() as *mut v128, x15);
    }

    let mut block_index = 0usize;
    while block_index < BLOCKS_PER_BATCH {
      let mut word_index = 0usize;
      while word_index < 16 {
        let offset = block_index.strict_mul(BLOCK_SIZE).strict_add(word_index.strict_mul(4));
        let keystream = words[word_index][block_index].to_le_bytes();
        chunk[offset..offset.strict_add(4)]
          .iter_mut()
          .zip(keystream)
          .for_each(|(dst, src)| *dst ^= src);
        word_index = word_index.strict_add(1);
      }
      block_index = block_index.strict_add(1);
    }

    counter = counter.wrapping_add(BLOCKS_PER_BATCH as u32);
  }

  let remainder = batches.into_remainder();
  if !remainder.is_empty() {
    xor_keystream_portable(key, counter, nonce, remainder);
  }
}

#[inline(always)]
fn quarter_round(a: &mut v128, b: &mut v128, c: &mut v128, d: &mut v128) {
  *a = u32x4_add(*a, *b);
  *d = rotl(v128_xor(*d, *a), 16);
  *c = u32x4_add(*c, *d);
  *b = rotl(v128_xor(*b, *c), 12);
  *a = u32x4_add(*a, *b);
  *d = rotl(v128_xor(*d, *a), 8);
  *c = u32x4_add(*c, *d);
  *b = rotl(v128_xor(*b, *c), 7);
}

#[inline(always)]
fn rotl(value: v128, bits: u32) -> v128 {
  v128_or(u32x4_shl(value, bits), u32x4_shr(value, 32u32.wrapping_sub(bits)))
}