egsl.h

Go to the documentation of this file.

#ifndef H_EASY_GSL
#define H_EASY_GSL
 
#include <gsl/gsl_vector.h>
#include <gsl/gsl_matrix.h>
 
 
#ifdef __cplusplus
extern "C" {
#endif
 
struct egsl_val {
        gsl_matrix * gslm;
        int cid : 16;
        int index : 16;
};
 
typedef struct egsl_val val;
 
/* Core functions */
 
/* Push a new context. */
void egsl_push();
void egsl_push_named(const char*name);
/* Pops a context */
void egsl_pop();
void egsl_pop_named(const char*name);
void egsl_free(void);
 
double* egsl_atmp(val v, size_t i, size_t j);
val egsl_alloc(size_t rows, size_t columns);
val egsl_alloc_in_context(int cid, size_t rows, size_t cols);
gsl_matrix * egsl_gslm(val v);
val egsl_promote(val v);
 
val egsl_scale(double, val);
val egsl_sum(val, val);
val egsl_sum3(val, val, val);
val egsl_mult(val, val);
val egsl_transpose(val);
val egsl_inverse(val);
val egsl_sub(val,val);
val egsl_sum(val v1,val v2);
val egsl_compose_col(val v1, val v2);
val egsl_compose_row(val v1, val v2);
void egsl_add_to(val v1, val v2);
void egsl_add_to_col(val v1, size_t j, val v2);
 
double egsl_norm(val);
 
void egsl_symm_eig(val v, double* eigenvalues, val* eigenvectors);
 
double egsl_atv(val, size_t i);
double egsl_atm(val, size_t i, size_t j);
 
/* File: egsl_conversions.c
  Conversions */
 
val egsl_vFa(size_t rows, const double*);
val egsl_vFda(size_t rows, size_t columns, const double*);
 
void egsl_v2a(val, double*);
void egsl_v2da(val, double*);
void egsl_v2vec(val, gsl_vector*);
 
val egsl_vFgslv(const gsl_vector*);
val egsl_vFgslm(const gsl_matrix*);
 
gsl_matrix* egsl_v2gslm(val);
 
/*/ File: egsl_misc.c
    Miscellaneous useful matrixes. */
val egsl_zeros(size_t rows, size_t columns);
val egsl_ones(size_t rows, size_t columns);
val egsl_vers(double theta);
val egsl_rot(double theta);
 
 
/* Misc */
void egsl_print(const char*str, val);   
void egsl_print_spectrum(const char*s, val v);
void egsl_print_stats(void);
 
 
        
void egsl_expect_size(val v, size_t rows, size_t cols);
void egsl_error(void);
 
 
#ifdef __cplusplus
}
#endif
 
#endif