replace old q direction generator of FFTIndexing

The current algorithm used for generating q direction over the unit sphere prior FFT indexing is not really clear and introduces a bias along z direction. I introduced an algorithm based on Fibonacci algorithm (see here for more information)