I posted to the ctypes-users mailing list about this
problem, but SourceForge failed to make that post and
its replies available in the ctypes-users archive, so
I'm repeating it here.

I'm using NumPy with ctypes. I would like to create a
callback function that calls into Python, allocates a
NumPy array and returns a suitable pointer to the C code.

NumPy has a function called ndpointer that allows one
to specify some details of NumPy arrays when dealing
with argtypes. This function also takes care of
extracting the array's data pointer.

Details here:

http://projects.scipy.org/scipy/numpy/browser/trunk/numpy/ctypeslib.py
http://projects.scipy.org/scipy/numpy/browser/trunk/numpy/core/_internal.py

Creating a callback function as follows works:

stuff = []
import numpy
def allocator1():
    x = numpy.array([...], dtype='f4')
    stuff.append(x)
    return x.ctypes.data_as(POINTER(c_float))
CFUNCTYPE(POINTER(c_float))(allocator1)

However, if one adds ndpointer to the mix, an error occurs:

stuff = []
import numpy
def allocator2():
    x = numpy.array([...], dtype='f4')
    stuff.append(x)
    return x
CFUNCTYPE(numpy.ctypeslib.ndpointer('f4'))(allocator2)

The error is:

SystemError: NULL result without error in PyObject_Call

Thomas Heller has a patch for this issue.

Logged In: YES 
user_id=11105

At least the 'SystemError: NULL result without error in
PyObject_Call' is fixed now in SVN revision 52367
(release25-maint), and revision 52365 (trunk).

You will now get a 'TypeError: unsupported result type for
callback function'.

That's all I can do for Python 2.5, anyway.

I suggest that possible workarounds should be discussed on
the ctypes-users list.