Initial commit

1 files changed, 49 insertions, 0 deletions

diff --git a/matrix.h b/matrix.h
new file mode 100644
index 0000000..e066777
--- /dev/null
+++ b/matrix.h
@@ -0,0 +1,49 @@
+#ifndef MATRIX_H
+#define MATRIX_H
+
+#include <stddef.h>
+
+#define matrix_at(mat, i, j) ((mat)->xs[(i) * (mat)->n + (j)])
+#define TRUE 1
+#define FALSE 0
+
+typedef enum _MatrixType {
+    MATRIX_DIAG,
+    MATRIX_TRIAG_UPPER,
+    MATRIX_TRIAG_LOWER,
+    MATRIX_TRIAG_SUPPER, /* S for strict */
+    MATRIX_TRIAG_SLOWER,
+    MATRIX_NONE
+} MatrixType;
+
+typedef struct _Matrix {
+    double *xs;
+    size_t m;
+    size_t n;
+} Matrix;
+
+Matrix *matrix_alloc(const size_t m, const size_t n);
+void matrix_free(Matrix *mat);
+void matrix_freen(Matrix **mats);
+void matrix_print(const Matrix *mat);
+Matrix *matrix_from_str(char *str);
+Matrix *matrix_from_arr(double arr[], size_t m, size_t n);
+char *matrix_to_str(const Matrix *mat);
+Matrix *matrix_id(const size_t n);
+Matrix *matrix_const(const size_t m, const size_t n, const double x);
+Matrix *matrix_copy (const Matrix *mat);
+double matrix_trace(const Matrix *mat);
+Matrix *matrix_transpose(const Matrix *mat);
+Matrix *matrix_add(const Matrix *A, const Matrix *B);
+Matrix *matrix_scale(const double x, const Matrix *A);
+Matrix *matrix_sub(const Matrix *A, const Matrix *B);
+Matrix *matrix_mult(const Matrix *A, const Matrix *B);
+Matrix *matrix_rand(const size_t m, const size_t n, const int bound_l, const int bound_u, MatrixType type);
+int matrix_is_square(const Matrix *mat);
+int matrix_eq(const Matrix *A, const Matrix *B, const double tol);
+double matrix_norm_frob(const Matrix *mat);
+Matrix **matrix_LR(const Matrix *A, const Matrix *b);
+Matrix *matrix_forwardel(const Matrix *L, const Matrix *b);
+Matrix *matrix_backsubst(const Matrix *R, const Matrix *y);
+
+#endif