Hi@all.
Ich habe ein großes "Array of Array of Array of double";
In diesem möchte ich per Finite-Differenzen-Methode die Lösung einer Poission Gleichung iterieren. Um eine Verbesserte Lösung an der Stelle (x,y,z) zu berechnen, muss auf alle Nachbarpunkte lesend, und auf die Stelle x,y,z schreibend zugegriffen werden.
Um das Verfahren zu beschleunigen, habe ich n
Threads eingeführt.
Jeder
Thread iteriert zunächst eine x-y-Ebene in diesem Array, erhöht dann z um 1 und iteriert die nächste x-y Ebene. Ist ein
Thread beim Maximalwert von z angekommen, wird z wieder auf 0 gesetzt.
Damit ein
Thread den anderen nicht "überholen" kann (Konsistenz der Lösung), fragt der i+1
Thread, ob die z Position des i-ten
Thread mindestens um 1 größer ist, als die Position des i+1-ten
Threads:
Delphi-Quellcode:
while (not Iterrationsgenauigkeit_erreicht) do
begin
inc(j);
for z := 0 to NodeCount_z-1 do
begin
current_z:=z;
i:=Thread_id-1;
if i=-1 then i:=Num_Thread-1;
if assigned(Threads[i]) then
while (Threads[i].Current_z>Current_z) and (Threads[i].Current_z-Current_z<=1) do
sleep(10);
for x := 0 to NodeCount_x-1 do
for y := 0 to NodeCount_y-1 do
begin
//Führe FDM-Iteration aus
end;
end;
end;
Zur kurzen und knappen Frage:
Benötige ich irgendwo eine Critical-Section?
Da ich ja nicht auf den selben Array-Eintrag mit zwei verschiedenen
Threads zugreife, gibts hier eig. kein Problem oder?
Benötige ich vielleicht für die variable current_z eine Critical-Section?
Besten Dank für eure Hilfe!
Euer Michael