gcc/testsuite/gcc.target/mips/lxc1-sxc1-1.c

   1 /* { dg-options "(HAS_LXC1) -mhard-float -mno-lxc1-sxc1" } */
   2 /* { dg-final { scan-assembler-not "\tldxc1\t" } } */
   3 /* { dg-final { scan-assembler-not "\tsdxc1\t" } } */
   4
   5 #ifndef __mips_no_lxc1_sxc1
   6 #error missing definition of __mips_no_lxc1_sxc1
   7 #endif
   8
   9 double ldexp(double x, int exp);
  10
  11 typedef struct
  12 {
  13   double** rows;
  14 } d_mat_struct;
  15
  16 typedef d_mat_struct d_mat_t[1];
  17
  18 #define d_mat_entry(mat,i,j) (*((mat)->rows[i] + (j)))
  19
  20 double __attribute__((noinline))
  21 ldxc1_test (int kappa, int zeros, double ctt, int* expo, d_mat_t r, double* s)
  22 {
  23   int kappa2 = kappa;
  24   double tmp = 0.0;
  25
  26   do
  27     {
  28       kappa--;
  29       if (kappa > zeros + 1)
  30         {
  31           tmp = d_mat_entry(r, kappa - 1, kappa - 1) * ctt;
  32           tmp = ldexp(tmp, (expo[kappa - 1] - expo[kappa2]));
  33         }
  34     }
  35   while ((kappa >= zeros + 2) && (s[kappa - 1] <= tmp));
  36
  37   return tmp;
  38 }
  39
  40 #define SIZE 20
  41
  42 int main(void)
  43 {
  44   int kappa = SIZE - 1;
  45   int zeros = 1;
  46   double ctt = 2;
  47
  48   int expo[SIZE] = {0};
  49   double s[SIZE] = {0};
  50   double rows_data[SIZE][SIZE] = {0};
  51   double* rows[SIZE];
  52
  53   for (int i = 0; i < SIZE; i++)
  54     rows[i] = rows_data[i];
  55
  56   d_mat_t r = { rows };
  57
  58   ldxc1_test(kappa, zeros, ctt, expo, r, s);
  59   return 0;
  60 }