тройной цикл

с ужасом понял что CUDA загибается от тройного цикла ,по крайней мере на длинных циклах ...
т.е. если i ,j ,k пустить до n =200 это ещё нормально ,но если 500 и больше то остановка видео драйвера и т.п.
как решить проблему?

файлик с проблемным участком прилагается )

AttachmentSize
Plain text icon cuda.txt507 bytes

Forums: 

вот собственно то место где начинаются страдания..

  1. for (int i=1;i<n;i++)
  2.         {
  3.                 for (int j=i;j<n;j++)
  4.                 {
  5.                         float sum=0;
  6.                 for (int k=0;k<i;k++)
  7.                         {
  8.                                 int idik=n*i + k ;
  9.                         int idjk= n*k + j ;
  10.                     sum=sum+L[idik]*U[idjk];
  11.                         }
  12.                 int id = n*i + j;U[id]=A[id]-sum;
  13.                 }
  14.          
  15.                 for (int j=i+1;j<n;j++)
  16.                 {
  17.                         float sum=0;
  18.                 for (int k=0;k<i;k++)
  19.                         {
  20.                                 int idik=n*j + k ;
  21.                         int idjk= n*k + i;
  22.                         sum+=L[idik]*U[idjk];
  23.                         }
  24.                         int idd= n*j + i;
  25.                         L[idd]=(A[idd]-sum)/U[i*n+i];
  26.                 }
  27.         }