Simple Hillis-Steele inclusive scan algorithm implementation.

gistfile1.d

V hillis_steele(alias op, V)(V array) {
    for(size_t i = 1; i < array.length; i *= 2) {
        for(size_t j = array.length; j > 0; --j) {
            if(j > i) {
                array[j-1] = op(array[j-1-i], array[j-1]);
            }
        }
    }
    return array;
}