Test with batch FFT's over entire grid
Towards faster y-direction FFTs.
All y-direction FFTs in separate loops, then we only have to call threadSynchronize once.
Single call for the complex->double and double->complex transformation
Cleanup of pres_2.cu. Overall speedup compared to last tag ~ 64% for 64^3 drycblles
Removed the unused complex arrays
Merge remote-tracking branch 'upstream/master' into cuda_fastpres
Added the faster FFTs to the 4th order solver
Normalization over entire domain instead of per slice
Also added normalization per domain to 2nd order version
Quick and dirty test with transposing the data for allowing batch FFT over entire domain in y-direction.
Cleanup, and added transpose for non-square grids
y-direction FFT with transpose on non-square grids working
Unfortunately the speedup by doing a transpose before the y-direction FFT only works (very well) for small grid sized, tipping point is somewhere around 128-256 grid points
Calculate dimensions thread block from grid dimensions. For small grid (an my slow home gpu) another factor ~2 faster
Quick check if cufftplan failed.
minor cleanup
Minor cleanup of unused FFT plans
Moved forward and backward FFTs to their own routines to cleanup the pressure exec function
As temporary performance fix, calculate hydrostatic pressure on host and copy back to device