Делается полный перебор по 4-х мерному массиву — выполняется очень долго, а нужно быстрее

 // Начинаем перебор
 for(i1=0;i1<=2*CellsNum;i1++)
 {tx = i1;
  for(i2=0;i2<=2*CellsNum;i2++)
  {ty = i2;
   for(i3=-Vg;i3<=Vg;i3++)
    {x = tx + i3*kof;
     for(i4=-Vg;i4<=Vg;i4++)
      {y = ty + i4*kof;
       Q1[i1][i2][i3+Vg][i4+Vg].dh += cl[x][y].h;

       if (Q1[i1][i2][i3][i4].dh<MinH)
        { // Сохраняем лучший вариант
         MinH = Q1[i1][i2][i3+Vg][i4+Vg].dh;
         r.x = i1;
         r.y = i2;
         r1.x = i3+Vg;
         r1.y = i4+Vg;
        }
      } // for i4
    } //for i3
  } // for i2
 } //for i1

аналогичная прога на паскале работает намного быстрее, но исходников нет (((

дайте пару советов по оптимизации такого кода