[pyar] Programación científica en python

Nico Chiaraviglio nico.chiaraviglio en gmail.com
Dom Ago 22 04:26:47 ART 2010


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.

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.

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!
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <http://listas.python.org.ar/pipermail/pyar/attachments/20100822/5c8490a6/attachment.html>


More information about the pyar mailing list