{"action":"created","issue":{"url":"https://api.github.com/repos/hgomersall/pyFFTW/issues/56","labels_url":"https://api.github.com/repos/hgomersall/pyFFTW/issues/56/labels{/name}","comments_url":"https://api.github.com/repos/hgomersall/pyFFTW/issues/56/comments","events_url":"https://api.github.com/repos/hgomersall/pyFFTW/issues/56/events","html_url":"https://github.com/hgomersall/pyFFTW/issues/56","id":56285610,"number":56,"title":"`Datatype not supported` with scipy","user":{"login":"caspervdw","id":8937919,"avatar_url":"https://avatars.githubusercontent.com/u/8937919?v=3","gravatar_id":"","url":"https://api.github.com/users/caspervdw","html_url":"https://github.com/caspervdw","followers_url":"https://api.github.com/users/caspervdw/followers","following_url":"https://api.github.com/users/caspervdw/following{/other_user}","gists_url":"https://api.github.com/users/caspervdw/gists{/gist_id}","starred_url":"https://api.github.com/users/caspervdw/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/caspervdw/subscriptions","organizations_url":"https://api.github.com/users/caspervdw/orgs","repos_url":"https://api.github.com/users/caspervdw/repos","events_url":"https://api.github.com/users/caspervdw/events{/privacy}","received_events_url":"https://api.github.com/users/caspervdw/received_events","type":"User","site_admin":false},"labels":[],"state":"open","locked":false,"assignee":null,"milestone":null,"comments":24,"created_at":"2015-02-02T20:03:00Z","updated_at":"2015-06-05T08:03:12Z","closed_at":null,"body":"Hi,\r\n\r\nI am working with trackpy for particle tracking and it has an optionality to use pyfftw for bandpass filters. On my windows (64-bit) box, I get a `RuntimeError: datatype not supported`, when I run a gaussian filter on the FFTd image. I can reproduce the problem as follows:\r\n\r\n    import numpy as np\r\n    import pyfftw\r\n    from scipy.ndimage.fourier import fourier_gaussian\r\n    \r\n    a = np.random.randint(0, 256, (50,50))\r\n    a = pyfftw.n_byte_align(a, a.dtype.alignment)\r\n    afftw = pyfftw.interfaces.numpy_fft.fftn(a)\r\n    \r\n    fourier_gaussian(afftw, (1., 1.))  # RuntimeError: data type not supported\r\n\r\nCopying the afftw to a new empty numpy array solves the problem. The difference I can see between `afftw` and an FFT which was produced by `numpy.fft` is that pyfftw has `strides == (800L, 16L)` whereas the numpy one has `strides == (16L, 800L)`. But both should work with scipy, and transposing or asfortranarray or ascontiguousarray doesn't help.\r\n\r\nSee also https://github.com/soft-matter/trackpy/issues/198"},"comment":{"url":"https://api.github.com/repos/hgomersall/pyFFTW/issues/comments/109197501","html_url":"https://github.com/hgomersall/pyFFTW/issues/56#issuecomment-109197501","issue_url":"https://api.github.com/repos/hgomersall/pyFFTW/issues/56","id":109197501,"user":{"login":"hgomersall","id":701883,"avatar_url":"https://avatars.githubusercontent.com/u/701883?v=3","gravatar_id":"","url":"https://api.github.com/users/hgomersall","html_url":"https://github.com/hgomersall","followers_url":"https://api.github.com/users/hgomersall/followers","following_url":"https://api.github.com/users/hgomersall/following{/other_user}","gists_url":"https://api.github.com/users/hgomersall/gists{/gist_id}","starred_url":"https://api.github.com/users/hgomersall/starred{/owner}{/repo}","subscriptions_url":"https://api.github.com/users/hgomersall/subscriptions","organizations_url":"https://api.github.com/users/hgomersall/orgs","repos_url":"https://api.github.com/users/hgomersall/repos","events_url":"https://api.github.com/users/hgomersall/events{/privacy}","received_events_url":"https://api.github.com/users/hgomersall/received_events","type":"User","site_admin":false},"created_at":"2015-06-05T08:03:12Z","updated_at":"2015-06-05T08:03:12Z","body":"So does dtype('longdouble') on windows show as dtype('float64')? I mean, is it possible to have an array that is dtype('longdouble')?"}}