Velocidad en programa

[rhervad]

Hola, estoy diseñando un programa que esta formado por multiples bucles en los que se hacen operaciones con arrays tridimensionales. Unas de las operaciones que se llevan a cabo dentro de los bucles es:

for ( i=1; i < ia; i++ ) {
for ( j=1; j < JE; j++ ) {
for ( k=1; k < KE; k++ ) {

curl_h = ( hz[i][j][k] - hz[i][j-1][k] - hy[i][j][k] + hy[i][j][k-1]);
}}}

Cambie el programa para utilizar arrays unidimensionales, pero acceder como si fueran matrices. El bucle anterior con el nuevo diseño seria:

for ( i=1; i < ia; i++ ) {
for ( j=1; j < JE; j++ ) {
for ( k=1; k < KE; k++ ) {

curl_h = hz.array[i+(j*hz.X)+(k*hz.X*hz.Y)] - hz.array[i+((j-1)*hz.X)+(k*hz.X*hz.Y)] - hy.array[i+(j*hy.X)+(k*hy.X*hy.Y)] + hy.array[i+(j*hy.X)+((k-1)*hy.X*hy.Y)];
}}}

Mi pregunta es si es normal que la velocidad de la ejecucion del programa llegue a multiplicarse entre 4 y 5 veces. Gracias.

[aloqui]

Yo diría que no es normal.
Internamente los compiladores transforman los arrays multidimensionales en unidimensionales ya que la memoria es unidimensional.

La única explicación que se me ocurre es que estes compilando con flags de depuración (o sin flags de optimización) y que el compilador incluya para cada indice de un array una funcion de comprobación de límites (para cercionarse de que no accedes a un índice fuera del array). De ése modo en los arrays multidimensionales dicha comprobación se hace 3 veces y en los arrays unidimensionales solo 1 vez.