Mailman 3 March 2000 - NumPy-Discussion

Binary kit for windows
by Janne Sinkkonen March 31, 2000

March 31, 2000

Could somebody give me a pointer to a easily installable Windows binary kit for NumPy? Even a bit older version (such as 11) would do fine. I couldn't find any reference to such a kit from sourceforge or the LLNL site. -- Janne

1 0

single precision patch for arrayfnsmodule.c :interp()
by Joe Van Andel March 29, 2000

March 29, 2000

As I've mentioned earlier, I'm using interp() with very large datasets, and I can't afford to use double precision arrays. Here's a patch that lets interp() accept an optional typecode argument. Passing 'f' calls the new single precision version, and returns a single precision array. Passing no argument or 'd' uses the previous, double precision version. (No other array types are supported - an error is returned.) I hope this can be added to the CVS version of Numeric. Thanks! -- Joe VanAndel National Center for Atmospheric Research http://www.atd.ucar.edu/~vanandel/ Internet: vanandel(a)ucar.edu *** arrayfnsmodule.c 1999/04/14 22:58:32 1.1 --- arrayfnsmodule.c 1999/08/12 23:27:28 *************** *** 6,11 **** --- 6,13 ---- #include <stdio.h> #include <stdlib.h> + #define MAX_INTERP_DIMS 6 + static PyObject *ErrorObject; /* Define 2 macros for error handling: *************** *** 34,39 **** --- 36,43 ---- #define A_DIM(a,i) (((PyArrayObject *)a)->dimensions[i]) #define GET_ARR(ap,op,type,dim) \ Py_Try(ap=(PyArrayObject *)PyArray_ContiguousFromObject(op,type,dim,dim)) + #define GET_ARR2(ap,op,type,min,max) \ + Py_Try(ap=(PyArrayObject *)PyArray_ContiguousFromObject(op,type,min,max)) #define ERRSS(s) ((PyObject *)(PyErr_SetString(ErrorObject,s),0)) #define SETERR(s) if(!PyErr_Occurred()) ERRSS(errstr ? errstr : s) #define DECREF_AND_ZERO(p) do{Py_XDECREF(p);p=0;}while(0) *************** *** 571,576 **** --- 575,613 ---- return result ; } + static int + binary_searchf(float dval, float dlist [], int len) + { + /* binary_search accepts three arguments: a numeric value and + * a numeric array and its length. It assumes that the array is sorted in + * increasing order. It returns the index of the array's + * largest element which is <= the value. It will return -1 if + * the value is less than the least element of the array. */ + /* self is not used */ + int bottom , top , middle, result; + + if (dval < dlist [0]) + result = -1 ; + else { + bottom = 0; + top = len - 1; + while (bottom < top) { + middle = (top + bottom) / 2 ; + if (dlist [middle] < dval) + bottom = middle + 1 ; + else if (dlist [middle] > dval) + top = middle - 1 ; + else + return middle ; + } + if (dlist [bottom] > dval) + result = bottom - 1 ; + else + result = bottom ; + } + + return result ; + } static char arr_interp__doc__[] = "" ; *************** *** 597,609 **** Py_DECREF(ay); Py_DECREF(ax); return NULL ;} ! GET_ARR(az,oz,PyArray_DOUBLE,1); lenz = A_SIZE (az); dy = (double *) A_DATA (ay); dx = (double *) A_DATA (ax); dz = (double *) A_DATA (az); ! Py_Try (_interp = (PyArrayObject *) PyArray_FromDims (1, &lenz, ! PyArray_DOUBLE)); dres = (double *) A_DATA (_interp) ; slopes = (double *) malloc ( (leny - 1) * sizeof (double)) ; for (i = 0 ; i < leny - 1; i++) { --- 634,647 ---- Py_DECREF(ay); Py_DECREF(ax); return NULL ;} ! GET_ARR2(az,oz,PyArray_DOUBLE,1,MAX_INTERP_DIMS); lenz = A_SIZE (az); dy = (double *) A_DATA (ay); dx = (double *) A_DATA (ax); dz = (double *) A_DATA (az); ! /* create output array with same size as 'Z' input array */ ! Py_Try (_interp = (PyArrayObject *) PyArray_FromDims ! (A_NDIM(az), az->dimensions, PyArray_DOUBLE)); dres = (double *) A_DATA (_interp) ; slopes = (double *) malloc ( (leny - 1) * sizeof (double)) ; for (i = 0 ; i < leny - 1; i++) { *************** *** 627,632 **** --- 665,724 ---- return PyArray_Return (_interp); } + /* return float, rather than double */ + + static PyObject * + arr_interpf(PyObject *self, PyObject *args) + { + /* interp (y, x, z) treats (x, y) as a piecewise linear function + * whose value is y [0] for x < x [0] and y [len (y) -1] for x > + * x [len (y) -1]. An array of floats the same length as z is + * returned, whose values are ordinates for the corresponding z + * abscissae interpolated into the piecewise linear function. */ + /* self is not used */ + PyObject * oy, * ox, * oz ; + PyArrayObject * ay, * ax, * az , * _interp; + float * dy, * dx, * dz , * dres, * slopes; + int leny, lenz, i, left ; + + Py_Try(PyArg_ParseTuple(args, "OOO", &oy, &ox, &oz)); + GET_ARR(ay,oy,PyArray_FLOAT,1); + GET_ARR(ax,ox,PyArray_FLOAT,1); + if ( (leny = A_SIZE (ay)) != A_SIZE (ax)) { + SETERR ("interp: x and y are not the same length."); + Py_DECREF(ay); + Py_DECREF(ax); + return NULL ;} + GET_ARR2(az,oz,PyArray_FLOAT,1,MAX_INTERP_DIMS); + lenz = A_SIZE (az); + dy = (float *) A_DATA (ay); + dx = (float *) A_DATA (ax); + dz = (float *) A_DATA (az); + /* create output array with same size as 'Z' input array */ + Py_Try (_interp = (PyArrayObject *) PyArray_FromDims + (A_NDIM(az), az->dimensions, PyArray_FLOAT)); + dres = (float *) A_DATA (_interp) ; + slopes = (float *) malloc ( (leny - 1) * sizeof (float)) ; + for (i = 0 ; i < leny - 1; i++) { + slopes [i] = (dy [i + 1] - dy [i]) / (dx [i + 1] - dx [i]) ; + } + for ( i = 0 ; i < lenz ; i ++ ) + { + left = binary_searchf (dz [i], dx, leny) ; + if (left < 0) + dres [i] = dy [0] ; + else if (left >= leny - 1) + dres [i] = dy [leny - 1] ; + else + dres [i] = slopes [left] * (dz [i] - dx [left]) + dy [left]; + } + + free (slopes); + Py_DECREF(ay); + Py_DECREF(ax); + Py_DECREF(az); + return PyArray_Return (_interp); + } static int incr_slot_ (float x, double *bins, int lbins) { int i ; *************** *** 1295,1300 **** --- 1387,1393 ---- {"array_set", arr_array_set, 1, arr_array_set__doc__}, {"index_sort", arr_index_sort, 1, arr_index_sort__doc__}, {"interp", arr_interp, 1, arr_interp__doc__}, + {"interpf", arr_interpf, 1, arr_interp__doc__}, {"digitize", arr_digitize, 1, arr_digitize__doc__}, {"zmin_zmax", arr_zmin_zmax, 1, arr_zmin_zmax__doc__}, {"reverse", arr_reverse, 1, arr_reverse__doc__},

2 2

del( distutils )
by Les Schaffer March 29, 2000

March 29, 2000

I am just trying to compile the latest CVS updates on NumPy, and am getting very cranky about the dependence of build process on distutils. let me say it politely, using distutils at this point is for the birds. maybe for distributing released versions to the general public, where the distutils setup.py script has been tested on many platforms, is a long term solution. but for right now, i miss the Makefile stuff. an example: (gustav)~/system/numpy/Numerical/: python setup.py build running build [snip] /usr/bin/egcc -c -I/usr/include/python1.5 -IInclude -O3 -mpentium -fpic Src/_numpymodule.c Src/arrayobject.c Src/ufuncobject.c Src/ufuncobject.c:783: conflicting types for `PyUFunc_FromFuncAndData' /usr/include/python1.5/ufuncobject.h:100: previous declaration of `PyUFunc_FromFuncAndData' Src/ufuncobject.c: In function `PyUFunc_FromFuncAndData': Src/ufuncobject.c:805: structure has no member named `doc' Src/ufuncobject.c: In function `ufunc_getattr': Src/ufuncobject.c:955: structure has no member named `doc' [blah blah blah] if this was from running make, i could be inside xemacs and click on the damn error and go right to the line in question. now, i gotta horse around, its a waste of time. another example, i posed to c.l.p a week or so ago. where it turned out that distutils doesnt know enough that, for example, arrayobject.h in distribution is newer than the one in /usr/include/python1.5 so breaks the build when API changes have been made. distutils should at leaast get the order of the -I's correct. but we're not a distutils SIG, we're NumPy, right? can we please go back, at least for now, to using -- or at a minimum distributing -- the Makefile's? les schaffer

2 3

Release of Pyclimate0.0
by Jon Saenz March 28, 2000

March 28, 2000

<P><A HREF="http://lcdx00.wm.lc.ehu.es/~jsaenz/pyclimate">Pyclimate 0.0</A> - Climate variability analysis using Numeric Python (28-Mar-00) Tuesday, 03/28/2000 Hello, all. We are making the first announce of a pre-alpha release (version 0.0) of our package pyclimate, which presents some tools used for climate variability analysis and which make extensive use of Numerical Python. It is released under the GNU Public License. We call them a pre-alpha release. Even though the routines are quite debugged, they are yet growing and we are thinking in making a stable release shortly after receiving some feedback from users. The package contains: IO functions ------------ -ASCII files (simple, but useful) -ncstruct.py: netCDF structure copier. From a COARDS compliant netCDF file, this module creates a COARDS compliant file, copying the needed attributes, dimensions, auxiliary variables, comments, and so on in one call. Time handling routines ---------------------- * JDTime.py -> Some C/Python functions to convert from date to Scaliger's Julian Day and from Julian Day to date. We are not trying to replace mxDate, but addressing a different problem. In particular, this module contains a routine especially suited to handling monthly time steps for climatological use. * JDTimeHandler.py -> Python module which parses the units attribute of the time variable in a COARDS file and which offsets and scales adequately the time values to read/save date fields. Interface to DCDFLIB.C ---------------------- A C/Python interface to the free DCDFLIB.C library is provided. This library allows direct and inverse computations of parameters for several probability distribution functions like Chi^2, normal, binomial, F, noncentral F, and many many more. EOF analysis ------------ Empirical Orthogonal Function analysis based on the SVD decomposition of the data matrix and related functions to test the reliability/degeneracy of eigenvalues (truncation rules). Monte Carlo test of the stability of eigenvectors to temporal subsampling. SVD decomposition ----------------- SVD decomposition of the correlation matrix of two datasets, functions to compute the expansion coefficients, the squared cumulative covariance fraction and the homogeneous and heterogeneous correlation maps. Monte Carlo test of the stability of singular vectors to temporal subsampling. Multivariate digital filter --------------------------- Multivariate digital filter (high and low pass) based on the Kolmogorov-Zurbenko filter Differential operators on the sphere ------------------------------------ Some classes to compute differential operators (gradient and divergence) on a regular latitude/longitude grid. PREREQUISITES ============= To be able to use it, you will need: 1. Python ;-) 2. netCDF library 3.4 or later 3. Scientific Python, by Konrad Hinsen 4. DCDFLIB.C version 1.1 IF AND ONLY IF you really want to change the C code (JDTime.[hc] and pycdf.[hc]), then, you will also need SWIG. COMPILATION =========== There is no a automatic compilation/installation procedure, but the Makefile is quite straightforward. After manually editing the Makefile for different platforms, the commands make make test -> Runs a (not infalible) regression test make install will do it. SORRY, we don't use it under Windows, only UNIX. Volunteers that generate a Windows installation file would be appreciated, but we will not do it. DOCUMENTATION ============= LaTeX, Postscript and PDF versions of the manual are included in the distribution. However, we are preparing a new set of documentation according to PSA rules. AVAILABILITY ============ http://lcdx00.wm.lc.ehu.es/~jsaenz/pyclimate (Europe) http://pyclimate.zubi.net/ (USA) http://starship.python.net/crew/~jsaenz (USA) Any feedback from the users of the package will be really appreciated by the authors. We will try to incorporate new developments, in case we are able to do so. Our time availability is scarce. Enjoy. Jon Saenz, jsaenz(a)wm.lc.ehu.es Juan Zubillaga, wmpzuesj(a)lg.ehu.es

1 0

Recent CVS NumPy Version with recent CVS distutil and Linux
by bhoel＠server.python.net March 25, 2000

March 25, 2000

Hello, I tried to install a NumPy version downloaded this morning using distutils also downloaded then. This failed on my Linux box. First the setup.py script in not compatible with the newer distutil API. Second the glibc math.h includes a mathtypes.h which defines a function named gamma. This leads to a problem with ranlibmodule.c. I tried to solve both problems. This is documented with the applied patch. Cheers Berthold -- bhoel(a)starship.python.net / http://starship.python.net/crew/bhoel/ It is unlawful to use this email address for unsolicited ads (USC Title 47 Sec.227). I will assess a US$500 charge for reviewing and deleting each unsolicited ad.

1 0

Doc strings for ufuncs
by Charles G Waldman March 25, 2000

March 25, 2000

I really, really, really like the Cephes module and a lot of the work Travis Oliphant has been doing on Numeric Python. Nice to see Numeric getting more and more powerful all the time. However, as more and more functions get added to the libraries, their names become more and more obscure, and it sure would be nice to have doc-strings for them. For the stock Numeric ufuncs like "add", the meanings are self-evident, but things like "pdtri" and "incbet" are a bit less obvious. (One of my pet peeves with Python is that there are so many functions and classes with empty doc-strings. Bit by bit, I'm trying to add them in). However, one little problem with the current Numeric implementation is that ufunc objects don't have support for a doc string, the doc string is hardwired in the type object, and comes up as: "Optimizied FUNCtions make it possible to implement arithmetic with matrices efficiently" This is not only non-helpful, it's misspelled ("Optimizied"?) Well, according to the charter, a well-behaved Nummie will "Fix bugs at will, without prior consultation." But a well-behaved Nummie wil also "Add new features only after consultation with other Nummies" So, I hereby declare doing something about the useless doc-strings to be fixing a bug and not adding a feature ;-) The patch below adds doc strings to all the ufuncs in the Cephes module. They are automatically extracted from the HTML documentation for the Cephes module. (In the process of doing this, I also added some missing items to said HTML documentation). This patch depends on another patch, which I am submitting via the SourceForge, which allows ufunc objects to have doc strings. With these patches, you get this: >>> import cephes >>> print cephes.incbet.__doc__ incbet(a,b,x) returns the incomplete beta integral of the arguments, evaluated from zero to x: gamma(a+b) / (gamma(a)*gamma(b)) * integral(t**(a-1) (1-t)**(b-1), t=0..x). instead of this: >>> print cephes.incbet.__doc__ Optimizied FUNCtions make it possible to implement arithmetic with matrices efficiently Isn't that nicer? "Ni-Ni-Numpy!" Here's the "gendoc.py" script to pull the docstrings out of included_functions.h: #!/usr/bin/env python """generate cephes_doc.h from included_functions.html""" import string def parse(infile): d={} key=None val='' for line in infile.readlines(): if not string.strip(line): continue if line[0]=='<': if key and val: d[key]=string.strip(val) key,val=None,None if line[:4]=='<DT>': tok=string.split(line) tok=string.split(tok[-1],'(') key=tok[0] elif line[:4]=='<DD>' and key: val=line[4:] else: if val: val=val+line return d if __name__=="__main__": d = parse(open("docs/included_functions.html",'r')) keys = d.keys() keys.sort() ofile=open("cephes_doc.h",'w') for key in keys: ofile.write('#define %s_doc "%s"\n'%(key,repr(d[key])[1:-1])) ofile.close() And here are the patches to cephes: --- cephes-1.2/Makefile Tue Dec 28 15:53:11 1999 +++ cephes-1.2.cgw/Makefile Thu Mar 23 14:40:02 2000 @@ -15,8 +15,11 @@ $(FC) -shared $(FFLAGS) -o cephesmodule.so $(OBJS) \ -L./cephes -lcephes -L./amos -lamos -L./toms -ltoms -cephesmodule.o: cephesmodule.c +cephesmodule.o: cephesmodule.c cephes_doc.h $(CC) -c $(CFLAGS) -fPIC cephesmodule.c + +cephes_doc.h: docs/included_functions.html + python gendoc.py amos/libamos.a: cd amos && make CC=$(CC) CFLAGS="$(CFLAGS)" FC=$(FC) FFLAGS="$(FFLAGS)" --- cephes-1.2/docs/included_functions.html Sat Mar 27 00:27:18 1999 +++ cephes-1.2.cgw/docs/included_functions.html Thu Mar 23 14:39:03 2000 @@ -7,7 +7,7 @@ <DT>(Ai,Aip,Bi,Bip) = airy(z) <DD>airy(z) calculates the Airy functions and their derivatives evaluated at real or complex number z. The Airy functions Ai and Bi -are two independent solutions of y"(x)=xy. Aip and Bip are the first derivatives +are two independent solutions of y''(x)=xy. Aip and Bip are the first derivatives evaluated at x of Ai and Bi respectively. <DT>(Aie,Aipe,Bie,Bipe) = airye(z) @@ -19,7 +19,7 @@ <h1> Elliptic Functions and Integrals </h1> <DL> -<DT>(sn,cn,dn,ph) = ellpj_sn(u,m) +<DT>(sn,cn,dn,ph) = ellpj(u,m) <DD>ellpj calculates the Jacobian elliptic functions of parameter m between 0 and 1, and real u. The returned functions are often written sn(u|m), cn(u|m), and dn(u|m). The value of ph is such @@ -430,11 +430,12 @@ <DT>y = cosm1(x) <DD>calculates cos(x) - 1 for use when x is near zero. -</DL> - - - +<DT>y = round(x) +<DD>Returns the nearest integer to x as a double precision +floating point result. If x ends in 0.5 exactly, the +nearest even integer is chosen. +</DL> --- cephes-1.2/cephesmodule.c Tue Dec 28 15:32:26 1999 +++ cephes-1.2.cgw/cephesmodule.c Thu Mar 23 14:24:37 2000 @@ -23,6 +23,8 @@ #include <math.h> #endif +#include "cephes_doc.h" + static PyUFuncGenericFunction cephes1_functions[] = { NULL, NULL, }; static PyUFuncGenericFunction cephes1_2_functions[] = { NULL, NULL, }; static PyUFuncGenericFunction cephes1c_4_functions[] = { NULL, NULL, NULL, NULL }; @@ -220,335 +222,435 @@ /* Create function objects for each function call and insert them in the dictionary */ - f = PyUFunc_FromFuncAndData(cephes3a_functions, bdtrc_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "bdtrc", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, bdtrc_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "bdtrc", bdtrc_doc, 1); PyDict_SetItemString(dictionary, "bdtrc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, bdtr_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "bdtr", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, bdtr_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "bdtr", bdtr_doc, 1); PyDict_SetItemString(dictionary, "bdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, bdtri_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "bdtri", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, bdtri_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "bdtri", bdtri_doc, 1); PyDict_SetItemString(dictionary, "bdtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, btdtr_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "btdtr", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, btdtr_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "btdtr", btdtr_doc, 1); PyDict_SetItemString(dictionary, "btdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, fdtrc_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "fdtrc", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, fdtrc_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "fdtrc", fdtrc_doc, 1); PyDict_SetItemString(dictionary, "fdtrc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, fdtr_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "fdtr", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, fdtr_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "fdtr", fdtr_doc, 1); PyDict_SetItemString(dictionary, "fdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, fdtri_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "fdtri", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, fdtri_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "fdtri", fdtri_doc, 1); PyDict_SetItemString(dictionary, "fdtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, gdtrc_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "gdtrc", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, gdtrc_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "gdtrc", gdtrc_doc, 1); PyDict_SetItemString(dictionary, "gdtrc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, gdtr_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "gdtr", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, gdtr_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "gdtr", gdtr_doc, 1); PyDict_SetItemString(dictionary, "gdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes4_functions, hyp2f1_data, cephes_5_types, 2, 4, 1, PyUFunc_None, "hyp2f1", 1); + f = PyUFunc_FromFuncAndData(cephes4_functions, hyp2f1_data, cephes_5_types, + 2, 4, 1, PyUFunc_None, "hyp2f1", hyp2f1_doc, 1); PyDict_SetItemString(dictionary, "hyp2f1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, hyperg_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "hyperg", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, hyperg_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "hyperg", hyperg_doc, 1); PyDict_SetItemString(dictionary, "hyperg", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes4a_2_functions, hyp2f0_data, cephes_6_types, 2, 4, 2, PyUFunc_None, "hyp2f0", 1); + f = PyUFunc_FromFuncAndData(cephes4a_2_functions, hyp2f0_data, cephes_6_types, + 2, 4, 2, PyUFunc_None, "hyp2f0", hyp2f0_doc, 1); PyDict_SetItemString(dictionary, "hyp2f0", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes4_2_functions, onef2_data, cephes_6_types, 2, 4, 2, PyUFunc_None, "onef2", 1); + f = PyUFunc_FromFuncAndData(cephes4_2_functions, onef2_data, cephes_6_types, + 2, 4, 2, PyUFunc_None, "onef2", onef2_doc, 1); PyDict_SetItemString(dictionary, "onef2", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes4_2_functions, threef0_data, cephes_6_types, 2, 4, 2, PyUFunc_None, "threef0", 1); + f = PyUFunc_FromFuncAndData(cephes4_2_functions, threef0_data, cephes_6_types, + 2, 4, 2, PyUFunc_None, "threef0", threef0_doc, 1); PyDict_SetItemString(dictionary, "threef0", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, incbet_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "incbet", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, incbet_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "incbet", incbet_doc, 1); PyDict_SetItemString(dictionary, "incbet", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, incbi_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "incbi", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, incbi_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "incbi", incbi_doc, 1); PyDict_SetItemString(dictionary, "incbi", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, nbdtrc_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "nbdtrc", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, nbdtrc_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "nbdtrc", nbdtrc_doc, 1); PyDict_SetItemString(dictionary, "nbdtrc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, nbdtr_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "nbdtr", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, nbdtr_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "nbdtr", nbdtr_doc, 1); PyDict_SetItemString(dictionary, "nbdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3a_functions, nbdtri_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "nbdtri", 1); + f = PyUFunc_FromFuncAndData(cephes3a_functions, nbdtri_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "nbdtri", nbdtri_doc, 1); PyDict_SetItemString(dictionary, "nbdtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, beta_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "beta", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, beta_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "beta", beta_doc, 1); PyDict_SetItemString(dictionary, "beta", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, lbeta_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "lbeta", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, lbeta_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "lbeta", lbeta_doc, 1); PyDict_SetItemString(dictionary, "lbeta", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, cbrt_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "cbrt", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, cbrt_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "cbrt", cbrt_doc, 1); PyDict_SetItemString(dictionary, "cbrt", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, chdtrc_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "chdtrc", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, chdtrc_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "chdtrc", chdtrc_doc, 1); PyDict_SetItemString(dictionary, "chdtrc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, chdtr_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "chdtr", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, chdtr_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "chdtr", chdtr_doc, 1); PyDict_SetItemString(dictionary, "chdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, chdtri_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "chdtri", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, chdtri_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "chdtri", chdtri_doc, 1); PyDict_SetItemString(dictionary, "chdtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, dawsn_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "dawsn", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, dawsn_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "dawsn", dawsn_doc, 1); PyDict_SetItemString(dictionary, "dawsn", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, ellie_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "ellie", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, ellie_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "ellie", ellie_doc, 1); PyDict_SetItemString(dictionary, "ellie", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, ellik_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "ellik", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, ellik_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "ellik", ellik_doc, 1); PyDict_SetItemString(dictionary, "ellik", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, ellpe_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "ellpe", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, ellpe_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "ellpe", ellpe_doc, 1); PyDict_SetItemString(dictionary, "ellpe", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, ellpk_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "ellpk", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, ellpk_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "ellpk", ellpk_doc, 1); PyDict_SetItemString(dictionary, "ellpk", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, exp10_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "exp10", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, exp10_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "exp10", exp10_doc, 1); PyDict_SetItemString(dictionary, "exp10", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, exp2_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "exp2", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, exp2_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "exp2", exp2_doc, 1); PyDict_SetItemString(dictionary, "exp2", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, Gamma_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "gamma", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, Gamma_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "gamma", gamma_doc, 1); PyDict_SetItemString(dictionary, "gamma", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, lgam_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "lgam", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, lgam_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "lgam", lgam_doc, 1); PyDict_SetItemString(dictionary, "lgam", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, i0_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "i0", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, i0_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "i0", i0_doc, 1); PyDict_SetItemString(dictionary, "i0", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, i0e_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "i0e", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, i0e_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "i0e", i0e_doc, 1); PyDict_SetItemString(dictionary, "i0e", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, i1_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "i1", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, i1_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "i1", i1_doc, 1); PyDict_SetItemString(dictionary, "i1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, i1e_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "i1e", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, i1e_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "i1e", i1e_doc, 1); PyDict_SetItemString(dictionary, "i1e", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, igamc_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "igamc", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, igamc_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "igamc", igamc_doc, 1); PyDict_SetItemString(dictionary, "igamc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, igam_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "igam", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, igam_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "igam", igam_doc, 1); PyDict_SetItemString(dictionary, "igam", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, igami_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "igami", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, igami_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "igami", igami_doc, 1); PyDict_SetItemString(dictionary, "igami", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2c_functions, iv_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "iv", 1); + f = PyUFunc_FromFuncAndData(cephes2c_functions, iv_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "iv", iv_doc, 1); PyDict_SetItemString(dictionary, "iv", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cp_functions, ive_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "ive", 1); + f = PyUFunc_FromFuncAndData(cephes2cp_functions, ive_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "ive", ive_doc, 1); PyDict_SetItemString(dictionary, "ive", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_4_functions, ellpj_data, cephes_6_types, 2, 2, 4, PyUFunc_None, "ellpj", 1); + f = PyUFunc_FromFuncAndData(cephes2_4_functions, ellpj_data, cephes_6_types, + 2, 2, 4, PyUFunc_None, "ellpj", ellpj_doc, 1); PyDict_SetItemString(dictionary, "ellpj", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, expn_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "expn", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, expn_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "expn", expn_doc, 1); PyDict_SetItemString(dictionary, "expn", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, jn_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "jn", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, jn_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "jn", jn_doc, 1); PyDict_SetItemString(dictionary, "jn", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, kn_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "kn", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, kn_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "kn", kn_doc, 1); PyDict_SetItemString(dictionary, "kn", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, pdtrc_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "pdtrc", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, pdtrc_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "pdtrc", pdtrc_doc, 1); PyDict_SetItemString(dictionary, "pdtrc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, pdtr_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "pdtr", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, pdtr_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "pdtr", pdtr_doc, 1); PyDict_SetItemString(dictionary, "pdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, pdtri_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "pdtri", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, pdtri_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "pdtri", pdtri_doc, 1); PyDict_SetItemString(dictionary, "pdtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, stdtr_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "stdtr", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, stdtr_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "stdtr", stdtr_doc, 1); PyDict_SetItemString(dictionary, "stdtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, stdtri_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "stdtri", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, stdtri_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "stdtri", stdtri_doc, 1); PyDict_SetItemString(dictionary, "stdtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, yn_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "yn", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, yn_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "yn", yn_doc, 1); PyDict_SetItemString(dictionary, "yn", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, smirnov_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "smirnov", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, smirnov_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "smirnov", smirnov_doc, 1); PyDict_SetItemString(dictionary, "smirnov", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2a_functions, smirnovi_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "smirnovi", 1); + f = PyUFunc_FromFuncAndData(cephes2a_functions, smirnovi_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "smirnovi", smirnovi_doc, 1); PyDict_SetItemString(dictionary, "smirnovi", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1c_4_functions, airy_data, cephes_5c_types, 4, 1, 4, PyUFunc_None, "airy", 1); + f = PyUFunc_FromFuncAndData(cephes1c_4_functions, airy_data, cephes_5c_types, + 4, 1, 4, PyUFunc_None, "airy", airy_doc, 1); PyDict_SetItemString(dictionary, "airy", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1cp_4_functions, airye_data, cephes_5c_types, 4, 1, 4, PyUFunc_None, "airye", 1); + f = PyUFunc_FromFuncAndData(cephes1cp_4_functions, airye_data, cephes_5c_types, + 4, 1, 4, PyUFunc_None, "airye", airye_doc, 1); PyDict_SetItemString(dictionary, "airye", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_2_functions, fresnl_data, cephes_3_types, 2, 1, 2, PyUFunc_None, "fresnl", 1); + f = PyUFunc_FromFuncAndData(cephes1_2_functions, fresnl_data, cephes_3_types, + 2, 1, 2, PyUFunc_None, "fresnl", fresnl_doc, 1); PyDict_SetItemString(dictionary, "fresnl", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_2_functions, shichi_data, cephes_3_types, 2, 1, 2, PyUFunc_None, "shichi", 1); + f = PyUFunc_FromFuncAndData(cephes1_2_functions, shichi_data, cephes_3_types, + 2, 1, 2, PyUFunc_None, "shichi", shichi_doc, 1); PyDict_SetItemString(dictionary, "shichi", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_2_functions, sici_data, cephes_3_types, 2, 1, 2, PyUFunc_None, "sici", 1); + f = PyUFunc_FromFuncAndData(cephes1_2_functions, sici_data, cephes_3_types, + 2, 1, 2, PyUFunc_None, "sici", sici_doc, 1); PyDict_SetItemString(dictionary, "sici", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, j0_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "j0", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, j0_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "j0", j0_doc, 1); PyDict_SetItemString(dictionary, "j0", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, y0_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "y0", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, y0_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "y0", y0_doc, 1); PyDict_SetItemString(dictionary, "y0", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, j1_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "j1", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, j1_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "j1", j1_doc, 1); PyDict_SetItemString(dictionary, "j1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, y1_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "y1", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, y1_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "y1", y1_doc, 1); PyDict_SetItemString(dictionary, "y1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2c_functions, jv_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "jv", 1); + f = PyUFunc_FromFuncAndData(cephes2c_functions, jv_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "jv", jv_doc, 1); PyDict_SetItemString(dictionary, "jv", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cp_functions, jve_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "jve", 1); + f = PyUFunc_FromFuncAndData(cephes2cp_functions, jve_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "jve", jve_doc, 1); PyDict_SetItemString(dictionary, "jve", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2c_functions, yv_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "yv", 1); + f = PyUFunc_FromFuncAndData(cephes2c_functions, yv_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "yv", yv_doc, 1); PyDict_SetItemString(dictionary, "yv", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cp_functions, yve_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "yve", 1); + f = PyUFunc_FromFuncAndData(cephes2cp_functions, yve_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "yve", yve_doc, 1); PyDict_SetItemString(dictionary, "yve", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, k0_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "k0", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, k0_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "k0", k0_doc, 1); PyDict_SetItemString(dictionary, "k0", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, k0e_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "k0e", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, k0e_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "k0e", k0e_doc, 1); PyDict_SetItemString(dictionary, "k0e", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, k1_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "k1", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, k1_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "k1", k1_doc, 1); PyDict_SetItemString(dictionary, "k1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, k1e_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "k1e", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, k1e_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "k1e", k1e_doc, 1); PyDict_SetItemString(dictionary, "k1e", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cp_functions, kv_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "kv", 1); + f = PyUFunc_FromFuncAndData(cephes2cp_functions, kv_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "kv", kv_doc, 1); PyDict_SetItemString(dictionary, "kv", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cp_functions, kve_data, cephes_3c_types, 4, 2, 1, PyUFunc_None, "kve", 1); + f = PyUFunc_FromFuncAndData(cephes2cp_functions, kve_data, cephes_3c_types, + 4, 2, 1, PyUFunc_None, "kve", kve_doc, 1); PyDict_SetItemString(dictionary, "kve", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel1_data, cephes_3cp_types, 2, 2, 1, PyUFunc_None, "hankel1", 1); + f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel1_data, cephes_3cp_types, + 2, 2, 1, PyUFunc_None, "hankel1", hankel1_doc, 1); PyDict_SetItemString(dictionary, "hankel1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel1e_data, cephes_3cp_types, 2, 2, 1, PyUFunc_None, "hankel1e", 1); + f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel1e_data, cephes_3cp_types, + 2, 2, 1, PyUFunc_None, "hankel1e", hankel1e_doc, 1); PyDict_SetItemString(dictionary, "hankel1e", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel2_data, cephes_3cp_types, 2, 2, 1, PyUFunc_None, "hankel2", 1); + f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel2_data, cephes_3cp_types, + 2, 2, 1, PyUFunc_None, "hankel2", hankel2_doc, 1); PyDict_SetItemString(dictionary, "hankel2", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel2e_data, cephes_3cp_types, 2, 2, 1, PyUFunc_None, "hankel2e", 1); + f = PyUFunc_FromFuncAndData(cephes2cpp_functions, hankel2e_data, cephes_3cp_types, + 2, 2, 1, PyUFunc_None, "hankel2e", hankel2e_doc, 1); PyDict_SetItemString(dictionary, "hankel2e", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, ndtr_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "ndtr", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, ndtr_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "ndtr", ndtr_doc, 1); PyDict_SetItemString(dictionary, "ndtr", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, erfc_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "erfc", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, erfc_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "erfc", erfc_doc, 1); PyDict_SetItemString(dictionary, "erfc", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, erf_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "erf", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, erf_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "erf", erf_doc, 1); PyDict_SetItemString(dictionary, "erf", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, ndtri_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "ndtri", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, ndtri_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "ndtri", ndtri_doc, 1); PyDict_SetItemString(dictionary, "ndtri", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, psi_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "psi", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, psi_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "psi", psi_doc, 1); PyDict_SetItemString(dictionary, "psi", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, rgamma_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "rgamma", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, rgamma_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "rgamma", rgamma_doc, 1); PyDict_SetItemString(dictionary, "rgamma", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, round_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "round", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, round_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "round", round_doc, 1); PyDict_SetItemString(dictionary, "round", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, sindg_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "sindg", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, sindg_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "sindg", sindg_doc, 1); PyDict_SetItemString(dictionary, "sindg", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, cosdg_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "cosdg", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, cosdg_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "cosdg", cosdg_doc, 1); PyDict_SetItemString(dictionary, "cosdg", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes3_functions, radian_data, cephes_4_types, 2, 3, 1, PyUFunc_None, "radian", 1); + f = PyUFunc_FromFuncAndData(cephes3_functions, radian_data, cephes_4_types, + 2, 3, 1, PyUFunc_None, "radian", radian_doc, 1); PyDict_SetItemString(dictionary, "radian", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, tandg_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "tandg", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, tandg_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "tandg", tandg_doc, 1); PyDict_SetItemString(dictionary, "tandg", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, cotdg_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "cotdg", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, cotdg_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "cotdg", cotdg_doc, 1); PyDict_SetItemString(dictionary, "cotdg", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, log1p_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "log1p", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, log1p_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "log1p", log1p_doc, 1); PyDict_SetItemString(dictionary, "log1p", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, expm1_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "expm1", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, expm1_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "expm1", expm1_doc, 1); PyDict_SetItemString(dictionary, "expm1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, cosm1_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "cosm1", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, cosm1_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "cosm1", cosm1_doc, 1); PyDict_SetItemString(dictionary, "cosm1", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, spence_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "spence", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, spence_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "spence", spence_doc, 1); PyDict_SetItemString(dictionary, "spence", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, zetac_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "zetac", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, zetac_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "zetac", zetac_doc, 1); PyDict_SetItemString(dictionary, "zetac", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, struve_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "struve", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, struve_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "struve", struve_doc, 1); PyDict_SetItemString(dictionary, "struve", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes2_functions, zeta_data, cephes_3_types, 2, 2, 1, PyUFunc_None, "zeta", 1); + f = PyUFunc_FromFuncAndData(cephes2_functions, zeta_data, cephes_3_types, + 2, 2, 1, PyUFunc_None, "zeta", zeta_doc, 1); PyDict_SetItemString(dictionary, "zeta", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, kolmogorov_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "kolmogorov", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, kolmogorov_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "kolmogorov", kolmogorov_doc, 1); PyDict_SetItemString(dictionary, "kolmogorov", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1_functions, kolmogi_data, cephes_2_types, 2, 1, 1, PyUFunc_None, "kolmogorovi", 1); + f = PyUFunc_FromFuncAndData(cephes1_functions, kolmogi_data, cephes_2_types, + 2, 1, 1, PyUFunc_None, "kolmogorovi", kolmogorovi_doc, 1); PyDict_SetItemString(dictionary, "kolmogorovi", f); Py_DECREF(f); - f = PyUFunc_FromFuncAndData(cephes1c_functions, wofz_data, cephes_1c_types, 2, 1, 1, PyUFunc_None, "wofz", 1); + f = PyUFunc_FromFuncAndData(cephes1c_functions, wofz_data, cephes_1c_types, + 2, 1, 1, PyUFunc_None, "wofz", wofz_doc, 1); PyDict_SetItemString(dictionary, "wofz", f); Py_DECREF(f); } NB: This stuff won't work at all without the corresponding patches to the Numeric core, to be posted shortly to the sourceforge site.

4 3

quadrature.py vs Multipack
by Jean-Bernard Addor March 23, 2000

March 23, 2000

Hey Numpy people! I have to integrate functions like: Int_0^1 (t*(t-1))**-(H/2) dt or Int_-1^1 1/abs(t)**(1-H) dt, with H around .3 I just tried quadrature on the 1st one: it needs very order of quadrature to be precise and is in that case slow. Would it work better with Multipack? (I have to upgrade to python 1.5.2 to try Multipack!) Thank you for your help. Jean-Bernard

2 2

how to compute a Gamma function with Numpy?
by Jean-Bernard Addor March 23, 2000

March 23, 2000

Hey Numpy people! I need to compute the Gamma function (the one related to factorial) with Numpy. It seems no to be included in the module, Am I right? Is any code for computing it available? Where could I find an algorythme to adapt ? Jean-Bernard

2 1

Documentation
by Paul Barrett March 14, 2000

March 14, 2000

I have a couple of question about the Numpy documentation: 1. Is there a recent version of the Numerical Python manual available anywhere? I can't find it at the SourceForge and I've tried xfiles.llnl.org, but can't get through. (But from what Paul Dubois has said recently about the LLNL site, I shouldn't expect to either.) 2. Have any changes been made to the documentation since about Q1 1999? I think my current version dates from about this period. -- Paul

5 4

single precision version of interp
by Joe Van Andel March 14, 2000

March 14, 2000

As I've previously mentioned to Paul, I need a single precision version of 'interp()', so I can use it on large single precision arrays, without returning a double precision array. In my own copy of NumPy, I've written such a routine, and added it to 'arrayfns.c '. Naturally, I want to see this functionality built into the official release, so I do not have to apply my own patches to new releases. Can we decide how such single precision needs are accomodated? Should there be a keyword argument on the 'interp()' call, that calls the single precision version? Or should the caller invoke 'interpf()', rather than 'interp()?' I don't much care what the solution looks like, as long as people agree that: 1) we need such functionality in NumPy. 2) we can establish a precedent on how single precision vs double precision methods are invoked. Please let me know your opinions on how this should be resolved. Thanks much! -- Joe VanAndel National Center for Atmospheric Research http://www.atd.ucar.edu/~vanandel/ Internet: vanandel(a)ucar.edu

2 1