From eda5bc26f44ee9a6f83dcf8c91f17296d7fc509d Mon Sep 17 00:00:00 2001 From: Nao Pross Date: Mon, 12 Feb 2024 14:52:43 +0100 Subject: Move into version control --- .../include/armadillo_bits/translate_fftw3.hpp | 106 +++++++++++++++++++++ 1 file changed, 106 insertions(+) create mode 100644 src/armadillo/include/armadillo_bits/translate_fftw3.hpp (limited to 'src/armadillo/include/armadillo_bits/translate_fftw3.hpp') diff --git a/src/armadillo/include/armadillo_bits/translate_fftw3.hpp b/src/armadillo/include/armadillo_bits/translate_fftw3.hpp new file mode 100644 index 0000000..1edd727 --- /dev/null +++ b/src/armadillo/include/armadillo_bits/translate_fftw3.hpp @@ -0,0 +1,106 @@ +// SPDX-License-Identifier: Apache-2.0 +// +// Copyright 2008-2016 Conrad Sanderson (http://conradsanderson.id.au) +// Copyright 2008-2016 National ICT Australia (NICTA) +// +// Licensed under the Apache License, Version 2.0 (the "License"); +// you may not use this file except in compliance with the License. +// You may obtain a copy of the License at +// http://www.apache.org/licenses/LICENSE-2.0 +// +// Unless required by applicable law or agreed to in writing, software +// distributed under the License is distributed on an "AS IS" BASIS, +// WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied. +// See the License for the specific language governing permissions and +// limitations under the License. +// ------------------------------------------------------------------------ + + +#if defined(ARMA_USE_FFTW3) + + +namespace fftw3 + { + template + arma_inline + void_ptr + plan_dft_1d(int N, eT* input, eT* output, int fftw3_sign, unsigned int fftw3_flags) + { + arma_type_check((is_cx::value == false)); + + if(is_cx_float::value) + { + return fftwf_plan_dft_1d(N, (cx_float*)input, (cx_float*)output, fftw3_sign, fftw3_flags); + } + else + if(is_cx_double::value) + { + return fftw_plan_dft_1d(N, (cx_double*)input, (cx_double*)output, fftw3_sign, fftw3_flags); + } + + return nullptr; + } + + + + template + arma_inline + void + execute(void_ptr plan) + { + arma_type_check((is_cx::value == false)); + + if(is_cx_float::value) + { + fftwf_execute(plan); + } + else + if(is_cx_double::value) + { + fftw_execute(plan); + } + } + + + + template + arma_inline + void + destroy_plan(void_ptr plan) + { + arma_type_check((is_cx::value == false)); + + if(is_cx_float::value) + { + fftwf_destroy_plan(plan); + } + else + if(is_cx_double::value) + { + fftw_destroy_plan(plan); + } + } + + + + template + arma_inline + void + cleanup() + { + arma_type_check((is_cx::value == false)); + + if(is_cx_float::value) + { + fftwf_cleanup(); + } + else + if(is_cx_double::value) + { + fftw_cleanup(); + } + } + } + + +#endif -- cgit v1.2.1