include/it/mat.h File Reference

#include <stddef.h>
#include <stdlib.h>
#include <string.h>
#include <it/types.h>
#include <it/vec.h>

Go to the source code of this file.

Basic vector allocation, initialization and deletion
#define	Mat_new_alloc(type_t, h, w, hmax, wmax)   ((type_t **) __Mat_new_alloc( sizeof(type_t), h, w, hmax, wmax ))
#define	Mat_new(type_t, h, w)   ((type_t **) __Mat_new( sizeof(type_t), h, w ))
#define	Mat_delete(m)   __Mat_delete((Mat) m)
Mat	__Mat_new_alloc (size_t elem_size, idx_t h, idx_t w, idx_t hmax, idx_t wmax)
mat	mat_new_alloc (idx_t h, idx_t w, idx_t hmax, idx_t wmax)
imat	imat_new_alloc (idx_t h, idx_t w, idx_t hmax, idx_t wmax)
bmat	bmat_new_alloc (idx_t h, idx_t w, idx_t hmax, idx_t wmax)
cmat	cmat_new_alloc (idx_t h, idx_t w, idx_t hmax, idx_t wmax)
Mat	__Mat_new (size_t elem_size, idx_t h, idx_t w)
mat	mat_new (idx_t h, idx_t w)
imat	imat_new (idx_t h, idx_t w)
bmat	bmat_new (idx_t h, idx_t w)
cmat	cmat_new (idx_t h, idx_t w)
void	__Mat_delete (Mat m)
void	mat_delete (mat m)
void	imat_delete (imat m)
void	bmat_delete (bmat m)
void	cmat_delete (cmat m)
Width and Height functions
#define	Mat_height(m)   Vec_length(m)
#define	Mat_height_max(m)   Vec_length_max(m)
#define	Mat_width(m)   Vec_length(m[0])
#define	Mat_width_max(m)   Vec_length_max(m[0])
#define	Mat_set_height_max(m, hmax)
#define	mat_set_height_max(m, h)   do { m = _mat_set_height_max(m, h); } while(0)
#define	imat_set_height_max(m, h)   do { m = _imat_set_height_max(m, h); } while(0)
#define	bmat_set_height_max(m, h)   do { m = _bmat_set_height_max(m, h); } while(0)
#define	cmat_set_height_max(m, h)   do { m = _cmat_set_height_max(m, h); } while(0)
#define	Mat_set_height(m, h)
#define	mat_set_height(m, h)   do { m = _mat_set_height(m, h); } while(0)
#define	imat_set_height(m, h)   do { m = _imat_set_height(m, h); } while(0)
#define	bmat_set_height(m, h)   do { m = _bmat_set_height(m, h); } while(0)
#define	cmat_set_height(m, h)   do { m = _cmat_set_height(m, h); } while(0)
#define	MAT_END_COL_PARAM(v, param)   do { if( param == end ) param = Mat_width( v ) - 1; } while(0)
#define	MAT_END_ROW_PARAM(v, param)   do { if( param == end ) param = Mat_height( v ) - 1; } while(0)
idx_t	mat_height (mat m)
idx_t	imat_height (imat m)
idx_t	bmat_height (bmat m)
idx_t	cmat_height (cmat m)
idx_t	mat_height_max (mat m)
idx_t	imat_height_max (imat m)
idx_t	bmat_height_max (bmat m)
idx_t	cmat_height_max (cmat m)
idx_t	mat_width (mat m)
idx_t	imat_width (imat m)
idx_t	bmat_width (bmat m)
idx_t	cmat_width (cmat m)
idx_t	mat_width_max (mat m)
idx_t	imat_width_max (imat m)
idx_t	bmat_width_max (bmat m)
idx_t	cmat_width_max (cmat m)
mat	_mat_set_height_max (mat m, idx_t hmax)
imat	_imat_set_height_max (imat m, idx_t hmax)
bmat	_bmat_set_height_max (bmat m, idx_t hmax)
cmat	_cmat_set_height_max (cmat m, idx_t hmax)
mat	_mat_set_height (mat m, idx_t h)
imat	_imat_set_height (imat m, idx_t h)
bmat	_bmat_set_height (bmat m, idx_t h)
cmat	_cmat_set_height (cmat m, idx_t h)
Basic matrix initialization/modification
#define	Mat_init(m, buf, w, h)
#define	mat_init(m, buf, w, h)   _mat_init(m, buf, w, h)
#define	imat_init(m, buf, w, h)   _imat_init(m, buf, w, h)
#define	bmat_init(m, buf, w, h)   _bmat_init(m, buf, w, h)
#define	cmat_init(m, buf, w, h)   _cmat_init(m, buf, w, h)
#define	Mat_set(m, val)
#define	Mat_set_between(m, r1, c1, r2, c2, val)
mat	_mat_init (mat m, double *buf, idx_t w, idx_t h)
imat	_imat_init (imat m, int *buf, idx_t w, idx_t h)
bmat	_bmat_init (bmat m, byte *buf, idx_t w, idx_t h)
cmat	_cmat_init (cmat m, cplx *buf, idx_t w, idx_t h)
mat	mat_set (mat m, double val)
imat	imat_set (imat m, int val)
bmat	bmat_set (bmat m, byte val)
cmat	cmat_set (cmat m, cplx val)
mat	mat_new_set (double val, idx_t h, idx_t w)
imat	imat_new_set (int val, idx_t h, idx_t w)
bmat	bmat_new_set (byte val, idx_t h, idx_t w)
cmat	cmat_new_set (cplx val, idx_t h, idx_t w)
mat	mat_set_between (mat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2, double val)
imat	imat_set_between (imat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2, int val)
bmat	bmat_set_between (bmat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2, byte val)
cmat	cmat_set_between (cmat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2, cplx val)
void	mat_void (mat m)
void	imat_void (imat m)
void	bmat_void (bmat m)
void	cmat_void (cmat m)
mat	mat_new_void ()
imat	imat_new_void ()
bmat	bmat_new_void ()
cmat	cmat_new_void ()
void	mat_zeros (mat m)
void	imat_zeros (imat m)
void	bmat_zeros (bmat m)
void	cmat_zeros (cmat m)
mat	mat_new_zeros (idx_t h, idx_t w)
imat	imat_new_zeros (idx_t h, idx_t w)
bmat	bmat_new_zeros (idx_t h, idx_t w)
cmat	cmat_new_zeros (idx_t h, idx_t w)
void	mat_eye (mat m)
void	imat_eye (imat m)
void	bmat_eye (bmat m)
void	cmat_eye (cmat m)
mat	mat_new_eye (idx_t n)
imat	imat_new_eye (idx_t n)
bmat	bmat_new_eye (idx_t n)
cmat	cmat_new_eye (idx_t n)
void	mat_ones (mat m)
void	imat_ones (imat m)
void	bmat_ones (bmat m)
void	cmat_ones (cmat m)
mat	mat_new_ones (idx_t h, idx_t w)
imat	imat_new_ones (idx_t h, idx_t w)
bmat	bmat_new_ones (idx_t h, idx_t w)
cmat	cmat_new_ones (idx_t h, idx_t w)
mat	mat_new_diag (vec v)
imat	imat_new_diag (ivec v)
bmat	bmat_new_diag (bvec v)
Copy and Conversions Functions
#define	Mat_copy(m1, m2)   __Mat_copy( (Mat) m1, (Mat) m2 )
Mat	__Mat_copy (Mat m1, Mat m2)
void	mat_copy (mat dest, mat orig)
void	imat_copy (imat dest, imat orig)
void	bmat_copy (bmat dest, bmat orig)
void	cmat_copy (cmat dest, cmat orig)
Mat	Mat_clone (Mat m)
mat	mat_clone (mat m)
imat	imat_clone (imat m)
bmat	bmat_clone (bmat m)
cmat	cmat_clone (cmat m)
vec	mat_to_vec (mat m)
ivec	imat_to_ivec (imat m)
bvec	bmat_to_bvec (bmat m)
cvec	cmat_to_cvec (cmat m)
mat	vec_to_mat (vec v, idx_t width)
imat	ivec_to_imat (ivec v, idx_t width)
bmat	bvec_to_bmat (bvec v, idx_t width)
cmat	cvec_to_cmat (cvec v, idx_t width)
mat	imat_to_mat (imat m)
mat	bmat_to_mat (bmat m)
imat	mat_to_imat (mat m)
imat	bmat_to_imat (bmat m)
Row/column manipulation
#define	Mat_del_row(m, pos)
#define	Mat_ins_row(m, pos, v)
#define	Mat_push_row(m, v)   Mat_ins_row( m, Mat_height(m), v )
#define	Mat_pop_row(m)   Mat_del_row( m, Mat_height(m)-1 )
#define	Mat_last_row(m)   ((m)[Mat_height(m)-1])
#define	mat_del_row(m, pos)   do { m = __mat_del_row( m, pos ); } while (0)
#define	imat_del_row(m, pos)   do { m = __imat_del_row( m, pos ); } while (0)
#define	bmat_del_row(m, pos)   do { m = __bmat_del_row( m, pos ); } while (0)
#define	cmat_del_row(m, pos)   do { m = __cmat_del_row( m, pos ); } while (0)
#define	mat_ins_row(m, pos, v)   do { m = __mat_ins_row( m, pos, v ); } while (0)
#define	imat_ins_row(m, pos, v)   do { m = __imat_ins_row( m, pos, v ); } while (0)
#define	bmat_ins_row(m, pos, v)   do { m = __bmat_ins_row( m, pos, v ); } while (0)
#define	cmat_ins_row(m, pos, v)   do { m = __cmat_ins_row( m, pos, v ); } while (0)
#define	mat_push_row(m, v)   do { m = __mat_push_row( m, v ); } while (0)
#define	imat_push_row(m, v)   do { m = __imat_push_row( m, v ); } while (0)
#define	bmat_push_row(m, v)   do { m = __bmat_push_row( m, v ); } while (0)
#define	cmat_push_row(m, v)   do { m = __cmat_push_row( m, v ); } while (0)
#define	mat_pop_row(m)   do { m = __mat_pop_row( m ); } while (0)
#define	imat_pop_row(m)   do { m = __imat_pop_row( m ); } while (0)
#define	bmat_pop_row(m)   do { m = __bmat_pop_row( m ); } while (0)
#define	cmat_pop_row(m)   do { m = __cmat_pop_row( m ); } while (0)
#define	mat_transpose(m)
#define	imat_transpose(m)
#define	bmat_transpose(m)
mat	__mat_del_row (mat m, idx_t pos)
imat	__imat_del_row (imat m, idx_t pos)
bmat	__bmat_del_row (bmat m, idx_t pos)
cmat	__cmat_del_row (cmat m, idx_t pos)
mat	__mat_ins_row (mat m, idx_t pos, vec v)
imat	__imat_ins_row (imat m, idx_t pos, ivec v)
bmat	__bmat_ins_row (bmat m, idx_t pos, bvec v)
cmat	__cmat_ins_row (cmat m, idx_t pos, cvec v)
mat	__mat_push_row (mat m, vec v)
imat	__imat_push_row (imat m, ivec v)
bmat	__bmat_push_row (bmat m, bvec v)
cmat	__cmat_push_row (cmat m, cvec v)
mat	__mat_pop_row (mat m)
imat	__imat_pop_row (imat m)
bmat	__bmat_pop_row (bmat m)
cmat	__cmat_pop_row (cmat m)
void	mat_swap_rows (mat m, idx_t i, idx_t j)
void	mat_swap_cols (mat m, idx_t i, idx_t j)
void	imat_swap_rows (imat m, idx_t i, idx_t j)
void	imat_swap_cols (imat m, idx_t i, idx_t j)
void	bmat_swap_rows (bmat m, idx_t i, idx_t j)
void	bmat_swap_cols (bmat m, idx_t i, idx_t j)
void	cmat_swap_rows (cmat m, idx_t i, idx_t j)
void	cmat_swap_cols (cmat m, idx_t i, idx_t j)
mat	_mat_transpose (mat m)
imat	_imat_transpose (imat m)
bmat	_bmat_transpose (bmat m)
mat	mat_new_transpose (mat m)
imat	imat_new_transpose (imat m)
bmat	bmat_new_transpose (bmat m)
Submatrices operations
#define	Mat_set_row(m, r, v)   Vec_copy((m)[r], (v))
#define	mat_set_row(m, r, v)   vec_copy((m)[r], (v))
#define	imat_set_row(m, r, v)   ivec_copy((m)[r], (v))
#define	bmat_set_row(m, r, v)   bvec_copy((m)[r], (v))
#define	cmat_set_row(m, r, v)   cvec_copy((m)[r], (v))
mat	mat_get_submatrix (mat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2)
imat	imat_get_submatrix (imat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2)
bmat	bmat_get_submatrix (bmat m, idx_t r1, idx_t c1, idx_t r2, idx_t c2)
mat	mat_set_submatrix (mat m, mat s, idx_t r, idx_t c)
imat	imat_set_submatrix (imat m, imat s, idx_t r, idx_t c)
bmat	bmat_set_submatrix (bmat m, bmat s, idx_t r, idx_t c)
void	mat_copy_row (vec v, mat m, idx_t r)
void	imat_copy_row (ivec v, imat m, idx_t r)
void	bmat_copy_row (bvec v, bmat m, idx_t r)
void	cmat_copy_row (cvec v, cmat m, idx_t r)
void	mat_copy_col (vec v, mat m, idx_t c)
void	imat_copy_col (ivec v, imat m, idx_t c)
void	bmat_copy_col (bvec v, bmat m, idx_t c)
void	cmat_copy_col (cvec v, cmat m, idx_t c)
vec	mat_get_row (mat m, idx_t r)
ivec	imat_get_row (imat m, idx_t r)
bvec	bmat_get_row (bmat m, idx_t r)
cvec	cmat_get_row (cmat m, idx_t r)
vec	mat_get_col (mat m, idx_t c)
ivec	imat_get_col (imat m, idx_t c)
bvec	bmat_get_col (bmat m, idx_t c)
cvec	cmat_get_col (cmat m, idx_t c)
void	Mat_set_col (Mat m, idx_t c, Vec v)
void	mat_set_col (mat m, idx_t c, vec v)
void	imat_set_col (imat m, idx_t c, ivec v)
void	bmat_set_col (bmat m, idx_t c, bvec v)
void	cmat_set_col (cmat m, idx_t c, cvec v)
vec	mat_get_diag (mat m)
ivec	imat_get_diag (imat m)
bvec	bmat_get_diag (bmat m)
cvec	cmat_get_diag (cmat m)
void	mat_diag (mat m, vec v)
void	imat_diag (imat m, ivec v)
void	bmat_diag (bmat m, bvec v)
void	cmat_diag (cmat m, cvec v)
void	mat_set_diag (mat m, vec v)
void	imat_set_diag (imat m, ivec v)
void	bmat_set_diag (bmat m, bvec v)
void	cmat_set_diag (cmat m, cvec v)
Component per component arbitrar function
#define	mat_eval(v, f, a)   mat_apply_function(v, f, a)
#define	mat_new_eval(v, f, a)   mat_new_apply_function(v, f, a)
#define	imat_eval(v, f, a)   imat_apply_function(v, f, a)
#define	imat_new_eval(v, f, a)   imat_new_apply_function(v, f, a)
void	mat_apply_function (mat v, it_function_t function, it_args_t args)
mat	mat_new_apply_function (mat v, it_function_t function, it_args_t args)
void	imat_apply_function (imat v, it_ifunction_t function, it_args_t args)
imat	imat_new_apply_function (imat v, it_ifunction_t function, it_args_t args)
Defines
#define	Mat_element_size(m)   Vec_element_size(m[0])
Typedefs
typedef Vec *	Mat
typedef vec *	mat
typedef ivec *	imat
typedef bvec *	bmat
typedef cvec *	cmat
Functions
Comparison functions
int	Mat_eq (Mat m1, Mat m2)
int	mat_eq (mat m1, mat m2)
int	imat_eq (imat m, imat m2)
int	bmat_eq (bmat m, bmat m2)
int	cmat_eq (cmat m, cmat m2)
int	mat_is_void (mat m)
int	imat_is_void (imat m)
int	bmat_is_void (bmat m)
int	cmat_is_void (cmat m)
Mean, median, min, max, sum and elementary statistic indicators
double	mat_sum (mat m)
int	imat_sum (imat m)
cplx	cmat_sum (cmat m)
double	mat_row_sum (mat m, idx_t c)
int	imat_row_sum (imat m, idx_t c)
cplx	cmat_row_sum (cmat m, idx_t c)
double	mat_col_sum (mat m, idx_t c)
int	imat_col_sum (imat m, idx_t c)
cplx	cmat_col_sum (cmat m, idx_t c)
vec	mat_rows_sum (mat m)
ivec	imat_rows_sum (imat m)
cvec	cmat_rows_sum (cmat m)
vec	mat_cols_sum (mat m)
ivec	imat_cols_sum (imat m)
cvec	cmat_cols_sum (cmat m)
double	mat_diag_sum (mat m)
int	imat_diag_sum (imat m)
int	bmat_diag_sum (bmat m)
cplx	cmat_diag_sum (cmat m)
double	mat_mean (mat m)
double	imat_mean (imat m)
cplx	cmat_mean (cmat m)
double	mat_max_index_submatrix (mat m, idx_t rmin, idx_t rmax, idx_t cmin, idx_t cmax, idx_t *r, idx_t *c)
double	mat_min_index_submatrix (mat m, idx_t rmin, idx_t rmax, idx_t cmin, idx_t cmax, idx_t *r, idx_t *c)
int	imat_max_index_submatrix (imat m, idx_t rmin, idx_t rmax, idx_t cmin, idx_t cmax, idx_t *r, idx_t *c)
int	imat_min_index_submatrix (imat m, idx_t rmin, idx_t rmax, idx_t cmin, idx_t cmax, idx_t *r, idx_t *c)
double	mat_max (mat m)
int	imat_max (imat m)
double	mat_min (mat m)
int	imat_min (imat m)
double	mat_max_index (mat m, idx_t *r, idx_t *c)
double	mat_max_col_index (mat m, idx_t c, idx_t *r)
double	mat_max_row_index (mat m, idx_t r, idx_t *c)
double	mat_min_index (mat m, idx_t *r, idx_t *c)
double	mat_min_col_index (mat m, idx_t c, idx_t *r)
double	mat_min_row_index (mat m, idx_t r, idx_t *c)
int	imat_max_index (imat m, idx_t *r, idx_t *c)
int	imat_max_col_index (imat m, idx_t c, idx_t *r)
int	imat_max_row_index (imat m, idx_t r, idx_t *c)
int	imat_min_index (imat m, idx_t *r, idx_t *c)
int	imat_min_col_index (imat m, idx_t c, idx_t *r)
int	imat_min_row_index (imat m, idx_t r, idx_t *c)
double	mat_variance (mat m)
void	mat_normalize (mat m)
mat	mat_cov (mat m)
void	mat_rows_normalize (mat m, double nr)
void	mat_cols_normalize (mat m, double nr)
double	mat_norm_1 (mat m)
double	mat_norm_inf (mat m)
Basic arithmetic operations
void	mat_incr (mat m, double a)
void	mat_decr (mat m, double a)
void	mat_mul_by (mat m, double a)
void	mat_div_by (mat m, double a)
void	imat_incr (imat m, int a)
void	imat_decr (imat m, int a)
void	imat_mul_by (imat m, int a)
void	imat_div_by (imat m, int a)
void	cmat_incr (cmat m, cplx a)
void	cmat_decr (cmat m, cplx a)
void	cmat_mul_by (cmat m, cplx a)
void	cmat_div_by (cmat m, cplx a)
Basic arithmetic operations on rows/columns
void	mat_col_set (mat m, idx_t col, double a)
void	mat_col_incr (mat m, idx_t col, double a)
void	mat_col_decr (mat m, idx_t col, double a)
void	mat_col_mul_by (mat m, idx_t col, double a)
void	mat_col_div_by (mat m, idx_t col, double a)
void	imat_col_set (imat m, idx_t col, int a)
void	imat_col_incr (imat m, idx_t col, int a)
void	imat_col_decr (imat m, idx_t col, int a)
void	imat_col_mul_by (imat m, idx_t col, int a)
void	imat_col_div_by (imat m, idx_t col, int a)
void	mat_row_incr (mat m, idx_t row, double a)
void	mat_row_decr (mat m, idx_t row, double a)
void	mat_row_mul_by (mat m, idx_t row, double a)
void	mat_row_div_by (mat m, idx_t row, double a)
void	imat_row_incr (imat m, idx_t row, int a)
void	imat_row_decr (imat m, idx_t row, int a)
void	imat_row_mul_by (imat m, idx_t row, int a)
void	imat_row_div_by (imat m, idx_t row, int a)
Component per component operations
void	mat_elem_add (mat m1, mat m2)
void	mat_elem_sub (mat m1, mat m2)
void	mat_elem_mul (mat m1, mat m2)
void	mat_elem_div (mat m1, mat m2)
void	mat_elem_sqr (mat m)
void	mat_elem_abs (mat m)
void	mat_elem_exp (mat m)
void	mat_elem_sqrt (mat m)
void	mat_elem_pow (mat m, double p)
void	imat_elem_add (imat m1, imat m2)
void	imat_elem_sub (imat m1, imat m2)
void	imat_elem_mul (imat m1, imat m2)
void	imat_elem_div (imat m1, imat m2)
void	imat_elem_sqr (imat m)
void	cmat_elem_add (cmat m1, cmat m2)
void	cmat_elem_sub (cmat m1, cmat m2)
void	cmat_elem_mul (cmat m1, cmat m2)
void	cmat_elem_div (cmat m1, cmat m2)
Matrix arithmetic
void	mat_add (mat m1, mat m2)
void	imat_add (imat m1, imat m2)
void	cmat_add (cmat m1, cmat m2)
void	mat_sub (mat m1, mat m2)
void	imat_sub (imat m1, imat m2)
void	cmat_sub (cmat m1, cmat m2)
mat	mat_new_add (mat m1, mat m2)
imat	imat_new_add (imat m1, imat m2)
cmat	cmat_new_add (cmat m1, cmat m2)
mat	mat_new_sub (mat m1, mat m2)
imat	imat_new_sub (imat m1, imat m2)
cmat	cmat_new_sub (cmat m1, cmat m2)
mat	mat_new_mul (mat m1, mat m2)
imat	imat_new_mul (imat m1, imat m2)
cmat	cmat_new_mul (cmat m1, cmat m2)
void	mat_mul (mat out, mat m1, mat m2)
void	imat_mul (imat out, imat m1, imat m2)
void	mat_mul_transpose_left (mat out, mat m1, mat m2)
void	mat_mul_transpose_right (mat out, mat m1, mat m2)
void	mat_mul_transpose_leftright (mat out, mat m1, mat m2)
vec	mat_vec_mul (mat m, vec v)
vec	mat_ivec_mul (mat m, ivec v)
vec	imat_vec_mul (imat m, vec v)
ivec	imat_bvec_mul (imat m, bvec v)
ivec	imat_ivec_mul (imat m, ivec v)
cvec	cmat_vec_mul (cmat m, vec v)
cvec	cmat_cvec_mul (cmat m, cvec v)
vec	vec_mat_mul (vec v, mat m)
ivec	ivec_imat_mul (ivec v, imat m)
Random matrices
void	mat_rand (mat m)
void	mat_randn (mat m)
mat	mat_new_rand (idx_t h, idx_t w)
mat	mat_new_randn (idx_t h, idx_t w)