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