FFTSS Library Version 3.0 User's Guide

Size: px

Start display at page:

Download "FFTSS Library Version 3.0 User's Guide"

かずまさこのえ
5 years ago
Views:

1 : FFTSS 3.0 Copyright (C) The Scalable Software Infrastructure Project, (CREST),,.

2 Contents (DFT) UNIX Microsoft Windows Visual Studio C/C MinGW UNIX Visual Studio fftss malloc fftss free fftss plan dft 1d fftss plan dft 2d fftss plan dft 3d fftss execute fftss execute dft fftss destroy plan fftss get wtime fftss init threads fftss plan with nthreads fftss cleanup threads MPI pfftss plan dft 2d pfftss execute pfftss execute dft pfftss destroy plan

3 8 FFTW List of FFT

4 1 FFTSS (FFT). (JST) (CREST).,. FFT.,., FFTW 3. FFTW. FFTW,. 2 (DFT) FFTSS n 1 DFT. Y k = Y k = n X j e 2πjk 1/n j=0 n X j e 2πjk 1/n. j=0 X, Y. FFTSS FFTW, 1/N. DFT 1 DFT. 3.,. 3.1 UNIX UNIX configure. 2. make. 3. make install. ( ) configure.

5 --without-simd --without-asm --with-bg --with-bg-compat --with-recommended --enable-openmp --enable-mpi SIMD.. IBM Blue Gene. ( ) Blue Gene FFT. CC CFLAGS. OpenMP. MPI. CC CFLAGS. $./configure CC=gcc CFLAGS= -O3 -msse2 3.2 Microsoft Windows Microsoft Windows Visual Studio win32 fftss.sln. Visual Studio.NET Visual Studio Visual Studio Microsoft SDK ( Platform SDK). win32 pfftss.sln MPI. Microsoft Compute Cluster Pack SDK MPI C/C++ win32 C/C++.. win32.. icl-x86.bat IA-32. (32 ) icl-amd64 EM64T (Intel 64) AMD64. (64 ) icl-ia64 IA-64. (64bit) Visual Studio,, MinGW MinGW UNIX configure UNIX. CFLAGS -I/path/to/header/file

6 LDFLAGS -L/path/to/library/file -lfftss LDFLAGS -lpfftss (MPI ) /path/to/header/file, /path/to/library/file fftss.h, libfftss.a Visual Studio Visual Studio, win32/release ( win32/debug) fftss.lib. MPI pfftss.lib. fftss.h pfftss.h include. FFTSS. 5 FFTSS. 2. FFT Stockham out-of-place transform. in-place transform out-of-place transform fftss malloc void *fftss malloc(long size); fftss malloc() size, fftss free void fftss free(void *ptr); fftss free() ptr. fftss malloc() fftss plan dft 1d fftss plan fftss plan dft 1d(long n, double *in, double *out, long sign, long flags );

7 fftss plan dft 1d() n 1 FFT. i in[i*2], in[i*2+1] fftss plan dft 2d fftss plan fftss plan dft 2d(long nx, long ny, long py, double *in, double *out, long sign, long flags ); fftss plan dft 2d() nx ny 2 FFT. (x,y) in[x*2+y*py*2], in[x*2+y*py*2+1] fftss plan dft 3d fftss plan fftss plan dft 3d(long nx, long ny, long nz, long py, long pz, double *in, double *out, long sign, long flags ); fftss plan dft 3d() nx ny nz 3 FFT. (x,y,z) in[x*2+y*py*2+z*pz*2], in[x*2+y*py*2+z*pz*2+1] , 2, 3 FFT. FFTSS VERBOSE.. FFT. FFTSS MEASURE FFT,. ( ) FFTSS ESTIMATE FFT. FFT,. FFTSS PATIENT FFTSS MEASURE. FFTSS EXHAUSTIVE FFTSS MEASURE. FFTSS NO SIMD SIMD ( SIMOMD). FFT. FFTSS UNALIGNED 16..

8 ,. fftss execute dft(),., fftss malloc(). FFTSS DESTROY INPUT. ( ) FFTSS PRESERVE INPUT.. FFTSS INOUT, in. out fftss execute void fftss execute(fftss plan p ); fftss execute() p fftss execute dft void fftss execute dft(fftss plan p, double *in, double *out); fftss execute() p. in out p fftss destroy plan void fftss destroy plan(fftss plan p ); fftss destroy plan() p, fftss get wtime double fftss get wtime(void);

9 fftss get wtime() fftss init threads int fftss init threads(void); fftss init threads(). FFTW fftss plan with nthreads void fftss plan with nthreads(int nthreads); fftss plan with nthreads(). FFTSS OpenMP, omp set num threads() OpenMP fftss cleanup threads void fftss cleanup threads(void); fftss cleanup threads(). FFTW. 6.7 MPI FFTSS 3 pfftss MPI pfftss plan dft 2d pfftss plan pfftss plan dft 2d(long nx, long ny, long py, long oy, long ly, double *inout, long sign, long flags, MPI Comm comm); pfftss plan dft 2d() MPI nx ny 2 FFT. inout,. oy ly, py ly 2. ( 1 double inout[2*py*ly])., comm. nx. ly pfftss execute void pfftss execute(pfftss plan p);

10 pfftss execute() p pfftss execute dft void pfftss execute dft(pfftss plan p, double *inout); pfftss execute dft() p. inout pfftss destroy plan void pfftss destroy plan(pfftss plan p); pfftss destroy plan() p. 7 FFTSS OpenMP., OpenMP, OpenMP. C/C++. -openmp, -xp SSE3. $./configure CC=icc CFLAGS= -O3 -openmp -xp OMP NUM THREADS., CPU 1. omp set num threads(). fftss plan with nthreads() FFTW, omp set num threads()...,. max_threads = omp_get_num_procs(); fftss_plan_with_nthreads(max_threads); plan = fftss_plan_dft_2d(nx, ny, py, vin, vout, FFTSS_FORWARD, FFTSS_MEASURE); { /*. */ } for (nthreads = 1; nthreads <= max_threads; nthreads ++) { fftss_plan_with_nthreads(nthreads); t = fftss_get_wtime(); fftss_execute(plan); t = fftss_get_wtime() - t; printf("%d %lf.\n", nthreads, t); }

11 MPI OpenMP. enable-mpi enable-openmp MPI OpenMP. 8 FFTW FFTW. FFTW FFTSS. FFTW fftw3.h, fftw3compat.h FFTSS. fftw3compat.h FFTW FFTSS, fftss.h. FFTW fftw3compat.h. MPI. 8.1 fftw malloc() fftw free() fftw plan dft 1d() fftw plan dft 2d() fftw plan dft 3d() fftw execute() fftw execute dft() fftw destroy plan() fftw init threads() fftw plan with nthreads() fftw cleanup threads() 8.2 FFTW MEASURE FFTW ESTIMATE FFTW PATIENT FFTW EXHAUSTIVE FFTW NO SIMD FFTW PRESERVE INPUT

12 FFTW DESTROY INPUT FFTW FORWARD FFTW BACKWARD 9 List of FFT FFT., FFTSS VERBOSE. normal. FMA FFT. SSE2 (1) Intel SSE2. SSE2 (2) Intel SSE2 (UNPCKHPD/UNPCKLPD). SSE3 Intel SSE2 (ADDSUBPD). SSE3 (H) Intel SSE2 (HADDPD/HSUBPD). C99 Complex C99. Blue Gene IBM Blue Gene. Blue Gene (PL) IBM Blue Gene ( ). Blue Gene asm IBM Blue Gene ( ). IA-64 asm Intel IA-64 ( ).

13 10 FFTSS. UltraSPARC III Sun Solaris 9 Sun ONE Studio 11 Itanium 2 Linux Intel C/C++ Compiler 9.1, gcc PowerPC G5 Mac OS X 10.4 IBM XL C Compiler 6.0, gcc 4.0 POWER5 Linux IBM XL C Compiler 7.0, gcc POWER4 AIX IBM XL C Compiler 6.0 PA-RISC HP-UX 11 Bundled C Compiler PPC440FP2 Blue Gene CNK IBM XL C Compiler 7.0/8.0 Opteron Linux gcc 3.3.3, gcc Pentium 4 Solaris 9 (IA-32) Sun ONE Studio 11, gcc Xeon Linux Intel C/C++ Compiler 8.1/9.0/9.1, gcc IA-32 Windows XP SP2 Visual Studio.NET 2003 IA-32 Windows XP SP2 Visual Studio 2005 IA-32 Windows XP SP2 Intel C/C++ Compiler 9.1 x64 Windows XP, 2003 Visual Studio.NET 2003 x64 Windows XP, 2003 Intel C/C++ Compiler for EM64T 9.1

CPU Levels in the memory hierarchy Level 1 Level 2... Increasing distance from the CPU in access time Level n Size of the memory at each level 1: 2.2

CPU Levels in the memory hierarchy Level 1 Level 2... Increasing distance from the CPU in access time Level n Size of the memory at each level 1: 2.2 FFT 1 Fourier fast Fourier transform FFT FFT FFT 1 FFT FFT 2 Fourier 2.1 Fourier FFT Fourier discrete Fourier transform DFT DFT n 1 y k = j=0 x j ω jk n, 0 k n 1 (1) x j y k ω n = e 2πi/n i = 1 (1) n DFT