in src/main/java/org/apache/datasketches/theta/BitPacking.java [2837:2909]
static void packBits57(final long[] values, final int i, final byte[] buf, int off) {
buf[off++] = (byte) (values[i + 0] >>> 49);
buf[off++] = (byte) (values[i + 0] >>> 41);
buf[off++] = (byte) (values[i + 0] >>> 33);
buf[off++] = (byte) (values[i + 0] >>> 25);
buf[off++] = (byte) (values[i + 0] >>> 17);
buf[off++] = (byte) (values[i + 0] >>> 9);
buf[off++] = (byte) (values[i + 0] >>> 1);
buf[off] = (byte) (values[i + 0] << 7);
buf[off++] |= values[i + 1] >>> 50;
buf[off++] = (byte) (values[i + 1] >>> 42);
buf[off++] = (byte) (values[i + 1] >>> 34);
buf[off++] = (byte) (values[i + 1] >>> 26);
buf[off++] = (byte) (values[i + 1] >>> 18);
buf[off++] = (byte) (values[i + 1] >>> 10);
buf[off++] = (byte) (values[i + 1] >>> 2);
buf[off] = (byte) (values[i + 1] << 6);
buf[off++] |= values[i + 2] >>> 51;
buf[off++] = (byte) (values[i + 2] >>> 43);
buf[off++] = (byte) (values[i + 2] >>> 35);
buf[off++] = (byte) (values[i + 2] >>> 27);
buf[off++] = (byte) (values[i + 2] >>> 19);
buf[off++] = (byte) (values[i + 2] >>> 11);
buf[off++] = (byte) (values[i + 2] >>> 3);
buf[off] = (byte) (values[i + 2] << 5);
buf[off++] |= values[i + 3] >>> 52;
buf[off++] = (byte) (values[i + 3] >>> 44);
buf[off++] = (byte) (values[i + 3] >>> 36);
buf[off++] = (byte) (values[i + 3] >>> 28);
buf[off++] = (byte) (values[i + 3] >>> 20);
buf[off++] = (byte) (values[i + 3] >>> 12);
buf[off++] = (byte) (values[i + 3] >>> 4);
buf[off] = (byte) (values[i + 3] << 4);
buf[off++] |= values[i + 4] >>> 53;
buf[off++] = (byte) (values[i + 4] >>> 45);
buf[off++] = (byte) (values[i + 4] >>> 37);
buf[off++] = (byte) (values[i + 4] >>> 29);
buf[off++] = (byte) (values[i + 4] >>> 21);
buf[off++] = (byte) (values[i + 4] >>> 13);
buf[off++] = (byte) (values[i + 4] >>> 5);
buf[off] = (byte) (values[i + 4] << 3);
buf[off++] |= values[i + 5] >>> 54;
buf[off++] = (byte) (values[i + 5] >>> 46);
buf[off++] = (byte) (values[i + 5] >>> 38);
buf[off++] = (byte) (values[i + 5] >>> 30);
buf[off++] = (byte) (values[i + 5] >>> 22);
buf[off++] = (byte) (values[i + 5] >>> 14);
buf[off++] = (byte) (values[i + 5] >>> 6);
buf[off] = (byte) (values[i + 5] << 2);
buf[off++] |= values[i + 6] >>> 55;
buf[off++] = (byte) (values[i + 6] >>> 47);
buf[off++] = (byte) (values[i + 6] >>> 39);
buf[off++] = (byte) (values[i + 6] >>> 31);
buf[off++] = (byte) (values[i + 6] >>> 23);
buf[off++] = (byte) (values[i + 6] >>> 15);
buf[off++] = (byte) (values[i + 6] >>> 7);
buf[off] = (byte) (values[i + 6] << 1);
buf[off++] |= values[i + 7] >>> 56;
buf[off++] = (byte) (values[i + 7] >>> 48);
buf[off++] = (byte) (values[i + 7] >>> 40);
buf[off++] = (byte) (values[i + 7] >>> 32);
buf[off++] = (byte) (values[i + 7] >>> 24);
buf[off++] = (byte) (values[i + 7] >>> 16);
buf[off++] = (byte) (values[i + 7] >>> 8);
buf[off] = (byte) (values[i + 7]);
}