la.h

Go to the documentation of this file.

 
#ifndef LA_H
#define LA_H
 
#include <stddef.h>
#include <assert.h>
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
 
typedef struct la_matrix {
  size_t rows;   
  size_t cols;   
  size_t stride; 
  double* data;  
} la_matrix;
#define la_matrix_element(matrix, row, col) \
  (matrix).data[(row) * (matrix).stride + (col)]
 
la_matrix la_matrix_alloc(size_t rows, size_t cols);
 
void la_matrix_free(la_matrix matrix);
 
void la_matrix_fill(la_matrix matrix, double value);
 
void la_matrix_rand(la_matrix matrix, double minvalue, double maxvalue);
 
void la_matrix_add(la_matrix a, la_matrix b, la_matrix result);
 
void la_matrix_subtract(la_matrix a, la_matrix b, la_matrix result);
void la_matrix_multiply(la_matrix a, la_matrix b, la_matrix result);
 
void la_matrix_copy(la_matrix dest, la_matrix src);
 
void la_matrix_T_copy(la_matrix dest, la_matrix src);
 
la_matrix la_submatrix(la_matrix src, size_t row_start, size_t col_start,
                       size_t new_rows, size_t new_cols);
 
void la_matrix_fprint(la_matrix matrix, char* name, int precision,
                      int digits_before_dot, FILE* fp);
 
#define la_matrix_print(matrix) \
  la_matrix_fprint((matrix), #matrix, 6, 3, stdout)
 
void la_matrix_fprint_simul(la_matrix* matrices, char** names, size_t n,
                            int precision, int digits_before_dot, FILE* fp);
 
void la_matrix_apply_function(la_matrix dest, la_matrix src,
                              double (*f)(double));
 
double la_matrix_sum(la_matrix src);
 
void la_matrix_sum_rows(la_matrix dest, la_matrix src);
 
void la_matrix_multiply_by_const(la_matrix src, const double x);
 
void la_matrix_get_max_index(la_matrix src, size_t* x, size_t* y);
 
void la_matrix_get_min_index(la_matrix src, size_t* x, size_t* y);
 
#endif