Hi, I came across the following problem: I was using cufft (using the ctypes wrapper provided in scikits.cuda) which returns results in form of pycuda.gpuarrays. However, since the arrays on the GPU seem to have Fortran (row-major) ordering and the default ordering for numpy arrays is C (column-major), the array entries appear to be mixed up.
The problem could be fixed by changing the get method of the gpuarray class (line 111 in gpuarray.py) from ary = numpy.empty(self.shape, self.dtype) to ary = numpy.empty(self.shape, self.dtype, order="F") However, I am new to using pycuda, and I am not sure if C ordering is assumed in other places. Maybe, somebody can comment on how to best assure compatibility with the cufft library. Thanks, Kilian _______________________________________________ PyCUDA mailing list [email protected] http://lists.tiker.net/listinfo/pycuda
