assemblyscript/std/assembly/internal/sort.ts

import {
  LOAD,
  STORE
} from "./arraybuffer";

import {
  compareUnsafe
} from "./string";

/** Obtains the default comparator for the specified value type. */
@inline
export function COMPARATOR<T>(): (a: T, b: T) => i32 {
  if (isInteger<T>()) {
    if (isSigned<T>() && sizeof<T>() <= 4) {
      return (a: T, b: T): i32 => (<i32>(a - b));
    } else {
      return (a: T, b: T): i32 => (<i32>(a > b) - <i32>(a < b));
    }
  } else if (isFloat<T>()) {
    if (sizeof<T>() == 4) {
      return (a: T, b: T): i32 => {
        var ia = reinterpret<i32>(a);
        var ib = reinterpret<i32>(b);
        ia ^= (ia >> 31) >>> 1;
        ib ^= (ib >> 31) >>> 1;
        return <i32>(ia > ib) - <i32>(ia < ib);
      };
    } else {
      return (a: T, b: T): i32 => {
        var ia = reinterpret<i64>(a);
        var ib = reinterpret<i64>(b);
        ia ^= (ia >> 63) >>> 1;
        ib ^= (ib >> 63) >>> 1;
        return <i32>(ia > ib) - <i32>(ia < ib);
      };
    }
  } else if (isString<T>()) {
    return (a: T, b: T): i32 => {
      if (a === b || a === null || b === null) return 0;
      var alen = (<string>a).length;
      var blen = (<string>b).length;
      if (!alen && !blen) return 0;
      if (!alen) return -1;
      if (!blen) return  1;
      return compareUnsafe(<string>a, 0, <string>b, 0, <usize>min(alen, blen));
    };
  } else {
    return (a: T, b: T): i32 => (<i32>(a > b) - <i32>(a < b));
  }
}

@inline
export function SORT<T>(
  buffer: ArrayBuffer,
  byteOffset: i32,
  length: i32,
  comparator: (a: T, b: T) => i32
): void {
  if (isReference<T>()) {
    // TODO replace this to faster stable sort (TimSort) when it implemented
    insertionSort<T>(buffer, byteOffset, length, comparator);
  } else {
    if (length < 256) {
      insertionSort<T>(buffer, byteOffset, length, comparator);
    } else {
      weakHeapSort<T>(buffer, byteOffset, length, comparator);
    }
  }
}

/** Sorts an Array with the 'Insertion Sort' algorithm. */
function insertionSort<T>(
  buffer: ArrayBuffer,
  byteOffset: i32,
  length: i32,
  comparator: (a: T, b: T) => i32
): void {
  for (let i = 0; i < length; i++) {
    let a = LOAD<T>(buffer, i, byteOffset);       // a = arr[i]
    let j = i - 1;
    while (j >= 0) {
      let b = LOAD<T>(buffer, j, byteOffset);     // b = arr[j]
      if (comparator(a, b) < 0) {
        STORE<T>(buffer, j-- + 1, b, byteOffset); // arr[j + 1] = b
      } else break;
    }
    STORE<T>(buffer, j + 1, a, byteOffset);       // arr[j + 1] = a
  }
}

/** Sorts an Array with the 'Weak Heap Sort' algorithm. */
function weakHeapSort<T>(
  buffer: ArrayBuffer,
  byteOffset: i32,
  length: i32,
  comparator: (a: T, b: T) => i32
): void {
  const shift32 = alignof<u32>();

  var bitsetSize = (length + 31) >> 5 << shift32;
  var bitset = memory.allocate(bitsetSize); // indexed in 32-bit chunks below
  memory.fill(bitset, 0, bitsetSize);

  // see: http://citeseerx.ist.psu.edu/viewdoc/download?doi=10.1.1.21.1863&rep=rep1&type=pdf

  for (let i = length - 1; i > 0; i--) {
    let j = i;
    while ((j & 1) == (load<u32>(bitset + (j >> 6 << shift32)) >> (j >> 1 & 31) & 1)) j >>= 1;

    let p = j >> 1;
    let a = LOAD<T>(buffer, p, byteOffset);  // a = arr[p]
    let b = LOAD<T>(buffer, i, byteOffset);  // b = arr[i]
    if (comparator(a, b) < 0) {
      store<u32>(
        bitset + (i >> 5 << shift32),
        load<u32>(bitset + (i >> 5 << shift32)) ^ (1 << (i & 31))
      );
      STORE<T>(buffer, i, a, byteOffset);    // arr[i] = a
      STORE<T>(buffer, p, b, byteOffset);    // arr[p] = b
    }
  }

  for (let i = length - 1; i >= 2; i--) {
    let a = LOAD<T>(buffer, 0, byteOffset);
    STORE<T>(buffer, 0, LOAD<T>(buffer, i, byteOffset), byteOffset);
    STORE<T>(buffer, i, a, byteOffset);

    let x = 1, y: i32;
    while ((y = (x << 1) + ((load<u32>(bitset + (x >> 5 << shift32)) >> (x & 31)) & 1)) < i) x = y;

    while (x > 0) {
      a = LOAD<T>(buffer, 0, byteOffset);     // a = arr[0]
      let b = LOAD<T>(buffer, x, byteOffset); // b = arr[x]

      if (comparator(a, b) < 0) {
        store<u32>(
          bitset + (x >> 5 << shift32),
          load<u32>(bitset + (x >> 5 << shift32)) ^ (1 << (x & 31))
        );
        STORE<T>(buffer, x, a, byteOffset);    // arr[x] = a
        STORE<T>(buffer, 0, b, byteOffset);    // arr[0] = b
      }
      x >>= 1;
    }
  }

  memory.free(bitset);

  var t = LOAD<T>(buffer, 1, byteOffset); // t = arr[1]
  STORE<T>(buffer, 1, LOAD<T>(buffer, 0, byteOffset), byteOffset);
  STORE<T>(buffer, 0, t, byteOffset); // arr[0] = t
}