core/vnl/algo/vnl_fft.cxx
Go to the documentation of this file.
00001 // This is core/vnl/algo/vnl_fft.cxx
00002 #ifdef VCL_NEEDS_PRAGMA_INTERFACE
00003 #pragma implementation
00004 #endif
00005 //:
00006 // \file
00007 // \author fsm
00008 
00009 #include "vnl_fft.h"
00010 
00011 #include <vnl/algo/vnl_netlib.h> // dgpfa_()
00012 
00013 void vnl_fft_setgpfa(float *triggs, long size, long pqr[3], long *info)
00014 {
00015   v3p_netlib_setgpfa_(triggs, &size, pqr, info);
00016 }
00017 
00018 void vnl_fft_setgpfa(double *triggs, long size, long pqr[3], long *info)
00019 {
00020   v3p_netlib_setdgpfa_(triggs, &size, pqr, info);
00021 }
00022 
00023 //----------------------------------------------------------------------
00024 
00025 void vnl_fft_gpfa(float  *a, float  *b, float const  *triggs,
00026                   long inc, long jump, long n,
00027                   long lot, long isign, long const pqr[3], long *info)
00028 {
00029   v3p_netlib_gpfa_(a, b, triggs, &inc, &jump, &n, &lot, &isign, pqr);
00030   *info = 0;
00031 }
00032 
00033 void vnl_fft_gpfa(double *a, double *b, double const *triggs,
00034                   long inc, long jump, long n,
00035                   long lot, long isign, long const pqr[3], long *info)
00036 {
00037   v3p_netlib_dgpfa_(a, b, triggs, &inc, &jump, &n, &lot, &isign, pqr);
00038   *info = 0;
00039 }