; ٮG@cm@sdkZdkZdkZeid\ZZZeideeeZdZ dZ dZ dZ dZ dZdZhd d fd <d d fd <d d fd <d d fd<d?d@<dAdB<dCdD<dEdFduZ?dvZ@dwdxdydzgd{d|d}dzgd~dddzgddddzgddddzgddddgddddgddddzgddddzgddddzgddddzgg ZAeAiBddddzge?oeAiBddddzgne@o6eAiBddddzgeAiBddddzgndfdYZCdd{gdd~gddgddgddggZDdfdYZEeEddZFdfdYZGeGddZHeGddZIdwd{fd{d{fd~d~fddfddfddfddfddfddfg ZJdwdfd{dfd~dfddfddfddfddfddfddfg ZKdwdwfd{dwfd~dwfddwfddwfddwfddwfddwfddwfg ZLdwdwfd{d{fd~d~fddfddfddfddfgZMdwdfd{dfd~dfddfddfddfddfgZNeJiBddfeKiBddfeLiBddwfeMiBddfeNiBddfe?oceJiBddfeKiBddfeLiBddwfeMiBddfeNiBddfnddfddfgZOddfddfgZPddfddfgZQddwfddwfgZRdZSdZTdZUdZVeiWeXe>dd{e4ZYe>dde5ZZe>dd~e6Z[e>dde7Z\e>dde9Z]e>dde8Z^e>dde:Z_e>dde;Z`ddde0dfeJdzgddde0dfeJdzgddde0dfeJdzgdddeYdfd{d{fgeSgdddeZdfddfgeSgddde[dfd~d~fgeSgddde\dfddfgeSgddde^dfddfgeTgddde]dfddfgeTgddde_dfddfgeUgddde`dfddfgeUgddde0dfeOdzgddde1dfeMdzgddde0dfeOdzgddde3dfeMdzgddde-dfeOdzgddde2dfeNdzgddde0dfeOdzgddde1dfeMdzgddde.dfeOdzgddde.dfeJdzgddde.dfeJdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde.dfeKdzgddde/dfdddwfgdze'dggddde/dfdddwfgdze'dggddde.dfeKdzgddde0dfeLdzgddde0dfeLdzgddde0dfeLdzgddde0dfeLdzgddde0dfeLdzgddde0dfeLdzgddde.dfeLdzgddde.dfeLdzgddde.dfeLdzgddde0dfeLdzgddde0dfeMdzgddde0dfeMdzgddde0dfeMdzgddde0dfeMdzgddde0dfeMdzgddde0dfeMdzgddde.dzfeMdzgddde.dfeOdzgddde.dzfeMdzgddde.dfeOdzgddde.dfeJeVgd dde.d feJeVgddde.dfeKdzgd dde.d feKdzgd dde.d feKdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde.dfeQdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-dfePdzgddde-d fePdzgddde-d!fePdzgddde-d"fePdzgddde-d#fePdzgddde-d$fePdzgddde-d%fePdzgddde.d#feRdzgddde.d&feRdzgddde.d'feRdzgddde.d(feRdzgddde.d)feRdzgddde.d*feRdzgddde-d+fePdzgddde-d,fePdzgddde.d-fePdzgd dde.d.fePdzgddde-d/fePdzgd dde-d0fePdzgd dde-d1fePdzgghZad2fd3YZbd4fd5YZcd6ecfd7YZdd8ecfd9YZed:ecfd;YZfd<ecfd=YZgedZheeeZiefeZjege ed>Zkeld?Zmeld@ZndAZoeoZpeqdBjo eondS(CNstemplates /* W W AAA RRRR N N III N N GGG !!! ** W W A A R R NN N I NN N G G !!! ** W W W AAAAA RRRR N N N I N N N G ! ** W W W A A R R N NN I N NN G GG ** W W A A R R N N III N N GGG !!! ** ** WARNING: This file is program generated by codegenerator.py. ** ** DO NOT EDIT THIS FILE! Any changes made to this file will be lost! */ #include #include #include "libnumarray.h" static PyObject *_Error; s5 /****************** *******************/ s static PyObject *asPyValue(void *data) { return Py_BuildValue(, *(( *) data)); } SELF_CHECKED_CFUNC_DESCR(asPyValue, CFUNC_AS_PY_VALUE); s static PyObject *asPyValue(void *data) { return PyComplex_FromDoubles((( *) data)->r, (( *) data)->i); } SELF_CHECKED_CFUNC_DESCR(asPyValue, CFUNC_AS_PY_VALUE); s static int fromPyValue(PyObject *value, void *dataptr) { *data = ( *) dataptr; if (!PyNumber_Check(value)) return 0; else { if (PyLong_Check(value)) { *data = (PyLong_AsLongLong(value)); } else if (PyInt_Check(value)) { *data = (PyInt_AsLong(value)); } else if (PyFloat_Check(value)) { *data = (PyFloat_AsDouble(value)); } else if (PyComplex_Check(value)) { *data = (PyComplex_RealAsDouble(value)); } else { return 0; } if (PyErr_Occurred()) return 0; else return 1; } } SELF_CHECKED_CFUNC_DESCR(fromPyValue, CFUNC_FROM_PY_VALUE); sb static int fromPyValue(PyObject *value, void *dataptr) { *data = ( *) dataptr; if (!PyNumber_Check(value)) return 0; else { if (PyLong_Check(value)) { data->r = PyLong_AsLong(value); data->i = 0; } else if (PyInt_Check(value)) { data->r = PyInt_AsLong(value); data->i = 0; } else if (PyFloat_Check(value)) { data->r = PyFloat_AsDouble(value); data->i = 0; } else if (PyComplex_Check(value)) { data->r = PyComplex_RealAsDouble(value); data->i = PyComplex_ImagAsDouble(value); } else { return 0; } if (PyErr_Occurred()) return 0; else return 1; } } SELF_CHECKED_CFUNC_DESCR(fromPyValue, CFUNC_FROM_PY_VALUE); s static int as(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; *tin = ( *) buffers[0]; *tout = ( *) buffers[1]; for (i=0; i } return 0; } UFUNC_DESCR2(as, sizeof(), sizeof()); sComplexs#tout->r = tin->r; tout->i = tin->i;sdefaults*tout = (tin->r);stout->r = *tin; tout->i = 0;s*tout = (*tin);cCs#|ddgjodSndSdS(Ns Complex32s Complex64sComplexsdefault(st(st((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pys _complex_typescCstt|t|fSdS(N(s_numarray_to_numarrays _complex_typest1st2(st1st2((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pysnumarray_to_numarrayss #define swap(i, j) { temp=v[i]; v[i]=v[j]; v[j]=temp; } #define wswap(i, j) { temp=v[i]; v[i]=v[j]; v[j]=temp; wtemp=w[i]; w[i]=w[j]; w[j]=wtemp; } s> #define STDC_LT(a,b) ((a) < (b)) #define STDC_LE(a,b) ((a) <= (b)) #define STDC_EQ(a,b) ((a) == (b)) void sort0( *v, long left, long right) { temp; long i, last, lastl, lastr, diff; diff = right - left; if (diff <= 0) return; /* choose the pivot randomly. */ i = left + (int) (((right-left)*1.0*rand())/(RAND_MAX+1.0)); /* Split array into values < pivot, and values > pivot. */ swap(left, i); last = left; for(i=left+1; i<=right; i++) if ((v[i], v[left])) { ++ last; swap(last, i); } swap(left, last); lastl = last-1; lastr = last+1; /* Exclude values equal to pivot from recursion */ while((left < lastl) && (v[last],v[lastl])) --lastl; while((lastr < right) && (v[last],v[lastr])) ++lastr; sort0(v, left, lastl); sort0(v, lastr, right); } static int sort(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { sort0( ( *) buffers[0], 0, niter-1); return 0; } UFUNC_DESCR1(sort, sizeof()); void asort0( *v, long *w, long left, long right) { temp; long wtemp, i, last, lastl, lastr; if (left >= right) return; /* choose the pivot randomly. */ i = left + (int) (((right-left)*1.0*rand())/(RAND_MAX+1.0)); wswap(left, i); last = left; for(i=left+1; i<=right; i++) if ((v[i], v[left])) { ++ last; wswap(last, i); } wswap(left, last); lastl = last-1; lastr = last+1; /* Exclude values equal to pivot from recursion */ while((left < lastl) && (v[last],v[lastl])) --lastl; while((lastr < right) && (v[last],v[lastr])) ++lastr; asort0(v, w, left, lastl); asort0(v, w, lastr, right); } static int asort(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { asort0( ( *) buffers[0], (long *) buffers[1], 0, niter-1); return 0; } CFUNC_DESCR(asort, CFUNC_UFUNC, CHECK_ALIGN, 0, 2, sizeof(), sizeof(long), 0, 0, 0, 0); static long search(long na, *a, v) { long i = 0; *b; while(na > 10) { long mid = na>>1; if ((v, a[mid])) na = mid; else { i += mid + 1; a += mid + 1; na = mid - 1 + (na&1); } } b = a; while(((*b, v)) && na--) ++ b; return i+(b-a); } static int searchsorted(int niter, int ninargs, int noutargs, void **buffers, long *bsizes) { maybelong nbins; *bins; *values; long *indices; maybelong i; if (NA_checkOneCBuffer("searchsorted", 1, buffers[0], bsizes[0], sizeof(maybelong))) return -1; nbins = *(maybelong *) buffers[0]; if (NA_checkOneCBuffer("searchsorted", nbins, buffers[1], bsizes[1], sizeof())) return -1; bins = ( *) buffers[1]; if (NA_checkOneCBuffer("searchsorted", niter, buffers[2], bsizes[2], sizeof())) return -1; values = ( *) buffers[2]; if (NA_checkOneCBuffer("searchsorted", niter, buffers[3], bsizes[3], sizeof(long))) return -1; indices = (long *) buffers[3]; if (NA_checkIo("searchsorted", 3, 1, ninargs, noutargs)) return -1; for(i=0; i(nbins, bins, values[i]); return 0; } SELF_CHECKED_CFUNC_DESCR(searchsorted, CFUNC_UFUNC); static int fillarray (long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; start, delta; *tin; *tparams; Int8 itemsizes[2] = { sizeof(), sizeof() }; Int8 iters[2] = { 0, 2 }; if (NA_checkIo("fillarray", 1, 1, ninargs, noutargs) || NA_checkNCBuffers("fillarray", 2, niter, buffers, bsizes, itemsizes, iters)) return -1; tin = ( *) buffers[0]; tparams = ( *) buffers[1]; start = tparams[0]; delta = tparams[1]; for (i=0; i } return 0; } SELF_CHECKED_CFUNC_DESCR(fillarray, CFUNC_UFUNC); static int mxmultiply( long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { *a, *b; maybelong maxi, maxj, maxk; *r; maybelong i, j, k; if (NA_checkIo("mxmultiply", 3, 1, ninargs, noutargs)) return -1; if (NA_checkOneCBuffer("mxmultiply", 1, buffers[2], bsizes[2], sizeof(maybelong))) return -1; else { maxi = ((maybelong *)buffers[2])[0]; maxj = ((maybelong *)buffers[2])[1]; maxk = ((maybelong *)buffers[2])[2]; } if (NA_checkOneCBuffer("mxmultiply", maxi*maxk, buffers[0], bsizes[0], sizeof())) return -1; else a = ( *) buffers[0]; if (NA_checkOneCBuffer("mxmultiply", maxj*maxk, buffers[1], bsizes[1], sizeof())) return -1; else b = ( *) buffers[1]; if (NA_checkOneCBuffer("mxmultiply", maxi*maxj, buffers[3], bsizes[3], sizeof())) return -1; else r = ( *) buffers[3]; for(i=0; i s; *ap = a + i*maxk; *bp = b + j*maxk; for(k=0; k } *r = s; } return 0; } SELF_CHECKED_CFUNC_DESCR(mxmultiply, CFUNC_UFUNC); static int nonzeroCoords(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { *source; maybelong *sstride; long **output; long i, minbsize, next = 0; if (NA_checkIo("nonzeroCoords", 2, noutargs, ninargs, noutargs)) return -1; if (NA_checkOneCBuffer("nonzeroCoords", niter, buffers[0], bsizes[0], sizeof())) return -1; if (NA_checkOneCBuffer("nonzeroCoords", noutargs, buffers[1], bsizes[1], sizeof(maybelong))) return -1; /* Check alignment only. Storage optimization makes numarray shorter than niter likely, possibly even 0! */ for (i=0; i", 0, buffers[2+i], bsizes[2+i], sizeof(long))) return -1; source = ( *) buffers[0]; sstride = (maybelong *) buffers[1]; output = (long **) &buffers[2]; minbsize = bsizes[0]; for(i=0; i: bad stride[%%ld].\n", i); return -1; } else if (bsizes[i+2] < minbsize) minbsize = bsizes[i+2]; next = 0; for(i=0; i) { long j, index = i; if (next >= minbsize) { PyErr_Format(_Error, "nonzeroCoords: insufficient index space.\n"); return -1; } for(j=0; j, CFUNC_UFUNC); slessthansSTDC_LTs lessequalsSTDC_LEsequalssSTDC_EQsclearss = 0;sdotstepss += *ap * *bp;sfillsteps*tin = start; start += delta;snonzerossource[i] != 0sNUM_CLTsNUM_CLEsNUM_CEQss.r = 0; s.i=0;s6Complex64 t; NUM_CMUL(*ap, *bp, t); NUM_CADD(s, t, s);s5NUM_CASS(*tin, start); NUM_CADD(start, delta, start);s)((source[i].r !=0) || (source[i].i != 0))s* #include #define NA_ACOPYN(i, o) memcpy(o, i, N) /* The following is used to copy nbytes of data for each element. ** ** As such it can be used to align any sort of data provided the ** ** output pointers used are aligned */ static int copyNbytes(long dim, long nbytes, maybelong *niters, void *input, long inboffset, maybelong *inbstrides, void *output, long outboffset, maybelong *outbstrides) { long i, j; char *tin = (char *) input + inboffset; char *tout = (char *) output + outboffset; if (dim == 0) { for (i=0; i=0; i--) { otemp = PySequence_GetItem(nitersObj, i); if (PyInt_Check(otemp)) ltemp = PyInt_AsLong(otemp); else if (PyLong_Check(otemp)) ltemp = PyLong_AsLong(otemp); else return PyErr_Format(_Error, "copyToString: non-integer shape element"); nelements *= ltemp; niters[nniters-i-1] = ltemp; Py_DECREF(otemp); otemp = PySequence_GetItem(inbstridesObj, i); if (PyInt_Check(otemp)) inbstrides[nniters-i-1] = PyInt_AsLong(otemp); else if (PyLong_Check(otemp)) inbstrides[nniters-i-1] = PyLong_AsLong(otemp); else return PyErr_Format(_Error, "copyToString: non-integer stride element"); Py_DECREF(otemp); } if (!nelements) return PyString_FromStringAndSize("", 0); outbstrides[0] = nbytes; for (i=1; i= max) i -= max; static int takeNbytes(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { maybelong i, cMode, N; maybelong *scatteredstrides, *scatteredshape, **indices; char *gathered, *scattered; maybelong nindices = ninargs-4, outi = ninargs+noutargs-1; if (NA_checkIo("takeNbytes", 4, 1, MIN(ninargs, 4), noutargs)) return -1; if (nindices == 0) return 0; if (NA_checkOneCBuffer("takeNbytes", 2, buffers[0], bsizes[0], sizeof(maybelong))) return -1; else { cMode = ((maybelong *) buffers[0])[0]; N = ((maybelong *) buffers[0])[1]; } if (NA_checkOneCBuffer("takeNbytes", nindices, buffers[2], bsizes[2], sizeof(maybelong))) return -1; else { scatteredstrides = (maybelong *) buffers[2]; } if (NA_checkOneCBuffer("takeNbytes", nindices, buffers[3], bsizes[3], sizeof(maybelong))) return -1; else { scatteredshape = (maybelong *) buffers[3]; } if (NA_checkOneStriding("takeNBytes", nindices, scatteredshape, 0, scatteredstrides, bsizes[1], N, 0)) return -1; else scattered = (char *) buffers[1]; for(i=4; i= scatteredshape[j]) k = scatteredshape[j]-1; index += scatteredstrides[j]*k; } memcpy( &gathered[i*N], scattered+index, N); } break; case RAISE: for(i=0; i= scatteredshape[j]) { PyErr_Format(PyExc_IndexError, "Index[" ML_DEC "," ML_DEC "]=" ML_DEC " out of range[" ML_DEC "]", i, j, k, scatteredshape[j]); return -1; } index += scatteredstrides[j]*k; } memcpy( &gathered[i*N], scattered+index, N); } break; } return 0; } SELF_CHECKED_CFUNC_DESCR(takeNbytes, CFUNC_UFUNC); static int putNbytes(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { maybelong i, cMode, N; maybelong *scatteredstrides, *scatteredshape, **indices; char *gathered, *scattered; long nindices = ninargs-4, outi = ninargs+noutargs-1; if (nindices == 0) return 0; if (NA_checkIo("putNbytes", 4, 1, MIN(ninargs, 4), noutargs)) return -1; if (NA_checkOneCBuffer("putNbytes", 2, buffers[0], bsizes[0], sizeof(maybelong))) return -1; else { cMode = ((maybelong *) buffers[0])[0]; N = ((maybelong *) buffers[0])[1]; } if (NA_checkOneCBuffer("putNbytes", niter*N, buffers[1], bsizes[1], 1)) return -1; else gathered = (char *) buffers[1]; if (NA_checkOneCBuffer("putNbytes", nindices, buffers[2], bsizes[2], sizeof(maybelong))) return -1; else { scatteredstrides = (maybelong *) buffers[2]; } if (NA_checkOneCBuffer("putNbytes", nindices, buffers[3], bsizes[3], sizeof(maybelong))) return -1; else { scatteredshape = (maybelong *) buffers[3]; } for(i=4; i= scatteredshape[j]) k = scatteredshape[j]-1; index += scatteredstrides[j]*k; } memcpy( scattered+index, &gathered[i*N], N); } break; case RAISE: for(i=0; i= scatteredshape[j]) { PyErr_Format(PyExc_IndexError, "Index[" ML_DEC "," ML_DEC "]=" ML_DEC " out of range[" ML_DEC "]", i, j, k, scatteredshape[j]); return -1; } index += scatteredstrides[j]*k; } memcpy( scattered+index, &gathered[i*N], N); } break; } return 0; } SELF_CHECKED_CFUNC_DESCR(putNbytes, CFUNC_UFUNC); s /******************************************* * * * These copy data to a contiguous buffer. * * They do not handle non-aligned data. * * Offsets and Strides are in byte units * * * *******************************************/ static int copybytes(long dim, long dummy, maybelong *niters, void *input, long inboffset, maybelong *inbstrides, void *output, long outboffset, maybelong *outbstrides) { long i; char *tin = (char *) input + inboffset; char *tout = (char *) output + outboffset; if (dim == 0) { for (i=0; i(tin, tout); tin += inbstrides[dim]; tout += outbstrides[dim]; } } else { for (i=0; ibytes(dim-1, dummy, niters, input, inboffset + i*inbstrides[dim], inbstrides, output, outboffset + i*outbstrides[dim], outbstrides); } } return 0; } STRIDING_DESCR2(copybytes, CHECK_ALIGN, , ); s static int alignbytes(long dim, long dummy, maybelong *niters, void *input, long inboffset, maybelong *inbstrides, void *output, long outboffset, maybelong *outbstrides) { return copyNbytes(dim, , niters, input, inboffset, inbstrides, output, outboffset, outbstrides); } STRIDING_DESCR2(alignbytes, !CHECK_ALIGN, , ); s /******* byteswap *****/ static int byteswap(long dim, long dummy, maybelong *niters, void *input, long inboffset, maybelong *inbstrides, void *output, long outboffset, maybelong *outbstrides) { long i; char *tin = (char *) input + inboffset; char *tout = (char *) output + outboffset; if (dim == 0) { for (i=0; i]; NA_COPY(tin, t); (t, tout); tin += inbstrides[dim]; tout += outbstrides[dim]; } } else { for (i=0; i(dim-1, dummy, niters, input, inboffset + i*inbstrides[dim], inbstrides, output, outboffset + i*outbstrides[dim], outbstrides); } } return 0; } STRIDING_DESCR2(byteswap, !CHECK_ALIGN, , ); s static int choosebytes(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { maybelong i, cMode, maxP, N, *selector; char **population, *output; int outi = ninargs + noutargs - 1; if (NA_checkIo("choosebytes", 2, 1, MIN(ninargs,2), noutargs)) return -1; if (NA_checkOneCBuffer("choosebytes", 2, buffers[0], bsizes[0], sizeof(maybelong))) return -1; else { cMode = ((maybelong *) buffers[0])[0]; N = ((maybelong *) buffers[0])[1]; } if (NA_checkOneCBuffer("choosebytes", niter, buffers[1], bsizes[1], sizeof(maybelong))) return -1; else selector = (maybelong *) buffers[1]; if (ninargs-2 == 0) return 0; else maxP = ninargs-2; for(i=2; ibytes", niter, buffers[i], bsizes[i], )) return -1; population = (char **) &buffers[2]; if (NA_checkOneCBuffer("choosebytes", niter, buffers[outi], bsizes[outi], )) return -1; else output = (char *) buffers[outi]; if (maxP == 0) return 0; switch(cMode) { case WRAPPED: for(i=0; i(&population[j][i*], &output[i*]); } break; default: case CLIPPED: for(i=0; i= maxP) j = maxP-1; NA_ACOPY(&population[j][i*], &output[i*]); } break; } return 0; } SELF_CHECKED_CFUNC_DESCR(choosebytes, CFUNC_UFUNC); sh static PyMethodDef _Methods[] = { {NULL, NULL} /* Sentinel */ }; s { "%s", %s, METH_VARARGS}, s /* platform independent*/ #ifdef MS_WIN32 __declspec(dllexport) #endif void init_(void) { PyObject *m, *d; m = Py_InitModule("_", _Methods); d = PyModule_GetDict(m); _Error = PyErr_NewException("_.error", NULL, NULL); PyDict_SetItemString(d, "error", _Error); import_libnumarray(); PyDict_SetItemString(d, "functionDict", init_funcDict()); } s #include #include #include #include "libnumarray.h" #ifdef MS_WIN32 #pragma warning(once : 4244) #endif #define logical_and(arg1, arg2) (arg1 != 0) & (arg2 != 0) #define logical_or(arg1, arg2) (arg1 != 0) | (arg2 != 0) #define logical_xor(arg1, arg2) ((arg1 != 0) ^ (arg2 != 0)) & 1 #define ufmaximum(arg1, arg2) (((temp1=arg1) > (temp2=arg2)) ? temp1 : temp2) #define ufminimum(arg1, arg2) (((temp1=arg1) < (temp2=arg2)) ? temp1 : temp2) static PyObject *_Error; s static int __vector(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; *tin = ( *) buffers[0]; *tout = ( *) buffers[1]; for (i=0; i } return 0; } UFUNC_DESCR2(__vector, sizeof(), sizeof()); s static int __vector_vector(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; *tin1 = ( *) buffers[0]; *tin2 = ( *) buffers[1]; *tout = ( *) buffers[2]; for (i=0; i } return 0; } UFUNC_DESCR3(__vector_vector, sizeof(), sizeof(), sizeof()); s static int __vector_scalar(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; *tin1 = ( *) buffers[0]; tscalar = *( *) buffers[1]; *tout = ( *) buffers[2]; for (i=0; i } return 0; } UFUNC_DESCR3vs(__vector_scalar, sizeof(), sizeof(), sizeof()); s static int __scalar_vector(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; tscalar = *( *) buffers[0]; *tin2 = ( *) buffers[1]; *tout = ( *) buffers[2]; for (i=0; i } return 0; } UFUNC_DESCR3sv(__scalar_vector, sizeof(), sizeof(), sizeof()); s static int __reduce(long dim, long dummy, maybelong *niters, void *input, long inboffset, maybelong *inbstrides, void *output, long outboffset, maybelong *outbstrides) { long i; *tin = ( *) ((char *) input + inboffset); *tout = ( *) ((char *) output + outboffset); net; if (dim == 0) { net = *tout; for (i=1; i *) ((char *) tin + inbstrides[dim]); } *tout = net; } else { for (i=0; i__reduce(dim-1, dummy, niters, input, inboffset + i*inbstrides[dim], inbstrides, output, outboffset + i*outbstrides[dim], outbstrides); } } return 0; } STRIDING_DESCR2(__reduce, CHECK_ALIGN, sizeof(), sizeof()); sF static int __accumulate( long dim, long dummy, maybelong *niters, void *input, long inboffset, maybelong *inbstrides, void *output, long outboffset, maybelong *outbstrides) { long i; *tin = ( *) ((char *) input + inboffset); *tout = ( *) ((char *) output + outboffset); lastval; if (dim == 0) { for (i=1; i *) ((char *) tin + inbstrides[dim]); tout = ( *) ((char *) tout + outbstrides[dim]); } } else { for (i=0; i__accumulate(dim-1, dummy, niters, input, inboffset + i*inbstrides[dim], inbstrides, output, outboffset + i*outbstrides[dim], outbstrides); } } return 0; } STRIDING_DESCR2(__accumulate, CHECK_ALIGN, sizeof(), sizeof()); s static int __vector_vector(long niter, long ninargs, long noutargs, void **buffers, long *bsizes) { long i; *tin1 = ( *) buffers[0]; *tin2 = ( *) buffers[1]; *tout = ( *) buffers[2]; for (i=0; i } return 0; } UFUNC_DESCR3(__vector_vector, sizeof(), sizeof(), sizeof()); sU static PyObject *init_funcDict(void) { PyObject *dict; dict = PyDict_New(); s2 NA_add_cfunc(dict, "%s", (void *) &%s_descr); s return dict; } s operator1s(*tin, *tout);s operator_vvs (*tin1, *tin2, *tout);s operator_vss"(*tin1, tscalar, *tout);s operator_svs"(tscalar, *tin2, *tout);soperator_reduces(net, *tin, net);soperator_accumulates!(lastval ,*tin, *tout);s*tout = (*tin);s!*tout = (*tin1, *tin2);s#*tout = (*tin1, tscalar);s#*tout = (tscalar, *tin2);snet = (net, *tin);s"*tout = (lastval ,*tin);s operator_pps*tout = (tin1, tin2);s*tout = *tin;s*tout = *tin1 *tin2;s!*tout = *tin1 tscalar;s!*tout = tscalar *tin2;snet = net *tin;s *tout = lastval *tin;sU*tout = ((*tin2==0) ? int_dividebyzero_error(*tin2, *tin1) : *tin1 *tin2);s[*tout = ((tscalar==0) ? int_dividebyzero_error(tscalar, *tin1) : *tin1 tscalar);sS*tout = ((*tin2==0) ? int_dividebyzero_error(*tin2, 0) : tscalar *tin2);sJnet = ((*tin==0) ? int_dividebyzero_error(*tin, 0) : net *tin);sP*tout = ((*tin==0) ? int_dividebyzero_error(*tin, 0) : lastval *tin);s^*tout = ((*tin2==0) ? int_dividebyzero_error(*tin2, *tin1) : *tin1 (double) *tin2);sd*tout = ((tscalar==0) ? int_dividebyzero_error(tscalar, *tin1) : *tin1 (double) tscalar);s\*tout = ((*tin2==0) ? int_dividebyzero_error(*tin2, 0) : tscalar (double) *tin2);sSnet = ((*tin==0) ? int_dividebyzero_error(*tin, 0) : net (double) *tin);sY*tout = ((*tin==0) ? int_dividebyzero_error(*tin, 0) : lastval (double) *tin);se*tout = ((*tin2==0) ? int_dividebyzero_error(*tin2, *tin1) : floor(*tin1 (double) *tin2));sk*tout = ((tscalar==0) ? int_dividebyzero_error(tscalar, *tin1) : floor(*tin1 (double) tscalar));sc*tout = ((*tin2==0) ? int_dividebyzero_error(*tin2, 0) : floor(tscalar (double) *tin2));sZnet = ((*tin==0) ? int_dividebyzero_error(*tin, 0) : floor(net (double) *tin));s`*tout = ((*tin==0) ? int_dividebyzero_error(*tin, 0) : floor(lastval (double) *tin));scif (temp > 127) temp = int_overflow_error(127.); if (temp < -128) temp = int_overflow_error(-128.);s0if (temp > 255) temp = int_overflow_error(255.);sjif (temp > 32767) temp = int_overflow_error(32767.); if (temp < -32768) temp = int_overflow_error(-32768);s4if (temp > 65535) temp = int_overflow_error(65535.);sif (temp > 2147483647) temp = (Float64) int_overflow_error(2147483647.); if (temp < (-2147483648.)) temp = (Float64) int_overflow_error(-2147483648.);sIif (temp > 4294967295U) temp = (Float64) int_overflow_error(4294967295.);s\if (smult64_overflow(, )) temp = (Float64) int_overflow_error(4611686018427387903.);s\if (umult64_overflow(, )) temp = (Float64) int_overflow_error(9223372036854775807.);s*tin1s*tin2s*toutstscalarsnets*tinslastvals`temp = (() ) % (() ); = () temp;c Csh}xtiD]}thd|<d|hdt|d<dt|d<<d|<dt|d<dt|d<dt|d <|| temp1, temp2;sminusis-saddis+ssubtractsmultiplys*sdivides/s floor_dividesNUM_FLOORDIVIDEs true_divides remainders % sfmodspowersnum_powsabssfabsssinscosstansarcsinsasinsarccossacossarctansatansarctan2satan2slogsnum_logslog10s num_log10sexpssinhscoshstanhsarcsinhs num_asinhsarccoshs num_acoshsarctanhs num_atanhsieeemasksNA_IeeeSpecial32s vector_vectorsNA_IeeeSpecial64ssqrtsequals==s not_equals!=sgreaters>s greater_equals>=slesss>slshifts< 2 inputssnames_skeyN(#sTypesselfsmodules ufuncconfigsscfgsapplys UFuncParamssufoscodelistsappends separatorspynamessigsssigsninputssnoutputsslensrangesissetattrsintype0stemplatesNonesUFUNC_UNARY_TEMPLATEsvslistsUFUNC_BINARY_TEMPLATEsUFUNC_COMPARE_TEMPLATEs ValueErrors opt_declss__dict__soptional_declarationssformssig2saddcfuncsrepr(sselfsTypesformsufosiscfgssig2ssig((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pys bodygeneratorsL  *#    "(s__name__s __module__s__init__s bodygenerator(((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pysUFuncCodeGenerator|s s3/********************* %s *********************/ cCsF|tjod||f}nd||||f}t|dS(Ns#generate_%s_code('Src/_%smodule.c')s)generate_%s_code('Src/_%smodule.c', '%s')(sTypesNonesNamessseval(sNamesTypess((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pysgenerates cCs|tjod|}nd||}tiid}tii|otii|}nd}||jp dt i jodG|GHt ||ndS(NsSrc/_%smodule.csLib/codegenerator.pyis --check-datessgenerating new version of( sTypesNonesNamesSourcesosspathsgetmtimes timemodgensexistsstimemodssyssargvsgenerate(sNamesTypesSourcestimemods timemodgen((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pys PseudoMakes  c Cstdtdtddddddd d d d d dg }|idto|idnx|D]}td|qtWdS(NsconvssortsbytessBoolsInt8sUInt8sInt16sUInt16sInt32sUInt32sFloat32sFloat64s Complex32s Complex64sInt64sUInt64sufunc(s PseudoMakestypessappends hasUInt64sT(sTstypes((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pysGenAlls   '  s__main__(rssyssossimps find_modulesfilespathsdescrs load_modulestemplates_HEADERs_TYPE_SEPARATORsAS_PyVALUE_TEMPLATEsCOMPLEX_AS_PyVALUE_TEMPLATEsFROM_PyVALUE_TEMPLATEsCOMPLEX_FROM_PyVALUE_TEMPLATEsCONVERSION_TEMPLATEs_numarray_to_numarrays _complex_typesnumarray_to_numarrays SORT_HEADERs SORT_TEMPLATEssort_operator_tdssort_complex_tds BYTES_HEADERs COPY_TEMPLATEsALIGN_TEMPLATEsBYTESWAP_TEMPLATEsCHOOSE_TEMPLATEsBYTES_TEMPLATEs _METHODS_DEFsADDMETHOD_TEMPLATEs_TAILs UFUNC_HEADERsUFUNC_UNARY_TEMPLATEsUFUNC_VECTOR_VECTORsUFUNC_VECTOR_SCALARsUFUNC_SCALAR_VECTORs UFUNC_REDUCEsUFUNC_ACCUMULATEsUFUNC_PVECTOR_PVECTORsUFUNC_COMPARE_TEMPLATEsUFUNC_BINARY_TEMPLATEs_CFUNCDICT_HEADERs_ADDCFUNC_TEMPLATEs_CFUNCDICT_TAILsmacro_tds function_tdsmask_tds operator_tds int_divide_tdsint_truedivide_tdsint_floordivide_tds int8_checks uint8_checks int16_checks uint16_checks int32_checks uint32_checks int64_checks uint64_checksintX_tdsintX_mult_templatesmake_int_template_dicts hasUInt64s hasFloat128s typeconfigsappends ConvParamss bytesconfigs BytesParamss NBytesParamssComplexBytesParamssComplex32BytesCfgsComplex64BytesCfgsgeneric_operator_sigssgeneric_mathfunction_sigssgeneric_comparisonfunction_sigssgeneric_int_sigsstruedivide_int_sigssgeneric_float_sigsscomplex_complex_sigsscomplex_real_sigsscomplex_bool_sigss opt_mult_declsopt_mult32_declsopt_mult64_declsopt_minmax_decls sugar_dictsglobalssmultiply_int8_tdsmultiply_uint8_tdsmultiply_int16_tdsmultiply_uint16_tdsmultiply_uint32_tdsmultiply_int32_tdsmultiply_int64_tdsmultiply_uint64_tds ufuncconfigss UFuncParamss CodeGeneratorsConvCodeGeneratorsSortCodeGeneratorsBytesCodeGeneratorsUFuncCodeGeneratorsgenerate_conv_codesgenerate_sort_codesgenerate_bytes_codesgenerate_ufunc_codesNonesgenerates PseudoMakesGenAllsmains__name__(kscomplex_real_sigss uint32_checksConvCodeGeneratorsgeneric_operator_sigssSortCodeGeneratorsgenerate_sort_codes_TYPE_SEPARATORsmultiply_int16_tdssort_operator_tdsgeneric_float_sigss typeconfigsUFUNC_ACCUMULATEs ConvParamss UFUNC_HEADERsmultiply_uint32_tdsmains _complex_typesUFUNC_BINARY_TEMPLATEsADDMETHOD_TEMPLATEstruedivide_int_sigss hasFloat128sgeneric_mathfunction_sigssAS_PyVALUE_TEMPLATEsCONVERSION_TEMPLATEsBYTES_TEMPLATEsmultiply_uint64_tds function_tds BytesParamss ufuncconfigssCOMPLEX_AS_PyVALUE_TEMPLATEs _METHODS_DEFsBYTESWAP_TEMPLATEstemplatesmultiply_uint16_tdsComplex64BytesCfgs SORT_TEMPLATEsUFUNC_VECTOR_SCALARsUFUNC_VECTOR_VECTORs_HEADERs int32_checks COPY_TEMPLATEs uint8_checkspathsgeneratesUFUNC_PVECTOR_PVECTORsnumarray_to_numarraysgeneric_comparisonfunction_sigss_CFUNCDICT_TAILs hasUInt64soss_CFUNCDICT_HEADERsgeneric_int_sigssmultiply_int8_tdsUFuncCodeGeneratorsgenerate_ufunc_codes_TAILsmask_tdsgenerate_bytes_codesALIGN_TEMPLATEsimpsmultiply_int32_tdsgenerate_conv_codesGenAlls SORT_HEADERs int8_checksComplex32BytesCfgs UFuncParamssUFUNC_SCALAR_VECTORsmake_int_template_dictsintX_tds int64_checksmultiply_int64_tds NBytesParamss CodeGenerators uint64_checksopt_minmax_declsComplexBytesParamssmacro_tds int16_checks opt_mult_decls uint16_checksfilessort_complex_tdsint_truedivide_tdsCHOOSE_TEMPLATEsdescrsopt_mult32_declsUFUNC_UNARY_TEMPLATEsUFUNC_COMPARE_TEMPLATEscomplex_bool_sigss_ADDCFUNC_TEMPLATEsFROM_PyVALUE_TEMPLATEsint_floordivide_tds operator_tds PseudoMakes bytesconfigs_numarray_to_numarrays int_divide_tds BYTES_HEADERssyssCOMPLEX_FROM_PyVALUE_TEMPLATEsopt_mult64_decls UFUNC_REDUCEsmultiply_uint8_tdsBytesCodeGeneratorsintX_mult_templatescomplex_complex_sigs((sO/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/codegenerator.pys?s   B   4E E e$ "I < < < 333`     3W W W E E   {K%'!9