Opened 16 years ago

Closed 16 years ago

#4 closed defect (fixed)

on powerpc, fft(zeros) sometimes gives phase = pi

Reported by: Paul Brossier Owned by: Paul Brossier
Priority: normal Milestone: 0.3.3
Component: corelib Version: 0.3.2
Severity: blocker Keywords:


when computing the fft on a vector of zeroes, the phase elements can be found equal to zero or pi on powerpc, whereas they are always zero on x86.

this is due to the way real and imaginary parts are computed in fftw on both platforms. on x86. the result is always +0.0000, whereas on ppc, it can be -0.0000; atan2f(-0.00,0.00) = pi

see test in source:tests/python/src/spectral/, test_do_zeroes

Change History (1)

comment:1 Changed 16 years ago by Paul Brossier <>

Resolution: fixed
Status: newclosed

(In ,543) closes #4, fft.c: if real part is zero, have phase = 0

Note: See TracTickets for help on using tickets.