3.0 API Reference

util.h

Go to the documentation of this file.

#ifndef utils_h_inc
#define utils_h_inc

#include <complex.h>

/*###########################################################################*/
/*###########################################################################*/
/*###########################################################################*/

#define NFFT_SWAP_complex(x,y) {double complex* temp; temp=(x); (x)=(y); (y)=temp;}

#define NFFT_SWAP_double(x,y) {double* temp; temp=(x); (x)=(y); (y)=temp;}

#define PI 3.1415926535897932384

#define NFFT_MAX(a,b) ((a)>(b)? (a) : (b))

#define NFFT_MIN(a,b) ((a)<(b)? (a) : (b))

/* ######################################################################### */
/* ########## Little helpers ############################################### */
/* ######################################################################### */

double nfft_second();

int nfft_total_used_memory();

int nfft_ld(int m);

int nfft_int_2_pow(int a);

int nfft_next_power_of_2(int N);

void nfft_next_power_of_2_exp(int N, int *N2, int *t);

/* ######################################################################### */
/* ########## Window function related ###################################### */
/* ######################################################################### */

double nfft_sinc(double x);

double nfft_bspline_old(int k,double x,double *A);

double nfft_bspline(int k, double x, double *scratch);

double nfft_i0(double x);

/* ######################################################################### */
/* ########## Vector routines ############################################## */
/* ######################################################################### */

int nfft_prod_int(int *vec, int d);

int nfct_prod_int(int *vec, int d);

int nfst_prod_minus_a_int(int *vec, int a, int d);

int nfft_plain_loop(int *idx,int *N,int d);

double nfft_prod_real(double *vec,int d);

double nfft_dot_complex(double complex* x, int n);

double nfft_dot_double( double*  x, int n);

double nfft_dot_w_complex(double complex* x, double* w, int n);

double nfft_dot_w_double( double*  x, double* w, int n);

double nfft_dot_w_w2_complex(double complex* x, double* w, double* w2, int n);

double nfft_dot_w2_complex(double complex* x, double* w2, int n);

void nfft_cp_complex(double complex* x, double complex* y, int n);

void nfft_cp_double( double*  x, double*  y, int n);

void nfft_cp_a_complex(double complex* x, double a, double complex* y, int n);

void nfft_cp_w_complex(double complex* x, double* w, double complex* y, int n);

void nfft_cp_w_double( double*  x, double* w, double*  y, int n);

void nfft_upd_axpy_complex(double complex* x, double a, double complex* y, int n);

void nfft_upd_axpy_double( double*  x, double a, double*  y, int n);

void nfft_upd_xpay_complex(double complex* x, double a, double complex* y, int n);

void nfft_upd_xpay_double( double*  x, double a, double*  y, int n);

void nfft_upd_axpby_complex(double complex* x, double a, double complex* y, double b, int n);

void nfft_upd_axpby_double(  double* x, double a, double*  y, double b, int n);

void nfft_upd_xpawy_complex(double complex* x, double a, double* w, double complex* y, int n);

void nfft_upd_xpawy_double( double*  x, double a, double* w, double*  y, int n);

void nfft_upd_axpwy_complex(double complex* x, double a, double* w, double complex* y, int n);

void nfft_upd_axpwy_double( double*  x, double a, double* w, double*  y, int n);

void nfft_fftshift_complex(double complex *x, int d, int* N);

double nfft_error_l_infty_complex(double complex *x, double complex *y, int n);

double nfft_error_l_infty_double(double *x, double *y, int n);

double nfft_error_l_infty_1_complex(double complex *x, double complex *y, int n, double complex *z,
                               int m);

double nfft_error_l_infty_1_double(double *x, double *y, int n, double *z,
            int m);

double nfft_error_l_2_complex(double complex *x, double complex *y, int n);

double  nfft_error_l_2_double(double *x, double *y, int n);

void nfft_vpr_int(int *x, int n, char *text);

void nfft_vpr_double(double *x, int n, char *text);

void nfft_vpr_complex(double complex *x, int n, char *text);

void nfft_vrand_unit_complex(double complex *x, int n);

void nfft_vrand_shifted_unit_double(double *x, int n);

/* ######################################################################### */
/* ########## Helpers for inverse transforms ############################### */
/* ######################################################################### */

void nfft_voronoi_weights_1d(double *w, double *x, int M);

double nfft_modified_fejer(int N,int kk);

double nfft_modified_jackson2(int N,int kk);

double nfft_modified_jackson4(int N,int kk);

double nfft_modified_sobolev(double mu,int kk);

double nfft_modified_multiquadric(double mu,double c,int kk);

#endif

---

Generated on 1 Nov 2006 by Doxygen 1.4.4