Skip to content

Commit

Permalink
feat: initial implementation of simd
Browse files Browse the repository at this point in the history
  • Loading branch information
@twoeths
twoeths committed Mar 25, 2024
1 parent 9cf6991 commit d5e5495
Show file tree
Hide file tree
Showing 15 changed files with 13,586 additions and 105 deletions.
105 changes: 66 additions & 39 deletions packages/as-sha256/assembly/index.ts
View file
Original file line number Diff line number Diff line change
@@ -1,45 +1,14 @@
import {digest64V128} from "./simd";
import {K, W64} from "./utils/const";

// As of Mar 2024, simd support v128, and each block is 32 bits, so we can hash 4 inputs at once
export const PARALLEL_FACTOR = 4;
export const INPUT_LENGTH = 512;
//https://github.com/dchest/fast-sha256-js/blob/master/src/sha256.ts
const DIGEST_LENGTH = 32;
export const INPUT_LENGTH = 512;

// constants used in the SHA256 compression function
const K: u32[] = [
0x428a2f98, 0x71374491, 0xb5c0fbcf, 0xe9b5dba5, 0x3956c25b,
0x59f111f1, 0x923f82a4, 0xab1c5ed5, 0xd807aa98, 0x12835b01,
0x243185be, 0x550c7dc3, 0x72be5d74, 0x80deb1fe, 0x9bdc06a7,
0xc19bf174, 0xe49b69c1, 0xefbe4786, 0x0fc19dc6, 0x240ca1cc,
0x2de92c6f, 0x4a7484aa, 0x5cb0a9dc, 0x76f988da, 0x983e5152,
0xa831c66d, 0xb00327c8, 0xbf597fc7, 0xc6e00bf3, 0xd5a79147,
0x06ca6351, 0x14292967, 0x27b70a85, 0x2e1b2138, 0x4d2c6dfc,
0x53380d13, 0x650a7354, 0x766a0abb, 0x81c2c92e, 0x92722c85,
0xa2bfe8a1, 0xa81a664b, 0xc24b8b70, 0xc76c51a3, 0xd192e819,
0xd6990624, 0xf40e3585, 0x106aa070, 0x19a4c116, 0x1e376c08,
0x2748774c, 0x34b0bcb5, 0x391c0cb3, 0x4ed8aa4a, 0x5b9cca4f,
0x682e6ff3, 0x748f82ee, 0x78a5636f, 0x84c87814, 0x8cc70208,
0x90befffa, 0xa4506ceb, 0xbef9a3f7, 0xc67178f2
];
const kPtr = K.dataStart;

//precomputed W + K for message block representing length 64 bytes for fixed input of 64 bytes for digest64
const W64: u32[] = [
0xc28a2f98, 0x71374491, 0xb5c0fbcf, 0xe9b5dba5,
0x3956c25b, 0x59f111f1, 0x923f82a4, 0xab1c5ed5,
0xd807aa98, 0x12835b01, 0x243185be, 0x550c7dc3,
0x72be5d74, 0x80deb1fe, 0x9bdc06a7, 0xc19bf374,
0x649b69c1, 0xf0fe4786, 0x0fe1edc6, 0x240cf254,
0x4fe9346f, 0x6cc984be, 0x61b9411e, 0x16f988fa,
0xf2c65152, 0xa88e5a6d, 0xb019fc65, 0xb9d99ec7,
0x9a1231c3, 0xe70eeaa0, 0xfdb1232b, 0xc7353eb0,
0x3069bad5, 0xcb976d5f, 0x5a0f118f, 0xdc1eeefd,
0x0a35b689, 0xde0b7a04, 0x58f4ca9d, 0xe15d5b16,
0x007f3e86, 0x37088980, 0xa507ea32, 0x6fab9537,
0x17406110, 0x0d8cd6f1, 0xcdaa3b6d, 0xc0bbbe37,
0x83613bda, 0xdb48a363, 0x0b02e931, 0x6fd15ca7,
0x521afaca, 0x31338431, 0x6ed41a95, 0x6d437890,
0xc39c91f2, 0x9eccabbd, 0xb5c9a0e6, 0x532fb63c,
0xd2c741c6, 0x07237ea3, 0xa4954b68, 0x4c191d76,
];
const w64Ptr = W64.dataStart;

// intermediate hash values stored in H0-H7
Expand All @@ -53,8 +22,8 @@ var a: u32, b: u32, c: u32, d: u32, e: u32, f: u32, g: u32, h: u32, i: u32, t1:
const M = new ArrayBuffer(64);
const mPtr = changetype<usize>(M);

// 64 32bit extended message blocks
const W = new ArrayBuffer(256);
// 64 32bit extended message blocks = 64 * 4 = 256
const W = new ArrayBuffer(256 * PARALLEL_FACTOR);
const wPtr = changetype<usize>(W);

// input buffer
Expand Down Expand Up @@ -316,7 +285,7 @@ export function final(outPtr: usize): void {
export function digest(length: i32): void {
init();
update(inputPtr, length);
final(outputPtr)
final(outputPtr);
}

export function digest64(inPtr: usize, outPtr: usize): void {
Expand All @@ -332,3 +301,61 @@ export function digest64(inPtr: usize, outPtr: usize): void {
store32(outPtr, 6, bswap(H6));
store32(outPtr, 7, bswap(H7));
}

/**
* Hash 4 inputs of exactly 64 bytes each
* Input pointer is 256 bytes as below:
* byte u32
* input 0 0 1 2 ... 63 <===> 0 1 ... 15
* input 1 64 65 ... 127 <===> 16 17 ... 31
* input 2 128 ... 191 <===> 32 33 ... 47
* input 3 192 ... 255 <===> 48 49 ... 63
*
* we need to transfer it to expanded message blocks, with 16 first items like:
*
* w_v128_0 0 16 32 48
* w_v128_1 1 17 33 49
* ...
* w_v128_15 15 31 47 63
*
* remaining 48 items are computed inside hashBlocksV128 loop.
* @param outPtr
*/
export function hash4Inputs(outPtr: usize): void {
for (i = 0; i < 16; i++) {
store32(wPtr, PARALLEL_FACTOR * i, load32be(inputPtr, i));
store32(wPtr, PARALLEL_FACTOR * i + 1, load32be(inputPtr, i + 16));
store32(wPtr, PARALLEL_FACTOR * i + 2, load32be(inputPtr, i + 32));
store32(wPtr, PARALLEL_FACTOR * i + 3, load32be(inputPtr, i + 48));
}

digest64V128(wPtr, outPtr);
}

/*
* Hash 8 hash objects which are 4 inputs of 64 bytes each similar to hash4Inputs
*
* Input pointer is 64 u32 (256 bytes) as below:
* input 0 input 1 input 2 input 3
* h0 h0 h0 h0
* h1 h1 h1 h1
* ...
* h7 h7 h7 h7
* h0 h0 h0 h0
* h1 h1 h1 h1
* ...
* h7 h7 h7 h7
*
* that's already the setup for wInputPtr, we only need to load be value of them to make
* it the first 16 items of expanded message blocks
*
* remaining 48 items are computed inside hashBlocksV128 loop.
*
*/
export function hash8HashObjects(outPtr: usize): void {
for (i = 0; i < 16 * PARALLEL_FACTOR; i++) {
store32(wPtr, i, load32be(inputPtr, i));
}

digest64V128(wPtr, outPtr);
}
238 changes: 238 additions & 0 deletions packages/as-sha256/assembly/simd.ts
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,238 @@
import {K, W64} from "./utils/const";
import {getV128, setV128} from "./utils/v128";

let H0V128: v128, H1V128: v128, H2V128: v128, H3V128: v128, H4V128: v128, H5V128: v128, H6V128: v128, H7V128: v128;
let aV128: v128, bV128: v128, cV128: v128, dV128: v128, eV128: v128, fV128: v128, gV128: v128, hV128: v128, iV128: v128, t1V128: v128, t2V128: v128;
let i: i32
let tmpW: v128;

const kV128ArrayBuffer = new ArrayBuffer(4 * 64 * 4);
const kV128Ptr = changetype<usize>(kV128ArrayBuffer);
for (let i = 0; i < 64; i++) {
setV128(kV128Ptr, i, i32x4.splat(K[i]));
}

const w64V12ArrayBuffer = new ArrayBuffer(4 * 64 * 4);
const w64V128Ptr = changetype<usize>(w64V12ArrayBuffer);
for (let i = 0; i < 64; i++) {
setV128(w64V128Ptr, i, i32x4.splat(W64[i]));
}

const DEFAULT_H0V128 = i32x4.splat(0x6a09e667);
const DEFAULT_H1V128 = i32x4.splat(0xbb67ae85);
const DEFAULT_H2V128 = i32x4.splat(0x3c6ef372);
const DEFAULT_H3V128 = i32x4.splat(0xa54ff53a);
const DEFAULT_H4V128 = i32x4.splat(0x510e527f);
const DEFAULT_H5V128 = i32x4.splat(0x9b05688c);
const DEFAULT_H6V128 = i32x4.splat(0x1f83d9ab);
const DEFAULT_H7V128 = i32x4.splat(0x5be0cd19);

function initV128(): void {
H0V128 = DEFAULT_H0V128;
H1V128 = DEFAULT_H1V128;
H2V128 = DEFAULT_H2V128;
H3V128 = DEFAULT_H3V128;
H4V128 = DEFAULT_H4V128;
H5V128 = DEFAULT_H5V128;
H6V128 = DEFAULT_H6V128;
H7V128 = DEFAULT_H7V128;
}

/**
* wInputPtr: pointer to extended blocks
* block 0 (4 bytes) block 1 (4 bytes) block 2 (4 bytes) block 3 (4 bytes)
* wV128_0 |--------------------|--------------------|--------------------|--------------------|
* wV128_1 |--------------------|--------------------|--------------------|--------------------|
* ... ...
* wV128_15 |--------------------|--------------------|--------------------|--------------------| ===> end of input data, below is extended area
* wV128_16 based on item 0 to 14
* wV128_17 based on item 1 to 15
* ...
* wV128_63 based on item 47 to 62
*
* outPtr is 32 bytes each x 4 (PARALLEL_FACTOR) = 128 bytes
*/
export function digest64V128(wInputPtr: usize, outPtr: usize): void {
initV128();
hashBlocksV128(wInputPtr);
hashPreCompWV128();

// extract lane manually otherwise get "Expression must be a compile-time constant.""
store32(outPtr, 0, bswap(i32x4.extract_lane(H0V128, 0)));
store32(outPtr, 1, bswap(i32x4.extract_lane(H1V128, 0)));
store32(outPtr, 2, bswap(i32x4.extract_lane(H2V128, 0)));
store32(outPtr, 3, bswap(i32x4.extract_lane(H3V128, 0)));
store32(outPtr, 4, bswap(i32x4.extract_lane(H4V128, 0)));
store32(outPtr, 5, bswap(i32x4.extract_lane(H5V128, 0)));
store32(outPtr, 6, bswap(i32x4.extract_lane(H6V128, 0)));
store32(outPtr, 7, bswap(i32x4.extract_lane(H7V128, 0)));

store32(outPtr, 8, bswap(i32x4.extract_lane(H0V128, 1)));
store32(outPtr, 9, bswap(i32x4.extract_lane(H1V128, 1)));
store32(outPtr, 10, bswap(i32x4.extract_lane(H2V128, 1)));
store32(outPtr, 11, bswap(i32x4.extract_lane(H3V128, 1)));
store32(outPtr, 12, bswap(i32x4.extract_lane(H4V128, 1)));
store32(outPtr, 13, bswap(i32x4.extract_lane(H5V128, 1)));
store32(outPtr, 14, bswap(i32x4.extract_lane(H6V128, 1)));
store32(outPtr, 15, bswap(i32x4.extract_lane(H7V128, 1)));

store32(outPtr, 16, bswap(i32x4.extract_lane(H0V128, 2)));
store32(outPtr, 17, bswap(i32x4.extract_lane(H1V128, 2)));
store32(outPtr, 18, bswap(i32x4.extract_lane(H2V128, 2)));
store32(outPtr, 19, bswap(i32x4.extract_lane(H3V128, 2)));
store32(outPtr, 20, bswap(i32x4.extract_lane(H4V128, 2)));
store32(outPtr, 21, bswap(i32x4.extract_lane(H5V128, 2)));
store32(outPtr, 22, bswap(i32x4.extract_lane(H6V128, 2)));
store32(outPtr, 23, bswap(i32x4.extract_lane(H7V128, 2)));

store32(outPtr, 24, bswap(i32x4.extract_lane(H0V128, 3)));
store32(outPtr, 25, bswap(i32x4.extract_lane(H1V128, 3)));
store32(outPtr, 26, bswap(i32x4.extract_lane(H2V128, 3)));
store32(outPtr, 27, bswap(i32x4.extract_lane(H3V128, 3)));
store32(outPtr, 28, bswap(i32x4.extract_lane(H4V128, 3)));
store32(outPtr, 29, bswap(i32x4.extract_lane(H5V128, 3)));
store32(outPtr, 30, bswap(i32x4.extract_lane(H6V128, 3)));
store32(outPtr, 31, bswap(i32x4.extract_lane(H7V128, 3)));
}

/**
* Expand message blocks (16 32bit blocks), into extended message blocks (64 32bit blocks)
* There are 4 inputs, each input is 64 bytes which is 16 v128 objects of wInputPtr
* The first 16 v128 objects are computed before this function
* The remaining 48 v128 objects are computed from the first 16 v128 objects in this function
* Apply SHA256 compression function on extended message blocks
* Update intermediate hash values
* @param WV128 64 v128 objects respective to 4 expanded message blocks memory
* @param mV12Arr 16 v128 objects respective to 4 message blocks memory
*
*
*/
@inline
function hashBlocksV128(wInputPtr: usize): void {
// this is a copy of data
aV128 = H0V128;
bV128 = H1V128;
cV128 = H2V128;
dV128 = H3V128;
eV128 = H4V128;
fV128 = H5V128;
gV128 = H6V128;
hV128 = H7V128;

// Apply SHA256 compression function on expanded message blocks
for (i = 0; i < 64; i++) {
tmpW = i < 16 ? getV128(wInputPtr, i) : i32x4.add(i32x4.add(i32x4.add(SIG1V128(getV128(wInputPtr, i - 2)), getV128(wInputPtr, i - 7)),
SIG0V128(getV128(wInputPtr, i - 15))), getV128(wInputPtr, i - 16));
if (i >= 16) {
setV128(wInputPtr, i, tmpW);
}
t1V128 = i32x4.add(i32x4.add(i32x4.add(i32x4.add(hV128, EP1V128(eV128)), CHV128(eV128, fV128, gV128)), getV128(kV128Ptr, i)), tmpW);
t2V128 = i32x4.add(EP0V128(aV128), MAJV128(aV128, bV128, cV128));
hV128 = gV128;
gV128 = fV128;
fV128 = eV128;
eV128 = i32x4.add(dV128, t1V128);
dV128 = cV128;
cV128 = bV128;
bV128 = aV128;
aV128 = i32x4.add(t1V128, t2V128);
}

H0V128 = i32x4.add(H0V128, aV128);
H1V128 = i32x4.add(H1V128, bV128);
H2V128 = i32x4.add(H2V128, cV128);
H3V128 = i32x4.add(H3V128, dV128);
H4V128 = i32x4.add(H4V128, eV128);
H5V128 = i32x4.add(H5V128, fV128);
H6V128 = i32x4.add(H6V128, gV128);
H7V128 = i32x4.add(H7V128, hV128);
}

function hashPreCompWV128(): void {
aV128 = H0V128;
bV128 = H1V128;
cV128 = H2V128;
dV128 = H3V128;
eV128 = H4V128;
fV128 = H5V128;
gV128 = H6V128;
hV128 = H7V128;

// Apply SHA256 compression function on expanded message blocks
for (i = 0; i < 64; i++) {
t1V128 = i32x4.add(i32x4.add(i32x4.add(hV128, EP1V128(eV128)), CHV128(eV128, fV128, gV128)), getV128(w64V128Ptr, i));
t2V128 = i32x4.add(EP0V128(aV128), MAJV128(aV128, bV128, cV128));
hV128 = gV128;
gV128 = fV128;
fV128 = eV128;
eV128 = i32x4.add(dV128, t1V128);
dV128 = cV128;
cV128 = bV128;
bV128 = aV128;
aV128 = i32x4.add(t1V128, t2V128);
}

H0V128 = i32x4.add(H0V128, aV128);
H1V128 = i32x4.add(H1V128, bV128);
H2V128 = i32x4.add(H2V128, cV128);
H3V128 = i32x4.add(H3V128, dV128);
H4V128 = i32x4.add(H4V128, eV128);
H5V128 = i32x4.add(H5V128, fV128);
H6V128 = i32x4.add(H6V128, gV128);
H7V128 = i32x4.add(H7V128, hV128);
}

@inline
function CHV128(x: v128, y: v128, z: v128): v128 {
return v128.xor(v128.and(x, y), v128.and(v128.not(x), z));
}

@inline
function MAJV128(x: v128, y: v128, z: v128): v128 {
return v128.xor(v128.xor(v128.and(x, y), v128.and(x, z)), v128.and(y, z));
}

@inline
function EP0V128(x: v128): v128 {
return v128.xor(v128.xor(rotrV128(x, 2), rotrV128(x, 13)), rotrV128(x, 22));
}

@inline
function EP1V128(x: v128): v128 {
return v128.xor(v128.xor(rotrV128(x, 6), rotrV128(x, 11)), rotrV128(x, 25));
}

@inline
function SIG0V128(x: v128): v128 {
return v128.xor(v128.xor(rotrV128(x, 7), rotrV128(x, 18)), i32x4.shr_u(x, 3));
}

@inline
function SIG1V128(x: v128): v128 {
return v128.xor(v128.xor(rotrV128(x, 17), rotrV128(x, 19)), i32x4.shr_u(x, 10));
}

/**
* rotr is not natively supported by v128 so we have to implement it manually
* @param value
* @param bits
* @returns
*/
@inline
function rotrV128(value: v128, bits: i32): v128 {
const maskBits = 32 - bits;

// Shift right (logical) each lane by 'bits'
const rightShifted = i32x4.shr_u(value, bits);

// Shift left each lane by (32 - bits) to handle the wrap-around part of rotation
const leftShifted = i32x4.shl(value, maskBits);

// Combine the shifted parts with bitwise OR to achieve rotation
return v128.or(rightShifted, leftShifted);
}

@inline
function store32(ptr: usize, offset: usize, u: u32): void {
store<u32>(ptr + (offset << alignof<u32>()), u);
}
Loading

0 comments on commit d5e5495

Please sign in to comment.