"julia-lang: Dieser Fall wird noch nicht unterstützt"

Question

Ich übte mit etwas Übung wie P3.jl. Wenn ich folgenden Code ausführe, hat es keinen Fehler und die Antwort ist korrekt.

N = 10 
A = rand(N,N) 
SA = convert(SharedArray, A) 
T = 100 
Jacobi_parallel(SA,N,T) 

Aber wenn ich mit N = 16386 laufen, hat es einen Fehler:

This case is not supported yet

Wie kann ich es beheben?

Vielen Dank.

Answer 1

In diesem Teil des Codes:

@everywhere function Jacobi_blockwise_parallel_one_time_step(A,x,y,k,threshold) 
    if k <= threshold 
    Jacobi_serial_one_block_one_time_step(A,x,y,k) 
    else 
    k1 = floor(Int64, k/2) 
    k2 = k - k1 
    if (k1 != k2) error("This case is not supported yet") end 
    lrefs = [@spawn Jacobi_blockwise_parallel_one_time_step(A,x,y,k1,threshold), 
      @spawn Jacobi_blockwise_parallel_one_time_step(A,x+k1,y,k1,threshold), 
      @spawn Jacobi_blockwise_parallel_one_time_step(A,x,y+k1,k1,threshold), 
      @spawn Jacobi_blockwise_parallel_one_time_step(A,x+k1,y+k1,k1,threshold)]; 
    pmap(fetch, lrefs) 
    end 
end 

Wenn Sie es das erste Mal mit k nennen = N-2 (dh für N = 2048 diese den Fehler für mich löst), k = 1023 und Wenn Sie zu den 'spawn'-Anweisungen kommen, wird im zweiten Aufruf' k '511 und k2 wird 512, so dass dieser Fehler ausgelöst wird.

Interessanterweise wird der Fehler nicht angezeigt, wenn Sie ausführen, ohne zuerst zusätzliche Arbeitsprozesse hinzuzufügen. Vermutlich verhält sich das @ spawn-Makro anders/unerwartet, wenn es keine definierten Worker gibt, bin ich mir nicht sicher. Wie auch immer, ich glaube nicht, dass das lrefs Bit tut, was der Autor denkt, dass es tut (vielleicht tat es in einer älteren Version von Julia?)