Glad to hear that it works. Note that dims is (optional) and different from d! Reading the docs again and looking at the code, I think it rather should be
d is the length of the transformed real array along the dims[1] dimension, which must satisfy d == (size(A,dims[1])-1)<<1 In my example above, size(a,1)==64 and size(a_from_rfft,1)==33 for instance - A.
