import std.traits : isFloatingPoint;

​

auto dot(T)(T[] x, T[] y) if (isFloatingPoint!T)

in (y.length == x.length)

{

    import std.range : iota;

    import std.parallelism : parallel, taskPool;

​

    auto sums = taskPool.workerLocalStorage(0.0L);

    foreach (i; parallel(iota(x.length)))

        sums.get += x[i] * y[i];

    T result = 0.0;

    foreach (threadResult; sums.toRange)

        result += threadResult;

    return result;

}

​

void main()

{

    double[] x = [1, 2, 3, 4, 5];

    double[] y = [6, 7, 8, 9, 10];

    assert(dot(x, y) == 130);

}