Julia function for calling threads

thread.jl

v = zeros(10)        # the "output" of the thread computation
index1 = [1,2,4,8]   # the elements handled by thread 1
index2 = [3,5,6]     # the elements handled by thread 2
indexwarmup = [7,10]    # elements handled on warmup

# The function run by each thread
function fillvals(out, index, val)
    out[index] = val
    return nothing
end

# This is a version that tries to avoid any type of allocation
# (without success). Is it the return statement?
function fillvals_forloop(out, index, val)
    for i = index
        arrayset(out, i, val)
    end
    return nothing
end

# Run the function once to make sure it is compiled
fillvals_forloop(v, indexwarmup, -1.0)

function runthreads(v, index1, index2)
    # Launch the threads
    thread1 = thread_create(fillvals_forloop, v, index1, 1.0)
    thread2 = thread_create(fillvals_forloop, v, index2, 2.0)
    
    # Wait until each thread has completed its task
    thread_join(thread1)
    thread_join(thread2)
    return v
end