matrix.h

 
#ifndef RC_MATRIX_H
#define RC_MATRIX_H
 
#ifdef __cplusplus
extern "C" {
#endif
 
#include <rc/math/vector.h>
 
typedef struct rc_matrix_t{
    int rows;   
    int cols;   
    double** d; 
    int initialized;
} rc_matrix_t;
 
#define RC_MATRIX_INITIALIZER {\
    .rows = 0,\
    .cols = 0,\
    .d = NULL,\
    .initialized = 0}
 
rc_matrix_t rc_matrix_empty(void);
 
int rc_matrix_alloc(rc_matrix_t* A, int rows, int cols);
 
int rc_matrix_free(rc_matrix_t* A);
 
int rc_matrix_zeros(rc_matrix_t* A, int rows, int cols);
 
int rc_matrix_identity(rc_matrix_t* A, int dim);
 
int rc_matrix_random(rc_matrix_t* A, int rows, int cols);
 
int rc_matrix_diagonal(rc_matrix_t* A, rc_vector_t v);
 
int rc_matrix_duplicate(rc_matrix_t A, rc_matrix_t* B);
 
int rc_matrix_print(rc_matrix_t A);
 
int rc_matrix_print_sci(rc_matrix_t A);
 
int rc_matrix_zero_out(rc_matrix_t* A);
 
int rc_matrix_times_scalar(rc_matrix_t* A, double s);
 
int rc_matrix_multiply(rc_matrix_t A, rc_matrix_t B, rc_matrix_t* C);
 
int rc_matrix_left_multiply_inplace(rc_matrix_t A, rc_matrix_t* B);
 
int rc_matrix_right_multiply_inplace(rc_matrix_t* A, rc_matrix_t B);
 
int rc_matrix_add(rc_matrix_t A, rc_matrix_t B, rc_matrix_t* C);
 
int rc_matrix_add_inplace(rc_matrix_t* A, rc_matrix_t B);
 
int rc_matrix_subtract_inplace(rc_matrix_t* A, rc_matrix_t B);
 
int rc_matrix_transpose(rc_matrix_t A, rc_matrix_t* T);
 
int rc_matrix_transpose_inplace(rc_matrix_t* A);
 
 
int   rc_matrix_times_col_vec(rc_matrix_t A, rc_vector_t v, rc_vector_t* c);
 
int   rc_matrix_row_vec_times_matrix(rc_vector_t v, rc_matrix_t A, rc_vector_t* c);
 
int rc_matrix_outer_product(rc_vector_t v1, rc_vector_t v2, rc_matrix_t* A);
 
 
double rc_matrix_determinant(rc_matrix_t A);
 
int rc_matrix_symmetrize(rc_matrix_t* P);
 
 
#ifdef __cplusplus
}
#endif
 
#endif // RC_MATRIX_H