1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
|
$NetBSD: patch-ad,v 1.3 1999/10/12 04:43:12 simonb Exp $
--- dct64.c.orig Sat Sep 20 21:16:04 1997
+++ dct64.c Thu Aug 12 23:28:33 1999
@@ -29,7 +29,7 @@
for(i=15;i>=0;i--)
*bs++ = (*b1++ + *--b2);
for(i=15;i>=0;i--)
- *bs++ = (*--b2 - *b1++) * *--costab;
+ *bs++ = REAL_MUL((*--b2 - *b1++), *--costab);
b1 = bufs;
costab = pnts[1]+8;
@@ -39,13 +39,13 @@
for(i=7;i>=0;i--)
*bs++ = (*b1++ + *--b2);
for(i=7;i>=0;i--)
- *bs++ = (*--b2 - *b1++) * *--costab;
+ *bs++ = REAL_MUL((*--b2 - *b1++), *--costab);
b2 += 32;
costab += 8;
for(i=7;i>=0;i--)
*bs++ = (*b1++ + *--b2);
for(i=7;i>=0;i--)
- *bs++ = (*b1++ - *--b2) * *--costab;
+ *bs++ = REAL_MUL((*b1++ - *--b2), *--costab);
b2 += 32;
}
@@ -58,12 +58,12 @@
for(i=3;i>=0;i--)
*bs++ = (*b1++ + *--b2);
for(i=3;i>=0;i--)
- *bs++ = (*--b2 - *b1++) * costab[i];
+ *bs++ = REAL_MUL((*--b2 - *b1++), costab[i]);
b2 += 16;
for(i=3;i>=0;i--)
*bs++ = (*b1++ + *--b2);
for(i=3;i>=0;i--)
- *bs++ = (*b1++ - *--b2) * costab[i];
+ *bs++ = REAL_MUL((*b1++ - *--b2), costab[i]);
b2 += 16;
}
@@ -75,13 +75,13 @@
{
*bs++ = (*b1++ + *--b2);
*bs++ = (*b1++ + *--b2);
- *bs++ = (*--b2 - *b1++) * costab[1];
- *bs++ = (*--b2 - *b1++) * costab[0];
+ *bs++ = REAL_MUL((*--b2 - *b1++), costab[1]);
+ *bs++ = REAL_MUL((*--b2 - *b1++), costab[0]);
b2 += 8;
*bs++ = (*b1++ + *--b2);
*bs++ = (*b1++ + *--b2);
- *bs++ = (*b1++ - *--b2) * costab[1];
- *bs++ = (*b1++ - *--b2) * costab[0];
+ *bs++ = REAL_MUL((*b1++ - *--b2), costab[1]);
+ *bs++ = REAL_MUL((*b1++ - *--b2), costab[0]);
b2 += 8;
}
bs = bufs;
@@ -92,10 +92,10 @@
real v0,v1;
v0=*b1++; v1 = *b1++;
*bs++ = (v0 + v1);
- *bs++ = (v0 - v1) * (*costab);
+ *bs++ = REAL_MUL((v0 - v1), (*costab));
v0=*b1++; v1 = *b1++;
*bs++ = (v0 + v1);
- *bs++ = (v1 - v0) * (*costab);
+ *bs++ = REAL_MUL((v1 - v0), (*costab));
}
}
|
