Модель частиц

[gusluk]

Т.е. вы пытались уравнение движения через полиномы Эрмита вывести?

[Pavia]

А как вы делали экстраполяцию через Эрмита? Ведь это метод для интерполяции и для его работу надо указать правую и левую границу. А правая нам неизвестна.

[type_Oleg]

Я давно это делал. Частности уже забыл. Помню, делал в разных системах - с центром в центре одного из тел ( обычно Солнца, или Земли ) или с центром в центре масс системы.
Вот, по этому документу - *.rar
Там , кстати очепятки нашел. Зеленым исправил, на рисунке.
Само вычисление по Эрмиту - на стр. 7, через промежуточные значения.

PS Точкой там производные обозначены. Например a с точкой - прозводная ускорения - рывок ( это только у Эрмита используется ).
r,r',a,a' - векторы конечно.

[gusluk]

Время от времени в программе вылазит ошибка доступа к памяти. Полагаю что неправильно работаю с потоками. Реализовано так:

Procedure ThreadEx_0(); stdcall;
var i,j,k:integer;
begin
k:=round(N/4);
for i:=1 to k do
    begin
      if (Atom[i].Mass>0) then
      begin
           for j:=1 to N do
             begin
             if (i<>j)and(Atom[j].Mass>0) then interactAtom(i,j);
             end;
      end;
    end;
end;
Procedure ThreadEx_1(); stdcall;
var i,j,k:integer;
begin
k:=round(N/4);
for i:=k+1 to 2*k do
    begin
      if (Atom[i].Mass>0) then
      begin
           for j:=1 to N do
             begin
             if (i<>j)and(Atom[j].Mass>0) then interactAtom(i,j);
             end;
      end;
    end;
end;
Procedure ThreadEx_2(); stdcall;
var i,j,k:integer;
begin
k:=round(N/4);
for i:=2*k+1 to 3*k do
    begin
      if (Atom[i].Mass>0) then
      begin
           for j:=1 to N do
             begin
             if (i<>j)and(Atom[j].Mass>0) then interactAtom(i,j);
             end;
      end;
    end;
end;
Procedure ThreadEx_3(); stdcall;
var i,j,k:integer;
begin
k:=round(N/4);
for i:=3*k+1 to 4*k do
    begin
      if (Atom[i].Mass>0) then
      begin
           for j:=1 to N do
             begin
             if (i<>j)and(Atom[j].Mass>0) then interactAtom(i,j);
             end;
      end;
    end;
end;

Procedure Step;
var  i,j:integer;  tid: cardinal; R:real;
begin
for i:=1 to N-1 do
      begin
      Atom[i].X:=Atom[i].X+Atom[i].u;
      Atom[i].Y:=Atom[i].Y+Atom[i].v;
      Atom[i].Z:=Atom[i].Z+Atom[i].w;
      end;
BMOut;

h[0] := CreateThread(nil, 0, @ThreadEx_0,nil, 0, tid);
h[1] := CreateThread(nil, 0, @ThreadEx_1,nil, 0, tid);
h[2] := CreateThread(nil, 0, @ThreadEx_2,nil, 0, tid);
h[3] := CreateThread(nil, 0, @ThreadEx_3,nil, 0, tid);
WaitForMultipleObjects(4, @h, true, INFINITE);
for i := 0 to 3 do CloseHandle(h[i]);

end;

Сам понимаю что криво, потому прошу пнуть в сторону правильного решения.
interactAtom приведен в первом сообщении.

[Pavia]

        for j:=1 to N do // Тут ошибка при J=N выход за приделы массива.
             begin
             if (i<>j)and(Atom[j].Mass>0) then interactAtom(i,j);

Почему не используете TThead?
Тут правильнее использовать ParallelFor в любой его реализации.
Какая версия Delphi? Для старых есть библиотека Omni
В новых смотри:
http://docwiki.embarcadero.com/RADSt...amming_Library
http://docwiki.embarcadero.com/RADSt...amming_Library

[gusluk]

Так как постоянно в разездах, а программирование для меня просто способ отвлечься от работы, пользуюсь portable версией Borland® Delphi® Version 10.0
TThead не стал использовать потому что решил для моей простой задачи достаточно и такого варианта. А чем он хуже использования TThead я не знаю.
Буду копать Parallel.For