linalg.hpp

Go to the documentation of this file.

/*#####################################################################
#
#  LINEAR ALGEBRA INTERFACE:
#
#  FOR A MOMENT, I STILL USE THOSE FROM THE EIGEN LIBRARY,
#  BUT THY MAY BE REPLACED BY DIRECT BLAS/LAPACK CALLS
#  AT A LATER TIME
#
#####################################################################*/

#ifndef _LINALG_H
#define _LINALG_H 1

#include <stdio.h>
#include <complex>
#include <cassert>
#include <cstring>

namespace linalg {

// interface to simple linear algebra routines:
// implementation (using EIGEN library) defined in linalg_eigen.cpp

//void cplx_sq_solve(void *a,void  *b,void *x,int n,int d);

void eigen_hermv(int size,std::complex<double> *A,double *eval,std::complex<double> *evec);


void real_sq_solve(double *ad,double *bd,double *xd,int dim);


void cplx_sq_solve(void *a,void *b,void *x,int dim);

void cplx_sq_solve_many(void *a,void *b,void *x,int dim,int d);

//void cplx_tri_solve(double *a,double *b,double *x,int dim);

void cplx_matrix_inverse(void *a,void *x,int n);

void linalg_matrix_inverse(double *a,double *x,int n);
//void QR_decomposition(double *aa,double *qq,double *rr,int n,int m);


} // namespace



#endif  // linalg