"julia-lang: Ce cas n'est pas encore supporté"

Question

Je m'entraînais avec un peu d'exercice comme P3.jl. Quand je cours le code suivant, il n'y a aucune erreur et la réponse est correcte.

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

Mais quand je lance avec N = 16386, il a une erreur:

This case is not supported yet

Comment puis-je résoudre?

Merci.

Answer 1

Dans cette partie du code:

@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 

Lorsque vous appelez la première fois avec k = N-2 (ie N = 2048 ce qui déclenche l'erreur pour moi), k = 1023, et lorsque vous arrivez aux instructions 'spawn', dans le deuxième appel 'k' devient 511 et k2 devient 512, donc cette erreur est déclenchée.

Fait intéressant, l'erreur n'apparaît pas si vous exécutez sans ajouter d'autres processus de travail en premier. On peut supposer que la macro @spawn se comporte différemment/de façon inattendue quand il n'y a pas de travailleurs définis, je ne suis pas sûr. Quoi qu'il en soit, je ne pense pas que le bit lrefs fait ce que l'auteur pense faire (peut-être l'a-t-il fait dans une ancienne version de julia?)