xlnscpp/xlns16_32test.cpp at main · GaneshPatil7517/xlnscpp · GitHub

1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51

//some numeric tests of xlns16.cpp vs xlns32lpvip
//xlns32_aicasb needs to be defined for accurate 32-bit sums
//xlns16_aicasb allows same hardware to compute 16-bit result
//this program proves this produces identical results
#define xlns16_aicasb
#define xlns32_aicasb

#include <stdio.h>
#include <stdlib.h>
#include <iostream>

#define xlns32_ideal
#include "xlns32.cpp"
#include "xlns32lpvip.cpp"

#define xlns16_alt
#define xlns16_altopt
#include "xlns16.cpp"

void testlns16()
{
	xlns16 z,s,one;
	xlns16 s32;
        xlns32 one32;
	int i;
        one = fp2xlns16(1.0);
        one32 = fp2xlns32(1.0);
	for (z=0x4000-0x400; z<=0x4000; z++)
	{
		s = xlns16_add(one,z);
                s32 = xlns32_add_lpvip(one32,z<<16)>>16;
                if (s != s32)
                {
	            printf("z=%0x %0x %0x\n",z, s, s32);
	            printf("z=%f %f %f\n", xlns162fp(z), xlns162fp(s), xlns162fp(s32));
                }

	}
}


int main()
{
     testlns16();
     return 1;

}