diff --git a/.gitattributes b/.gitattributes
new file mode 100644
index 0000000..7b28786
--- /dev/null
+++ b/.gitattributes
@@ -0,0 +1,13 @@
+csv/xpa2.csv filter=lfs diff=lfs merge=lfs -text
+csv/acars.csv filter=lfs diff=lfs merge=lfs -text
+csv/ads-b.csv filter=lfs diff=lfs merge=lfs -text
+csv/flex.csv filter=lfs diff=lfs merge=lfs -text
+csv/jorn.csv filter=lfs diff=lfs merge=lfs -text
+csv/lofar-20m-x.csv filter=lfs diff=lfs merge=lfs -text
+csv/test-reverse.csv filter=lfs diff=lfs merge=lfs -text
+csv/xpa1.csv filter=lfs diff=lfs merge=lfs -text
+wav/lofar-20m-x.wav filter=lfs diff=lfs merge=lfs -text
+wav/xpa1-8000.wav filter=lfs diff=lfs merge=lfs -text
+wav/acars1-iq-48000.wav filter=lfs diff=lfs merge=lfs -text
+wav/dvb-t-iq-5333333.wav filter=lfs diff=lfs merge=lfs -text
+wav/flex-iq-39062.wav filter=lfs diff=lfs merge=lfs -text
diff --git a/.gitlab-ci.yml b/.gitlab-ci.yml
index 7f1d1b8..ba2e807 100644
--- a/.gitlab-ci.yml
+++ b/.gitlab-ci.yml
@@ -1,23 +1,77 @@
+image: $CI_REGISTRY_IMAGE/ci-build-runner:$CI_COMMIT_REF_SLUG
+variables:
+ GIT_SUBMODULE_STRATEGY: recursive
+
stages:
+ - prepare
- build-lib
- build
- test
+ - release
+
+include:
+ - template: Security/SAST.gitlab-ci.yml
+ - template: Security/Dependency-Scanning.gitlab-ci.yml
+ - template: Security/Secret-Detection.gitlab-ci.yml
+
+default:
+ before_script:
+ - mkdir -p build
+ - cd build || exit
+ cache:
+ paths:
+ - "build/*"
+
+trigger-prepare:
+ stage: prepare
+ trigger:
+ strategy: depend
+ include: .prepare.gitlab-ci.yml
+
+prepare-release:
+ stage: prepare
+ image: alpine:latest
+ rules:
+ - if: '$CI_COMMIT_TAG && $CI_COMMIT_REF_PROTECTED == "true"'
+ before_script:
+ - apk add curl jq
+ script:
+ - 'curl -H "PRIVATE-TOKEN: $CI_API_TOKEN" "$CI_API_V4_URL/projects/$CI_PROJECT_ID/repository/changelog?version=$CI_COMMIT_TAG" | jq -r .notes > release_notes.md'
+ - cat release_notes.md
+ artifacts:
+ paths:
+ - release_notes.md
+
build-library:
stage: build-lib
+ needs:
+ - trigger-prepare
script:
- - cd $CI_PROJECT_DIR
- - mkdir -p build
- - cd build || exit
- cmake ..
- make oCLFFTsup -j $(nproc)
artifacts:
paths:
- build/lib/liboCLFFTsup.so
- cache:
+
+build-lookup:
+ stage: build-lib
+ needs:
+ - trigger-prepare
+ script:
+ - cmake ..
+ - make play-lookup play-lookup-fp32
+ - ./playground/play-lookup > ../oclfft/ard-ocl/src/lookup.cl
+ - ./playground/play-lookup > ../oclfft/bit-ocl/src/lookup.cl
+ - ./playground/play-lookup > ../oclfft/flk-ocl/src/lookup.cl
+ - ./playground/play-lookup-fp32 > ../oclfft/flk-ocl-fp32/src/lookup.cl
+ artifacts:
paths:
- - "build/*"
+ - oclfft/**/lookup.cl
+
build-documents:
stage: build
+ needs:
+ - trigger-prepare
script:
- cd $CI_PROJECT_DIR/proposal || exit
- texi2pdf proposal.tex
@@ -26,40 +80,70 @@ build-documents:
artifacts:
paths:
- "*/*.pdf"
+
build-documentation:
stage: build
+ needs:
+ - trigger-prepare
script:
- - cd $CI_PROJECT_DIR
- - mkdir -p build
- - cd build || exit
- cmake -DENABLE_DOXYGEN=on ..
- make doxygen
artifacts:
paths:
- "docs/*"
+
build-binaries:
stage: build
+ needs:
+ - trigger-prepare
+ - build-lookup
+ - build-library
script:
- - cd $CI_PROJECT_DIR || exit
- - mkdir -p build
- - cd build || exit
- cmake -DENABLE_DOXYGEN=off ..
- - make -j $(nproc)
+ - make ard-ocl
+ - make ard-omp
+ - make ard-seq
+ - make bit-ocl
+ - make dft-seq
+ - make dft-seq-fp32
+ - make flk-ocl
+ - make flk-ocl-fp32
+ - make numerical-analysis
+ - make play-objcopy
+ - make play-example
+ - make play-reversegen
artifacts:
paths:
- "*/**/play-objcopy"
- "*/**/play-example"
- - "*/**/testexample"
- cache:
- paths:
- - "build/*"
+
test:
stage: test
+ needs:
+ - trigger-prepare
+ - build-lookup
+ - build-library
+ - build-binaries
+ variables:
+ CTEST_OUTPUT_ON_FAILURE: 1
script:
- - cd $CI_PROJECT_DIR || exit
- - mkdir -p build
- - cd build || exit
- - make test
- cache:
- paths:
- - "build/*"
\ No newline at end of file
+ - make run_tests
+
+release_job:
+ stage: release
+ image: registry.gitlab.com/gitlab-org/release-cli:latest
+ environment: gitlab-releases
+ needs:
+ - job: prepare-release
+ artifacts: true
+ rules:
+ - if: '$CI_COMMIT_TAG && $CI_COMMIT_REF_PROTECTED == "true"'
+ before_script:
+ - echo "Creating release ${CI_COMMIT_TAG}"
+ script:
+ - cat release_notes.md
+ release:
+ name: '${CI_COMMIT_TAG}'
+ description: release_notes.md
+ tag_name: '${CI_COMMIT_TAG}'
+ ref: '$CI_COMMIT_SHA'
diff --git a/.prepare.gitlab-ci.yml b/.prepare.gitlab-ci.yml
new file mode 100644
index 0000000..87b71e3
--- /dev/null
+++ b/.prepare.gitlab-ci.yml
@@ -0,0 +1,21 @@
+stages:
+ - build
+
+build_ci_runner_image:
+ stage: build
+ image: docker:latest
+ tags:
+ - dind
+ before_script:
+ - apk add yq
+ script:
+ - docker login -u $CI_REGISTRY_USER -p $CI_REGISTRY_PASSWORD $CI_REGISTRY
+ - docker context create my-builder
+ - docker buildx create my-builder --driver docker-container --use
+ - docker buildx build --push -t $CI_REGISTRY_IMAGE/ci-build-runner:$CI_COMMIT_REF_SLUG --cache-to type=registry,ref=$CI_REGISTRY_IMAGE/ci-build-runner-cache:$CI_COMMIT_REF_SLUG,mode=max --cache-from type=registry,ref=$CI_REGISTRY_IMAGE/ci-build-runner-cache:$CI_COMMIT_REF_SLUG docker/
+
+
+
+workflow:
+ rules:
+ - when: always
diff --git a/CMakeLists.txt b/CMakeLists.txt
index 9150bd5..d3fb2ff 100644
--- a/CMakeLists.txt
+++ b/CMakeLists.txt
@@ -14,7 +14,7 @@
# You should have received a copy of the GNU General Public License
# along with this program. If not, see .
-CMAKE_MINIMUM_REQUIRED(VERSION 3.0)
+CMAKE_MINIMUM_REQUIRED(VERSION 3.10)
PROJECT(oCLFFT)
# 3.0 is minimum version as it is the first version to include FindOpenCL
@@ -43,16 +43,16 @@ use_cxx_warning_pedantic()
# Configurable options and defaults
OPTION(ENABLE_TESTS "Enable unit tests" ON)
OPTION(ENABLE_CODECOV "Measure code coverage" OFF)
-OPTION(ENABLE_DOXYGEN "Generate documentation" ON)
+OPTION(ENABLE_DOXYGEN "Generate documentation" OFF)
OPTION(ENABLE_CI "Enable continuous integration" OFF)
-include_directories(${CMAKE_SOURCE_DIR}/lib/include)
+#include_directories(${CMAKE_SOURCE_DIR}/lib/include)
# Required packages
set(CMAKE_THREAD_PREFER_PTHREAD TRUE)
set(THREADS_PREFER_PTHREAD_FLAG TRUE)
find_package(Threads REQUIRED)
-find_package(OpenCL 2.0 REQUIRED)
+find_package(OpenCL 2.2 REQUIRED)
# Include Commonly required headers
# include_directories(${ZeroMQ_INCLUDE_DIR})
@@ -89,8 +89,8 @@ if(CMAKE_BUILD_TYPE STREQUAL "Debug")
set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -g3 -O0")
elseif(CMAKE_BUILD_TYPE STREQUAL "Release")
message("${PRJ_PRX}: Enabling extensive optimizations")
- set(CMAKE_CXX_FLAGS
- "${CMAKE_CXX_FLAGS} -Ofast -ffast-math -ftree-vectorize")
+# set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -march=native -Ofast -ffast-math -ftree-vectorize")
+ set(CMAKE_CXX_FLAGS "${CMAKE_CXX_FLAGS} -march=native -O3 -ftree-vectorize")
endif()
# Commonly shared library
diff --git a/README.md b/README.md
index 648509b..8d245b8 100644
--- a/README.md
+++ b/README.md
@@ -5,6 +5,10 @@ Project for University of Amsterdam (UvA) performance engineering course,
boilerplate CMake files are from
[a previous project](https://github.com/hexoxide/O2-Balancer2)
+
+
+
+
### Directory structure
* cmake - small cmake snippets to enable various features
@@ -13,6 +17,7 @@ boilerplate CMake files are from
* lib - support library with functions and definitions
* measurements - performance evaluation scripts as well as measured values
* midterm - LaTeX files for the presentation at the first half of the course
+* numerical-analysis - Results and measurements to determine numerical error
* playground - small toy examples or other experiments
* [python](python/README.md) - python scripts to aid in visualization or measurements
* proposal - LaTeX source for project proposal
@@ -23,31 +28,44 @@ boilerplate CMake files are from
### Implementations
-| Executable | Hardware Target | Base Algorithm | Properties |
-|------------|-------------------|----------------|-------------------------------|
-| ard-seq | CPU Single Thread | Cooley-Tukey | Bit-reversal in-place radix-2 |
-| dft-seq | CPU Single Thread | DFT | |
-| ard-omp | CPU OpenMP | Cooley-Tukey | Bit-reversal in-place radix-2 |
-| ard-ocl | GPU | Cooley-Tukey | Slow bit-reversal, fft lookup |
-| bit-ocl | GPU | Cooley-Tukey | 2D bit-reversal, fft lookup |
-| flk-ocl | GPU | Cooley-Tukey | 2D bit-reversal, fast lookup |
+| Executable | Hardware Target | Base Algorithm | Precision | Properties |
+|--------------------|-------------------|----------------|--------------|------------------------------------------------------------|
+| ard-seq | CPU Single Thread | Cooley-Tukey | 64bit double | Bit-reversal in-place radix-2 |
+| dft-seq | CPU Single Thread | DFT | 64bit double | |
+| dft-seq-fp32 | CPU Single Thread | DFT | 32bit float | |
+| ard-omp | CPU OpenMP | Cooley-Tukey | 64bit double | Bit-reversal in-place radix-2 |
+| ard-ocl | GPU | Cooley-Tukey | 64bit double | Slow bit-reversal, fft lookup |
+| bit-ocl | GPU | Cooley-Tukey | 64bit double | 2D bit-reversal, fft lookup |
+| flk-ocl | GPU | Cooley-Tukey | 64bit double | 2D bit-reversal, fast lookup |
+| flk-ocl-fp32 | GPU | Cooley-Tukey | 32bit float | 2D bit-reversal, fast lookup, optimized local groups |
+| numerical-analysis | CPU/GPU | N.A | arbitrary | Arbitrary precision MPFR FFR, compute L1,L2 and Linf norms |
#### Dependencies
-* cmake 3.0 or higher
-* cmake 3.9 or higher (ard-omp)
-* doxygen (documentation)
+* cmake 3.10 or higher
* boost 1.32 or higher
* boost 1.53 or higher (unit tests)
+* OpenCL 2.0 or higher
+* OpenCL C++ headers
+* MPFR
+* MPFR C++, MPREAL (rms-error)
+* doxygen (documentation)
* python 3.x (python)
* virtualenv (python)
#### Setup
+Initialize lfs:
+
+```bash
+git lfs install
+git lfs checkout
+```
+
Python environment:
```bash
-virtualenv -p python3 python/e
+virtualenv -p python3 python/
cd python
source bin/activate
pip install -r requirements.txt
@@ -58,11 +76,14 @@ Generating lookup tables:
```bash
mkdir build
cd build
-cmake ..
+cmake -DCMAKE_BUILD_TYPE=Release ..
make play-lookup
./playground/play-lookup > ../oclfft/ard-ocl/src/lookup.cl
./playground/play-lookup > ../oclfft/bit-ocl/src/lookup.cl
./playground/play-lookup > ../oclfft/flk-ocl/src/lookup.cl
+make play-lookup-fp32
+./playground/play-lookup-fp32 > ../oclfft/flk-ocl-fp32/src/lookup.cl
+make
```
#### Licensing
@@ -81,13 +102,24 @@ A few examples of files licensed under different authors:
* [Emebedding source files into binaries](https://www.linuxjournal.com/content/embedding-file-executable-aka-hello-world-version-5967)
* [OpenCL, SyCL and SPIR-V progress 2016](https://www.youtube.com/watch?v=TYp1d6yzHUQ)
+* [RDNA whitepaper](https://www.amd.com/system/files/documents/rdna-whitepaper.pdf)
+* Arbitrary Precision Arithmetic (FFTs)
+ * [Python Arbitrary Precision FFT](https://github.com/dbstein/apfft/tree/master)
+ * [C++ Arbitrary Precision Arithmetic](https://github.com/anxieuse/Arbitrary-precision-arithmetic/tree/main)
+ * [BenchFFT](https://fftw.org/benchfft/)
+ * [RRMSE](https://stats.stackexchange.com/questions/413209/is-there-something-like-a-root-mean-square-relative-error-rmsre-or-what-is-t#413217)
#### Snippets
```bash
-../cmake-build-debug/oclfft/bit-ocl/bit-ocl -f ../csv/ads-b.csv -s 262144 -o real | python evaluate-error.py
+cd python
+../cmake-build-debug/oclfft/flk-ocl-fp32/flk-ocl-fp32 -w 256 -f ../csv/lofar-20m-x.csv -s 4194304 -o real | python evaluate-error.py
```
```bash
- /opt/rocm/bin/rocprof --obj-tracking on --hsa-trace binary
-```
\ No newline at end of file
+/opt/rocm/bin/rocprofv3 --obj-tracking on --hsa-trace binary
+rocprofv2 --sys-trace --kernel-trace --plugin perfetto -o results.json binary
+rocprofv3 -r -s --output-format pftrace -i ../counters.txt -- oclfft/flk-ocl-fp32/flk-ocl-fp32 -w 256 -f ../csv/lofar-20m-x.csv -s 4194304 -o TIME
+```
+
+Open pftrace in https://ui.perfetto.dev
\ No newline at end of file
diff --git a/cmake/FindFFTW.cmake b/cmake/FindFFTW.cmake
index 68d97fc..f21b8d8 100644
--- a/cmake/FindFFTW.cmake
+++ b/cmake/FindFFTW.cmake
@@ -1,46 +1,419 @@
-# Copyright jedbrown, ptooley, mdave ,johnfettig, gruenich, ckhroulev, ihnorton
-# All rights reserved.
+# - Find the FFTW library
#
-# Redistribution and use in source and binary forms, with or without modification,
-# are permitted provided that the following conditions are met:
+# Original version of this file:
+# Copyright (c) 2015, Wenzel Jakob
+# https://github.com/wjakob/layerlab/blob/master/cmake/FindFFTW.cmake, commit 4d58bfdc28891b4f9373dfe46239dda5a0b561c6
+# Modifications:
+# Copyright (c) 2017, Patrick Bos
#
-# * Redistributions of source code must retain the above copyright notice, this
-# list of conditions and the following disclaimer.
+# Usage:
+# find_package(FFTW [REQUIRED] [QUIET] [COMPONENTS component1 ... componentX] )
#
-# * Redistributions in binary form must reproduce the above copyright notice, this
-# list of conditions and the following disclaimer in the documentation and/or
-# other materials provided with the distribution.
+# It sets the following variables:
+# FFTW_FOUND ... true if fftw is found on the system
+# FFTW_[component]_LIB_FOUND ... true if the component is found on the system (see components below)
+# FFTW_LIBRARIES ... full paths to all found fftw libraries
+# FFTW_[component]_LIB ... full path to one of the components (see below)
+# FFTW_INCLUDE_DIRS ... fftw include directory paths
#
-# THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS" AND
-# ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE IMPLIED
-# WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
-# DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE FOR
-# ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES
-# (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR SERVICES;
-# LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON
-# ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT
-# (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE OF THIS
-# SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
-
-# - Find FFTW
-# Find the native FFTW includes and library
+# The following variables will be checked by the function
+# FFTW_USE_STATIC_LIBS ... if true, only static libraries are found, otherwise both static and shared.
+# FFTW_ROOT ... if set, the libraries are exclusively searched
+# under this path
#
-# FFTW_INCLUDES - where to find fftw3.h
-# FFTW_LIBRARIES - List of libraries when using FFTW.
-# FFTW_FOUND - True if FFTW found.
+# This package supports the following components:
+# FLOAT_LIB
+# DOUBLE_LIB
+# LONGDOUBLE_LIB
+# FLOAT_THREADS_LIB
+# DOUBLE_THREADS_LIB
+# LONGDOUBLE_THREADS_LIB
+# FLOAT_OPENMP_LIB
+# DOUBLE_OPENMP_LIB
+# LONGDOUBLE_OPENMP_LIB
+#
+
+# TODO (maybe): extend with ExternalProject download + build option
+# TODO: put on conda-forge
+
+
+if( NOT FFTW_ROOT AND DEFINED ENV{FFTWDIR} )
+ set( FFTW_ROOT $ENV{FFTWDIR} )
+endif()
+
+# Check if we can use PkgConfig
+find_package(PkgConfig)
+
+#Determine from PKG
+if( PKG_CONFIG_FOUND AND NOT FFTW_ROOT )
+ pkg_check_modules( PKG_FFTW QUIET "fftw3" )
+endif()
+
+#Check whether to search static or dynamic libs
+set( CMAKE_FIND_LIBRARY_SUFFIXES_SAV ${CMAKE_FIND_LIBRARY_SUFFIXES} )
+
+if( ${FFTW_USE_STATIC_LIBS} )
+ set( CMAKE_FIND_LIBRARY_SUFFIXES ${CMAKE_STATIC_LIBRARY_SUFFIX} )
+else()
+ set( CMAKE_FIND_LIBRARY_SUFFIXES ${CMAKE_FIND_LIBRARY_SUFFIXES_SAV} )
+endif()
+
+if( FFTW_ROOT )
+ # find libs
+
+ find_library(
+ FFTW_DOUBLE_LIB
+ NAMES "fftw3" libfftw3-3
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_DOUBLE_THREADS_LIB
+ NAMES "fftw3_threads"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_DOUBLE_OPENMP_LIB
+ NAMES "fftw3_omp"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_DOUBLE_MPI_LIB
+ NAMES "fftw3_mpi"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_FLOAT_LIB
+ NAMES "fftw3f" libfftw3f-3
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_FLOAT_THREADS_LIB
+ NAMES "fftw3f_threads"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_FLOAT_OPENMP_LIB
+ NAMES "fftw3f_omp"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_FLOAT_MPI_LIB
+ NAMES "fftw3f_mpi"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_LONGDOUBLE_LIB
+ NAMES "fftw3l" libfftw3l-3
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_LONGDOUBLE_THREADS_LIB
+ NAMES "fftw3l_threads"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_LONGDOUBLE_OPENMP_LIB
+ NAMES "fftw3l_omp"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ find_library(
+ FFTW_LONGDOUBLE_MPI_LIB
+ NAMES "fftw3l_mpi"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "lib" "lib64"
+ NO_DEFAULT_PATH
+ )
+
+ #find includes
+ find_path(FFTW_INCLUDE_DIRS
+ NAMES "fftw3.h"
+ PATHS ${FFTW_ROOT}
+ PATH_SUFFIXES "include"
+ NO_DEFAULT_PATH
+ )
+
+else()
+
+ find_library(
+ FFTW_DOUBLE_LIB
+ NAMES "fftw3"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_DOUBLE_THREADS_LIB
+ NAMES "fftw3_threads"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_DOUBLE_OPENMP_LIB
+ NAMES "fftw3_omp"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_DOUBLE_MPI_LIB
+ NAMES "fftw3_mpi"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_FLOAT_LIB
+ NAMES "fftw3f"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_FLOAT_THREADS_LIB
+ NAMES "fftw3f_threads"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_FLOAT_OPENMP_LIB
+ NAMES "fftw3f_omp"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_FLOAT_MPI_LIB
+ NAMES "fftw3f_mpi"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_LONGDOUBLE_LIB
+ NAMES "fftw3l"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(
+ FFTW_LONGDOUBLE_THREADS_LIB
+ NAMES "fftw3l_threads"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(FFTW_LONGDOUBLE_OPENMP_LIB
+ NAMES "fftw3l_omp"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_library(FFTW_LONGDOUBLE_MPI_LIB
+ NAMES "fftw3l_mpi"
+ PATHS ${PKG_FFTW_LIBRARY_DIRS} ${LIB_INSTALL_DIR}
+ )
+
+ find_path(FFTW_INCLUDE_DIRS
+ NAMES "fftw3.h"
+ PATHS ${PKG_FFTW_INCLUDE_DIRS} ${INCLUDE_INSTALL_DIR}
+ )
+
+endif( FFTW_ROOT )
+
+#--------------------------------------- components
+
+if (FFTW_DOUBLE_LIB)
+ set(FFTW_DOUBLE_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_DOUBLE_LIB})
+ add_library(FFTW::Double INTERFACE IMPORTED)
+ set_target_properties(FFTW::Double
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_DOUBLE_LIB}"
+ )
+else()
+ set(FFTW_DOUBLE_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_FLOAT_LIB)
+ set(FFTW_FLOAT_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_FLOAT_LIB})
+ add_library(FFTW::Float INTERFACE IMPORTED)
+ set_target_properties(FFTW::Float
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_FLOAT_LIB}"
+ )
+else()
+ set(FFTW_FLOAT_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_LONGDOUBLE_LIB)
+ set(FFTW_LONGDOUBLE_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_LONGDOUBLE_LIB})
+ add_library(FFTW::LongDouble INTERFACE IMPORTED)
+ set_target_properties(FFTW::LongDouble
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_LONGDOUBLE_LIB}"
+ )
+else()
+ set(FFTW_LONGDOUBLE_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_DOUBLE_THREADS_LIB)
+ set(FFTW_DOUBLE_THREADS_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_DOUBLE_THREADS_LIB})
+ add_library(FFTW::DoubleThreads INTERFACE IMPORTED)
+ set_target_properties(FFTW::DoubleThreads
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_DOUBLE_THREADS_LIB}"
+ )
+else()
+ set(FFTW_DOUBLE_THREADS_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_FLOAT_THREADS_LIB)
+ set(FFTW_FLOAT_THREADS_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_FLOAT_THREADS_LIB})
+ add_library(FFTW::FloatThreads INTERFACE IMPORTED)
+ set_target_properties(FFTW::FloatThreads
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_FLOAT_THREADS_LIB}"
+ )
+else()
+ set(FFTW_FLOAT_THREADS_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_LONGDOUBLE_THREADS_LIB)
+ set(FFTW_LONGDOUBLE_THREADS_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_LONGDOUBLE_THREADS_LIB})
+ add_library(FFTW::LongDoubleThreads INTERFACE IMPORTED)
+ set_target_properties(FFTW::LongDoubleThreads
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_LONGDOUBLE_THREADS_LIB}"
+ )
+else()
+ set(FFTW_LONGDOUBLE_THREADS_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_DOUBLE_OPENMP_LIB)
+ set(FFTW_DOUBLE_OPENMP_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_DOUBLE_OPENMP_LIB})
+ add_library(FFTW::DoubleOpenMP INTERFACE IMPORTED)
+ set_target_properties(FFTW::DoubleOpenMP
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_DOUBLE_OPENMP_LIB}"
+ )
+else()
+ set(FFTW_DOUBLE_OPENMP_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_FLOAT_OPENMP_LIB)
+ set(FFTW_FLOAT_OPENMP_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_FLOAT_OPENMP_LIB})
+ add_library(FFTW::FloatOpenMP INTERFACE IMPORTED)
+ set_target_properties(FFTW::FloatOpenMP
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_FLOAT_OPENMP_LIB}"
+ )
+else()
+ set(FFTW_FLOAT_OPENMP_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_LONGDOUBLE_OPENMP_LIB)
+ set(FFTW_LONGDOUBLE_OPENMP_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_LONGDOUBLE_OPENMP_LIB})
+ add_library(FFTW::LongDoubleOpenMP INTERFACE IMPORTED)
+ set_target_properties(FFTW::LongDoubleOpenMP
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_LONGDOUBLE_OPENMP_LIB}"
+ )
+else()
+ set(FFTW_LONGDOUBLE_OPENMP_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_DOUBLE_MPI_LIB)
+ set(FFTW_DOUBLE_MPI_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_DOUBLE_MPI_LIB})
+ add_library(FFTW::DoubleMPI INTERFACE IMPORTED)
+ set_target_properties(FFTW::DoubleMPI
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_DOUBLE_MPI_LIB}"
+ )
+else()
+ set(FFTW_DOUBLE_MPI_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_FLOAT_MPI_LIB)
+ set(FFTW_FLOAT_MPI_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_FLOAT_MPI_LIB})
+ add_library(FFTW::FloatMPI INTERFACE IMPORTED)
+ set_target_properties(FFTW::FloatMPI
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_FLOAT_MPI_LIB}"
+ )
+else()
+ set(FFTW_FLOAT_MPI_LIB_FOUND FALSE)
+endif()
+
+if (FFTW_LONGDOUBLE_MPI_LIB)
+ set(FFTW_LONGDOUBLE_MPI_LIB_FOUND TRUE)
+ set(FFTW_LIBRARIES ${FFTW_LIBRARIES} ${FFTW_LONGDOUBLE_MPI_LIB})
+ add_library(FFTW::LongDoubleMPI INTERFACE IMPORTED)
+ set_target_properties(FFTW::LongDoubleMPI
+ PROPERTIES INTERFACE_INCLUDE_DIRECTORIES "${FFTW_INCLUDE_DIRS}"
+ INTERFACE_LINK_LIBRARIES "${FFTW_LONGDOUBLE_MPI_LIB}"
+ )
+else()
+ set(FFTW_LONGDOUBLE_MPI_LIB_FOUND FALSE)
+endif()
-if (FFTW_INCLUDES)
- # Already in cache, be silent
- set (FFTW_FIND_QUIETLY TRUE)
-endif (FFTW_INCLUDES)
+#--------------------------------------- end components
-find_path (FFTW_INCLUDES fftw3.h)
+set( CMAKE_FIND_LIBRARY_SUFFIXES ${CMAKE_FIND_LIBRARY_SUFFIXES_SAV} )
-find_library (FFTW_LIBRARIES NAMES fftw3)
+include(FindPackageHandleStandardArgs)
-# handle the QUIETLY and REQUIRED arguments and set FFTW_FOUND to TRUE if
-# all listed variables are TRUE
-include (FindPackageHandleStandardArgs)
-find_package_handle_standard_args (FFTW DEFAULT_MSG FFTW_LIBRARIES FFTW_INCLUDES)
+find_package_handle_standard_args(FFTW
+ REQUIRED_VARS FFTW_INCLUDE_DIRS
+ HANDLE_COMPONENTS
+)
-mark_as_advanced (FFTW_LIBRARIES FFTW_INCLUDES)
\ No newline at end of file
+mark_as_advanced(
+ FFTW_INCLUDE_DIRS
+ FFTW_LIBRARIES
+ FFTW_FLOAT_LIB
+ FFTW_DOUBLE_LIB
+ FFTW_LONGDOUBLE_LIB
+ FFTW_FLOAT_THREADS_LIB
+ FFTW_DOUBLE_THREADS_LIB
+ FFTW_LONGDOUBLE_THREADS_LIB
+ FFTW_FLOAT_OPENMP_LIB
+ FFTW_DOUBLE_OPENMP_LIB
+ FFTW_LONGDOUBLE_OPENMP_LIB
+ FFTW_FLOAT_MPI_LIB
+ FFTW_DOUBLE_MPI_LIB
+ FFTW_LONGDOUBLE_MPI_LIB
+)
diff --git a/cmake/FindMPFR.cmake b/cmake/FindMPFR.cmake
new file mode 100644
index 0000000..d2dedea
--- /dev/null
+++ b/cmake/FindMPFR.cmake
@@ -0,0 +1,83 @@
+# Try to find the MPFR library
+# See http://www.mpfr.org/
+#
+# This module supports requiring a minimum version, e.g. you can do
+# find_package(MPFR 2.3.0)
+# to require version 2.3.0 to newer of MPFR.
+#
+# Once done this will define
+#
+# MPFR_FOUND - system has MPFR lib with correct version
+# MPFR_INCLUDES - the MPFR include directory
+# MPFR_LIBRARIES - the MPFR library
+# MPFR_VERSION - MPFR version
+
+# Copyright (c) 2006, 2007 Montel Laurent,
+# Copyright (c) 2008, 2009 Gael Guennebaud,
+# Copyright (c) 2010 Jitse Niesen,