[Scipy-svn] r4517 - trunk/scipy/sandbox/mkufunc/examples

scipy-svn@scip... scipy-svn@scip...
Tue Jul 1 16:38:29 CDT 2008


Author: ilan
Date: 2008-07-01 16:38:28 -0500 (Tue, 01 Jul 2008)
New Revision: 4517

Modified:
   trunk/scipy/sandbox/mkufunc/examples/benchmark.py
Log:
simplified benchmark C code

Modified: trunk/scipy/sandbox/mkufunc/examples/benchmark.py
===================================================================
--- trunk/scipy/sandbox/mkufunc/examples/benchmark.py	2008-07-01 16:24:20 UTC (rev 4516)
+++ trunk/scipy/sandbox/mkufunc/examples/benchmark.py	2008-07-01 21:38:28 UTC (rev 4517)
@@ -1,6 +1,6 @@
 #!/usr/bin/env python
 from math import sin, cos
-import time
+import time, hashlib
 
 from numpy import arange, vectorize, allclose
 from scipy import weave
@@ -18,13 +18,6 @@
 
 #####################################################################
 support_code = '''
-double foo(double x)
-{
-        return x;
-}
-
-typedef double Func_0(double a);
-
 static void
 PyUFunc_0(char **args, npy_intp *dimensions, npy_intp *steps, void *func)
 {
@@ -33,13 +26,13 @@
         npy_intp os = steps[1];
         char *ip0 = args[0];
         char *op = args[1];
-        Func_0 *f = (Func_0 *) func;
         n = dimensions[0];
+        double x;
         
         for(i = 0; i < n; i++) {
-                double x = (double *)ip0;
+                x = *(double *)ip0;
                 double *out = (double *)op;
-
+                
                 *out = 4.2 * x*x + 3.7 * x + 1.5;
 
                 ip0 += is0;
@@ -51,23 +44,19 @@
         PyUFunc_0,
 };
 
-static void *f_data[] = {
-        (void *) foo,
-};
-
 static char f_types[] = {
-        NPY_DOUBLE, NPY_DOUBLE,   /* 0 */
+        NPY_DOUBLE, NPY_DOUBLE,
 };
 '''
 ufunc_info = weave.base_info.custom_info()
 ufunc_info.add_header('"numpy/ufuncobject.h"')
 
-ufunc = weave.inline('''
+ufunc = weave.inline('/*' + hashlib.md5(support_code).hexdigest() + '''*/
 import_ufunc();
 
 return_val = PyUFunc_FromFuncAndData(
             f_functions,
-            f_data,
+            NULL,
             f_types,
             1,             /* ntypes */
             1,             /* nin */



More information about the Scipy-svn mailing list