fpt.h

/*
 * Copyright (c) 2002, 2017 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.
 */

#ifndef _FPT_H_
#define _FPT_H_

#include <stdbool.h>

void fpt_precompute_1(fpt_set set, const int m, int k_start);
void fpt_precompute_2(fpt_set set, const int m, double *alpha, double *beta, double *gam, int k_start, const double threshold);

typedef struct fpt_step_
{
  bool stable;                            
  int Ns;                                 
  int ts;                                 
  double *a;                              
//  double *a11,*a12,*a21,*a22;         /**< The matrix components          */
  double g;                              
} fpt_step;

typedef struct fpt_data_
{
  fpt_step **steps;                       
  int k_start;                            
  double *alphaN;                         
  double *betaN;                          
  double *gammaN;                         
  double alpha_0;                         
  double beta_0;                          
  double gamma_m1;                        
  /* Data for direct transform. */        
  double *_alpha;                         
  double *_beta;                          
  double *_gamma;                         
  bool precomputed;
} fpt_data;

typedef struct fpt_set_s_
{
  unsigned int flags;                     
  int M;                                  
  int N;                                  
  int t;                                  
  fpt_data *dpt;                          
  double **xcvecs;                        
  double *xc;                             
  double _Complex *temp;                          
  double _Complex *work;                          
  double _Complex *result;                        
  double _Complex *vec3;
  double _Complex *vec4;
  double _Complex *z;
  fftw_plan *plans_dct3;                  
  fftw_plan *plans_dct2;                  
  fftw_r2r_kind *kinds;                   
  fftw_r2r_kind *kindsr;                  
  /* Data for slow transforms. */
  double *xc_slow;
} fpt_set_s;

#endif /*_FPT_H_*/