import newlib-2000-02-17 snapshot

1 files changed, 291 insertions, 0 deletions

diff --git a/newlib/libm/test/test.c b/newlib/libm/test/test.c
new file mode 100644
index 000000000..615c54aa5
--- /dev/null
+++ b/newlib/libm/test/test.c
@@ -0,0 +1,291 @@
+#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);
+
+}