Lorenzo, you can indeed use f2py to write extensions around some C code: http://cens.ioc.ee/projects/f2py2e/usersguide/index.html http://www.scipy.org/Cookbook/f2py_and_NumPy I think you should also be able to find some actual examples in the scipy sources...