[pyar] Programación científica en python

[Octavio Bruzzone]

On Dom 22 Ago 2010 04:26:47 Nico Chiaraviglio escribió:
> Nico la pregunta es qué hacés dentro del loop?
> 
> Porque ciertamente lo mejor es tratar de resolverlo con objetos de
> NumPy que encapsulan la iteracion (lo que llaman 'vectorizacion' y
> 'broadcast'). La idea es evitar el costo de loopear sobre vectores  a
> nivel Python y  en cambio que eso se haga sólo mediantes las rutinas
> compiladas en C de numpy,
> 
> Fijate la explicación que da en la intro:
> http://docs.scipy.org/doc/numpy/user/whatisnumpy.html
> 
> saludos
> 
> 
> En realidad lo que voy a loopear va a ser variable de acuerdo al problema.
> Siempre van a ser simulaciones neutrónicas y para que se den una idea
> cantidades típicas en un reactor son del orden de 10^9, 10^10 n/cm3,
> entonces simular 10^8 historias en un montecarlo no es ni cerca
> descabellado.

En ese caso no te serviría tampoco numpy.

> 
> La consulta que hice fue por una especie de decepción que tuve ya que hasta
> convencí al profesor (un físico que usa fortran desde antes que exista la
> computadora) a que pruebe python y cuando voy a ver los tiempos de
> ejecución encontré resultados no muy agradables.

Bueno, eso es común cuando te encuentras con alguien que programa en Fortran 
(yo también estoy tratando de convencer a uno de esos de usar python y la 
velocidad es un tema clave para esa gente), en general python es muy superior 
a Fortran para la parte de E/S de datos, para graficar y todo eso, pero cuando 
se trata de "triturar números", no hay nada que hacerle, ni siquiera con C 
puedes ganarle al viejo Fortran. Debes convencer al profesor de usar f2py para 
reusar sus viejas rutinas, y dejar las cosas de mas alto nivel a Python.

> 
> La idea tampoco es hacer las cosas demasiado complicadas y correr parte en
> un lenguaje y parte en otro, venía contento con python y creo que lo voy a
> usar igual (no estoy haciendo un desarrollo de alguna aplicación no
> perecedera, solo busco aprobar una materia) pero fue una sorpresa no
> agradable.
> 
> Saludos!

Deja esas cosas para cuando hagas la tesis, en las materias te corre mucho el 
tiempo y no puedes experimentar (ademas de exasperar al docente).



-- 
Octavio A. Bruzzone
Investigador Adjunto (CONICET).
Ecología Numérica y Dinámica de Poblaciones de Insectos,
Laboratorio de Ecologia de Insectos.
INTA EEA Bariloche.
__________________________________________________
Correo Yahoo!
Espacio para todos tus mensajes, antivirus y antispam ¡gratis! 
¡Abrí tu cuenta ya! - http://correo.yahoo.com.ar