#if !defined(invertor_h)
#define invertor_h


typedef struct {
    double d_max;
    double *x;
    double *y;
    double *err;
    int npoints;    
    int ny;    
    int nerr;    
} Invertor_params;

void invertor_dealloc(Invertor_params *pars);

void invertor_init(Invertor_params *pars);


double pr_sphere(double R, double r);
double ortho(double d_max, int n, double r);
double ortho_transformed(double d_max, int n, double q);
double ortho_derived(double d_max, int n, double r);
double iq(double *c, double d_max, int n_c, double q);
double pr(double *c, double d_max, int n_c, double r);

#endif