#include <gsl/gsl_matrix.h>

#ifndef SOLVE_H
#define SOLVE_H

typedef int (*invert_t)(gsl_matrix*, gsl_matrix*);

int
inv_lu
( gsl_matrix *A, gsl_matrix *inverse );

int
inv_chol
( gsl_matrix *A, gsl_matrix *inverse );

#endif // !SOLVE_H