int main() { int32_t i, j, n; alloc_constants(); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("iszero(a) = %d\n", bvconst_is_zero(a, 4)); bvconst_add_one(a, 4); printf("a+1 = "); bvconst_print(stdout, a, 128); printf("\n"); printf("iszero(a+1) = %d\n", bvconst_is_zero(a, 4)); bvconst_sub_one(a, 4); printf("a+1-1 = "); bvconst_print(stdout, a, 128); printf("\n"); printf("iszero(a+1-1) = %d\n", bvconst_is_zero(a, 4)); bvconst_sub_one(a, 4); printf("a+1-2 = "); bvconst_print(stdout, a, 128); printf("\n"); printf("iszero(a+1-2) = %d\n", bvconst_is_zero(a, 4)); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); printf("iszero(b) = %d\n", bvconst_is_zero(b, 4)); bvconst_add_one(b, 4); printf("b+1 = "); bvconst_print(stdout, b, 128); printf("\n"); printf("iszero(b+1) = %d\n", bvconst_is_zero(b, 4)); bvconst_sub_one(b, 4); bvconst_sub_one(b, 4); printf("b-1 = "); bvconst_print(stdout, b, 128); printf("\n"); printf("iszero(b-1) = %d\n", bvconst_is_zero(b, 4)); printf("c = "); bvconst_print(stdout, c, 128); printf("\n"); printf("iszero(c) = %d\n", bvconst_is_zero(c, 4)); printf("d = "); bvconst_print(stdout, d, 128); printf("\n"); printf("iszero(d) = %d\n", bvconst_is_zero(d, 4)); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); printf("iszero(e) = %d\n", bvconst_is_zero(e, 4)); printf("e := a * b\n"); bvconst_mul2(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n\n"); random_vector(vector, 32); bvconst_set_array(a, vector, 32); printf("a = "); bvconst_print(stdout, a, 32); printf("\n"); random_vector(vector, 73); bvconst_set_array(b, vector, 73); printf("b = "); bvconst_print(stdout, b, 73); printf("\n"); bvconst_concat(c, a, 32, b, 73); printf("conc(a, b) = "); bvconst_print(stdout, c, 105); printf("\n"); bvconst_concat(d, b, 73, a, 32); printf("conc(b, a) = "); bvconst_print(stdout, d, 105); printf("\n\n"); // random_vector(vector, 63); bvconst_set_minus_one(a, 4); // random_vector(vector, 63); bvconst_clear(b, 4); for (n=2; n<=128; n++) { for (i=1; i<n; i++) { printf("--- n = %"PRId32", i = %"PRId32" ---\n", n, i); printf("a = "); bvconst_print(stdout, a, i); printf("\n"); printf("b = "); bvconst_print(stdout, b, n - i); printf("\n"); bvconst_concat(c, a, i, b, n - i); bvconst_concat(d, b, n - i, a, i); printf("conc(a, b) = "); bvconst_print(stdout, c, n); printf("\n"); printf("conc(b, a) = "); bvconst_print(stdout, d, n); printf("\n\n"); } } // exit(0); for (n=20; n>0; n--) { random_vector(vector, 20); bvconst_set_array(a, vector, 20); printf("a = "); bvconst_print(stdout, a, 20); printf("\n"); for (i=0; i<=20-n; i++) { bvconst_extract(b, a, i, i+n); printf("a[%2"PRId32", %2"PRId32") = ", i, i+n); bvconst_print(stdout, b, n); printf("\n"); } printf("\n"); } for (n=50; n>20; n--) { random_vector(vector, 62); bvconst_set_array(a, vector, 62); printf("a = "); bvconst_print(stdout, a, 62); printf("\n"); for (i=0; i<=62-n; i++) { bvconst_extract(b, a, i, i+n); printf("a[%2"PRId32", %2"PRId32") = ", i, i+n); bvconst_print(stdout, b, n); printf("\n"); } printf("\n"); } random_vector(vector, 20); bvconst_set_array(a, vector, 20); printf("a = "); bvconst_print(stdout, a, 20); printf("\n"); bvconst_set_extend(b, 31, a, 20, 0); printf("ext(a, 31, 0) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, 1); printf("ext(a, 31, 1) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, -1); printf("sgnext(a, 31) = "); bvconst_print(stdout, b, 31); printf("\n\n"); random_vector(vector, 20); bvconst_set_array(a, vector, 20); printf("a = "); bvconst_print(stdout, a, 20); printf("\n"); bvconst_set_extend(b, 31, a, 20, 0); printf("ext(a, 31, 0) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, 1); printf("ext(a, 31, 1) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, -1); printf("sgnext(a, 31) = "); bvconst_print(stdout, b, 31); printf("\n\n"); random_vector(vector, 20); bvconst_set_array(a, vector, 20); printf("a = "); bvconst_print(stdout, a, 20); printf("\n"); bvconst_set_extend(b, 31, a, 20, 0); printf("ext(a, 31, 0) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, 1); printf("ext(a, 31, 1) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, -1); printf("sgnext(a, 31) = "); bvconst_print(stdout, b, 31); printf("\n\n"); random_vector(vector, 20); bvconst_set_array(a, vector, 20); printf("a = "); bvconst_print(stdout, a, 20); printf("\n"); bvconst_set_extend(b, 31, a, 20, 0); printf("ext(a, 31, 0) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, 1); printf("ext(a, 31, 1) = "); bvconst_print(stdout, b, 31); printf("\n"); bvconst_set_extend(b, 31, a, 20, -1); printf("sgnext(a, 31) = "); bvconst_print(stdout, b, 31); printf("\n\n"); for (i=20; i<=99; i++) { random_vector(vector, 20); bvconst_set_array(a, vector, 20); printf("a = "); bvconst_print(stdout, a, 20); printf("\n"); bvconst_set_extend(b, i, a, 20, 0); printf("ext(a, %2"PRId32", 0) = ", i); bvconst_print(stdout, b, i); printf("\n"); bvconst_set_extend(b, i, a, 20, 1); printf("ext(a, %2"PRId32", 1) = ", i); bvconst_print(stdout, b, i); printf("\n"); bvconst_set_extend(b, i, a, 20, -1); printf("sgnext(a, %2"PRId32") = ", i); bvconst_print(stdout, b, i); printf("\n\n"); } for (i=32; i<=99; i++) { random_vector(vector, 32); bvconst_set_array(a, vector, 32); printf("a = "); bvconst_print(stdout, a, 32); printf("\n"); bvconst_set_extend(b, i, a, 32, 0); printf("ext(a, %2"PRId32", 0) = ", i); bvconst_print(stdout, b, i); printf("\n"); bvconst_set_extend(b, i, a, 32, 1); printf("ext(a, %2"PRId32", 1) = ", i); bvconst_print(stdout, b, i); printf("\n"); bvconst_set_extend(b, i, a, 32, -1); printf("sgnext(a, %2"PRId32") = ", i); bvconst_print(stdout, b, i); printf("\n\n"); } random_vector(vector, 120); bvconst_set_array(a, vector, 120); printf("a = "); bvconst_print(stdout, a, 120); printf("\n"); for (i=0; i<=120; i++) { printf("rshift %3"PRId32": ", i); bvconst_shift_right(a, 120, i, 0); bvconst_print(stdout, a, 120); printf("\n"); bvconst_set_array(a, vector, 120); } printf("\n"); random_vector(vector, 120); bvconst_set_array(a, vector, 120); printf("a = "); bvconst_print(stdout, a, 120); printf("\n"); for (i=0; i<=120; i++) { printf("rshift %3"PRId32": ", i); bvconst_shift_right(a, 120, i, 1); bvconst_print(stdout, a, 120); printf("\n"); bvconst_set_array(a, vector, 120); } printf("\n"); random_vector(vector, 120); bvconst_set_array(a, vector, 120); printf("a = "); bvconst_print(stdout, a, 120); printf("\n"); for (i=0; i<=120; i++) { printf("lshift %3"PRId32": ", i); bvconst_shift_left(a, 120, i, 0); bvconst_print(stdout, a, 120); printf("\n"); bvconst_set_array(a, vector, 120); } printf("\n"); random_vector(vector, 120); bvconst_set_array(a, vector, 120); printf("a = "); bvconst_print(stdout, a, 120); printf("\n"); for (i=0; i<=120; i++) { printf("lshift %3"PRId32": ", i); bvconst_shift_left(a, 120, i, 1); bvconst_print(stdout, a, 120); printf("\n"); bvconst_set_array(a, vector, 120); } printf("\n"); // exit(0); printf("c = "); bvconst_print(stdout, c, 128); printf("\n"); printf("d = "); bvconst_print(stdout, d, 128); printf("\n"); printf("e := c * d\n"); bvconst_mul2(e, 4, c, d); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); printf("e := d * c\n"); bvconst_mul2(e, 4, d, c); printf("e = "); bvconst_print(stdout, e, 128); printf("\n\n"); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); printf("d = "); bvconst_print(stdout, d, 128); printf("\n"); printf("e := b * d\n"); bvconst_mul2(e, 4, b, d); printf("e = "); bvconst_print(stdout, e, 128); printf("\n\n"); bvconst_set64(a, 4, 1372919719782793ULL); bvconst_set32(b, 4, 12670371); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); printf("e := a * b\n"); bvconst_set(c, 4, a); printf("c = "); bvconst_print(stdout, c, 128); printf("\n"); bvconst_mul(c, 4, b); printf("c = "); bvconst_print(stdout, c, 128); printf("\n"); bvconst_mul2(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n\n"); mpz_init(z0); mpz_init(z1); mpz_init(z2); bvconst_get_mpz(a, 4, z0); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(b, 4, z0); printf("--> conversion b to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(e, 4, z0); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); printf("e -= a * b\n"); bvconst_submul(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_get_mpz(e, 4, z0); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); for (i=0; i<10; i++) { printf("test %"PRId32": e *= a\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); random_vector(vector, 128); bvconst_set_array(e, vector, 128); bvconst_get_mpz(e, 4, z1); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_mul(e, 4, a); printf("e' = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_get_mpz(a, 4, z0); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); bvconst_get_mpz(e, 4, z2); printf("--> conversion e' to mpz: "); mpz_out_str(stdout, 10, z2); printf("\n"); mpz_mul(z0, z0, z1); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z2) == 0); } for (i=0; i<10; i++) { printf("test %"PRId32", a := - a\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); bvconst_get_mpz(a, 4, z0); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); bvconst_negate(a, 4); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(a, 4, z1); printf("--> conversion -a to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); mpz_neg(z0, z0); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z1) == 0); } for (i=0; i<10; i++) { printf("test %"PRId32": e := - a\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); bvconst_negate2(e, 4, a); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); printf("--> conversion a to mpz: "); bvconst_get_mpz(a, 4, z0); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(e, 4, z1); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); mpz_neg(z0, z0); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z1) == 0); } for (i=0; i<10; i++) { printf("test %"PRId32", e = a * b\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); random_vector(vector, 128); bvconst_set_array(b, vector, 128); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); bvconst_mul2(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_get_mpz(a, 4, z0); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(b, 4, z1); printf("--> conversion b to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); bvconst_get_mpz(e, 4, z2); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z2); printf("\n"); mpz_mul(z0, z0, z1); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z2) == 0); } for (i=0; i<10; i++) { printf("test %"PRId32": e = - (a * b)\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); random_vector(vector, 128); bvconst_set_array(b, vector, 128); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); bvconst_clear(e, 4); bvconst_submul(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_get_mpz(a, 4, z0); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(b, 4, z1); printf("--> conversion b to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); bvconst_get_mpz(d, 4, z2); printf("--> conversion d to mpz: "); mpz_out_str(stdout, 10, z2); printf("\n"); bvconst_get_mpz(e, 4, z2); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z2); printf("\n"); mpz_mul(z0, z0, z1); mpz_neg(z0, z0); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z2) == 0); } for (i=0; i<10; i++) { printf("test %"PRId32": e = a + b\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); random_vector(vector, 128); bvconst_set_array(b, vector, 128); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); bvconst_add2(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_get_mpz(a, 4, z0); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(b, 4, z1); printf("--> conversion b to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); bvconst_get_mpz(e, 4, z2); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z2); printf("\n"); mpz_add(z0, z0, z1); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z2) == 0); } for (i=0; i<10; i++) { printf("test %"PRId32": e = a - b\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); random_vector(vector, 128); bvconst_set_array(b, vector, 128); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); printf("b = "); bvconst_print(stdout, b, 128); printf("\n"); bvconst_sub2(e, 4, a, b); printf("e = "); bvconst_print(stdout, e, 128); printf("\n"); bvconst_get_mpz(a, 4, z0); printf("--> conversion a to mpz: "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_get_mpz(b, 4, z1); printf("--> conversion b to mpz: "); mpz_out_str(stdout, 10, z1); printf("\n"); bvconst_get_mpz(e, 4, z2); printf("--> conversion e to mpz: "); mpz_out_str(stdout, 10, z2); printf("\n"); mpz_sub(z0, z0, z1); mpz_fdiv_r_2exp(z0, z0, 128); printf("--> check: "); mpz_out_str(stdout, 10, z0); printf("\n\n"); assert(mpz_cmp(z0, z2) == 0); } for (i=0; i<10; i++) { printf("\ntest %"PRId32": comparisons\n", i); random_vector(vector, 25); bvconst_set_array(a, vector, 25); random_vector(vector, 25); bvconst_set_array(b, vector, 25); printf("a = "); bvconst_print(stdout, a, 25); printf("\n"); printf("b = "); bvconst_print(stdout, b, 25); printf("\n"); bvconst_normalize(a, 25); bvconst_normalize(b, 25); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 25)); printf(" a ge b: %d\n", bvconst_ge(a, b, 25)); printf(" a lt b: %d\n", bvconst_lt(a, b, 25)); printf(" a gt b: %d\n", bvconst_gt(a, b, 25)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 25)); printf(" a sge b: %d\n", bvconst_sge(a, b, 25)); printf(" a slt b: %d\n", bvconst_slt(a, b, 25)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 25)); } for (i=0; i<10; i++) { printf("\ntest %"PRId32": comparisons\n", i); random_vector(vector, 32); bvconst_set_array(a, vector, 32); random_vector(vector, 32); bvconst_set_array(b, vector, 32); printf("a = "); bvconst_print(stdout, a, 32); printf("\n"); printf("b = "); bvconst_print(stdout, b, 32); printf("\n"); bvconst_normalize(a, 32); bvconst_normalize(b, 32); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 32)); printf(" a ge b: %d\n", bvconst_ge(a, b, 32)); printf(" a lt b: %d\n", bvconst_lt(a, b, 32)); printf(" a gt b: %d\n", bvconst_gt(a, b, 32)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 32)); printf(" a sge b: %d\n", bvconst_sge(a, b, 32)); printf(" a slt b: %d\n", bvconst_slt(a, b, 32)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 32)); } for (i=0; i<10; i++) { printf("\ntest %"PRId32": comparisons\n", i); random_vector(vector, 33); bvconst_set_array(a, vector, 33); random_vector(vector, 33); bvconst_set_array(b, vector, 33); printf("a = "); bvconst_print(stdout, a, 33); printf("\n"); printf("b = "); bvconst_print(stdout, b, 33); printf("\n"); bvconst_normalize(a, 33); bvconst_normalize(b, 33); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 33)); printf(" a ge b: %d\n", bvconst_ge(a, b, 33)); printf(" a lt b: %d\n", bvconst_lt(a, b, 33)); printf(" a gt b: %d\n", bvconst_gt(a, b, 33)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 33)); printf(" a sge b: %d\n", bvconst_sge(a, b, 33)); printf(" a slt b: %d\n", bvconst_slt(a, b, 33)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 33)); } for (i=0; i<10; i++) { printf("\ntest %"PRId32": comparisons\n", i); random_vector(vector, 63); bvconst_set_array(a, vector, 63); random_vector(vector, 63); bvconst_set_array(b, vector, 63); printf("a = "); bvconst_print(stdout, a, 63); printf("\n"); printf("b = "); bvconst_print(stdout, b, 63); printf("\n"); bvconst_normalize(a, 63); bvconst_normalize(b, 63); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 63)); printf(" a ge b: %d\n", bvconst_ge(a, b, 63)); printf(" a lt b: %d\n", bvconst_lt(a, b, 63)); printf(" a gt b: %d\n", bvconst_gt(a, b, 63)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 63)); printf(" a sge b: %d\n", bvconst_sge(a, b, 63)); printf(" a slt b: %d\n", bvconst_slt(a, b, 63)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 63)); } for (i=0; i<10; i++) { printf("\ntest %"PRId32": comparisons\n", i); random_vector(vector, 64); bvconst_set_array(a, vector, 64); random_vector(vector, 64); bvconst_set_array(b, vector, 64); printf("a = "); bvconst_print(stdout, a, 64); printf("\n"); printf("b = "); bvconst_print(stdout, b, 64); printf("\n"); bvconst_normalize(a, 64); bvconst_normalize(b, 64); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 64)); printf(" a ge b: %d\n", bvconst_ge(a, b, 64)); printf(" a lt b: %d\n", bvconst_lt(a, b, 64)); printf(" a gt b: %d\n", bvconst_gt(a, b, 64)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 64)); printf(" a sge b: %d\n", bvconst_sge(a, b, 64)); printf(" a slt b: %d\n", bvconst_slt(a, b, 64)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 64)); } printf("\nMore comparisons\n"); random_vector(vector, 64); bvconst_set_array(a, vector, 64); bvconst_set_array(b, vector, 64); printf("a = "); bvconst_print(stdout, a, 64); printf("\n"); printf("b = "); bvconst_print(stdout, b, 64); printf("\n"); bvconst_normalize(a, 64); bvconst_normalize(b, 64); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 64)); printf(" a ge b: %d\n", bvconst_ge(a, b, 64)); printf(" a lt b: %d\n", bvconst_lt(a, b, 64)); printf(" a gt b: %d\n", bvconst_gt(a, b, 64)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 64)); printf(" a sge b: %d\n", bvconst_sge(a, b, 64)); printf(" a slt b: %d\n", bvconst_slt(a, b, 64)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 64)); random_vector(vector, 65); bvconst_set_array(a, vector, 65); bvconst_set_array(b, vector, 65); printf("a = "); bvconst_print(stdout, a, 65); printf("\n"); printf("b = "); bvconst_print(stdout, b, 65); printf("\n"); bvconst_normalize(a, 65); bvconst_normalize(b, 65); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 65)); printf(" a ge b: %d\n", bvconst_ge(a, b, 65)); printf(" a lt b: %d\n", bvconst_lt(a, b, 65)); printf(" a gt b: %d\n", bvconst_gt(a, b, 65)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 65)); printf(" a sge b: %d\n", bvconst_sge(a, b, 65)); printf(" a slt b: %d\n", bvconst_slt(a, b, 65)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 65)); random_vector(vector, 63); bvconst_set_array(a, vector, 63); bvconst_set_array(b, vector, 63); printf("a = "); bvconst_print(stdout, a, 63); printf("\n"); printf("b = "); bvconst_print(stdout, b, 63); printf("\n"); bvconst_normalize(a, 63); bvconst_normalize(b, 63); printf("Unsigned tests\n"); printf(" a le b: %d\n", bvconst_le(a, b, 63)); printf(" a ge b: %d\n", bvconst_ge(a, b, 63)); printf(" a lt b: %d\n", bvconst_lt(a, b, 63)); printf(" a gt b: %d\n", bvconst_gt(a, b, 63)); printf("Signed tests\n"); printf(" a sle b: %d\n", bvconst_sle(a, b, 63)); printf(" a sge b: %d\n", bvconst_sge(a, b, 63)); printf(" a slt b: %d\n", bvconst_slt(a, b, 63)); printf(" a sgt b: %d\n", bvconst_sgt(a, b, 63)); printf("\nTest powers of two\n"); mpz_set_ui(z0, 1); for (i=0; i<128; i++) { printf("z0 = "); mpz_out_str(stdout, 10, z0); printf("\n"); bvconst_set_mpz(a, 4, z0); printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } mpz_add(z0, z0, z0); } a[0] = 0; a[1] = 1024; a[2] = 0; a[3] = 36136287; printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } a[0] = 0; a[1] = 1024; a[2] = 31209; a[3] = 0; printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } a[0] = 0; a[1] = 1024; a[2] = 0; a[3] = 0; printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } a[0] = 128; a[1] = 1024; a[2] = 0; a[3] = 0; printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } bvconst_clear(a, 4); for (i=0; i<256; i++) { printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } bvconst_add_one(a, 4); } for (i=0; i<10; i++) { printf("test %"PRId32": power of two\n", i); random_vector(vector, 128); bvconst_set_array(a, vector, 128); for (j=0; j<64; j++) { printf("a = "); bvconst_print(stdout, a, 128); printf("\n"); n = bvconst_is_power_of_two(a, 4); if (n < 0) { printf("--> not a power of 2\n\n"); } else { printf("--> a = 2^%"PRId32"\n\n", n); } bvconst_clr_bit(a, j); bvconst_clr_bit(a, 127 - j); } } // test mulpower for (i=0; i<130; i++) { bvconst_set_one(a, 4); // a = 1 bvconst_set32(b, 4, 2); // b = 2 // compute a * b^i = 2^i bvconst_mulpower(a, 4, b, i); printf("---> mulpower: a = b^%"PRIu32" = 2^%"PRIu32"\n", i, i); printf(" b = "); bvconst_print(stdout, b, 128); printf("\n"); printf(" a = "); bvconst_print(stdout, a, 128); printf("\n"); } // test mulpower for (i=0; i<40; i++) { bvconst_set32(a, 4, 4); // a = 0b0..00100 bvconst_set_minus_one(b, 4); // b = 0b11..1111 // compute a * b^i = -4 or +4 bvconst_mulpower(a, 4, b, i); printf("---> mulpower: a = 4 * b^%"PRIu32" = 4 * (-1)^%"PRIu32"\n", i, i); printf(" b = "); bvconst_print(stdout, b, 128); printf("\n"); printf(" a = "); bvconst_print(stdout, a, 128); printf("\n"); } // test mulpower for (i=0; i<40; i++) { bvconst_set32(a, 4, 4); // a = 0b0..00100 bvconst_clear(b, 4); // b = 0 // compute a * b^i = 0 bvconst_mulpower(a, 4, b, i); printf("---> mulpower: a = 4 * b^%"PRIu32"\n", i); printf(" b = "); bvconst_print(stdout, b, 128); printf("\n"); printf(" a = "); bvconst_print(stdout, a, 128); printf("\n"); } // test is_one and is_minus_one printf("\n\n"); for (i=1; i<=128; i++) { bvconst_set_one(a, 4); // a = 1 bvconst_set_minus_one(b, 4); // b = -1 bvconst_clear(c, 4); // c = 0 random_vector(vector, i); bvconst_set_array(d, vector, i); // d = random // make all i-bit length constant bvconst_normalize(a, i); bvconst_normalize(b, i); bvconst_normalize(c, i); bvconst_normalize(d, i); n = (i+31) >> 5; // number of words printf("---> %"PRIu32" bits\n", i); printf("a = "); bvconst_print(stdout, a, i); printf(": is_one = %s, is_minus_one = %s\n", b2str(bvconst_is_one(a, n)), b2str(bvconst_is_minus_one(a, i))); printf("b = "); bvconst_print(stdout, b, i); printf(": is_one = %s, is_minus_one = %s\n", b2str(bvconst_is_one(b, n)), b2str(bvconst_is_minus_one(b, i))); printf("c = "); bvconst_print(stdout, c, i); printf(": is_one = %s, is_minus_one = %s\n", b2str(bvconst_is_one(c, n)), b2str(bvconst_is_minus_one(c, i))); printf("d = "); bvconst_print(stdout, d, i); printf(": is_one = %s, is_minus_one = %s\n", b2str(bvconst_is_one(d, n)), b2str(bvconst_is_minus_one(d, i))); printf("\n"); fflush(stdout); } // test is_min_signed and is_max_signed printf("\n\n"); for (i=1; i<=128; i++) { bvconst_clear(a, 4); // a = 0 bvconst_set_minus_one(b, 4); // b = -1 bvconst_clear(c, 4); bvconst_set_bit(c, i-1); // c = 0b10000000 (smallest negative integer) bvconst_set_minus_one(d, 4); bvconst_clr_bit(d, i-1); // d = 0b01111111 (largest positive integer) random_vector(vector, i); bvconst_set_array(e, vector, i); // e = random // make all i-bit length constant bvconst_normalize(a, i); bvconst_normalize(b, i); bvconst_normalize(c, i); bvconst_normalize(d, i); bvconst_normalize(e, i); printf("---> %"PRIu32" bits\n", i); printf("a = "); bvconst_print(stdout, a, i); printf(": is_min_signed = %s, is_max_signed = %s\n", b2str(bvconst_is_min_signed(a, i)), b2str(bvconst_is_max_signed(a, i))); printf("b = "); bvconst_print(stdout, b, i); printf(": is_min_signed = %s, is_max_signed = %s\n", b2str(bvconst_is_min_signed(b, i)), b2str(bvconst_is_max_signed(b, i))); printf("c = "); bvconst_print(stdout, c, i); printf(": is_min_signed = %s, is_max_signed = %s\n", b2str(bvconst_is_min_signed(c, i)), b2str(bvconst_is_max_signed(c, i))); printf("d = "); bvconst_print(stdout, d, i); printf(": is_min_signed = %s, is_max_signed = %s\n", b2str(bvconst_is_min_signed(d, i)), b2str(bvconst_is_max_signed(d, i))); printf("e = "); bvconst_print(stdout, e, i); printf(": is_min_signed = %s, is_max_signed = %s\n", b2str(bvconst_is_min_signed(e, i)), b2str(bvconst_is_max_signed(e, i))); printf("\n"); fflush(stdout); } for (i=100; i<256; i += 10) { for (j=1; j<=i; j += 9) { test_set_extend(i, j); } } mpz_clear(z0); mpz_clear(z1); mpz_clear(z2); free_constants(); return 0; }
void separation(const char* filename, double background_integral, double* stream_integrals) { int q[ap->number_streams]; double nstars[ap->number_streams]; int total; double sprob[ap->number_streams]; double prob_s[ap->number_streams]; double prob_b; double pbx; double psg[ap->number_streams]; double d; int twoPanel; double** cmatrix; double dnormal[3]; double dortho[3]; double xsun[3]; double epsilon_s[ap->number_streams]; double epsilon_b; double star_coords[3]; double starxyz[3]; double starxyzTransform[3]; int s_ok = 0; int i, j, retval; FILE* file; double reff_xr_rp3, *qw_r3_N, *r_point; twoPanel = 1; for (j = 0; j < ap->number_streams; j++) { nstars[j] = 0; q[j] = 0; } total = 0; prob_ok_init(); printf("Integral complete.\n Beginning probability calculations...\n"); file = fopen(filename, "w"); if (ap->sgr_coordinates == 0) { stripe_normal(ap->wedge, dnormal); } else if (ap->sgr_coordinates == 1) { sgr_stripe_normal(ap->wedge, dnormal); } else { printf("Error: ap->sgr_coordinates not valid"); } free_star_points(sp); free(sp); sp = (STAR_POINTS*)malloc(sizeof(STAR_POINTS)); retval = read_star_points(star_points_file, sp); if (retval) { fprintf(stderr, "APP: error reading star points: %d\n", retval); exit(1); } printf("read %d stars.\n", sp->number_stars); cmatrix = (double**)malloc(sizeof(double*) * 3); for (i = 0; i < 3; i++) cmatrix[i] = (double*)malloc(sizeof(double) * 3); dortho[0] = 0.0; dortho[1] = 0.0; dortho[2] = 1.0; get_transform(dnormal, dortho, cmatrix); printf("\nTransformation matrix:\n"); printf("\t%lf %lf %lf\n", cmatrix[0][0], cmatrix[0][1], cmatrix[0][2]); printf("\t%lf %lf %lf\n", cmatrix[1][0], cmatrix[1][1], cmatrix[1][2]); printf("\t%lf %lf %lf\n\n", cmatrix[2][0], cmatrix[2][1], cmatrix[2][2]); xsun[0] = -8.5; xsun[1] = 0.0; xsun[2] = 0.0; d = dotp(dnormal, xsun); printf("==============================================\n"); printf("bint: %lf", background_integral); for (j = 0; j < ap->number_streams; j++) { printf(", "); printf("sint[%d]: %lf", j, stream_integrals[j]); } printf("\n"); /*get stream & background weight constants*/ double denom = 1.0; for (j = 0; j < ap->number_streams; j++) { denom += exp(ap->stream_weights[j]); } for (j = 0; j < ap->number_streams; j++) { epsilon_s[j] = exp(ap->stream_weights[j]) / denom; printf("epsilon_s[%d]: %lf\n", j, epsilon_s[j]); } epsilon_b = 1.0 / denom; printf("epsilon_b: %lf\n", epsilon_b); r_point = (double*)malloc(sizeof(double) * ap->convolve); qw_r3_N = (double*)malloc(sizeof(double) * ap->convolve); init_constants(ap); printf("initialized constants\n"); for (i = 0; i < sp->number_stars; i++) { MW_DEBUG("[%d/%d] setting star coords\n", i, sp->number_stars); star_coords[0] = sp->stars[i][0]; star_coords[1] = sp->stars[i][1]; star_coords[2] = sp->stars[i][2]; MW_DEBUG("star_coords: %g %g %g\n", star_coords[0], star_coords[1], star_coords[2]); MW_DEBUG("twoPanel: %d\n", twoPanel); if (twoPanel == 1) { MW_DEBUG("setting probability constants\n"); set_probability_constants(ap->convolve, star_coords[2], r_point, qw_r3_N, &reff_xr_rp3); MW_DEBUG("calculating probabilities\n"); calculate_probabilities(r_point, qw_r3_N, reff_xr_rp3, star_coords, ap, &prob_b, prob_s); MW_DEBUG("calculated probabilities\n"); MW_DEBUG("prob_s: %lf\n", prob_s[0]); MW_DEBU("prob_b: %lf\n", prob_b); pbx = epsilon_b * prob_b / background_integral; for (j = 0; j < ap->number_streams; j++) { psg[j] = epsilon_s[j] * prob_s[j] / stream_integrals[j]; } MW_DEBUG("pbx: %g\n", pbx); MW_DEBUG("psg: %g\n", psg[0]); double psgSum = 0; for (j = 0; j < ap->number_streams; j++) { psgSum += psg[j]; } for (j = 0; j < ap->number_streams; j++) { sprob[j] = psg[j] / (psgSum + pbx); } MW_DEBUG("sprob: %g\n", sprob[0]); for (j = 0; j < ap->number_streams; j++) { nstars[j] += sprob[j]; } MW_DEBUG("nstars: %g\n", nstars[0]); } else { for (j = 0; j < ap->number_streams; j++) { sprob[j] = 1.0; nstars[j] += 1.0; } } /*determine if star with sprob should be put into stream*/ //for(j = 0; j < ap->number_streams; j++) { s_ok = prob_ok(ap->number_streams, sprob); // if (s_ok == 1) { // s_ok += j; // break; // } //} MW_DEBUG("s_ok: %d\n", s_ok); if (s_ok >= 1) { q[s_ok-1]++; } lbr2xyz(star_coords, starxyz); transform_point(starxyz, cmatrix, xsun, starxyzTransform); fprintf(file, "%d %lf %lf %lf\n", s_ok, starxyzTransform[0], starxyzTransform[1], starxyzTransform[2]); //free(starxyz); //free(starxyzTransform); total += 1; if ( (total % 10000) == 0 ) printf("%d\n", total); } printf("%d total stars\n", total); for (j = 0; j < ap->number_streams; j++) { printf("%lf in stream[%d] (%lf%%)\n", nstars[j], j, (nstars[j] / total * 100)); } for (j = 0; j < ap->number_streams; j++) { printf("%d stars separated into stream\n", q[j]); } fclose(file); printf("Output written to: %s\n", filename); free(r_point); free(qw_r3_N); free_constants(ap); }