fourier.h

Go to the documentation of this file.

//
//

/*
    Copyright (C) 2003- 2011, Hammersmith Imanet Ltd
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/
#ifndef __stir_numerics_stir_fourier_h__
#define __stir_numerics_stir_fourier_h__
#include "stir/VectorWithOffset.h"
#include "stir/Array_complex_numbers.h"
START_NAMESPACE_STIR

template <typename T>
void fourier(T& c, const int sign = 1);
// fourier(VectorWithOffset<elemT>& c, const int sign = 1);

template <typename T>
inline void
inverse_fourier(T& c, const int sign = 1)
{
  fourier(c, -sign);
#ifdef _MSC_VER
// disable warning about conversion
#  pragma warning(disable : 4244)
#endif
  c /= c.size_all();
#ifdef _MSC_VER
// disable warning about conversion
#  pragma warning(default : 4244)
#endif
}

template <typename T>
void fourier_1d(T& c, const int sign);

template <typename T>
inline void
inverse_fourier_1d(T& c, const int sign = 1)
{
  fourier_1d(c, -sign);
#ifdef _MSC_VER
// disable warning about conversion
#  pragma warning(disable : 4244)
#endif
  c /= c.size();
#ifdef _MSC_VER
#  pragma warning(default : 4244)
#endif
}

template <typename T>
Array<1, std::complex<T>> fourier_1d_for_real_data(const Array<1, T>& c, const int sign = 1);

template <typename T>
Array<1, T> inverse_fourier_1d_for_real_data(const Array<1, std::complex<T>>& c, const int sign = 1);

template <typename T>
Array<1, T> inverse_fourier_1d_for_real_data_corrupting_input(Array<1, std::complex<T>>& c, const int sign);

template <int num_dimensions, typename T>
Array<num_dimensions, std::complex<T>> fourier_for_real_data(const Array<num_dimensions, T>& c, const int sign = 1);

template <int num_dimensions, typename T>
Array<num_dimensions, T> inverse_fourier_for_real_data(const Array<num_dimensions, std::complex<T>>& c, const int sign = 1);

template <int num_dimensions, typename T>
Array<num_dimensions, T> inverse_fourier_for_real_data_corrupting_input(Array<num_dimensions, std::complex<T>>& c,
                                                                        const int sign = 1);

template <int num_dimensions, typename T>
Array<num_dimensions, std::complex<T>> pos_frequencies_to_all(const Array<num_dimensions, std::complex<T>>& c);
 
END_NAMESPACE_STIR
 
#endif