8 typedef struct gs_matrix_t gs_matrix_t;
11 * Allocate a new matrix and init internal data for a matrix of size
12 * row_init X col_init. Matrix cannot grow beyond these init values.
13 * All elements are initially (implicitly) set to 0.
15 gs_matrix_t *gs_new_matrix(int n_init_rows, int n_init_cols);
18 * Free space used by matrix m
20 void gs_delete_matrix(gs_matrix_t *m);
22 void gs_matrix_assure_row_capacity(gs_matrix_t *m, int row, int min_capacity);
24 void gs_matrix_trim_row_capacities(gs_matrix_t *m);
26 void gs_matrix_delete_zero_entries(gs_matrix_t *m);
29 * Sets m[row, col] to val
30 * @p increase If non-zero @p val is added to the existing
32 void gs_matrix_set(gs_matrix_t *m, int row, int col, double val);
35 * Returns the value stored in m[row, col].
37 double gs_matrix_get(const gs_matrix_t *m, int row, int col);
40 * Performs one step of the Gauss-Seidel algorithm
41 * @p m The iteration matrix
42 * @p x The iteration vector
43 * @p a The dimension of the matrix (axa matrix)
45 double gs_matrix_gauss_seidel(const gs_matrix_t *m, double *x, int n);
47 unsigned gs_matrix_get_n_entries(const gs_matrix_t *m);
50 * Dumps the matrix factor*m to the stream @p out.
52 void gs_matrix_dump(const gs_matrix_t *m, int a, int b, FILE *out);
54 int gs_matrix_get_sizeof_allocated_memory(const gs_matrix_t *m);
56 void gs_matrix_export(const gs_matrix_t *m, double *nw, int size);