nfft3util.h

Go to the documentation of this file.

/*
 * Copyright (c) 2002, 2012 Jens Keiner, Stefan Kunis, Daniel Potts
 *
 * This program is free software; you can redistribute it and/or modify it under
 * the terms of the GNU General Public License as published by the Free Software
 * Foundation; either version 2 of the License, or (at your option) any later
 * version.
 *
 * This program is distributed in the hope that it will be useful, but WITHOUT
 * ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS
 * FOR A PARTICULAR PURPOSE.  See the GNU General Public License for more
 * details.
 *
 * You should have received a copy of the GNU General Public License along with
 * this program; if not, write to the Free Software Foundation, Inc., 51
 * Franklin Street, Fifth Floor, Boston, MA 02110-1301, USA.
 */

/* $Id: nfft3util.h 3896 2012-10-10 12:19:26Z tovo $ */

#ifndef __UTIL_H__
#define __UTIL_H__

#include <fftw3.h>

#ifdef __cplusplus
extern "C"
{
#endif /* __cplusplus */

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

#define NFFT_SWAP_complex(x,y) {fftw_complex* NFFT_SWAP_temp__; \
  NFFT_SWAP_temp__=(x); (x)=(y); (y)=NFFT_SWAP_temp__;}

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

#define PI 3.141592653589793238462643383279502884197169399375105820974944592
#define PI2 6.283185307179586476925286766559005768394338798750211641949889185
#define PI4 12.56637061435917295385057353311801153678867759750042328389977837

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

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

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

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

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

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

int nfft_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(fftw_complex* x, int n);

double nfft_dot_double( double*  x, int n);

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

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

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

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

void nfft_cp_complex(fftw_complex* x, fftw_complex* y, int n);

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

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

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

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

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

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

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

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

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

void nfft_upd_axpby_complex(fftw_complex* x, double a, fftw_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(fftw_complex* x, double a, double* w, fftw_complex* y, int n);

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

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

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

void nfft_fftshift_complex(fftw_complex *x, int d, int* N);

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

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

void nfft_vpr_complex(fftw_complex *x, int n, const char *text);

void nfft_vrand_unit_complex(fftw_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);

void nfft_voronoi_weights_S2(double *w, double *xi, 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);

int nfft_smbi(const double x, const double alpha, const int nb, const int ize,
  double *b);

double nfft_drand48(void);

void nfft_srand48(long int seed);

void nfft_sort_node_indices_radix_lsdf(int n, int *keys0, int *keys1, int rhigh);

void nfft_sort_node_indices_radix_msdf(int n, int *keys0, int *keys1, int rhigh);

int nfft_get_num_threads(void);

#ifdef _OPENMP
int nfft_get_omp_num_threads(void);
#endif

#ifdef __cplusplus
}  /* extern "C" */
#endif /* __cplusplus */

#endif