Как написать параллельный цикл в julia?

у меня есть следующий код Julia, и я хотел бы распараллелить его.

using DistributedArrays

function f(x)
    return x^2;
end
y = DArray[]
@parallel for i in 1:100
    y[i] = f(i)
end
println(y)

выход DistributedArrays.DArray[]. Я хотел бы иметь значение y следующим образом:y=[1,4,9,16,...,10000]

12
julia-lang parallel-processing

2 ответов

Привет Кира,

Я новичок в Джулии, но сталкиваюсь с той же проблемой. Попробуйте этот подход и посмотрите, соответствует ли он вашим потребностям.

function f(x)
  return x^2;
end

y=@parallel vcat for i= 1:100
 f(i);
end;

println(y)

С уважением, RN

3

вы можете использовать N-мерные представления распределенного массива:

сначала вам нужно добавить еще несколько процессов, локальных или удаленных:

julia> addprocs(CPU_CORES - 1);

затем вы должны использовать DistributedArrays на каждом из порожденных процессов:

julia> @everywhere using DistributedArrays

наконец, вы можете использовать @DArray макрос такой:

julia> x = @DArray [@show x^2 for x = 1:10];
        From worker 2:  x ^ 2 = 1
        From worker 2:  x ^ 2 = 4
        From worker 4:  x ^ 2 = 64
        From worker 2:  x ^ 2 = 9
        From worker 4:  x ^ 2 = 81
        From worker 4:  x ^ 2 = 100
        From worker 3:  x ^ 2 = 16
        From worker 3:  x ^ 2 = 25
        From worker 3:  x ^ 2 = 36
        From worker 3:  x ^ 2 = 49

вы можете видеть, что он делает то, что вы ожидаете:

julia> x
10-element DistributedArrays.DArray{Int64,1,Array{Int64,1}}:
   1
   4
   9
  16
  25
  36
  49
  64
  81
 100

помните, что он работает с произвольное количество измерений:

julia> y = @DArray [@show i + j for i = 1:3, j = 4:6];
        From worker 4:  i + j = 7
        From worker 4:  i + j = 8
        From worker 4:  i + j = 9
        From worker 2:  i + j = 5
        From worker 2:  i + j = 6
        From worker 2:  i + j = 7
        From worker 3:  i + j = 6
        From worker 3:  i + j = 7
        From worker 3:  i + j = 8

julia> y
3x3 DistributedArrays.DArray{Int64,2,Array{Int64,2}}:
 5  6  7
 6  7  8
 7  8  9

julia>

это самое Юлиан способ сделать то, что вы намеревались ИМХО.

мы можем посмотреть на macroexpand вывод, чтобы увидеть, что происходит:

Примечание: этот вывод был слегка отредактирован для удобства чтения,T означает:

DistributedArrays.Tuple{DistributedArrays.Vararg{DistributedArrays.UnitRange{DistributedArrays.Int}}}

julia> macroexpand(:(@DArray [i^2 for i = 1:10]))
  :(
    DistributedArrays.DArray(
      (
        #231#I::T -> begin
          [i ^ 2 for i = (1:10)[#231#I[1]]]
        end
      ),
      DistributedArrays.tuple(DistributedArrays.length(1:10))
    )
  )

что в основном то же самое, что и ручной ввод:

julia> n = 10; dims = (n,);

julia> DArray(x -> [i^2 for i = (1:n)[x[1]]], dims)
10-element DistributedArrays.DArray{Any,1,Array{Any,1}}:
   1
   4
   9
  16
  25
  36
  49
  64
  81
 100

julia>
12