mat.h

Go to the documentation of this file.

/*
 * Copyright (c) 2018 XLAB d.o.o.
 *
 * Licensed under the Apache License, Version 2.0 (the "License");
 * you may not use this file except in compliance with the License.
 * You may obtain a copy of the License at
 *
 * http://www.apache.org/licenses/LICENSE-2.0
 *
 * Unless required by applicable law or agreed to in writing, software
 * distributed under the License is distributed on an "AS IS" BASIS,
 * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
 * See the License for the specific language governing permissions and
 * limitations under the License.
 */
 
#ifndef CIFER_MATRIX_H
#define CIFER_MATRIX_H
 
#include <gmp.h>
 
#include "cifer/internal/errors.h"
#include "cifer/data/vec.h"
 
typedef struct cfe_mat {
    cfe_vec *mat; 
    size_t rows; 
    size_t cols; 
} cfe_mat;
 
void cfe_mat_init(cfe_mat *m, size_t rows, size_t cols);
 
void cfe_mat_inits(size_t rows, size_t cols, cfe_mat *m, ...);
 
void cfe_mat_set_const(cfe_mat *m, mpz_t c);
 
void cfe_mat_copy(cfe_mat *res, cfe_mat *m);
 
void cfe_mat_free(cfe_mat *m);
 
void cfe_mat_frees(cfe_mat *m, ...);
 
void cfe_mat_print(cfe_mat *m);
 
void cfe_mat_add(cfe_mat *res, cfe_mat *m1, cfe_mat *m2);
 
void cfe_mat_neg(cfe_mat *res, cfe_mat *m);
 
void cfe_mat_dot(mpz_t res, cfe_mat *m1, cfe_mat *m2);
 
void cfe_mat_mod(cfe_mat *res, cfe_mat *m, mpz_t modulo);
 
void cfe_mat_mul_vec(cfe_vec *res, cfe_mat *m, cfe_vec *v);
 
void cfe_mat_mul(cfe_mat *res, cfe_mat *m1, cfe_mat *m2);
 
bool cfe_mat_check_bound(cfe_mat *m, mpz_t bound);
 
void cfe_mat_get(mpz_t res, cfe_mat *m, size_t i, size_t j);
 
void cfe_mat_get_col(cfe_vec *res, cfe_mat *m, size_t i);
 
void cfe_mat_get_row(cfe_vec *res, cfe_mat *m, size_t i);
 
cfe_vec *cfe_mat_get_row_ptr(cfe_mat *m, size_t i);
 
void cfe_mat_set(cfe_mat *m, mpz_t el, size_t i, size_t j);
 
void cfe_mat_set_vec(cfe_mat *m, cfe_vec *v, size_t i);
 
void cfe_mat_to_vec(cfe_vec *res, cfe_mat *m);
 
void cfe_mat_from_vec(cfe_mat *m, cfe_vec *v);
 
void cfe_mat_transpose(cfe_mat *t, cfe_mat *m);
 
void cfe_mat_extract_submatrix(cfe_mat *m, cfe_mat *min, size_t i, size_t j);
 
void cfe_mat_determinant(mpz_t det, cfe_mat *m);
 
cfe_error cfe_mat_inverse_mod(cfe_mat *inverse_mat, cfe_mat *m, mpz_t mod);
 
void cfe_mat_mul_x_mat_y(mpz_t res, cfe_mat *mat, cfe_vec *x, cfe_vec *y);
 
void cfe_mat_mul_scalar(cfe_mat *res, cfe_mat *mat, mpz_t s);
 
void cfe_mat_gaussian_elimination(cfe_mat *res, cfe_mat *mat, mpz_t p);
 
cfe_error cfe_mat_inverse_mod_gauss(cfe_mat *res, mpz_t det, cfe_mat *m, mpz_t p);
 
void cfe_mat_determinant_gauss(mpz_t det, cfe_mat *m, mpz_t p);
 
cfe_error cfe_gaussian_elimination_solver(cfe_vec *res, cfe_mat *mat, cfe_vec *vec, mpz_t p);
 
#endif