fastsum.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: fastsum.h 3896 2012-10-10 12:19:26Z tovo $ */

#ifndef fastsum_h_inc
#define fastsum_h_inc

#include "config.h"

#ifdef HAVE_COMPLEX_H
#include <complex.h>
#endif

#include "nfft3util.h"
#include "nfft3.h"

#if !(defined(NF_LIN) || defined(NF_QUADR) || defined(NF_KUB))
  #define NF_KUB
#endif

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

typedef double _Complex (*kernel)(double , int , const double *);

#define EXACT_NEARFIELD  (1U<< 0)

#define NEARFIELD_BOXES (1U<< 1)

typedef struct fastsum_plan_
{
  int d;                                
  int N_total;                          
  int M_total;                          
  double _Complex *alpha;                       
  double _Complex *f;                           
  double *x;                            
  double *y;                            
  kernel k;  
  double *kernel_param;                 
  unsigned flags;                       
  double _Complex *pre_K;                       
  int n;                                
  fftw_complex *b;                      
  int p;                                
  double eps_I;                           /* fixed to p/n so far  */
  double eps_B;                           /* fixed to 1/16 so far */

  nfft_plan mv1;                        
  nfft_plan mv2;                        
  int Ad;                               
  double _Complex *Add;                 
  /* things for computing *b - are they used only once?? */
  fftw_plan fft_plan;

  int box_count;
  int box_count_per_dim;
  int *box_offset;
  double *box_x;
  double _Complex *box_alpha;

  double MEASURE_TIME_t[8]; 
} fastsum_plan;

void fastsum_init_guru(fastsum_plan *ths, int d, int N_total, int M_total, kernel k, double *param, unsigned flags, int nn, int m, int p, double eps_I, double eps_B);

void fastsum_finalize(fastsum_plan *ths);

void fastsum_exact(fastsum_plan *ths);

void fastsum_precompute(fastsum_plan *ths);

void fastsum_trafo(fastsum_plan *ths);
/* \} */

double _Complex regkern(kernel k, double xx, int p, const double *param, double a, double b);

double _Complex kubintkern(const double x, const double _Complex *Add,
  const int Ad, const double a);

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

#endif
/* fastsum.h */