#10058: Segfault in backward and inverse FFT for 2**n elements
-------------------------+--------------------------------------------------
Reporter: mpatel | Owner: jason, jkantor
Type: defect | Status: new
Priority: major | Milestone: sage-4.6.1
Component: numerical | Keywords:
Author: | Upstream: N/A
Reviewer: | Merged:
Work_issues: |
-------------------------+--------------------------------------------------
Comment(by mpatel):
This happens with `a.inverse_transform()`, too. From near the end of
`sage/gsl/fft.pyx`:
{{{
#!python
def backward_transform(self):
cdef gsl_fft_complex_wavetable * wt
cdef gsl_fft_complex_workspace * mem
N = Integer(self.n)
e = N.exact_log(2)
if N==2**e:
gsl_fft_complex_backward(self.data, self.stride, self.n, wt,
mem)
}}}
Should this be
{{{
#!python
gsl_fft_complex_radix2_backward(self.data, self.stride,
self.n)
}}}
(and similarly for `inverse_transform`)?
{{{
#!python
if N!=2**e:
mem = gsl_fft_complex_workspace_alloc(self.n)
wt = gsl_fft_complex_wavetable_alloc(self.n)
gsl_fft_complex_backward(self.data, self.stride, self.n, wt,
mem)
gsl_fft_complex_workspace_free(mem)
gsl_fft_complex_wavetable_free(wt)
}}}
--
Ticket URL: <http://trac.sagemath.org/sage_trac/ticket/10058#comment:1>
Sage <http://www.sagemath.org>
Sage: Creating a Viable Open Source Alternative to Magma, Maple, Mathematica,
and MATLAB
--
You received this message because you are subscribed to the Google Groups
"sage-trac" group.
To post to this group, send email to [email protected].
To unsubscribe from this group, send email to
[email protected].
For more options, visit this group at
http://groups.google.com/group/sage-trac?hl=en.
