How to compute gradients or derivatives.

Dear all, I was lookogn for u numerical extension to python and was therefore glad that i ran into numpy. I did all the calculations on mataln and I want to migrate to python. Numpy and numarray pretend to have the basic functionality of matlab and therefore i expected als a function that is able to compute gradients of high-order arrays which is a basic matlab thing. I could not find something like that. I found something in SciPy but I think that this one cannot compute high-order arrays. I'm thinkin of something as in matlab e.g. say that T=f(x,y,z,p) [gx,gy] = gradientT; Is something like that possible? Thanks in advance

take a look to some parts of the code in http://www.pyclimate.org. There are some functions to compute gradients, divergence and vertical component of the rotational of scalar and vector fields. BUT... we only support regular grids in a certain ordering, which is the one most used for climate data sets. Moreover, our derivatives are computed in spherical coordinates, like the ones used in geophysical data sets. Still, you can use the code as a starting point to develop your own functions. Jon Saenz. | Tfno: +34 946012445 Depto. Fisica Aplicada II | Fax: +34 944648500 Facultad de Ciencias. \\ Universidad del Pais Vasco \\ Apdo. 644 \\ 48080 - Bilbao \\ SPAIN On Thu, 24 Jul 2003, Maurice van de Rijzen wrote:
Dear all,
I was lookogn for u numerical extension to python and was therefore glad that i ran into numpy. I did all the calculations on mataln and I want to migrate to python. Numpy and numarray pretend to have the basic functionality of matlab and therefore i expected als a function that is able to compute gradients of high-order arrays which is a basic matlab thing. I could not find something like that. I found something in SciPy but I think that this one cannot compute high-order arrays. I'm thinkin of something as in matlab e.g. say that T=f(x,y,z,p)
[gx,gy] = gradientT;
Is something like that possible?
Thanks in advance
------------------------------------------------------- This SF.Net email sponsored by: Free pre-built ASP.NET sites including Data Reports, E-commerce, Portals, and Forums are available now. Download today and enter to win an XBOX or Visual Studio .NET. http://aspnet.click-url.com/go/psa00100003ave/direct;at.aspnet_072303_01/01 _______________________________________________ Numpy-discussion mailing list Numpy-discussion@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/numpy-discussion
participants (2)
-
Jon Saenz
-
Maurice van de Rijzen