[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: [ccp4bb]: Error running FFT in alpha

***  For details on how to be removed from this list visit the  ***
***    CCP4 home page http://www.dl.ac.uk/CCP/CCP4/main.html    ***


first off, what Alpha machine are you using.  Previously we have had
reports of problems using fft on the XP1000, divide by zero in that case.
That seemed to have been caused by various arrays being overwritten,
although the problem was not fully tracked down.

Several solutions were found though.

1.  use FFTBIG.  This worked.  It uses dynamic memory allocation.
2.  increase the size of the working array, X.  Set IDUM to 1000000.  This
    will work if the problem arises in the sorting routines by bypassing
    the multiple passes.
3.  turning off optimization seemed to also cure the problem.  Could be
    that it arose due code reordering.  Warning: this leads to a large
    performance hit.

Hope these suggestions help.

Charles Ballard

ps- anyone using a new Alpha have any further information?