URL
https://opencores.org/ocsvn/or1k/or1k/trunk
Subversion Repositories or1k
[/] [or1k/] [tags/] [alpha/] [newlib/] [newlib/] [libm/] [test/] [test.c] - Rev 39
Go to most recent revision | Compare with Previous | Blame | View Log
#include <signal.h> #include "test.h" #include <math.h> #include <ieeefp.h> #include <string.h> int verbose; static int count; int inacc; int _DEFUN(main,(ac, av), int ac _AND char **av) { int i; int math2 = 1; int string= 1; int is = 1; int math= 1; int cvt = 1; int ieee= 1; bt(); for (i = 1; i < ac; i++) { if (strcmp(av[i],"-v")==0) verbose ++; if (strcmp(av[i],"-nomath2") == 0) math2 = 0; if (strcmp(av[i],"-nostrin") == 0) string= 0; if (strcmp(av[i],"-nois") == 0) is = 0; if (strcmp(av[i],"-nomath") == 0) math= 0; if (strcmp(av[i],"-nocvt") == 0) cvt = 0; if (strcmp(av[i],"-noiee") == 0) ieee= 0; } if (cvt) test_cvt(); if (math2) test_math2(); if (string) test_string(); if (math) test_math(); if (is) test_is(); if (ieee) test_ieee(); printf("Tested %d functions, %d errors detected\n", count, inacc); return 0; } static _CONST char *iname = "foo"; void _DEFUN(newfunc,(string), _CONST char *string) { if (strcmp(iname, string)) { printf("testing %s\n", string); fflush(stdout); iname = string; } } static int theline; void line(li) int li; { if (verbose) { printf(" %d\n", li); } theline = li; count++; } int redo = 0; int reduce = 0; int strtod_vector = 0; int _DEFUN(bigger,(a,b), __ieee_double_shape_type *a _AND __ieee_double_shape_type *b) { if (a->parts.msw > b->parts.msw) { return 1; } else if (a->parts.msw == b->parts.msw) { if (a->parts.lsw > b->parts.lsw) { return 1; } } return 0; } /* Return the first bit different between two double numbers */ int _DEFUN(mag_of_error,(is, shouldbe), double is _AND double shouldbe) { __ieee_double_shape_type a,b; int i; int a_big; unsigned int mask; unsigned long int __x; unsigned long int msw, lsw; a.value = is; b.value = shouldbe; if (a.parts.msw == b.parts.msw && a.parts.lsw== b.parts.lsw) return 64; /* Subtract the larger from the smaller number */ a_big = bigger(&a, &b); if (!a_big) { int t; t = a.parts.msw; a.parts.msw = b.parts.msw; b.parts.msw = t; t = a.parts.lsw; a.parts.lsw = b.parts.lsw; b.parts.lsw = t; } __x = (a.parts.lsw) - (b.parts.lsw); msw = (a.parts.msw) - (b.parts.msw) - (__x > (a.parts.lsw)); lsw = __x; /* Find out which bit the difference is in */ mask = 0x80000000; for (i = 0; i < 32; i++) { if (((msw) & mask)!=0) return i; mask >>=1; } mask = 0x80000000; for (i = 0; i < 32; i++) { if (((lsw) & mask)!=0) return i+32; mask >>=1; } return 64; } int ok_mag; void _DEFUN(test_sok,(is, shouldbe), char *is _AND char *shouldbe) { if (strcmp(is,shouldbe)) { printf("%s:%d, inacurate answer: (%s should be %s)\n", iname, theline, is, shouldbe); inacc++; } } void _DEFUN(test_iok,(is, shouldbe), int is _AND int shouldbe) { if (is != shouldbe){ printf("%s:%d, inacurate answer: (%08x should be %08x)\n", iname, theline, is, shouldbe); inacc++; } } /* Compare counted strings upto a certain length - useful to test single prec float conversions against double results */ void _DEFUN(test_scok,(is, shouldbe, count), char *is _AND char *shouldbe _AND int count) { if (strncmp(is,shouldbe, count)) { printf("%s:%d, inacurate answer: (%s should be %s)\n", iname, theline, is, shouldbe); inacc++; } } void _DEFUN(test_eok,(is, shouldbe), int is _AND int shouldbe) { if (is != shouldbe){ printf("%s:%d, bad errno answer: (%d should be %d)\n", iname, theline, is, shouldbe); inacc++; } } void _DEFUN(test_mok,(value, shouldbe, okmag), double value _AND double shouldbe _AND int okmag) { __ieee_double_shape_type a,b; int mag = mag_of_error(value, shouldbe); if (mag == 0) { /* error in the first bit is ok if the numbers are both 0 */ if (value == 0.0 && shouldbe == 0.0) return; } a.value = shouldbe; b.value = value; if (mag < okmag) { printf("%s:%d, wrong answer: bit %d ", iname, theline, mag); printf("%08x%08x %08x%08x) ", a.parts.msw, a.parts.lsw, b.parts.msw, b.parts.lsw); printf("(%g %g)\n", a.value, b.value); inacc++; } } #ifdef __PCCNECV70__ kill() {} getpid() {} #endif bt(){ double f1,f2; f1 = 0.0; f2 = 0.0/f1; printf("(%g)\n", f2); }
Go to most recent revision | Compare with Previous | Blame | View Log