#include <gsl/gsl_matrix.h>

#ifndef KERNELS_H
#define KERNELS_H
typedef int(*kernel)(double, double, gsl_matrix* , gsl_matrix*, gsl_matrix*);

typedef struct {
  kernel kern;
  double sigmaf;
  double length_scale;
} kernel_info;

int 
squared_exp 
( double sigmaf, double ell, gsl_matrix *A, gsl_matrix *B, gsl_matrix *result );

int
matern32
( double sigmaf, double ell, gsl_matrix *A, gsl_matrix *B, gsl_matrix *result );

#endif // !KERNELS_H