|
111
|
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
|
|
|
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
|
|
|
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 }
|
