none/tests/ppc32/test_gx.c

   1
   2 #include <stdio.h>
   3 #include "tests/malloc.h"
   4
   5 /* This is a Marie Celeste instruction.  Some IBM documents think it
   6    exists, others don't.  The same appears to be true for
   7    implementations - ppc970 doesn't think it exists, but POWER5
   8    does. */
   9 double do_fre ( double x )
  10 {
  11   double block[2];
  12   block[0] = x;
  13   __asm__ __volatile__(
  14      "lfd %%f1, 0(%0)\n\t"
  15      ".long 0xfc200830\n\t" /* == fre %%f1,%%f1 */
  16      "stfd %%f1, 8(%0)"
  17      : /*out*/
  18      : /*in*/ "b" (&block[0])
  19      : /*trash*/ "memory", "fr1"
  20   );
  21   return block[1];
  22 }
  23
  24 double do_fres ( double x )
  25 {
  26   double block[2];
  27   block[0] = x;
  28   __asm__ __volatile__(
  29      "lfd %%f1, 0(%0)\n\t"
  30      "fres %%f1,%%f1\n\t"
  31      "stfd %%f1, 8(%0)"
  32      : /*out*/
  33      : /*in*/ "b" (&block[0])
  34      : /*trash*/ "memory", "fr1"
  35      );
  36   return block[1];
  37 }
  38
  39 double do_frsqrte ( double x )
  40 {
  41   double block[2];
  42   block[0] = x;
  43   __asm__ __volatile__(
  44      "lfd %%f1, 0(%0)\n\t"
  45      "frsqrte %%f1,%%f1\n\t"
  46      "stfd %%f1, 8(%0)"
  47      : /*out*/
  48      : /*in*/ "b" (&block[0])
  49      : /*trash*/ "memory", "fr1"
  50      );
  51   return block[1];
  52 }
  53
  54 /* Another Marie Celeste insn. */
  55 double do_frsqrtes ( double x )
  56 {
  57   double block[2];
  58   block[0] = x;
  59   __asm__ __volatile__(
  60      "lfd %%f1, 0(%0)\n\t"
  61      ".long 0xec200834\n\t" /* == frsqrtes %%f1,%%f1 */
  62      "stfd %%f1, 8(%0)"
  63      : /*out*/
  64      : /*in*/ "b" (&block[0])
  65      : /*trash*/ "memory", "fr1"
  66      );
  67   return block[1];
  68 }
  69
  70 ////////////////////////////////////////////////////////////
  71
  72 void do_one ( char* name,
  73               double(*f)(double),
  74               double* args, int nargs,
  75               char* argfmt, char* resfmt )
  76 {
  77   int i;
  78   double a, r;
  79   printf("\n");
  80
  81   for (i = 0; i < nargs; i++) {
  82     a = args[i];
  83     r = f(a);
  84     printf("%s ", name);
  85     printf(argfmt, a);
  86     printf(" -> ");
  87     printf(resfmt, r);
  88     printf("\n");
  89   }
  90 }
  91
  92 int main ( void )
  93 {
  94   int nargs = 19;
  95   double* args = malloc(nargs * sizeof(double));
  96   args[0]  =  0.0;
  97   args[1]  =  1.0 / 0.0; // inf
  98   args[2]  = -args[1]; //  -inf
  99   args[3]  = args[2]/args[2]; // nan
 100   args[4]  = -args[3]; // -nan
 101   args[5]  = -5e100;
 102   args[6]  = -5e20;
 103   args[7]  = -501.0;
 104   args[8]  = -6.0;
 105   args[9]  = -1.01;
 106   args[10] = -2e-20;
 107   args[11] = -2e-200;
 108   args[12] =  2e-200;
 109   args[13] =  2e-20;
 110   args[14] =  1.01;
 111   args[15] =  6.0;
 112   args[16] =  501.0;
 113   args[17] =  5e20;
 114   args[18] =  5e100;
 115
 116   do_one( "fre",  do_fre,  args, nargs, "%e", "%4.1e");
 117   do_one( "fres", do_fres, args, nargs, "%e", "%4.1e");
 118
 119   do_one( "frsqrte",  do_frsqrte,  args, nargs, "%e", "%4.1e");
 120   do_one( "frsqrtes", do_frsqrtes, args, nargs, "%e", "%4.1e");
 121
 122   free(args);
 123   return 0;
 124 }