Difference between revisions of "Development/FFTW"

From bwHPC Wiki
Jump to: navigation, search
m (Description)
m (S Richling moved page FFTW to Development/FFTW: Umzug in Development-Bereich)
 
(132 intermediate revisions by 6 users not shown)
Line 1: Line 1:
  +
{{Softwarepage|numlib/mkl}}
  +
 
{| width=600px class="wikitable"
 
{| width=600px class="wikitable"
 
|-
 
|-
Line 4: Line 6:
 
|-
 
|-
 
| module load
 
| module load
| numlib/fftw
+
| numlib/mkl
|-
 
| Availability
 
| [[bwUniCluster]] | [[BwForCluster_Chemistry]]
 
 
|-
 
|-
 
| License
 
| License
  +
| Commercial. See [https://software.intel.com/en-us/articles/end-user-license-agreement EULA].
| [http://www.fftw.org/faq/section1.html#isfftwfree GPL] | [http://www.fftw.org/fftw3_doc/License-and-Copyright.html License and Copyright]
 
 
|-
 
|-
 
|Citing
 
|Citing
  +
| n/a
| [https://github.com/FFTW/fftw3/blob/master/AUTHORS Authors]
 
 
|-
 
|-
 
| Links
 
| Links
| [http://www.fftw.org/ FFTW Homepage] | [http://www.fftw.org/#documentation Documentation]
+
| [https://software.intel.com/en-us/intel-mkl Intel MKL Homepage] | [http://www.fftw.org/ FFTW Homepage]
 
|-
 
|-
 
| Graphical Interface
 
| Graphical Interface
 
| No
 
| No
 
|}
 
|}
<br>
+
<br>
 
= Description =
 
= Description =
'''FFTW''' is a C subroutine library for computing the discrete Fourier transform (DFT) in one or more dimensions, of arbitrary input size, and of both real and complex data (as well as of even/odd data, i.e. the discrete cosine/sine transforms or DCT/DST). We believe that FFTW, which is free software, should become the FFT library of choice for most applications.
 
<br>
 
[http://www.fftw.org/ More infos about FFTW]
 
<br>
 
This package provides three versions of the fftw3 library depending on precision:
 
* libfft3
 
* libfftw3f
 
* libfftw3l <small>for double, single and long-double precision libraries</small>
 
<br>
 
 
= Documentation =
 
== Online ==
 
[http://www.fftw.org/fftw3_doc/ Online Documentation]
 
== Local ==
 
* info fftw3
 
* man fftw-wisdom
 
* man fftw-wisdom-to-conf
 
<br>
 
See also documentation folder pointed to by shell variable $FFTW_DOC_DIR
 
<br>
 
<br>
 
= Versions and Availability =
 
 
A list of versions currently available on all bwHPC-C5-Clusters can be obtained from the
 
<br>
 
<big>
 
 
[https://cis-hpc.uni-konstanz.de/prod.cis/ Cluster Information System CIS]
 
 
</big>
 
{{#widget:Iframe
 
|url=https://cis-hpc.uni-konstanz.de/prod.cis/bwUniCluster/numlib/fftw
 
|width=99%
 
|height=200
 
|border=0
 
}}
 
On the command line interface of any bwHPC cluster you'll get a list of available versions
 
by using the command 'module avail chem/gamess'.
 
<pre>
 
$ module avail chem/gamess
 
------------------------ /opt/bwhpc/common/modulefiles -------------------------
 
chem/gamess/12052014
 
</pre>
 
<br>
 
 
= Loading =
 
Load the fftw module, and, if needed, the corresponding openmpi module.
 
 
<br>
 
<br>
 
= Hints for compiling and linking =
 
== FFTW Specific Environments ==
 
To see a list of all FFTW environments set by the 'module load'-command use the command
 
'module display numlib/fftw'.
 
<br>
 
After having loaded the appropriate module(s), you can use several environment variables to compile and link your application.
 
<pre>
 
$ module display numlib/fftw
 
-------------------------------------------------------------------
 
/opt/bwhpc/common/modulefiles/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1:
 
FTW_VERSION = 3.3.3-impi-4.1.1-intel-13.1
 
FFTW_HOME = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1
 
FFTW_BIN_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/bin
 
FFTW_INC_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/include
 
FFTW_LIB_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/lib
 
FFTW_STA_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/lib
 
FFTW_SHA_DIR ??
 
FFTW_MAN_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/share/man
 
FFTW_INF_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/share/info
 
FFTW_DOC_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/doc
 
FFTW_EXA_DIR = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/examples
 
FFTW_WWW = http://www.fftw.org/
 
PATH = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/bin:$PATH
 
MANPATH = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/share/man:$MANPATH
 
INFOPATH = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/share/info:$INFOPATH
 
LD_LIBRARY_PATH = /opt/bwhpc/common/numlib/fftw/3.3.3-impi-4.1.1-intel-13.1/lib:$LD_LIBRARY_PATH
 
[...]
 
</pre>
 
== Compile a Serial Program ==
 
<pre>
 
$ gcc example.c -o example -I$FFTW_INC_DIR -L$FFTW_LIB_DIR -lfftw3 -lm
 
</pre>
 
=== with POSIX Threads ===
 
Compile program with support for POSIX threads:
 
<pre>
 
$ gcc example.c -o example -I$FFTW_INC_DIR -L$FFTW_LIB_DIR -lfftw3_threads -lfftw3 -lpthread -lm
 
</pre>
 
   
  +
The '''Fastest Fourier Transform in the West (FFTW)''' is a software library for computing discrete Fourier transforms in one or more dimensions, of arbitrary input size, and of both real and complex data (as well as of even/odd data, i.e. the discrete cosine/sine transforms or DCT/DST). FFTW was developed by Matteo Frigo and Steven G. Johnson at the Massachusetts Institute of Technology.
== Compile program with support for OpenMP threads ==
 
<pre>
 
$ gcc example.c -o example -fopenmp -I$FFTW_INC_DIR -L$FFTW_LIB_DIR -lfftw3_omp -lfftw3 -lm
 
</pre>
 
   
  +
The '''Intel Math Kernel Library (Intel MKL)''' offers FFTW2 (for version 2.x) and FFTW3 (for version 3.x) interfaces to the Intel MKL Fast Fourier Transform and Trigonometric Transform functionality. These interfaces enable applications using FFTW to gain performance with Intel MKL without changing the application source code. Therefore, it is highly recommended to use Intel MKL instead of a separate FFTW installation.
== Compile program with support for MPI ==
 
<pre>
 
$ mpicc example.c -o example -I$FFTW_INC_DIR -L$FFTW_LIB_DIR -lfftw3_mpi -lfftw3 -lm
 
</pre>
 
   
  +
= FAQ =
== Run program with MPI support ==
 
<pre>
 
$ mpirun -n <ncpu> ./example
 
</pre>
 
   
  +
[[File:comparison.png|right|border|300px|Copyright: KIZ (Ulm University)]]
(Replace <ncpu> by number of processor cores.)
 
   
  +
'''Q:''' Why is there no FFTW module on the cluster?
Replace -lfftw3, -lfftw3_threads, etc. by -lfftw3f, -lfftw3f_threads, etc. for single
 
precision and by -lfftw3l, -lfftw3l_threads etc. for long-double precision codes, respectively.
 
   
  +
'''A:''' MKL exhibits better performance than FFTS libraries (see Figure on the right). Therefore, we recommend to use MKL and do not offer a separate FFTW installation.
These commands will compile your program with dynamic fftw library versions in
 
which case you also have to have the fftw module loaded for running the program.
 
Alternatively, you may want to link your program with static fftw library versions.
 
With static fftw libraries it is only necessary to load the fftw module for compiling
 
but not for executing the program.
 
   
  +
'''Q:''' Why does my code complain about <span style="background:#edeae2;margin:2px;padding:1px;border:1px dotted #808080"> argument of type "long double *" is incompatible with parameter of type "double *" </span>?
== Compile program with static fftw library versions ==
 
Example for POSIX threads support
 
<pre>
 
$ gcc example.c -o example -I$FFTW_INC_DIR $FFTW_LIB_DIR/{libfftw3_threads.a,libfftw3.a} -lpthread -lm
 
</pre>
 
or:
 
<pre>
 
$ gcc example.c -o example -I$FFTW_INC_DIR -L$FFTW_LIB_DIR -Wl,-Bstatic -lfftw3 -lfftw3_threads \
 
-Wl,-Bdynamic -lpthread -lm
 
</pre>
 
   
  +
'''A:''' The interfaces do not support long double precision because Intel MKL FFT functions operate only on single- and double-precision floating point data types. For the very rare case that you need extended data types, please submit a support ticket at https://www.bwhpc.de/supportportal.
Environment variables $FFTW_INC_DIR, $FFTW_LIB_DIR etc. are available after
 
loading the module.
 
   
  +
= Useful links =
Sample code for various test cases is provided in folder pointed to by environment
 
variable $FFTW_EXA_DIR.
 
   
  +
* [https://software.intel.com/content/www/us/en/develop/articles/intel-math-kernel-library-documentation.html Documentation (english)]
<!--
 
  +
* [https://de.wikipedia.org/wiki/Math_Kernel_Library Wikipedia article (german)]
[[Category:Numerical_libraries]][[Category:bwUniCluster]]
 
  +
* [https://en.wikipedia.org/wiki/Math_Kernel_Library Wikipedia article (english)]
-->
 

Latest revision as of 00:04, 15 March 2023

The main documentation is available via module help numlib/mkl on the cluster. Most software modules for applications provide working example batch scripts.


Description Content
module load numlib/mkl
License Commercial. See EULA.
Citing n/a
Links Intel MKL Homepage | FFTW Homepage
Graphical Interface No


1 Description

The Fastest Fourier Transform in the West (FFTW) is a software library for computing discrete Fourier transforms in one or more dimensions, of arbitrary input size, and of both real and complex data (as well as of even/odd data, i.e. the discrete cosine/sine transforms or DCT/DST). FFTW was developed by Matteo Frigo and Steven G. Johnson at the Massachusetts Institute of Technology.

The Intel Math Kernel Library (Intel MKL) offers FFTW2 (for version 2.x) and FFTW3 (for version 3.x) interfaces to the Intel MKL Fast Fourier Transform and Trigonometric Transform functionality. These interfaces enable applications using FFTW to gain performance with Intel MKL without changing the application source code. Therefore, it is highly recommended to use Intel MKL instead of a separate FFTW installation.

2 FAQ

Copyright: KIZ (Ulm University)

Q: Why is there no FFTW module on the cluster?

A: MKL exhibits better performance than FFTS libraries (see Figure on the right). Therefore, we recommend to use MKL and do not offer a separate FFTW installation.

Q: Why does my code complain about argument of type "long double *" is incompatible with parameter of type "double *" ?

A: The interfaces do not support long double precision because Intel MKL FFT functions operate only on single- and double-precision floating point data types. For the very rare case that you need extended data types, please submit a support ticket at https://www.bwhpc.de/supportportal.

3 Useful links