none/tests/ppc32/test_fx.c

   1
   2 #include <stdio.h>
   3 #include <stdlib.h>
   4
   5 double do_fsqrt ( double x )
   6 {
   7   double block[2];
   8   block[0] = x;
   9   __asm__ __volatile__(
  10      "lfd %%f1, 0(%0)\n\t"
  11      "fsqrt %%f1,%%f1\n\t"
  12      "stfd %%f1, 8(%0)"
  13      : /*out*/
  14      : /*in*/ "b" (&block[0])
  15      : /*trash*/ "memory", "fr1"
  16   );
  17   return block[1];
  18 }
  19
  20 double do_fsqrts ( double x )
  21 {
  22   double block[2];
  23   block[0] = x;
  24   __asm__ __volatile__(
  25      "lfd %%f1, 0(%0)\n\t"
  26      "fsqrts %%f1,%%f1\n\t"
  27      "stfd %%f1, 8(%0)"
  28      : /*out*/
  29      : /*in*/ "b" (&block[0])
  30      : /*trash*/ "memory", "fr1"
  31      );
  32   return block[1];
  33 }
  34
  35 ////////////////////////////////////////////////////////////
  36
  37 void do_one ( char* name,
  38               double(*f)(double),
  39               double* args, int nargs,
  40               char* argfmt, char* resfmt )
  41 {
  42   int i;
  43   double a, r;
  44   printf("\n");
  45
  46   for (i = 0; i < nargs; i++) {
  47     a = args[i];
  48     r = f(a);
  49     printf("%s ", name);
  50     printf(argfmt, a);
  51     printf(" -> ");
  52     printf(resfmt, r);
  53     printf("\n");
  54   }
  55 }
  56
  57 int main ( void )
  58 {
  59   int nargs = 19;
  60   double* args = malloc(nargs * sizeof(double));
  61   args[0]  =  0.0;
  62   args[1]  =  1.0 / 0.0; // inf
  63   args[2]  = -args[1]; //  -inf
  64   args[3]  = args[2]/args[2]; // nan
  65   args[4]  = -args[3]; // -nan
  66   args[5]  = -5e100;
  67   args[6]  = -5e20;
  68   args[7]  = -501.0;
  69   args[8]  = -6.0;
  70   args[9]  = -1.0;
  71   args[10] = -2e-20;
  72   args[11] = -2e-200;
  73   args[12] =  2e-200;
  74   args[13] =  2e-20;
  75   args[14] =  1.0;
  76   args[15] =  6.0;
  77   args[16] =  501.0;
  78   args[17] =  5e20;
  79   args[18] =  5e100;
  80
  81   do_one( "fsqrt",  do_fsqrt,  args, nargs, "%e", "%20.14e");
  82   do_one( "fsqrts", do_fsqrts, args, nargs, "%e", "%e");
  83
  84   free(args);
  85   return 0;
  86 }