1 /* PR target/9732
2    This testcase segfaulted on PPC because PIC_OFFSET_TABLE_REGNUM was no
3    fixed register.
4    Distilled from the xvid sources by Guillaume Morin <guillaume@morinfr.org>
5    and Benjamin Herrenschmidt <benh@kernel.crashing.org>.  */
6 /* { dg-do run { target fpic } } */
7 /* { dg-options "-O2 -fPIC" } */
8 
9 extern void exit (int);
10 
11 #define W1 2841					/* 2048*sqrt(2)*cos(1*pi/16) */
12 #define W2 2676					/* 2048*sqrt(2)*cos(2*pi/16) */
13 #define W3 2408					/* 2048*sqrt(2)*cos(3*pi/16) */
14 #define W5 1609					/* 2048*sqrt(2)*cos(5*pi/16) */
15 #define W6 1108					/* 2048*sqrt(2)*cos(6*pi/16) */
16 #define W7 565					/* 2048*sqrt(2)*cos(7*pi/16) */
17 
18 
19 /* private data */
20 static short iclip[1024];		/* clipping table */
21 static short *iclp;
22 
23 void
idct_int32(short * const block)24 idct_int32(short *const block)
25 {
26   static short *blk;
27   static long i;
28   static long X0, X1, X2, X3, X4, X5, X6, X7, X8;
29 
30   for (i = 0; i < 8; i++)		/* idct columns */
31     {
32       blk = block + i;
33       /* shortcut  */
34       if (! ((X1 = (blk[8 * 4] << 8)) | (X2 = blk[8 * 6])
35 	     | (X3 = blk[8 * 2]) | (X4 = blk[8 * 1]) | (X5 = blk[8 * 7])
36 	     | (X6 = blk[8 * 5]) | (X7 = blk[8 * 3])))
37 	{
38 	  blk[8 * 0] = blk[8 * 1] = blk[8 * 2] =
39 	    blk[8 * 3] = blk[8 * 4] =
40 	    blk[8 * 5] = blk[8 * 6] = blk[8 * 7] =
41 	    iclp[(blk[8 * 0] + 32) >> 6];
42 	  continue;
43 	}
44       X0 = (blk[8 * 0] << 8) + 8192;
45 
46       /* first stage  */
47       X8 = W7 * (X4 + X5) + 4;
48       X4 = (X8 + (W1 - W7) * X4) >> 3;
49       X5 = (X8 - (W1 + W7) * X5) >> 3;
50       X8 = W3 * (X6 + X7) + 4;
51       X6 = (X8 - (W3 - W5) * X6) >> 3;
52       X7 = (X8 - (W3 + W5) * X7) >> 3;
53 
54       /* second stage  */
55       X8 = X0 + X1;
56       X0 -= X1;
57       X1 = W6 * (X3 + X2) + 4;
58       X2 = (X1 - (W2 + W6) * X2) >> 3;
59       X3 = (X1 + (W2 - W6) * X3) >> 3;
60       X1 = X4 + X6;
61       X4 -= X6;
62       X6 = X5 + X7;
63       X5 -= X7;
64 
65       /* third stage  */
66       X7 = X8 + X3;
67       X8 -= X3;
68       X3 = X0 + X2;
69       X0 -= X2;
70       X2 = (181 * (X4 + X5) + 128) >> 8;
71       X4 = (181 * (X4 - X5) + 128) >> 8;
72 
73       /* fourth stage  */
74       blk[8 * 0] = iclp[(X7 + X1) >> 14];
75       blk[8 * 1] = iclp[(X3 + X2) >> 14];
76       blk[8 * 2] = iclp[(X0 + X4) >> 14];
77       blk[8 * 3] = iclp[(X8 + X6) >> 14];
78       blk[8 * 4] = iclp[(X8 - X6) >> 14];
79       blk[8 * 5] = iclp[(X0 - X4) >> 14];
80       blk[8 * 6] = iclp[(X3 - X2) >> 14];
81       blk[8 * 7] = iclp[(X7 - X1) >> 14];
82     }
83 }								/* end function idct_int32(block) */
84 
85 
main(void)86 int main(void) {
87   int i;
88   unsigned short tab[64];
89 
90   for (i = 0 ; i < 64 ; ++i)
91     {
92       tab[i] = (1+(int) (65535)*8000/(2147483647+1.0));
93     }
94 
95   iclp = iclip + 512;
96   for (i = -512; i < 512; i++)
97     iclp[i] = (i < -256) ? -256 : ((i > 255) ? 255 : i);
98 
99   idct_int32((short *) tab);
100   exit (0);
101 }
102