🚀 go-pugleaf

On Sat, 26 Oct 2002, Recurrente wrote:

# �De verdad no creeis que usando rutinas en ensamblador se puede mejorar la
# velocidad del c�digo. mi rutina es as� de simple en C (por poner un
# ejemplo)
#
# v=array grande.
# for (i=0;i<Ngrande;i++) v[i]/=262144;

Podrias postear el codigo que genera tu compilador, para ese bucle.
En general cualquier buen compilador hara un trabajo mejor que tu,
a menos que seas un genio de la optimizacion :)

Asegurate de que tu compilador (y si no, hazlo tu) no DIVIDE, sino
que MULTIPLICA por 64^(-3), ya que ese numero lo puedes precalcular
antes de iniciar el bucle:

  double K = 1.0/262144.0;
  for(i = 0; i < Ngrande; i ++) v[i] *= K;

Tambien puedes hacer loop-unrolling, aunque eso lo hara el
compilador mejor que tu -- a menos que sepas el valor de Ngrande,
o alguna propidad de el --

Si, por ejemplo, Ngrande es multiplo de 4 *siempre*, puedes hacer:

for(i = 0; i < Ngrande; i += 4) {
	v[i] *= K;
	v[i + 1] *= K;
	v[i + 2] *= K;
	v[i + 3] *= K;
}

Mejor dejar el incremento de 'i' en el for(), por si el compilador
no optimiza bien el contenido del bucle, al haber dependencias:

  v[i ++] *= K;
      ^^ aqui cambia 'i', asi que la siguiente instruccion espera..
  v[i ++] *= K;

# 262144=64^3. Si os preguntais que porque necesito esto es que uso la
# libreria fftw para realizar transformadas r�pidas de Fourier (libreria
# excelente, os la recomiendo si necesitais alg�n d�a hacer FFT's) pero al
#  transformar-destransformar no se tienen los datos

No se realmente lo que haces, pero quiza fuera conveniente pre-escalar
el vector antes de hacer la fft(), de forma que al hacer la transformada
inversa tengas los datos con la escala que quieres.

Lepton.