TeaSpeak/tommath
1
0
Fork 0
Code Issues Pull Requests Projects Releases Wiki Activity

explicit condition

Browse Source
This commit is contained in:
Francois Perrad 2018-02-02 17:44:40 +01:00
parent d051d6ba58
commit e3598dc8b2
2 changed files with 19 additions and 19 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

32
demo/demo.c
View File

@ -696,7 +696,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
cmd[strlen(cmd) - 1] = 0; cmd[strlen(cmd) - 1] = 0;
printf("%-6s ]\r", cmd); printf("%-6s ]\r", cmd);
fflush(stdout); fflush(stdout);
if (!strcmp(cmd, "mul2d")) { if (strcmp(cmd, "mul2d") == 0) {
++mul2d_n; ++mul2d_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -713,7 +713,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&b); draw(&b);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "div2d")) { } else if (strcmp(cmd, "div2d") == 0) {
++div2d_n; ++div2d_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -733,7 +733,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&b); draw(&b);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "add")) { } else if (strcmp(cmd, "add") == 0) {
++add_n; ++add_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -777,7 +777,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "sub")) { } else if (strcmp(cmd, "sub") == 0) {
++sub_n; ++sub_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -795,7 +795,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&d); draw(&d);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "mul")) { } else if (strcmp(cmd, "mul") == 0) {
++mul_n; ++mul_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -813,7 +813,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&d); draw(&d);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "div")) { } else if (strcmp(cmd, "div") == 0) {
++div_n; ++div_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -837,7 +837,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "sqr")) { } else if (strcmp(cmd, "sqr") == 0) {
++sqr_n; ++sqr_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -852,7 +852,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&c); draw(&c);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "gcd")) { } else if (strcmp(cmd, "gcd") == 0) {
++gcd_n; ++gcd_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -871,7 +871,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&d); draw(&d);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "lcm")) { } else if (strcmp(cmd, "lcm") == 0) {
++lcm_n; ++lcm_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -890,7 +890,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&d); draw(&d);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "expt")) { } else if (strcmp(cmd, "expt") == 0) {
++expt_n; ++expt_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -911,7 +911,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&e); draw(&e);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "invmod")) { } else if (strcmp(cmd, "invmod") == 0) {
++inv_n; ++inv_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -933,7 +933,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "div2")) { } else if (strcmp(cmd, "div2") == 0) {
++div2_n; ++div2_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -947,7 +947,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&c); draw(&c);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "mul2")) { } else if (strcmp(cmd, "mul2") == 0) {
++mul2_n; ++mul2_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -961,7 +961,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
draw(&c); draw(&c);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "add_d")) { } else if (strcmp(cmd, "add_d") == 0) {
++add_d_n; ++add_d_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -978,7 +978,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
printf("d == %d\n", ix); printf("d == %d\n", ix);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "sub_d")) { } else if (strcmp(cmd, "sub_d") == 0) {
++sub_d_n; ++sub_d_n;
FGETS(buf, 4095, stdin); FGETS(buf, 4095, stdin);
mp_read_radix(&a, buf, 64); mp_read_radix(&a, buf, 64);
@ -995,7 +995,7 @@ printf("compare no compare!\n"); return EXIT_FAILURE;
printf("d == %d\n", ix); printf("d == %d\n", ix);
return EXIT_FAILURE; return EXIT_FAILURE;
} }
} else if (!strcmp(cmd, "exit")) { } else if (strcmp(cmd, "exit") == 0) {
printf("\nokay, exiting now\n"); printf("\nokay, exiting now\n");
break; break;
} }

6
demo/timing.c
View File

@ -36,7 +36,7 @@ static unsigned long lfsr = 0xAAAAAAAAUL;
static int lbit(void) static int lbit(void)
{ {
if (lfsr & 0x80000000UL) { if ((lfsr & 0x80000000UL) != 0UL) {
lfsr = ((lfsr << 1) ^ 0x8000001BUL) & 0xFFFFFFFFUL; lfsr = ((lfsr << 1) ^ 0x8000001BUL) & 0xFFFFFFFFUL;
return 1; return 1;
} else { } else {
@ -258,7 +258,7 @@ int main(void)
logb = FOPEN("logs/expt_dr.log", "w"); logb = FOPEN("logs/expt_dr.log", "w");
logc = FOPEN("logs/expt_2k.log", "w"); logc = FOPEN("logs/expt_2k.log", "w");
logd = FOPEN("logs/expt_2kl.log", "w"); logd = FOPEN("logs/expt_2kl.log", "w");
for (n = 0; primes[n]; n++) { for (n = 0; primes[n] != NULL; n++) {
SLEEP; SLEEP;
mp_read_radix(&a, primes[n], 10); mp_read_radix(&a, primes[n], 10);
mp_zero(&b); mp_zero(&b);
@ -283,7 +283,7 @@ int main(void)
mp_sub(&e, &b, &b); mp_sub(&e, &b, &b);
mp_exptmod(&c, &b, &a, &e); /* c^(p-1-b) mod a */ mp_exptmod(&c, &b, &a, &e); /* c^(p-1-b) mod a */
mp_mulmod(&e, &d, &a, &d); /* c^b * c^(p-1-b) == c^p-1 == 1 */ mp_mulmod(&e, &d, &a, &d); /* c^b * c^(p-1-b) == c^p-1 == 1 */
if (mp_cmp_d(&d, 1)) { if (mp_cmp_d(&d, 1) != MP_EQ) {
printf("Different (%d)!!!\n", mp_count_bits(&a)); printf("Different (%d)!!!\n", mp_count_bits(&a));
draw(&d); draw(&d);
exit(0); exit(0);