[drcham]

Добрый день всем. Требуется распараллелить вычисление интеграла. Алгоритм сам выглядит так:
Необходимо посчитать интеграл в диапазоне [a,b] f(t)=sin(x*t)dx (на деле формула в три раза длиннее).

Поскольку интервал очень большой и из-за специфики функции его расчет возможен только по расчету суммы на каждом периоде sin. Т.е. в программе сделан цикл по каждому периоду синуса и общая сумма потом складывается. Считает все это программа довольно медленно. Поэтому решено сделать расчет с использованием потоков.

Прочитал что есть класс TTask.
Сделал:

Код:

Для теста можно функцию integrator в любом виде представить.

var
     tasks: array of ITask;
     sums:array of double;
     i:integer;
     int_period_start,int_period_end:extended;
     total:extended;
begin
   Setlength (tasks ,norm_period_count); //число периодов синуса в [a,b]
   Setlength (sums ,norm_period_count);// массив куда суммы на каждом периоде записывать
   for i := 0 to trunc(norm_period_count)-1 do //цикл по периодам
   begin
     int_period_start:=i*period; //a 
     int_period_end:=(i+1)*period; //b
     // диапазон на котором интегрировать на каждом потоке
      Tasks[i] := TTask.Create (procedure()
        begin
          //в поток выносим вычисление.
          res:=Integrator(i*period, (i+1)*period,current_time);
          //current_time - это "t" или точка времени в которой надо посчитать интеграл. 
          //TInterlocked.Exchange(sums[i],res);
          //Result:=res;
//              TThread.Synchronize(nil,
//            procedure
//            begin
//                  sums[1]:=res;
//            end);
       end);
      Tasks[i].Start; //запускаем поток

    end;

  TTask.WaitForAll(tasks); //ждем завершения
  total:=0;
 for i := 0 to High(sums) do
  begin
    total:=total+sums[i];// считает сумму и выдаем результат.
  end;
   Result:=total;

К сожалению моих знаний катастрофически не хватает.
В текущем варианте функция даже работает но в массив sums ничего не записывает. Пробовал оба варианта. Они закомментированы сейчас.

Прошу помощи как отладить данный алгоритм. Может быть надо iFuture использовать... хотя я тоже пробовал, но не смог реализовать.

Повторюсь что задача состоит в том, чтобы запустить параллельно в отдельном потоке вычисление интегрирования на заданном участке. В ответ получить значение суммы и потом просто сложить их. Количество периодов будет расти и довольно сильно т.к. sin (2*pi/t), а t это время которое считает от 1мкс до 800мкс в среднем.

Заранее спасибо. Могу поблагодарить если потребуется.