gcc/testsuite/gcc.target/i386/pr59133.c

   1 /* { dg-do compile } */
   2 /* { dg-options "-O3 -ffast-math -march=core-avx2" } */
   3
   4 #define XX 0
   5 #define YY 1
   6 #define ZZ 2
   7 #define DIM 3
   8 typedef float matrix[DIM][DIM];
   9 typedef float rvec[DIM];
  10 extern int det (matrix);
  11 extern void foo(matrix);
  12
  13 void bar1 (int n,int *index,rvec x[],matrix trans)
  14 {
  15   float   xt,yt,zt;
  16   int    i,ii;
  17
  18   for(i=0; (i<n); i++) {
  19     ii=index ? index[i] : i;
  20     xt=x[ii][XX];
  21     yt=x[ii][YY];
  22     zt=x[ii][ZZ];
  23     x[ii][XX]=trans[XX][XX]*xt+trans[XX][YY]*yt+trans[XX][ZZ]*zt;
  24     x[ii][YY]=trans[YY][XX]*xt+trans[YY][YY]*yt+trans[YY][ZZ]*zt;
  25     x[ii][ZZ]=trans[ZZ][XX]*xt+trans[ZZ][YY]*yt+trans[ZZ][ZZ]*zt;
  26   }
  27 }
  28
  29
  30 void bar2 (int n, rvec x[])
  31 {
  32   int     m;
  33   matrix  trans;
  34
  35   foo (trans);
  36
  37   if (det (trans) < 0) {
  38     for(m=0; (m<DIM); m++)
  39       trans[ZZ][m] = -trans[ZZ][m];
  40   }
  41   bar1 (n,(int*) 0,x,trans);
  42 }