Нужен знающий человек.

[Dimedrolum]

Имеются функции на с++

double F(double me, double N, double t0, double j0, double lm, double x, double T)
{
double REZ=0, sqrtN=sqrt(N);
 for(double nx=1; nx<=sqrtN/2; nx++)
  for(double ny=1; ny<=sqrtN/2; ny++)
   {
	if ((me+lm*((1+x)*t0*(cos(2*M_PI*nx/sqrtN)+cos(2*M_PI*ny/sqrtN))-j0*(1-x)/2))/T<700)
	REZ+=1/(exp((me+lm*((1+x)*t0*(cos(2*M_PI*nx/sqrtN)+cos(2*M_PI*ny/sqrtN))-j0*(1-x)/2))/T)+1);
   }
return (4*(1+x)/N)*REZ+x-1;
}

double ds(double N, double t0, double j0, double x, double T, double lm_left, double lm_right, double lm_step)

{
double REZ=0e0, Meps0=T*log(2*x/(1-x)), sqrtN=sqrt(N),a,b;
 for(double lm=lm_left;lm<=lm_right;lm+=lm_step)
  {
   double Meps=Me(lm, x, T);
	for(double nx=1;nx<=sqrtN/2;nx++)
	  for(double ny=1;ny<=sqrtN/2;ny++)
	   {
		a=t0*(1+x)*(cos(2*M_PI*nx/sqrtN)+cos(2*M_PI*ny/sqrtN))-j0*(1-x)/2;
		b=Meps+lm*a;
		if (b/T<700) REZ+=b*(b-Meps0)*exp(b/T)/(lm*T*T*(exp(b/T)+1)*(exp(b/T)+1));
	   }
  }
    return 4*lm_step*REZ/N;
}

double sum1(double N, double t0, double j0, double lm, double x, double T)
{
 double REZ=0.0,sqrtN=sqrt(N),a,b,Meps=Me(lm,x,T);
  for(double nx=1;nx<=sqrtN/2;nx++)
   for(double ny=1;ny<=sqrtN/2;ny++)
	{
	 a=t0*(1+x)*(cos(2*M_PI*nx/sqrtN)+cos(2*M_PI*ny/sqrtN))+j0*(1-x)/2; b=Meps-lm*a;
	 if (b/T<200)
	 REZ+=(t0*(1+x)*(cos(2*M_PI*nx/sqrtN)+cos(2*M_PI*ny/sqrtN)))*exp(b/T)/((exp(b/T)+1)*(exp(b/T)+1));
	};
 return 4*REZ;
}

double sum2(double N, double t0, double j0, double lm, double x, double T)
{
 double REZ=0.0,sqrtN=sqrt(N),a,b,Meps=Me(lm,x,T);
  for(double nx=1;nx<=sqrtN/2;nx++)
   for(double ny=1;ny<=sqrtN/2;ny++)
	{
	 a=t0*(1+x)*(cos(2*M_PI*nx/sqrtN)+cos(2*M_PI*ny/sqrtN))+j0*(1-x)/2;
	 b=Meps-lm*a;
	 if (b/T<200) REZ+=exp(b/T)/((exp(b/T)+1)*(exp(b/T)+1));
	};
 return 4*REZ;
}

Собственно вопрос, будет ли заметен прирост производительности, если эти функции сделать вставками на ассемблере? (Интересен прирост в разы) Если да, то кто за какую сумму может взяться?

[Black Fregat]

Мое личное мнение: просто от ассемблера прироста в разы не будет. Хороший компилятор должен оптимизировать такой код настолько, что руками намного лучше не сделаешь.

Можно попытаться оптимизировать математику под расширенные инструкции процессора, но навскидку не очевидно, что будет польза.

Можно попытаться покрутить алгоритмы преобразований, но это снова чисто математика.

[nworm]

1. Лучше всего лезть в математику самой задачи (если повезёт, можно добиться серьезного успеха).

2. Если ничего не выйдет, далее разумно работать с гтовой программой, например попробовать сделать цикл по Im внутренним в функции ds и т.д. Здесь опять же нужно учитывать особенности реальных данных, смотреть, какие циклы длинные, какие - нет.

3. Далее, выигрыш может дать выбор хорошего компилятора (но не большой, не в несколько раз ) и параметров компиляции в этом компиляторе.

4. Потом можно переводить на ассемблер, здесь также можно получить прирост (особенно если всё это считается на каких-то специфических процессорах), но всё равно лучше всего смотреть математику самой задачи.

Математику могу посмотреть, если хотите, пишите (e-mail такой же, как сайт из подписи).