dropbear: tomsfastmath/demo/test.c annotate

annotate tomsfastmath/demo/test.c @ 646:f10335e5e42f dropbear-tfm

- More asm constraint fixes. Now seems to build OK on 32-bit OS X.

author	Matt Johnston <matt@ucc.asn.au>
date	Wed, 30 Nov 2011 23:03:47 +0800
parents	a362b62d38b2
children

rev	line source
643 a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	1 /* TFM demo program */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	2 #include <tfm.h>
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	3
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	4 void draw(fp_int *a)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	5 {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	6 int x;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	7 printf("%d, %d, ", a->used, a->sign);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	8 for (x = a->used - 1; x >= 0; x--) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	9 printf("%08lx ", a->dp[x]);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	10 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	11 printf("\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	12 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	13
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	14 int myrng(unsigned char dst, int len, void dat)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	15 {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	16 int x;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	17 for (x = 0; x < len; x++) dst[x] = rand() & 0xFF;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	18 return len;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	19 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	20
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	21 /* RDTSC from Scott Duplichan */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	22 static ulong64 TIMFUNC (void)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	23 {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	24 #if defined __GNUC__
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	25 #if defined(INTEL_CC)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	26 ulong64 a;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	27 asm ("rdtsc":"=A"(a));
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	28 return a;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	29 #elif defined(__i386__) \|\| defined(__x86_64__)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	30 ulong64 a;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	31 __asm__ __volatile__ ("rdtsc\nmovl %%eax,%0\nmovl %%edx,4+%0\n"::"m"(a):"%eax","%edx");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	32 return a;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	33 #elif defined(TFM_PPC32)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	34 unsigned long a, b;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	35 __asm__ __volatile__ ("mftbu %1 \nmftb %0\n":"=r"(a), "=r"(b));
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	36 return (((ulong64)b) << 32ULL) \| ((ulong64)a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	37 #elif defined(TFM_AVR32)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	38 FILE *in;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	39 char buf[20];
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	40 in = fopen("/sys/devices/system/cpu/cpu0/pccycles", "r");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	41 fgets(buf, 20, in);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	42 fclose(in);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	43 return strtoul(buf, NULL, 10);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	44 #else /* gcc-IA64 version */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	45 unsigned long result;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	46 __asm__ __volatile__("mov %0=ar.itc" : "=r"(result) :: "memory");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	47 while (__builtin_expect ((int) result == -1, 0))
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	48 __asm__ __volatile__("mov %0=ar.itc" : "=r"(result) :: "memory");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	49 return result;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	50 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	51
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	52 // Microsoft and Intel Windows compilers
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	53 #elif defined _M_IX86
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	54 __asm rdtsc
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	55 #elif defined _M_AMD64
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	56 return __rdtsc ();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	57 #elif defined _M_IA64
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	58 #if defined __INTEL_COMPILER
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	59 #include <ia64intrin.h>
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	60 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	61 return __getReg (3116);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	62 #else
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	63 #error need rdtsc function for this build
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	64 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	65 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	66
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	67 char cmd[4096], buf[4096];
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	68
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	69 int main(void)
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	70 {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	71 fp_int a,b,c,d,e,f;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	72 fp_digit fp;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	73 int n, err;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	74 unsigned long expt_n, add_n, sub_n, mul_n, div_n, sqr_n, mul2d_n, div2d_n, gcd_n, lcm_n, inv_n,
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	75 div2_n, mul2_n, add_d_n, sub_d_n, mul_d_n, t, cnt, rr, ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	76 ulong64 t1, t2;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	77
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	78 srand(time(NULL));
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	79 printf("TFM Ident string:\n%s\n\n", fp_ident());
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	80 fp_zero(&b); fp_zero(&c); fp_zero(&d); fp_zero(&e); fp_zero(&f);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	81 fp_zero(&a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	82
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	83 /* test set and simple shifts */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	84 printf("Testing mul/div 2\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	85 fp_set(&a, 1); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	86 for (n = 0; n <= DIGIT_BIT; n++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	87 fp_mul_2(&a, &a); printf("(%d) ", fp_count_bits(&a));
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	88 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	89
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	90 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	91 for (n = 0; n <= (DIGIT_BIT + 1); n++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	92 fp_div_2(&a, &a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	93 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	94 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	95 fp_set(&a, 1);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	96
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	97 /* test lshd/rshd */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	98 printf("testing lshd/rshd\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	99 fp_lshd(&a, 3); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	100 fp_rshd(&a, 3); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	101
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	102 /* test more complicated shifts */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	103 printf("Testing mul/div 2d\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	104 fp_mul_2d(&a, DIGIT_BIT/2, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	105 fp_div_2d(&a, DIGIT_BIT/2, &a, NULL); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	106
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	107 fp_mul_2d(&a, DIGIT_BIT + DIGIT_BIT/2, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	108 fp_div_2d(&a, DIGIT_BIT + DIGIT_BIT/2, &a, NULL); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	109
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	110 /* test neg/abs */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	111 printf("testing neg/abs\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	112 fp_neg(&a, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	113 fp_neg(&a, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	114 fp_neg(&a, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	115 fp_abs(&a, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	116
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	117 /* test comparisons */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	118 fp_set(&b, 3);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	119 fp_set(&c, 4); fp_neg(&c, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	120 fp_set(&d, 1);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	121 printf("Testing compares\n%d, %d, %d, %d\n", fp_cmp(&a, &b), fp_cmp(&a, &c), fp_cmp(&a, &d), fp_cmp(&b, &c));
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	122
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	123 /* test add/sub */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	124 printf("Testing add/sub \n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	125 fp_set(&a, ((fp_digit)1)<<(DIGIT_BIT-1)); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	126 fp_set(&b, ((fp_digit)1)<<(DIGIT_BIT-2));
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	127 fp_add(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	128 fp_add(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	129 fp_add(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	130 printf("sub...\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	131 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	132 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	133 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	134 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	135 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	136 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	137 fp_read_radix(&a, "FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFF000000000000000000000001", 16); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	138 fp_sub_d(&a, 3, &b); draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	139 fp_read_radix(&a, "FFFFFFFFFFFFFFFFFFFFFFFFFFFFFFFEFFFFFFFFFFFFFFFFFFFFFFFE", 16);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	140 printf("cmp returns: %d, ", fp_cmp(&a, &b)); fp_sub(&a, &b, &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	141
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	142 /* test mul_d */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	143 printf("Testing mul_d and div_d\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	144 fp_set(&a, 1);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	145 fp_mul_d(&a, ((fp_digit)1)<<(DIGIT_BIT/2), &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	146 fp_mul_d(&a, ((fp_digit)1)<<(DIGIT_BIT/2), &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	147 fp_mul_d(&a, ((fp_digit)1)<<(DIGIT_BIT/2), &a); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	148 printf("div_d\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	149 fp_div_d(&a, ((fp_digit)1)<<(DIGIT_BIT/2), &a, NULL); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	150 fp_div_d(&a, ((fp_digit)1)<<(DIGIT_BIT/2), &a, NULL); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	151 fp_div_d(&a, ((fp_digit)1)<<(DIGIT_BIT/2), &a, NULL); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	152
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	153 /* testing read radix */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	154 printf("Testing read_radix\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	155 fp_read_radix(&a, "123456789012345678901234567890", 16); draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	156
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	157 #if 0
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	158 /* test mont */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	159 printf("Montgomery test #1\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	160 fp_set(&a, 0x1234567ULL);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	161 fp_montgomery_setup(&a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	162 fp_montgomery_calc_normalization(&b, &a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	163
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	164 fp_read_radix(&d, "123456789123", 16);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	165 for (n = 0; n < 1000000; n++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	166 fp_add_d(&d, 1, &d); fp_sqrmod(&d, &a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	167 fp_mul(&d, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	168 fp_montgomery_reduce(&c, &a, fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	169 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	170 printf("Failed mont %d\n", n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	171 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	172 draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	173 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	174 return EXIT_FAILURE;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	175 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	176 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	177 printf("Passed.\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	178
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	179 printf("Montgomery test #2\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	180 fp_set(&a, 0x1234567ULL);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	181 fp_lshd(&a, 4);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	182 fp_add_d(&a, 1, &a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	183 fp_montgomery_setup(&a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	184 fp_montgomery_calc_normalization(&b, &a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	185
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	186 fp_read_radix(&d, "123456789123", 16);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	187 for (n = 0; n < 1000000; n++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	188 fp_add_d(&d, 1, &d); fp_sqrmod(&d, &a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	189 fp_mul(&d, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	190 fp_montgomery_reduce(&c, &a, fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	191 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	192 printf("Failed mont %d\n", n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	193 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	194 draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	195 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	196 return EXIT_FAILURE;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	197 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	198 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	199 printf("Passed.\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	200
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	201 /* test for size */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	202 for (ix = 8DIGIT_BIT; ix < 10DIGIT_BIT; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	203 printf("Testing (not safe-prime): %9lu bits \r", ix); fflush(stdout);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	204 err = fp_prime_random_ex(&a, 8, ix, (rand()&1)?TFM_PRIME_2MSB_OFF:TFM_PRIME_2MSB_ON, myrng, NULL);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	205 if (err != FP_OKAY) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	206 printf("failed with err code %d\n", err);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	207 return EXIT_FAILURE;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	208 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	209 if ((unsigned long)fp_count_bits(&a) != ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	210 printf("Prime is %d not %lu bits!!!\n", fp_count_bits(&a), ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	211 return EXIT_FAILURE;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	212 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	213 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	214 printf("\n\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	215 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	216
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	217 #ifdef TESTING
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	218 goto testing;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	219 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	220
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	221 #if 1
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	222
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	223 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	224 sleep(1);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	225 printf("Ticks per second: %llu\n", TIMFUNC() - t1);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	226
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	227 goto multtime;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	228 /* do some timings... */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	229 printf("Addition:\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	230 for (t = 2; t <= FP_SIZE/2; t += 2) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	231 fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	232 fp_zero(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	233 fp_zero(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	234 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	235 a.dp[ix] = ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	236 b.dp[ix] = ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	237 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	238 a.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	239 b.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	240 t2 = -1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	241 for (ix = 0; ix < 25000; ++ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	242 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	243 fp_add(&a, &b, &c); fp_add(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	244 fp_add(&a, &b, &c); fp_add(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	245 fp_add(&a, &b, &c); fp_add(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	246 fp_add(&a, &b, &c); fp_add(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	247 t2 = (TIMFUNC() - t1)>>3;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	248 if (t1<t2) { --ix; t2 = t1; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	249 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	250 printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	251 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	252 multtime:
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	253 printf("Multiplication:\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	254 for (t = 2; t < FP_SIZE/2; t += 2) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	255 fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	256 fp_zero(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	257 fp_zero(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	258 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	259 a.dp[ix] = ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	260 b.dp[ix] = ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	261 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	262 a.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	263 b.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	264 t2 = -1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	265 for (ix = 0; ix < 100; ++ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	266 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	267 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	268 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	269 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	270 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	271 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	272 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	273 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	274 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	275 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	276 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	277 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	278 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	279 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	280 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	281 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	282 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	283 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	284 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	285 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	286 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	287 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	288 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	289 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	290 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	291 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	292 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	293 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	294 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	295 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	296 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	297 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	298 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	299 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	300 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	301 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	302 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	303 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	304 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	305 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	306 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	307 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	308 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	309 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	310 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	311 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	312 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	313 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	314 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	315 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	316 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	317 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	318 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	319 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	320 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	321 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	322 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	323 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	324 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	325 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	326 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	327 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	328 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	329 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	330 fp_mul(&a, &b, &c); fp_mul(&a, &b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	331 t2 = (TIMFUNC() - t1)>>7;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	332 if (t1<t2) { --ix; t2 = t1; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	333 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	334 printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	335 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	336 //#else
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	337 sqrtime:
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	338 printf("Squaring:\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	339 for (t = 2; t < FP_SIZE/2; t += 2) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	340 fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	341 fp_zero(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	342 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	343 a.dp[ix] = ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	344 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	345 a.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	346 t2 = -1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	347 for (ix = 0; ix < 100; ++ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	348 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	349 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	350 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	351 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	352 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	353 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	354 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	355 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	356 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	357 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	358 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	359 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	360 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	361 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	362 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	363 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	364 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	365 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	366 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	367 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	368 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	369 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	370 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	371 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	372 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	373 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	374 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	375 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	376 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	377 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	378 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	379 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	380 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	381 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	382 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	383 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	384 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	385 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	386 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	387 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	388 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	389 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	390 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	391 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	392 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	393 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	394 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	395 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	396 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	397 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	398 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	399 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	400 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	401 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	402 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	403 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	404 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	405 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	406 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	407 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	408 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	409 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	410 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	411 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	412 fp_sqr(&a, &b); fp_sqr(&a, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	413 t2 = (TIMFUNC() - t1)>>7;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	414 if (t1<t2) { --ix; t2 = t1; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	415 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	416 printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	417 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	418 invmodtime:
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	419 printf("Invmod:\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	420 for (t = 2; t < FP_SIZE/2; t += 2) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	421 fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	422 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	423 a.dp[ix] = ix \| 1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	424 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	425 a.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	426 fp_zero(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	427 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	428 b.dp[ix] = rand();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	429 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	430 b.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	431 fp_clamp(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	432 fp_zero(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	433 t2 = -1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	434 for (ix = 0; ix < 100; ++ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	435 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	436 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	437 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	438 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	439 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	440 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	441 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	442 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	443 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	444 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	445 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	446 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	447 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	448 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	449 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	450 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	451 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	452 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	453 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	454 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	455 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	456 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	457 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	458 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	459 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	460 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	461 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	462 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	463 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	464 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	465 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	466 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	467 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	468 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	469 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	470 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	471 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	472 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	473 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	474 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	475 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	476 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	477 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	478 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	479 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	480 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	481 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	482 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	483 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	484 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	485 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	486 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	487 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	488 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	489 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	490 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	491 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	492 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	493 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	494 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	495 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	496 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	497 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	498 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	499 fp_invmod(&b, &a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	500 t2 = (TIMFUNC() - t1)>>6;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	501 if (t1<t2) { --ix; t2 = t1; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	502 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	503 printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	504 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	505 //#else
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	506 monttime:
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	507 printf("Montgomery:\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	508 for (t = 2; t <= (FP_SIZE/2)-4; t += 2) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	509 // printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	510 fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	511 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	512 a.dp[ix] = ix \| 1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	513 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	514 a.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	515
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	516 fp_montgomery_setup(&a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	517 fp_sub_d(&a, 3, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	518 fp_sqr(&b, &b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	519 fp_copy(&b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	520 fp_copy(&b, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	521
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	522 t2 = -1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	523 for (ix = 0; ix < 100; ++ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	524 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	525 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	526 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	527 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	528 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	529 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	530 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	531 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	532 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	533 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	534 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	535 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	536 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	537 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	538 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	539 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	540 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	541 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	542 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	543 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	544 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	545 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	546 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	547 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	548 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	549 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	550 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	551 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	552 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	553 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	554 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	555 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	556 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	557 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	558 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	559 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	560 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	561 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	562 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	563 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	564 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	565 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	566 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	567 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	568 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	569 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	570 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	571 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	572 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	573 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	574 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	575 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	576 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	577 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	578 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	579 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	580 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	581 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	582 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	583 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	584 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	585 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	586 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	587 fp_montgomery_reduce(&c, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	588 fp_montgomery_reduce(&d, &a, &fp);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	589 t2 = (TIMFUNC() - t1)>>6;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	590 fp_copy(&b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	591 fp_copy(&b, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	592 if (t1<t2) { --ix; t2 = t1; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	593 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	594 printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	595 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	596 //#else
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	597 expttime:
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	598 printf("Exptmod:\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	599
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	600 for (t = 512/DIGIT_BIT; t <= (FP_SIZE/2)-2; t += 256/DIGIT_BIT) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	601 fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	602 fp_zero(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	603 fp_zero(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	604 for (ix = 0; ix < t; ix++) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	605 a.dp[ix] = ix+1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	606 b.dp[ix] = (fp_digit)rand() * (fp_digit)rand();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	607 c.dp[ix] = ix;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	608 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	609 a.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	610 b.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	611 c.used = t;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	612
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	613 t2 = -1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	614 for (ix = 0; ix < 500; ++ix) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	615 t1 = TIMFUNC();
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	616 fp_exptmod(&c, &b, &a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	617 fp_exptmod(&c, &b, &a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	618 t2 = (TIMFUNC() - t1)>>1;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	619 fp_copy(&b, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	620 fp_copy(&b, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	621 if (t1<t2) { t2 = t1; --ix; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	622 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	623 printf("%5lu-bit: %9llu\n", t * DIGIT_BIT, t2);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	624 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	625 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	626 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	627
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	628 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	629 testing:
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	630
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	631 fp_zero(&b); fp_zero(&c); fp_zero(&d); fp_zero(&e); fp_zero(&f); fp_zero(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	632
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	633
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	634 div2_n = mul2_n = inv_n = expt_n = lcm_n = gcd_n = add_n =
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	635 sub_n = mul_n = div_n = sqr_n = mul2d_n = div2d_n = cnt = add_d_n = sub_d_n= mul_d_n = 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	636
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	637 for (;;) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	638 printf("%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu/%4lu ", add_n, sub_n, mul_n, div_n, sqr_n, mul2d_n, div2d_n, gcd_n, lcm_n, expt_n, inv_n, div2_n, mul2_n, add_d_n, sub_d_n, mul_d_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	639 fgets(cmd, 4095, stdin);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	640 cmd[strlen(cmd)-1] = 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	641 printf("%s ]\r",cmd); fflush(stdout);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	642 if (!strcmp(cmd, "mul2d")) { ++mul2d_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	643 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	644 fgets(buf, 4095, stdin); sscanf(buf, "%lu", &rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	645 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	646
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	647 fp_mul_2d(&a, rr, &a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	648 a.sign = b.sign;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	649 if (fp_cmp(&a, &b) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	650 printf("mul2d failed, rr == %lu\n",rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	651 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	652 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	653 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	654 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	655 } else if (!strcmp(cmd, "div2d")) { ++div2d_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	656 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	657 fgets(buf, 4095, stdin); sscanf(buf, "%lu", &rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	658 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	659
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	660 fp_div_2d(&a, rr, &a, &e);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	661 a.sign = b.sign;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	662 if (a.used == b.used && a.used == 0) { a.sign = b.sign = FP_ZPOS; }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	663 if (fp_cmp(&a, &b) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	664 printf("div2d failed, rr == %lu\n",rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	665 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	666 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	667 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	668 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	669 } else if (!strcmp(cmd, "add")) { ++add_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	670 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	671 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	672 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	673 fp_copy(&a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	674 fp_add(&d, &b, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	675 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	676 printf("add %lu failure!\n", add_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	677 draw(&a);draw(&b);draw(&c);draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	678 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	679 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	680
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	681 /* test the sign/unsigned storage functions */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	682
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	683 rr = fp_signed_bin_size(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	684 fp_to_signed_bin(&c, (unsigned char *)cmd);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	685 memset(cmd+rr, rand()&255, sizeof(cmd)-rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	686 fp_read_signed_bin(&d, (unsigned char *)cmd, rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	687 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	688 printf("fp_signed_bin failure!\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	689 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	690 draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	691 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	692 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	693
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	694 rr = fp_unsigned_bin_size(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	695 fp_to_unsigned_bin(&c, (unsigned char *)cmd);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	696 memset(cmd+rr, rand()&255, sizeof(cmd)-rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	697 fp_read_unsigned_bin(&d, (unsigned char *)cmd, rr);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	698 if (fp_cmp_mag(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	699 printf("fp_unsigned_bin failure!\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	700 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	701 draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	702 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	703 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	704
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	705 } else if (!strcmp(cmd, "sub")) { ++sub_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	706 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	707 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	708 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	709 fp_copy(&a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	710 fp_sub(&d, &b, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	711 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	712 printf("sub %lu failure!\n", sub_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	713 draw(&a);draw(&b);draw(&c);draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	714 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	715 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	716 } else if (!strcmp(cmd, "mul")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	717 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	718 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	719 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	720 //continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	721 fp_copy(&a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	722 fp_mul(&d, &b, &d); ++mul_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	723 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	724 printf("mul %lu failure!\n", mul_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	725 draw(&a);draw(&b);draw(&c);draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	726 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	727 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	728 } else if (!strcmp(cmd, "div")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	729 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	730 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	731 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	732 fgets(buf, 4095, stdin); fp_read_radix(&d, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	733 // continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	734 fp_div(&a, &b, &e, &f); ++div_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	735 if (fp_cmp(&c, &e) != FP_EQ \|\| fp_cmp(&d, &f) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	736 printf("div %lu failure!\n", div_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	737 draw(&a);draw(&b);draw(&c);draw(&d); draw(&e); draw(&f);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	738 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	739 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	740
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	741 } else if (!strcmp(cmd, "sqr")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	742 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	743 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	744 // continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	745 fp_copy(&a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	746 fp_sqr(&c, &c); ++sqr_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	747 if (fp_cmp(&b, &c) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	748 printf("sqr %lu failure!\n", sqr_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	749 draw(&a);draw(&b);draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	750 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	751 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	752 } else if (!strcmp(cmd, "gcd")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	753 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	754 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	755 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	756 // continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	757 fp_copy(&a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	758 fp_gcd(&d, &b, &d); ++gcd_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	759 d.sign = c.sign;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	760 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	761 printf("gcd %lu failure!\n", gcd_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	762 draw(&a);draw(&b);draw(&c);draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	763 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	764 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	765 } else if (!strcmp(cmd, "lcm")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	766 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	767 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	768 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	769 //continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	770 fp_copy(&a, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	771 fp_lcm(&d, &b, &d); ++lcm_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	772 d.sign = c.sign;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	773 if (fp_cmp(&c, &d) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	774 printf("lcm %lu failure!\n", lcm_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	775 draw(&a);draw(&b);draw(&c);draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	776 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	777 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	778 } else if (!strcmp(cmd, "expt")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	779 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	780 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	781 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	782 fgets(buf, 4095, stdin); fp_read_radix(&d, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	783 // continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	784 fp_copy(&a, &e);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	785 fp_exptmod(&e, &b, &c, &e); ++expt_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	786 if (fp_cmp(&d, &e) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	787 printf("expt %lu failure!\n", expt_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	788 draw(&a);draw(&b);draw(&c);draw(&d); draw(&e);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	789 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	790 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	791 } else if (!strcmp(cmd, "invmod")) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	792 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	793 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	794 fgets(buf, 4095, stdin); fp_read_radix(&c, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	795 //continue;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	796 fp_invmod(&a, &b, &d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	797 #if 1
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	798 fp_mulmod(&d,&a,&b,&e); ++inv_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	799 if (fp_cmp_d(&e, 1) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	800 #else
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	801 if (fp_cmp(&d, &c) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	802 #endif
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	803 printf("inv [wrong value from MPI?!] failure\n");
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	804 draw(&a);draw(&b);draw(&c);draw(&d);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	805 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	806 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	807
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	808 } else if (!strcmp(cmd, "div2")) { ++div2_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	809 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	810 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	811 fp_div_2(&a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	812 if (fp_cmp(&c, &b) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	813 printf("div_2 %lu failure\n", div2_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	814 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	815 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	816 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	817 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	818 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	819 } else if (!strcmp(cmd, "mul2")) { ++mul2_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	820 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	821 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	822 fp_mul_2(&a, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	823 if (fp_cmp(&c, &b) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	824 printf("mul_2 %lu failure\n", mul2_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	825 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	826 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	827 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	828 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	829 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	830 } else if (!strcmp(cmd, "add_d")) { ++add_d_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	831 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	832 fgets(buf, 4095, stdin); sscanf(buf, "%lu", &ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	833 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	834 fp_add_d(&a, ix, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	835 if (fp_cmp(&b, &c) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	836 printf("add_d %lu failure\n", add_d_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	837 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	838 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	839 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	840 printf("d == %lu\n", ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	841 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	842 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	843 } else if (!strcmp(cmd, "sub_d")) { ++sub_d_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	844 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	845 fgets(buf, 4095, stdin); sscanf(buf, "%lu", &ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	846 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	847 fp_sub_d(&a, ix, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	848 if (fp_cmp(&b, &c) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	849 printf("sub_d %lu failure\n", sub_d_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	850 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	851 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	852 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	853 printf("d == %lu\n", ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	854 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	855 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	856 } else if (!strcmp(cmd, "mul_d")) { ++mul_d_n;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	857 fgets(buf, 4095, stdin); fp_read_radix(&a, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	858 fgets(buf, 4095, stdin); sscanf(buf, "%lu", &ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	859 fgets(buf, 4095, stdin); fp_read_radix(&b, buf, 64);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	860 fp_mul_d(&a, ix, &c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	861 if (fp_cmp(&b, &c) != FP_EQ) {
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	862 printf("mul_d %lu failure\n", sub_d_n);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	863 draw(&a);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	864 draw(&b);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	865 draw(&c);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	866 printf("d == %lu\n", ix);
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	867 return 0;
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	868 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	869 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	870
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	871 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	872 }
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	873
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	874
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	875
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	876 /* $Source$ */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	877 /* $Revision$ */
a362b62d38b2 Add tomsfastmath from git rev bfa4582842bc3bab42e4be4aed5703437049502a Matt Johnston <matt@ucc.asn.au> parents: diff changeset	878 /* $Date$ */

Mercurial > dropbear

annotate tomsfastmath/demo/test.c @ 646:f10335e5e42f dropbear-tfm