test/reservoir/mg-psinv.good.c

   1 /* Generated from ../../../git/cloog/test/./reservoir/mg-psinv.cloog by CLooG 0.14.0-72-gefe2fc2 gmp bits in 0.05s. */
   2 extern void hash(int);
   3
   4 /* Useful macros. */
   5 #define floord(n,d) (((n)<0) ? -((-(n)+(d)-1)/(d)) : (n)/(d))
   6 #define ceild(n,d)  (((n)<0) ? -((-(n))/(d)) : ((n)+(d)-1)/(d))
   7 #define max(x,y)    ((x) > (y) ? (x) : (y))
   8 #define min(x,y)    ((x) < (y) ? (x) : (y))
   9
  10 #define S1(i,j,k) { hash(1); hash(i); hash(j); hash(k); }
  11 #define S2(i,j,k) { hash(2); hash(i); hash(j); hash(k); }
  12 #define S3(i,j,k) { hash(3); hash(i); hash(j); hash(k); }
  13
  14 void test(int M, int N, int O)
  15 {
  16   /* Scattering iterators. */
  17   int c2, c4, c6;
  18   /* Original iterators. */
  19   int i, j, k;
  20   if ((M >= 3) && (N >= 4)) {
  21     for (c2=2;c2<=O-1;c2++) {
  22       for (c6=1;c6<=M;c6++) {
  23         S1(c2,2,c6) ;
  24         S2(c2,2,c6) ;
  25       }
  26       for (c4=4;c4<=2*N-3;c4++) {
  27         for (c6=1;c6<=M;c6++) {
  28           if ((c4+1)%2 == 0) {
  29             j = (c4+1)/2 ;
  30             S1(c2,(c4+1)/2,c6) ;
  31             S2(c2,(c4+1)/2,c6) ;
  32           }
  33         }
  34         for (c6=2;c6<=M-1;c6++) {
  35           if (c4%2 == 0) {
  36             S3(c2,c4/2,c6) ;
  37           }
  38         }
  39       }
  40       c4 = 2*N-2 ;
  41       for (c6=2;c6<=M-1;c6++) {
  42         j = N-1 ;
  43         S3(c2,N-1,c6) ;
  44       }
  45     }
  46   }
  47   if ((M >= 3) && (N == 3)) {
  48     for (c2=2;c2<=O-1;c2++) {
  49       for (c6=1;c6<=M;c6++) {
  50         S1(c2,2,c6) ;
  51         S2(c2,2,c6) ;
  52       }
  53       for (c6=2;c6<=M-1;c6++) {
  54         S3(c2,2,c6) ;
  55       }
  56     }
  57   }
  58   if ((M >= 1) && (M <= 2) && (N >= 3)) {
  59     for (c2=2;c2<=O-1;c2++) {
  60       for (c4=3;c4<=2*N-3;c4++) {
  61         for (c6=1;c6<=M;c6++) {
  62           if ((c4+1)%2 == 0) {
  63             j = (c4+1)/2 ;
  64             S1(c2,(c4+1)/2,c6) ;
  65             S2(c2,(c4+1)/2,c6) ;
  66           }
  67         }
  68       }
  69     }
  70   }
  71 }