/dports/math/m4ri/m4ri-release-20200125/m4ri/ |
H A D | brilliantrussian.h | 56 void mzd_make_table(mzd_t const *M, rci_t r, rci_t c, int k, mzd_t *T, rci_t *L); 91 void mzd_process_rows2(mzd_t *M, rci_t startrow, rci_t endrow, rci_t startcol, int k, mzd_t const *… 111 mzd_t const *T0, rci_t const *L0, mzd_t const *T1, rci_t const *L1, 134 mzd_t const *T0, rci_t const *L0, mzd_t const *T1, rci_t const *L1, 135 mzd_t const *T2, rci_t const *L2, mzd_t const *T3, rci_t const *L3); 159 mzd_t const *T0, rci_t const *L0, mzd_t const *T1, rci_t const *L1, 160 mzd_t const *T2, rci_t const *L2, mzd_t const *T3, rci_t const *L3, 254 mzd_t *mzd_inv_m4ri(mzd_t *dst, const mzd_t* src, int k); 272 mzd_t *mzd_mul_m4rm(mzd_t *C, mzd_t const *A, mzd_t const *B, int k); 290 mzd_t *mzd_addmul_m4rm(mzd_t *C, mzd_t const *A, mzd_t const *B, int k); [all …]
|
H A D | strassen.c | 44 mzd_t *_mzd_mul_even(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in _mzd_mul_even() 157 mzd_free_window((mzd_t*)A11); mzd_free_window((mzd_t*)A12); in _mzd_mul_even() 158 mzd_free_window((mzd_t*)A21); mzd_free_window((mzd_t*)A22); in _mzd_mul_even() 160 mzd_free_window((mzd_t*)B11); mzd_free_window((mzd_t*)B12); in _mzd_mul_even() 208 mzd_t *_mzd_sqr_even(mzd_t *C, mzd_t const *A, int cutoff) { in _mzd_sqr_even() 260 mzd_t *Wmk; in _mzd_sqr_even() 341 mzd_t *mzd_mul(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in mzd_mul() 368 mzd_t *_mzd_addmul_even(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in _mzd_addmul_even() 524 mzd_t *_mzd_addsqr_even(mzd_t *C, mzd_t const *A, int cutoff) { in _mzd_addsqr_even() 660 mzd_t *_mzd_addmul(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in _mzd_addmul() [all …]
|
H A D | mp.c | 42 mzd_t *_mzd_addmul_mp4(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in _mzd_addmul_mp4() 142 mzd_free_window((mzd_t*)A00); mzd_free_window((mzd_t*)A01); in _mzd_addmul_mp4() 143 mzd_free_window((mzd_t*)A10); mzd_free_window((mzd_t*)A11); in _mzd_addmul_mp4() 145 mzd_free_window((mzd_t*)B00); mzd_free_window((mzd_t*)B01); in _mzd_addmul_mp4() 146 mzd_free_window((mzd_t*)B10); mzd_free_window((mzd_t*)B11); in _mzd_addmul_mp4() 155 mzd_t *_mzd_mul_mp4(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in _mzd_mul_mp4() 255 mzd_free_window((mzd_t*)A00); mzd_free_window((mzd_t*)A01); in _mzd_mul_mp4() 256 mzd_free_window((mzd_t*)A10); mzd_free_window((mzd_t*)A11); in _mzd_mul_mp4() 258 mzd_free_window((mzd_t*)B00); mzd_free_window((mzd_t*)B01); in _mzd_mul_mp4() 268 mzd_t *mzd_mul_mp(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff) { in mzd_mul_mp() [all …]
|
H A D | mzd.h | 740 mzd_t *mzd_transpose(mzd_t *DST, mzd_t const *A); 756 mzd_t *mzd_mul_naive(mzd_t *C, mzd_t const *A, mzd_t const *B); 772 mzd_t *mzd_addmul_naive(mzd_t *C, mzd_t const *A, mzd_t const *B); 900 mzd_t *mzd_copy(mzd_t *DST, mzd_t const *A); 920 mzd_t *mzd_concat(mzd_t *C, mzd_t const *A, mzd_t const *B); 939 mzd_t *mzd_stack(mzd_t *C, mzd_t const *A, mzd_t const *B); 967 mzd_t *mzd_invert_naive(mzd_t *INV, mzd_t const *A, mzd_t const *I); 980 mzd_t *mzd_add(mzd_t *C, mzd_t const *A, mzd_t const *B); 990 mzd_t *_mzd_add(mzd_t *C, mzd_t const *A, mzd_t const *B); 1317 mzd_t *mzd_extract_u(mzd_t *U, mzd_t const *A); [all …]
|
H A D | strassen.h | 52 mzd_t *mzd_mul(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 68 mzd_t *mzd_addmul(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 88 mzd_t *_mzd_mul_even(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 109 mzd_t *_mzd_addmul_even(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 126 mzd_t *_mzd_addmul(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff);
|
H A D | triangular.h | 50 void mzd_trsm_upper_right(mzd_t const *U, mzd_t *B, const int cutoff); 64 void _mzd_trsm_upper_right(mzd_t const *U, mzd_t *B, const int cutoff); 82 void mzd_trsm_lower_right(mzd_t const *L, mzd_t *B, const int cutoff); 100 void _mzd_trsm_lower_right(mzd_t const *L, mzd_t *B, const int cutoff); 115 void mzd_trsm_lower_left(mzd_t const *L, mzd_t *B, const int cutoff); 127 void _mzd_trsm_lower_left(mzd_t const *L, mzd_t *B, const int cutoff); 142 void mzd_trsm_upper_left(mzd_t const *U, mzd_t *B, const int cutoff); 153 void _mzd_trsm_upper_left (mzd_t const *U, mzd_t *B, const int cutoff); 163 mzd_t *mzd_trtri_upper(mzd_t *A);
|
H A D | mp.h | 47 mzd_t *mzd_mul_mp(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 62 mzd_t *mzd_addmul_mp(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 75 mzd_t *_mzd_addmul_mp4(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff); 87 mzd_t *_mzd_mul_mp4(mzd_t *C, mzd_t const *A, mzd_t const *B, int cutoff);
|
H A D | triangular.c | 40 void _mzd_trsm_upper_right_base(mzd_t const *U, mzd_t *B); 42 void mzd_trsm_upper_right(mzd_t const *U, mzd_t *B, const int cutoff) { in mzd_trsm_upper_right() 51 void _mzd_trsm_upper_right_trtri(mzd_t const *U, mzd_t *B) { in _mzd_trsm_upper_right_trtri() 60 void _mzd_trsm_upper_right(mzd_t const *U, mzd_t *B, const int cutoff) { in _mzd_trsm_upper_right() 112 void _mzd_trsm_upper_right_base(mzd_t const *U, mzd_t *B) { in _mzd_trsm_upper_right_base() 195 void _mzd_trsm_lower_right_base(mzd_t const *L, mzd_t *B); 197 void mzd_trsm_lower_right(mzd_t const *L, mzd_t *B, const int cutoff) { in mzd_trsm_lower_right() 255 void _mzd_trsm_lower_right_base(mzd_t const *L, mzd_t *B) { in _mzd_trsm_lower_right_base() 338 void mzd_trsm_lower_left(mzd_t const *L, mzd_t *B, const int cutoff) { in mzd_trsm_lower_left() 399 void mzd_trsm_upper_left(mzd_t const *U, mzd_t *B, const int cutoff) { in mzd_trsm_upper_left() [all …]
|
H A D | solve.h | 48 int mzd_solve_left(mzd_t *A, mzd_t *B, int const cutoff, int const inconsistency_check); 74 int mzd_pluq_solve_left (mzd_t const *A, rci_t rank, 76 mzd_t *B, int const cutoff, int const inconsistency_check); 102 int _mzd_pluq_solve_left(mzd_t const *A, rci_t rank, 104 mzd_t *B, int const cutoff, int const inconsistency_check); 123 int _mzd_solve_left(mzd_t *A, mzd_t *B, int const cutoff, int const inconsistency_check); 140 mzd_t *mzd_kernel_left_pluq(mzd_t *A, int const cutoff);
|
H A D | solve.c | 30 int mzd_solve_left(mzd_t *A, mzd_t *B, int const cutoff, int const inconsistency_check) { in mzd_solve_left() 37 int mzd_pluq_solve_left (mzd_t const *A, rci_t rank, in mzd_pluq_solve_left() 50 int _mzd_pluq_solve_left(mzd_t const *A, rci_t rank, in _mzd_pluq_solve_left() 72 mzd_t *Y1 = mzd_init_window(B, 0, 0, rank, B->ncols); in _mzd_pluq_solve_left() 92 mzd_free_window((mzd_t*)H); in _mzd_pluq_solve_left() 97 mzd_free_window((mzd_t*)LU); in _mzd_pluq_solve_left() 120 int _mzd_solve_left(mzd_t *A, mzd_t *B, int const cutoff, int const inconsistency_check) { in _mzd_solve_left() 146 mzd_t *mzd_kernel_left_pluq(mzd_t *A, int const cutoff) { in mzd_kernel_left_pluq() 159 mzd_t *U = mzd_init_window(A, 0, 0, r, r); in mzd_kernel_left_pluq() 161 mzd_t *R = mzd_init(A->ncols, A->ncols - r); in mzd_kernel_left_pluq() [all …]
|
H A D | ple_russian.h | 40 mzd_t *T; /*!< the actual table with 2^k entries */ 82 rci_t _mzd_ple_russian(mzd_t *A, mzp_t *P, mzp_t *Q, int k); 99 rci_t _mzd_pluq_russian(mzd_t *A, mzp_t *P, mzp_t *Q, int k); 123 int _mzd_ple_submatrix(mzd_t *A, 140 mzd_t *_mzd_ple_to_e(mzd_t *E, mzd_t const *A, rci_t r, rci_t c, int k, rci_t *offsets); 153 void _mzd_process_rows_ple_2(mzd_t *M, rci_t startrow, rci_t stoprow, rci_t startcol, int const *k,… 166 void _mzd_process_rows_ple_3(mzd_t *M, rci_t startrow, rci_t stoprow, rci_t startcol, int const *k,… 179 void _mzd_process_rows_ple_4(mzd_t *M, rci_t startrow, rci_t stoprow, rci_t startcol, int const *k,… 192 void _mzd_process_rows_ple_5(mzd_t *M, rci_t startrow, rci_t stoprow, rci_t startcol, int const *k,… 205 void _mzd_process_rows_ple_6(mzd_t *M, rci_t startrow, rci_t stoprow, rci_t startcol, int const *k,…
|
H A D | io.h | 44 void mzd_fprint_row(FILE* stream, mzd_t const *M, const rci_t i); 55 static inline void mzd_print_row(mzd_t const *M, const rci_t i) { in mzd_print_row() 68 static inline void mzd_fprint(FILE* stream, mzd_t const *M) { in mzd_fprint() 82 static inline void mzd_print(mzd_t const *M) { in mzd_print() 95 void mzd_info(const mzd_t *A, int do_rank); 109 mzd_t * mzd_from_png(const char *fn, int verbose); 135 int mzd_to_png(const mzd_t *A, const char *fn, int compression_level, const char *comment, int verb… 178 mzd_t *mzd_from_jcf(const char *fn, int verbose); 198 mzd_t *mzd_from_str(rci_t m, rci_t n, const char *str);
|
H A D | echelonform.c | 30 rci_t mzd_echelonize(mzd_t *A, int full) { in mzd_echelonize() 34 rci_t mzd_echelonize_m4ri(mzd_t *A, int full, int k) { in mzd_echelonize_m4ri() 38 rci_t mzd_echelonize_pluq(mzd_t *A, int full) { in mzd_echelonize_pluq() 46 mzd_t *U = mzd_init_window(A, 0, 0, r, r); in mzd_echelonize_pluq() 47 mzd_t *B = mzd_init_window(A, 0, r, r, A->ncols); in mzd_echelonize_pluq() 58 mzd_t *A0 = mzd_init_window(A, 0, 0, r, A->ncols); in mzd_echelonize_pluq() 67 mzd_t *U = mzd_init_window(A, 0, 0, r, r); in mzd_echelonize_pluq() 72 mzd_t *B = mzd_init_window(A, 0, r, r, A->ncols); in mzd_echelonize_pluq() 95 mzd_t *Bw = mzd_init_window(A, 0, r_radix, r, A->ncols); in mzd_echelonize_pluq() 109 mzd_t *A0 = mzd_init_window(A, 0, 0, r, A->ncols); in mzd_echelonize_pluq() [all …]
|
H A D | mzp.h | 121 void mzd_apply_p_left(mzd_t *A, mzp_t const *P); 132 void mzd_apply_p_left_trans(mzd_t *A, mzp_t const *P); 143 void mzd_apply_p_right(mzd_t *A, mzp_t const *P); 154 void mzd_apply_p_right_trans(mzd_t *A, mzp_t const *P); 167 void mzd_apply_p_right_even_capped(mzd_t *A, mzp_t const *P, rci_t start_row, rci_t start_col); 180 void mzd_apply_p_right_trans_even_capped(mzd_t *A, mzp_t const *P, rci_t start_row, rci_t start_col… 191 void mzd_apply_p_right_trans(mzd_t *A, mzp_t const *P); 202 void mzd_apply_p_right_trans_tri(mzd_t *A, mzp_t const *Q); 222 void _mzd_compress_l(mzd_t *A, rci_t r1, rci_t n1, rci_t r2);
|
/dports/math/m4rie/m4rie-release-20200115/src/ |
H A D | mzd_ptr.h | 64 static inline mzd_t *_mzd_ptr_add_to_all(const gf2e *ff, mzd_t *A, mzd_t **X, const int n, ...) { in _mzd_ptr_add_to_all() 83 static inline void _mzd_ptr_add(mzd_t **c, const mzd_t **a, const mzd_t **b, const deg_t length) { in _mzd_ptr_add() 148 void _mzd_ptr_addmul_karatsuba2(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 154 void _mzd_ptr_addmul_karatsuba3(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 160 void _mzd_ptr_addmul_karatsuba4(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 166 void _mzd_ptr_addmul_karatsuba5(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 172 void _mzd_ptr_addmul_karatsuba6(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 178 void _mzd_ptr_addmul_karatsuba7(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 184 void _mzd_ptr_addmul_karatsuba8(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); 190 void _mzd_ptr_addmul_karatsuba9(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B); [all …]
|
H A D | mzd_poly.h | 41 mzd_t **x; /**< Coefficients. */ 61 _mzd_ptr_add(C->x+offset, (const mzd_t**)A->x, (const mzd_t**)B->x, A->depth); in _mzd_poly_add() 92 A->x = (mzd_t**)m4ri_mm_malloc(sizeof(mzd_t*)*(d+1)); in mzd_poly_init() 184 case 2: _mzd_ptr_addmul_karatsuba2(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() 185 case 3: _mzd_ptr_addmul_karatsuba3(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() 186 case 4: _mzd_ptr_addmul_karatsuba4(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() 187 case 5: _mzd_ptr_addmul_karatsuba5(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() 188 case 6: _mzd_ptr_addmul_karatsuba6(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() 189 case 7: _mzd_ptr_addmul_karatsuba7(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() 190 case 8: _mzd_ptr_addmul_karatsuba8(NULL, C->x, (const mzd_t**)A->x, (const mzd_t**)B->x); break; in _mzd_poly_addmul_karatsubs_balanced() [all …]
|
H A D | blm.h | 52 mzd_t *H; /*!< final linear map H*/ 55 mzd_t *F; /*!< lineatr map on A */ 58 mzd_t *G; /*!< lineatr map on B */ 131 void _mzd_ptr_apply_blm_mzd(mzd_t **X, const mzd_t **A, const mzd_t **B, const blm_t *f); 142 void _mzd_ptr_apply_blm_djb(mzd_t **X, const mzd_t **A, const mzd_t **B, const blm_t *f); 153 static inline void _mzd_ptr_apply_blm(mzd_t **X, const mzd_t **A, const mzd_t **B, const blm_t *f) { in _mzd_ptr_apply_blm()
|
H A D | karatsuba.c | 4 void _mzd_ptr_addmul_karatsuba2(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba2() 25 void _mzd_ptr_addmul_karatsuba3(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba3() 50 void _mzd_ptr_addmul_karatsuba4(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba4() 94 void _mzd_ptr_addmul_karatsuba5(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba5() 162 void _mzd_ptr_addmul_karatsuba6(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba6() 235 void _mzd_ptr_addmul_karatsuba7(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba7() 353 void _mzd_ptr_addmul_karatsuba8(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba8() 457 void _mzd_ptr_addmul_karatsuba9(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba9() 580 void _mzd_ptr_addmul_karatsuba10(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba10() 720 void _mzd_ptr_addmul_karatsuba11(const gf2e *ff, mzd_t **X, const mzd_t **A, const mzd_t **B) { in _mzd_ptr_addmul_karatsuba11() [all …]
|
H A D | mzd_slice.h | 520 _mzd_ptr_add(C->x, (const mzd_t**)A->x, (const mzd_t**)B->x, A->depth); in _mzd_slice_add() 628 …case 2: _mzd_ptr_addmul_karatsuba2(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 629 …case 3: _mzd_ptr_addmul_karatsuba3(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 630 …case 4: _mzd_ptr_addmul_karatsuba4(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 631 …case 5: _mzd_ptr_addmul_karatsuba5(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 632 …case 6: _mzd_ptr_addmul_karatsuba6(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 633 …case 7: _mzd_ptr_addmul_karatsuba7(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 634 …case 8: _mzd_ptr_addmul_karatsuba8(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 635 …case 9: _mzd_ptr_addmul_karatsuba9(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() 636 …case 10: _mzd_ptr_addmul_karatsuba10(A->finite_field, C->x, (const mzd_t**)A->x, (const mzd_t**)B-… in _mzd_slice_addmul_karatsuba() [all …]
|
H A D | m4ri_functions.h | 40 static inline word __mzd_read_bits(const mzd_t *M, const rci_t x, const rci_t y, const rci_t n) { in __mzd_read_bits() 52 static inline void __mzd_xor_bits(const mzd_t *M, const rci_t x, const rci_t y, const rci_t n, word… in __mzd_xor_bits() 62 static inline void __mzd_clear_bits(const mzd_t *M, const rci_t x, const rci_t y, const rci_t n) { in __mzd_clear_bits() 81 static inline mzd_t *mzd_sum(mzd_t *A, const int n, ...) { in mzd_sum() 86 mzd_add(A, va_arg(b_list, mzd_t *), va_arg(b_list, mzd_t *)); in mzd_sum() 89 mzd_t *B = va_arg(b_list, mzd_t *); in mzd_sum()
|
/dports/math/m4ri/m4ri-release-20200125/tests/ |
H A D | test_transpose.c | 44 mzd_t* A = mzd_init(m, n); in test_transpose() 45 mzd_t* B = mzd_init(m, n); in test_transpose() 48 mzd_t* C = mzd_add(NULL, A, B); in test_transpose() 49 mzd_t* AT = mzd_init(n, m); in test_transpose() 82 mzd_t* A = mzd_init(m, n); in main() 83 mzd_t* B = mzd_init(m, n); in main() 86 mzd_t* C = mzd_add(NULL, A, B); in main() 87 mzd_t* AT = mzd_init(n, m); in main() 90 mzd_t* BT = mzd_transpose(NULL, B); in main() 91 mzd_t* CT = mzd_add(NULL, AT, BT); in main() [all …]
|
H A D | test_trsm.c | 9 mzd_t* Ubase = mzd_init (2048, 2048); in test_trsm_upper_right() 10 mzd_t* Bbase = mzd_init (2048, 2048); in test_trsm_upper_right() 13 mzd_t* Bbasecopy = mzd_copy (NULL, Bbase); in test_trsm_upper_right() 17 mzd_t* W = mzd_copy (NULL, B); in test_trsm_upper_right() 62 mzd_t* Lbase = mzd_init (2048, 2048); in test_trsm_lower_right() 63 mzd_t* Bbase = mzd_init (2048, 2048); in test_trsm_lower_right() 70 mzd_t* W = mzd_copy (NULL, B); in test_trsm_lower_right() 116 mzd_t* Lbase = mzd_init (2048, 2048); in test_trsm_lower_left() 117 mzd_t* Bbase = mzd_init (2048, 2048); in test_trsm_lower_left() 124 mzd_t* W = mzd_copy (NULL, B); in test_trsm_lower_left() [all …]
|
/dports/math/cado-nfs/cado-nfs-f4284e2391121b2bfb97bc4880b6273c7250dc2f/tests/linalg/bblas/ |
H A D | test_bblas_m4ri.cpp | 21 void my_mzd_randomize(mzd_t * A, gmp_randstate_t rstate)/*{{{*/ in my_mzd_randomize() 29 void test_bblas_base::mzd_set_mem(mzd_t * M, const uint64_t * s, unsigned int n)/*{{{*/ in mzd_set_mem() 41 void test_bblas_base::mzd_set_memT(mzd_t * M, const uint64_t * s, unsigned int n)/*{{{*/ in mzd_set_memT() 47 void test_bblas_base::mzd_check_mem(mzd_t * M, uint64_t * s, unsigned int n)/*{{{*/ in mzd_check_mem() 64 void mzd_check_memT(mzd_t * M, uint64_t * s, unsigned int n)/*{{{*/ 66 mzd_t * tmp = mzd_transpose(NULL, M); 82 static inline void mzd_mypluq(mzd_t * LU, mzd_t * M, mzp_t * P, mzp_t * Q, int c) in mzd_mypluq() 87 static inline void mzd_myechelonize_m4ri(mzd_t * E, mzd_t * M, int full, int k) in mzd_myechelonize_m4ri() 92 static inline void mzd_myechelonize_pluq(mzd_t * E, mzd_t * M, int full) in mzd_myechelonize_pluq() 103 mzd_t * M; in m4ri_plu_tests() [all …]
|
H A D | test_bblas_base.hpp | 39 mzd_t *R; 40 mzd_t *A; 41 mzd_t *R64; 42 mzd_t *A64; 43 mzd_t *W; 44 mzd_t *WT; 47 static void mzd_set_mem(mzd_t * M, const uint64_t * s, unsigned int n); 48 static void mzd_set_memT(mzd_t * M, const uint64_t * s, unsigned int n); 49 static void mzd_check_mem(mzd_t * M, uint64_t * s, unsigned int n);
|
/dports/math/m4ri/m4ri-release-20200125/bench/ |
H A D | bench_mzd.c | 125 mzd_t* volatile vA; 138 mzd_t* const A = mzd_init(64, 64); in BENCHMARK_PREFIX() 150 mzd_t* const A = mzd_init(p->m, p->n); in BENCHMARK_PREFIX() 168 mzd_t* const A = mzd_init(p->m, p->n); in BENCHMARK_PREFIX() 182 mzd_t* const A = mzd_init(p->m, p->n); in BENCHMARK_PREFIX() 364 mzd_t* const V = mzd_init(1, p->m); in BENCHMARK_PREFIX() 365 mzd_t* const C = mzd_init(1, p->n); in BENCHMARK_PREFIX() 381 mzd_t** A = malloc(sizeof(mzd_t) * (p->count + 1)); in BENCHMARK_PREFIX() 406 mzd_t** A = malloc(sizeof(mzd_t) * (p->count + 1)); in BENCHMARK_PREFIX() 432 mzd_t* const B = mzd_copy(NULL, A); in BENCHMARK_PREFIX() [all …]
|