FTransform.h

//-*-C++-*-
/***************************************************************************
 *
 *   Copyright (C) 2006 by Willem van Straten
 *   Licensed under the Academic Free License version 2.1
 *
 ***************************************************************************/
 
#ifndef __FTransform_h_
#define __FTransform_h_
 
#include "Reference.h"
 
#include <string>
#include <vector>

namespace FTransform {
 
  /* ////////////////////////////////////////////////////////////////////
     
  One-dimensional FFT library interface
  


  void frc1d (size_t nfft, float* into, const float* from);


  void bcr1d (size_t nfft, float* into, const float* from);

  void fcc1d (size_t nfft, float* into, const float* from);

  void bcc1d (size_t nfft, float* into, const float* from);

  extern bool optimize;

  extern bool simd;

  extern unsigned nthread;

  std::string get_library ();
  
  void set_library (const std::string& name);

  unsigned get_num_libraries ();

  std::string get_library_name (unsigned i);

  void clean_plans();

  enum normalization { normalized, unnormalized };

  enum direction { forward = 0x10, backward = 0x20 };

  enum dimension { real = 0x01, analytic = 0x02 };

  enum type {
    frc = forward  | real,
    bcr = backward | real,
    fcc = forward  | analytic,
    bcc = backward | analytic
  };

  normalization get_norm ();

  double get_scale (size_t nfft, type t);

  void shift (unsigned npts, float* arr, double shift);

  void derivative (unsigned npts, float* data);

  int inplace_frc1d(size_t nfft, float* srcdest);
  int inplace_fcc1d(size_t nfft, float* srcdest);
  int inplace_bcc1d(size_t nfft, float* srcdest);
  int inplace_bcr1d(size_t nfft, float* srcdest);
 
 
  /* ////////////////////////////////////////////////////////////////////
     
  Two-dimensional FFT library interface
  

  void fcc2d (size_t nx, size_t ny, float* into, const float* from);

  void bcc2d (size_t nx, size_t ny, float* into, const float* from);
 
}
 
#endif