algorithm
Algoritmi multithread
Ricerca…
introduzione
Esempi per alcuni algoritmi multithread.
Sintassi
- parallelo prima di un loop significa che ogni iterazione del loop è indipendente l'una dall'altra e può essere eseguita in parallelo.
- spawn indica la creazione di un nuovo thread.
- la sincronizzazione è sincronizzare tutti i thread creati.
- Array / matrice sono indicizzati da 1 a n negli esempi.
Multithread di moltiplicazione a matrice quadrata
multiply-square-matrix-parallel(A, B)
n = A.lines
C = Matrix(n,n) //create a new matrix n*n
parallel for i = 1 to n
parallel for j = 1 to n
C[i][j] = 0
pour k = 1 to n
C[i][j] = C[i][j] + A[i][k]*B[k][j]
return C
Multithread vettoriale a matrice di moltiplicazione
matrix-vector(A,x)
n = A.lines
y = Vector(n) //create a new vector of length n
parallel for i = 1 to n
y[i] = 0
parallel for i = 1 to n
for j = 1 to n
y[i] = y[i] + A[i][j]*x[j]
return y
unire il multithread
A è un array e indici p e q della matrice, come ad esempio ordinare il sub-array A [p..r] . B è un sotto-array che verrà popolato dall'ordinamento.
Una chiamata a p-merge-sort (A, p, r, B, s) ordina gli elementi da A [p..r] e li inserisce in B [s..s + rp] .
p-merge-sort(A,p,r,B,s)
n = r-p+1
if n==1
B[s] = A[p]
else
T = new Array(n) //create a new array T of size n
q = floor((p+r)/2))
q_prime = q-p+1
spawn p-merge-sort(A,p,q,T,1)
p-merge-sort(A,q+1,r,T,q_prime+1)
sync
p-merge(T,1,q_prime,q_prime+1,n,B,s)
Ecco la funzione ausiliaria che esegue l'unione in parallelo.
p-merge presuppone che le due sotto-matrici da unire si trovano nello stesso array ma non assumono che siano adiacenti nell'array. Ecco perché abbiamo bisogno di p1, r1, p2, r2 .
p-merge(T,p1,r1,p2,r2,A,p3)
n1 = r1-p1+1
n2 = r2-p2+1
if n1<n2 //check if n1>=n2
permute p1 and p2
permute r1 and r2
permute n1 and n2
if n1==0 //both empty?
return
else
q1 = floor((p1+r1)/2)
q2 = dichotomic-search(T[q1],T,p2,r2)
q3 = p3 + (q1-p1) + (q2-p2)
A[q3] = T[q1]
spawn p-merge(T,p1,q1-1,p2,q2-1,A,p3)
p-merge(T,q1+1,r1,q2,r2,A,q3+1)
sync
Ed ecco la funzione dicotomica di ricerca ausiliaria.
x è la chiave da cercare nell'array secondario T [p..r].
dichotomic-search(x,T,p,r)
inf = p
sup = max(p,r+1)
while inf<sup
half = floor((inf+sup)/2)
if x<=T[half]
sup = half
else
inf = half+1
return sup
Modified text is an extract of the original Stack Overflow Documentation
Autorizzato sotto CC BY-SA 3.0
Non affiliato con Stack Overflow