rpolator.h

Go to the documentation of this file.

/* ============================================================================

   Copyright (C) 2017, 2018, 2019  Konrad Bernloehr
   
   This file is part of sim_telarray (also known as sim_hessarray)
   and also part of the IACT/atmo package for CORSIKA.

   The software in this file is free software: you can redistribute it 
   and/or modify it under the terms of the GNU Lesser General Public
   License as published by the Free Software Foundation; either
   version 2.1 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 Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with this program. If not, see <http://www.gnu.org/licenses/>.

============================================================================ */

/* =============================================================== */
/* =============================================================== */

#ifndef RPOLATOR_H__LOADED
#define RPOLATOR_H__LOADED 1

#define WITH_RPOLATOR 1

#ifdef __cplusplus
extern "C" {
#endif

struct cubic_params
{
   double a, b, c, d;  
};
typedef struct cubic_params CsplinePar;

struct rpol_table
{
   int ndim;          
   size_t nx, ny;     
   double *x;         
   double *y;         
   double *z;         
   double *zxmax;     
   double *zxmin;     
   double aux;        
   double xmin, xmax; 
   double dx, dxi;    
   double xrise;      
   double ymin, ymax; 
   double dy, dyi;    
   double yrise;      
   double zmin, zmax; 
   char *fname;       
   char *options;     
   int equidistant;   
   int scheme;        
   int clipping;      
   int remapped;      
   int zxreq;         
   int logs, xlog, ylog, zlog; 
   CsplinePar *csp;   
   int use_count;     
};
typedef struct rpol_table RpolTable;

/* rpolator.c */
int read_table(const char *fname, int maxrow, double *col1, double *col2);
int read_table2(const char *fname, int maxrow, double *col1, double *col2);
int read_table3(const char *fname, int maxrow, double *col1, double *col2, double *col3);
int read_table4(const char *fname, int maxrow, double *col1, double *col2, double *col3, double *col4);
int read_table5(const char *fname, int maxrow, double *col1, double *col2, double *col3, double *col4, double *col5);

int read_table_v(const char *fname, FILE *fptr, size_t *nrow, size_t ncol, double ***col, const size_t *selcol);

struct rpol_table *read_rpol1d_table(const char *fn);
struct rpol_table *read_rpol2d_table(const char *fn, const char *ymarker);
struct rpol_table *read_rpol3d_table(const char *fn);
struct rpol_table *read_rpol_table(const char *fname, int nd, const char *ymarker, const char *options);

struct rpol_table *simple_rpol1d_table(const char *label, double *x, double *y, int n, int clip);

int rpol_is_verbose(void);
int rpol_set_verbose(int lvl);

void rpol_info(struct rpol_table *rpt);
void rpol_info_lvl (struct rpol_table *rpt, int lvl);
void rpol_free(struct rpol_table *rpt, int removing);
void rpol_check_equi_range(struct rpol_table *rpt);

double rpol(double *x, double *y, int n, double xp);
double rpol_nearest(double *x, double *y, int n, double xp, int eq, int clip);
double rpol_linear(double *x, double *y, int n, double xp, int eq, int clip);
double rpol_2nd_order(double *x, double *y, int n, double xp, int eq, int clip);
CsplinePar* set_1d_cubic_params (double *x, double *y, int n, int clamped);
double rpol_cspline ( double *x, double *y, const CsplinePar *csp, int n, double xp, int eq, int clip );
double rpol_2d_linear(double *x, double *y, double *z, int nx, int ny, double xp, double yp, int eq, int clip);

double rpolate_1d(struct rpol_table *rpt, double x, int scheme);
double rpolate_1d_lin(struct rpol_table *rpt, double x);
double rpolate_2d(struct rpol_table *rpt, double x, double y, int scheme);
double rpolate(struct rpol_table *rpt, double x, double y, int scheme);

#ifdef __cplusplus
}
#endif

#endif