matrix_tools.h

Go to the documentation of this file.

#ifndef _MAT_H
#define _MAT_H

#include <gsl/gsl_matrix.h>
#include <gsl/gsl_vector.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_linalg.h>

#include <vector>
using namespace std; 

//matrix-tools
gsl_matrix * colcat(const gsl_matrix *a,const gsl_matrix *b);
gsl_matrix * rowcat(const gsl_matrix *a,const gsl_matrix *b);
double gsl_det(gsl_matrix * m);
int gsl_inverse(gsl_matrix * m,gsl_matrix * im);
double gsl_inverse_det(gsl_matrix * m,gsl_matrix * im);
int gsl_linsolve_lu(gsl_matrix * A,gsl_vector * b,gsl_vector * x);
int gsl_linsolve_svd(gsl_matrix * A,gsl_vector * b,gsl_vector * x);

int encaps_gsl_linsolveSVD(vector<vector<double> > &A, vector<double> &b, vector<double> &x);

int encaps_gsl_linsolveLU(double A[], double b[], double x[], int size);  
void gsl_vector2double(gsl_vector *gsl_v, double *v, size_t size);
void gsl_vector2double_offset(gsl_vector *gsl_v, double *v, size_t size, size_t offset);
void sub_double(double *v, double *vi, size_t size, size_t offset);
void gsl_matrix2double(gsl_matrix *gsl_m, double *m, size_t sizei, size_t sizej);
int gsl_rank_svd(gsl_matrix* A);
int gsl_inverse_svd(gsl_matrix *A, gsl_matrix *invA);
int encaps_gsl_inverse_svd(double *A, double *invA, int noc, int nor);
double gsl_vector_norme(gsl_vector *v);
int gsl_vector_cross_product3D(gsl_vector *A, gsl_vector *B, gsl_vector *AB);
int gsl_vector_create_row(gsl_vector *row, gsl_vector * A, gsl_vector *B);
double norme_diff_vecteur(double *v1, double *v2, size_t size);

void vect2double(vector<double> &vecA, double *A);
void vectormat2double(vector< vector<double> > &matA, double *A);

void gsl_dbl2matrix(double *A, gsl_matrix *Matr);
int gsl_matrixprod(gsl_matrix *A, gsl_matrix *B, gsl_matrix *C);
int gsl_matrixvecprod(gsl_matrix *A, gsl_vector *B, gsl_vector *C);



#endif