GFourier

The implementation of the multiply operation in operations.c appears to be
violating the assumption
that the complex product is distinct from the operands (in the case where
new_main == main).
Either the product should be a temp, or the GFCOMPLEX_MUL operation should
incorporate a temp
itself to avoid other inadvertent problems. Here is a possible
definition:

#define GFCOMPLEX_MUL(a,b,c) \
{GFREAL _a_re; _a_re = (b).re * (c).re - (b).im * (c).im; (a).im =
(b).re
* (c).im + (b).im * (c).re; (a).re = _a_re;}