Linear_algebra.hpp

Go to the documentation of this file.

 
#pragma once
 
#include <vector>
#include <iostream> //for bit push << overriding
 
#include <gsl/gsl_matrix.h>
#include <gsl/gsl_blas.h>
#include <gsl/gsl_linalg.h>
 
class Matrix
{
public:
    gsl_matrix *mat;
 
    Matrix(const int &nrow, const int &ncol, const bool &setToZero);
 
 
 
    Matrix();
 
 
    Matrix(const int &nrow);
 
    Matrix(const int &nrow, const int &ncol);
 
    Matrix(const Matrix &mat_src);
 
    Matrix(const std::vector<std::vector<double>> &vec);
 
    Matrix& operator=(const Matrix &mat_src);
 
    ~Matrix();
 
    void free();
 
    int nrow() const;
 
    int ncol() const;
 
    double get(const int &i, const int &j) const;
 
    double operator()(const int &i, const int &j) const;
 
    void set(const int &i, const int &j, const double &v);
 
    void set_all(const double &v);
 
    void transpose();
 
    void set_identity();
 
    Matrix get_row(const int &irow) const;
 
    Matrix get_col(const int &icol) const;
 
    void set_row_from_col(const int &destRow, const Matrix &srcMat, const int &srcCol);
 
    Matrix get_transposed() const;
 
    void get_inverted(Matrix &matInverted);
 
    double calculate_norm() const;
 
    double calculate_quad() const;
 
};
 
Matrix mat_calloc(const int &nrow, const int &ncol);
 
std::ostream &operator<<(std::ostream &stream, const Matrix &matrix);
 
Matrix &operator+=(Matrix &m1, const Matrix &m2);
 
Matrix operator+(const Matrix &m1, const Matrix &m2);
 
Matrix &operator-=(Matrix &m1, const Matrix &m2);
 
Matrix operator-(const Matrix &m1, const Matrix &m2);
 
Matrix &operator*=(Matrix &mat, const double &scalar);
 
Matrix operator*(const Matrix &matrix, const double &scalar);
Matrix operator*(const double &scalar, const Matrix &matrix);
 
Matrix &operator*=(Matrix m1, const Matrix &m2);
 
Matrix operator*(const Matrix &m1, const Matrix &m2);
Matrix dot(const Matrix &m1, const Matrix &m2);
 
Matrix &operator/=(Matrix &mat, const double &scalar);
 
Matrix operator/(Matrix matrix, const double &scalar);
 
double dot_col(const Matrix &m1, const Matrix &m2);
 
double dot_row(const Matrix &m1, const Matrix &m2);
 
Matrix cross_col(const Matrix &m1, const Matrix &m2);
 
Matrix cross_row(const Matrix &m1, const Matrix &m2);
 
void negative(const Matrix &m1, Matrix &m_neg);
 
void get_unit_vector(Matrix &m);
 
void get_unit_vector(const Matrix &m1, Matrix &m_unit);
 
void cross(const Matrix &m1, const Matrix &m2, Matrix &temp);
 
void addition(const Matrix &m1, const Matrix &m2, Matrix &tmp);
 
void subtraction(const Matrix &m1, const Matrix &m2, Matrix &tmp);
 
void multiplication(const Matrix &m1, const Matrix &m2, Matrix &tmp);
 
void const_multiplication(const Matrix &m1, const double num, Matrix &tmp);
 
void const_division(const Matrix &m1, const double num, Matrix &tmp);
 
void colvec_matrix_multiplication(const Matrix &v1, const Matrix &m1, Matrix &tmp);
 
Matrix kron(const Matrix& v1, const Matrix& v2);
 
Matrix kron(const Matrix &v1, const Matrix &v2, Matrix &m_result);
 
void assign_rowVec_to_colVec(const Matrix &srcRowVec, Matrix &destColVec);
 
void a_cross_b_plus_c_cross_d(const Matrix &a,
                              const Matrix &b,
                              const Matrix &c,
                              const Matrix &d,
                              Matrix &tmp_f,
                              Matrix &tmp_l,
                              Matrix &v_result);