00001
00002 #ifdef VCL_NEEDS_PRAGMA_INTERFACE
00003 #pragma implementation
00004 #endif
00005
00006
00007
00008
00009 #include "vnl_fft.h"
00010
00011 #include <vnl/algo/vnl_netlib.h>
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 }