Foro - Perl en Español

Hola gente, escribo el post en esta sección para asegurarme de obtener una respuesta segura. Sucede que estoy programando una simulación "cosmológica" en Perl, y necesito conseguir la mejor eficiencia del código para que el script "vuele". Sucede que estoy trabajando con una matriz 3D de 1000x1000x3, y necesito saber cómo recorrerla con eficiencia... ya que también trabajo en ella otras 1000 veces más.

El programa me está demorando sus 6000 seg. aprox. en Ubuntu.

El ciclo es así:


Agradezco de antemano cualquier ayuda... eso.

No es posible mejorar lo que has escrito, salvo que algunas de las operaciones que has escrito entre medias pueda ser cacheada con la ayuda del módulo Memoize.

La respuesta del foro Experto es... use PDL; Aprenderlo no es sencillo; tanto como aprender cualquier otra librería matemática; pero las prestaciones son tremendas.

Usa el sistema de paquetes de Ubuntu para localizar perl-pdl o perlpdl para instalarle. Luego, únete a su lista de correo. Existen un par de manuales, uno de ellos de introducción de Xavier Calbet.

Aquí tienes un ejemplo de uso.

Hay más comentarios sobre PDL en estos foros.

De todas formas, si nos das más detalles de lo que quieres hacer, podremos darte más pistas.

He visto que has modificado el mensaje casi al mismo tiempo de que respondiera.

Por las operaciones que has puesto:

se ve que es un ejemplo muy bueno para el uso de PDL.

Mientras tanto, se puede acelerar un poco el programa, optimizando de la misma forma que lo haríamos en C:


Bueno, la idea es intentar reducir los accesos a estructuras de datos, cacheando operaciones y valores intermedios. En algunas ocasiones, se aprecia un aumento de la velocidad de procesado.

De todas maneras, como ya digo, esto es más bien trabajo de PDL. O incluso hasta se podría hacer en Inline::C.

¿De qué operaciones estamos hablando? Sería interesante verlas...

Un truco simple consiste en cambiar el orden de las dimensiones en la matriz:

Con tu esquema de 1000x1000x3 utilizas 1 millón de arrays para almacenar la matriz; si lo cambiases por 3x1000x1000, utilizarías solo 3000 arrays. Esto supondría un ahorro de memoria probablemente cercano al 50% que repercutiría en un aumento de la velocidad.

Otra solución es utilizar un único array para guardar toda la matriz y calcular el índice de cada elemento a mano (por ejemplo: $a[$i*3000+$j*3+$k]). El cálculo del índice, más o menos, tendrá un coste similar a las indirecciones, pero en la mayoría de los algoritmos se puede simplificar en función de los valores que tome en las iteraciones anteriores de los bucles con lo que podría ser más rápido. Una representación a medio camino con una matriz de 3x1000_000 también podría estar bien en algunos casos.

Bueno, no sé si se entiende lo que trato de decir demasiado bien... de todas formas mi consejo final también es que ¡¡¡uses PDL!!!