testsuite/ramp.c

   1
   2 #include <config.h>
   3
   4 #include <stdio.h>
   5
   6
   7 void ramp(int n);
   8 void test(int n);
   9
  10 int sum;
  11
  12 int
  13 main (int argc, char *argv[])
  14 {
  15   int i;
  16
  17   for(i=2;i<=32;i+=2){
  18     ramp(i);
  19   }
  20
  21   sum = 0;
  22   for(i=0;i<256;i++){
  23     test(i);
  24   }
  25   printf("sum %d\n", sum);
  26
  27   return 0;
  28 }
  29
  30
  31 void
  32 ramp(int n)
  33 {
  34   int i;
  35   int xoff = n/2;
  36
  37   printf("%d: ", n);
  38   for(i=0;i<n;i++){
  39     printf("%d ", 1 + (6*(i+1) + xoff)/(2*xoff + 1));
  40   }
  41   printf ("\n");
  42
  43   printf("%d: ", n);
  44   for(i=0;i<n;i++){
  45     printf("%g ", 1.0 + (6.0*(i+1) + xoff)/(2.0*xoff + 1));
  46   }
  47   printf ("\n");
  48
  49   printf ("\n");
  50 }
  51
  52
  53 void
  54 test (int n)
  55 {
  56   int i;
  57   int v, ref;
  58
  59   printf("%d: ", n);
  60   for(i=0;i<8;i++){
  61     //v = ((n*i + 4)>>3) + ((n*(8-i) + 4)>>3);
  62     v = ((n*i)>>3) + ((n*(8-i))>>3);
  63     ref = ((n*i + n*(8-i) + 4)>>3);
  64     printf("%d ", v);
  65     sum += (ref - v) * (ref - v);
  66   }
  67   printf ("\n");
  68 }
  69