excess_map/kernels.h at master · foreverstupid/excess_map · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
#ifndef KERNELS_MODULE_H
#define KERNELS_MODULE_H

#include <math.h>
#include <gsl/gsl_math.h>
#include <gsl/gsl_multiroots.h>

#include "vector.h"

/*
 * Params for calculation method
 */
struct params{
    double k;                   /* excess kurtosis value */
    double d;                   /* dispersion value */

    struct vector_func buffer;  /* calculation buffer */
};


/* Kurtic kernel */
int kurtic_f(const gsl_vector *x, void *params, gsl_vector *f);

int kurtic_df(const gsl_vector *x, void *params, gsl_matrix *J);

int kurtic_fdf(const gsl_vector *x, void *params, gsl_vector *f,
    gsl_matrix *J);


/* Exonent polynomial kernel */
int polyexp_f(const gsl_vector *x, void *params, gsl_vector *f);

/*int kurtic_df(const gsl_vector *x, void *params, gsl_matrix *J);

int kurtic_fdf(const gsl_vector *x, void *params, gsl_vector *f,
    gsl_matrix *J);
*/


/* Roughgarden kernel */
int rgarden_f(const gsl_vector *x, void *params, gsl_vector *f);

#endif