<#part sign=pgpmime> On Tue, 27 Mar 2012 14:30:44 -0400, Alex Nitz <[email protected]> wrote: > It seems that pycuda does not currently support the ** operator (pow) for > complex dtypes (complex64, complex128). Is this something that might be on > the horizon? > > I am currently working on a project that would benefit from the pycuda > gpuarray working the same as the numpy array. We would like to accelerate > portions of our code, but need to maintain full compatibility on machines > that do not have a GPU. How alike is the gpuarray intended to become to > the numpy array? > > I understand that this may not be on anyone else's priority list, so I may > try to make a stab at a patch for this if nobody else comes along. I don't > have much experience with the inner workings of pycuda yet, so that may > take a while, but if I get something together, can someone enlighten me on > the patch process, or point me to a link?
Huh, weird. pow(complex, complex) is implemented in the C++ header, so in principle this should work. Can you insert a keep=True in the call to get_elwise_kernel in get_pow_kernel in pycuda/elementwise.py? That should give you the code for the failing kernel, and a better shot at making sense of the nvcc error message. Andreas _______________________________________________ PyCUDA mailing list [email protected] http://lists.tiker.net/listinfo/pycuda
