void cauchy_improve_coding_matrix(int k, int m, int w, int *matrix)
{
int index, i, j, x;
int tmp;
int bno, tno, bno_index;
for (j = 0; j < k; j++)
{
if (matrix[j] != 1)
{
tmp = galois_single_divide(1, matrix[j], w);
index = j;
for (i = 0; i < m; i++)
{
matrix[index] = galois_single_multiply(matrix[index], tmp, w);
index += k;
}
}
}
for (i = 1; i < m; i++)
{
bno = 0;
index = i*k;
for (j = 0; j < k; j++) bno += cauchy_n_ones(matrix[index+j], w);
bno_index = -1;
for (j = 0; j < k; j++)
{
if (matrix[index+j] != 1)
{
tmp = galois_single_divide(1, matrix[index+j], w);
tno = 0;
for (x = 0; x < k; x++)
{
tno += cauchy_n_ones(galois_single_multiply(matrix[index+x], tmp, w), w);
}
if (tno < bno)
{
bno = tno;
bno_index = j;
}
}
}
if (bno_index != -1)
{
tmp = galois_single_divide(1, matrix[index+bno_index], w);
for (j = 0; j < k; j++)
{
matrix[index+j] = galois_single_multiply(matrix[index+j], tmp, w);
}
}
}
}
int main(int argc, char **argv)
{
int n;
int i, no, w;
int *bitmatrix;
if (argc != 3) usage(NULL);
if (sscanf(argv[1], "0x%x", &n) == 0) {
if (sscanf(argv[1], "%d", &n) == 0) usage("Bad n");
}
if (sscanf(argv[2], "%d", &w) == 0 || w <= 0 || w > 32) usage("Bad w");
if (w == 31) {
if (n & 0x80000000L) usage("Bad n/w combination (n not between 0 and 2^w-1)\n");
} else if (w < 31) {
if (n >= (1 << w)) usage("Bad n/w combination (n not between 0 and 2^w-1)\n");
}
bitmatrix = jerasure_matrix_to_bitmatrix(1, 1, w, &n);
printf("<HTML><title>cauchy_01 %u %d</title>\n", w, n);
printf("<HTML><h3>cauchy_01 %u %d</h3>\n", w, n);
printf("<pre>\n");
if (w == 32) {
printf("Converted the value 0x%x to the following bitmatrix:\n\n", n);
} else {
printf("Converted the value %d (0x%x) to the following bitmatrix:\n\n", n, n);
}
jerasure_print_bitmatrix(bitmatrix, w, w, w);
printf("\n");
no = 0;
for (i = 0; i < w*w; i++) no += bitmatrix[i];
if (no != cauchy_n_ones(n, w)) {
fprintf(stderr, "Jerasure error: # ones in the bitmatrix (%d) doesn't match cauchy_n_ones() (%d).\n",
no, cauchy_n_ones(n, w));
exit(1);
}
printf("# Ones: %d\n", cauchy_n_ones(n, w));
return 0;
}
int main(int argc, char **argv)
{
int n, i, no, w;
int *bitmatrix;
if (argc != 3) usage(NULL);
if (sscanf(argv[1], "%d", &n) == 0 || n <= 0) usage("Bad n");
if (sscanf(argv[2], "%d", &w) == 0 || w <= 0 || w > 32) usage("Bad w");
if (w < 30 && n >= (1 << w)) usage("n is too big");
bitmatrix = jerasure_matrix_to_bitmatrix(1, 1, w, &n);
no = 0;
for (i = 0; i < w*w; i++) no += bitmatrix[i];
printf("# Ones: %d\n", cauchy_n_ones(n, w));
printf("\n");
printf("Bitmatrix has %d ones\n", no);
printf("\n");
jerasure_print_bitmatrix(bitmatrix, w, w, w);
return 0;
}
/*
* Class: eu_vandertil_jerasure_jni_Cauchy
* Method: cauchy_n_ones
* Signature: (II)I
*/
JNIEXPORT jint JNICALL Java_eu_vandertil_jerasure_jni_Cauchy_cauchy_1n_1ones
(JNIEnv *env, jclass clazz, jint jn, jint jw)
{
return cauchy_n_ones(jn, jw);
}