int run_mzed_add(void *_p, unsigned long long *data, int *data_len) {
struct smallops_params *p = (struct smallops_params *)_p;
*data_len = 2;
gf2e *ff = gf2e_init(irreducible_polynomials[p->k][1]);
mzed_t *A = mzed_init(ff,p->m,p->n);
mzed_randomize(A);
mzed_t *B = mzed_init(ff,p->m,p->n);
mzed_randomize(B);
mzed_t *C = mzed_init(ff,p->m,p->n);
data[0] = walltime(0);
data[1] = cpucycles();
mzed_add(C, A, B);
data[1] = cpucycles() - data[1];
data[0] = walltime(data[0]);
mzed_free(A);
mzed_free(B);
mzed_free(C);
gf2e_free(ff);
return 0;
}
void generate_matrix_MS_inverse()
{
uint16_t int Irred[9]; /*Irreducible polynomials */
Irred[0] = 1; Irred[1] = 1; Irred[2] = 0; Irred[3] = 1;
Irred[4] = 1; Irred[5] = 0; Irred[6] = 0; Irred[7] = 0;
Irred[8] = 1;
gf2e *ff;
int k;
for(k = 2; k <= 8; k++)
{
ff = gf2e_init(Irred[k]);
}
const size_t mat_dim = m;
int i,j;
mzed_t *S = mzed_init(ff, ,mat_dim, mat_dim);
for(i = 0; i < m; i++)
{
for(j = 0; j < m; j++)
{
mzed_add_elem(S, i, j, matrix_MS[i][j]);
}
}
mzed_t *S_INV = mzed_init(ff, mat_dim, mat_dim);
mzed_invert_newton_john(S_INV, S);
mzed_free(S);
matrix_MS_inverse = malloc(m*sizeof(int*));
for(i = 0; i < m; i++)
{
matrix_MS_inverse[i] = malloc((m+1)*sizeof(int));
}
for(i = 0; i < m; i++)
{
for(j = 0; j < m; j++)
{
matrix_MS_inverse[i][j] = mzed_read_elem(S_INV, i, j);
}
matrix_MS_inverse[i][m] = matrix_MS[i][m];
}
mzed_free(S_INV);
gf2e_free(ff);
}
void generate_matrix_MT_inverse()
{
uint16_t int Irred[9];
Irred[0] = 1; Irred[1] = 1; Irred[2] = 0; Irred[3] = 1;
Irred[4] = 1; Irred[5] = 0; Irred[6] = 0; Irred[7] = 0;
Irred[8] = 1;
gf2e *ff;
int k;
for(k = 2; k <= 8; k++)
{
ff = gf2e_init(Irred[k]);
}
const size_t mat_dim = n;
mzed_t *T = mzed_init(ff, mat_dim, mat_dim);
int i, j;
for(i = 0; i < n; i++)
{
for(j = 0; j < n; j++)
{
mzed_add_elem(T, i, j, matrix_MT[i][j]);
}
}
mzed_t *T_INV = mzed_init(ff,mat_dim, mat_dim);
mzed_invert_newton_john(T_INV, T);
mzed_free(T);
matrix_MT_inverse = malloc(n*sizeof(int*));
for(i = 0; i < n; i++)
{
matrix_MT_inverse[i] = malloc((n+1)*int);
}
for(i = 0; i < n; i++)
{
for(j = 0; j < n; j++)
{
matrix_MT_inverse[i][j] = mzed_read_elem(T_INV, i, j);
}
matrix_MT_inverse[i][n] = matrix_MT[i][n];
}
mzed_free(T_INV);
gf2e_free(ff);
}
void generate_matrix_MS22_inverse()
{
uint16_t int Irred[9];
Irred[0] = 1; Irred[1] = 1; Irred[2] = 0; Irred[3] = 1;
Irred[4] = 1; Irred[5] = 0; Irred[6] = 0; Irred[7] = 0;
Irred[8] = 1;
gf2e *ff;
int k;
for(k = 2; k <= 8; k++)
{
ff = gf2e_init(Irred[k]);
}
const size_t mat_dim = o2;
mzed_t *S22 = mzed_init(ff, mat_dim, mat_dim);
int i, j;
for(i = 0; i < o2; i++)
{
for(j = 0; j < o2; j++)
{
/* o1 + i >= o2 || o1 + j >= o2 .*/
mzed_add_elem(S22, i, j, matrix_MS[o1+i][o1+j]);
}
}
mzed_t *S22_INV = mzed_init(ff, mat_dim, mat_dim);
mzed_invert_newton_john(S22_INV, S22);
mzed_free(S22);
matrix_MS22_inverse = malloc(o2*sizeof(int*));
for(i = 0; i < o2; i++)
{
matrix_MS22_inverse[i] = malloc(o2*sizeof(int));
}
for(i = 0; i < o2; i++)
{
for(j = 0; j < o2; j++)
{
matrix_MS22_inverse[i][j] = mzed_read_elem(S22_INV, i, j);
}
}
mzed_free(S22_INV);
gf2e_free(ff);
}
void generate_matrix_A_inverse()
{
uint16_t int Irred[9];
Irred[0] = 1; Irred[1] = 1; Irred[2] = 0; Irred[3] = 1;
Irred[4] = 1; Irred[5] = 0; Irred[6] = 0; Irred[7] = 0;
Irred[8] = 1;
int k;
gf2e *ff;
for(k = 2; k <= 8; k++)
{
ff = gf2e_init(Irred[k]);
}
const size_t mat_dim = D2;
mzed_t *A = mzed_init(ff, mat_dim, mat_dim);
int i,j;
for(i = 0; i < D2; i++)
{
for(j = 0; j < D2; j++)
{
mzed_add_elem(A, i, j, matrix_A[i][j]);
}
}
mzed_t *A_INV = mzed_init(ff, mat_dim, mat_dim);
mzed_invert_newton_john(A_INV, A);
matrix_A_inverse_T = malloc(D2*sizeof(int*));
mzed_free(A);
for(i = 0; i < D2; i++)
{
matrix_A_inverse_T[i] = malloc(D2*sizeof(int));
}
for(i = 0; i < D2; i++)
{
for(j = 0; j < D2; j++)
{
matrix_A_inverse_T[j][i] = mzed_read_elem(A_INV, i, j);
}
}
mzed_free(A_INV);
gf2e_free(ff);
}
int test_equality(gf2e *ff, rci_t m, rci_t n) {
int fail_ret = 0;
mzed_t *A0 = random_mzed_t(ff, m, n);
mzed_t *A1 = mzed_copy(NULL, A0);
mzed_t *A2 = mzed_copy(NULL, A0);
mzed_t *A3 = mzed_copy(NULL, A0);
mzed_set_canary(A1);
mzed_set_canary(A2);
mzed_set_canary(A3);
const rci_t r0 = mzed_echelonize_newton_john(A0,1);
const rci_t r1 = mzed_echelonize_naive(A1,1);
const rci_t r2 = mzed_echelonize(A2,1);
const rci_t r3 = mzed_echelonize_ple(A3,1);
m4rie_check( r0 == r1);
m4rie_check( mzed_cmp(A0, A1) == 0);
m4rie_check( r1 == r2);
m4rie_check( mzed_cmp(A1, A2) == 0);
m4rie_check( r2 == r3);
m4rie_check( mzed_cmp(A2, A3) == 0);
m4rie_check( r3 == r0);
m4rie_check( mzed_cmp(A3, A0) == 0);
m4rie_check( mzed_canary_is_alive(A0) );
m4rie_check( mzed_canary_is_alive(A1) );
m4rie_check( mzed_canary_is_alive(A2) );
m4rie_check( mzed_canary_is_alive(A3) );
mzed_free(A0);
mzed_free(A1);
mzed_free(A2);
mzed_free(A3);
return fail_ret;
}
MatrixFn::~MatrixFn()
{
mzed_free(matrix);
gf2e_free(field);
}
