1 /*****************************************************************************
2 * idct_sparse.h : Sparse IDCT functions (must be include at the end)
3 *****************************************************************************
4 * Copyright (C) 1999, 2000 VideoLAN
5 * $Id: idct_sparse.h,v 1.1 2002/08/04 17:23:42 sam Exp $
7 * Author: Gaƫl Hendryckx <jimmy@via.ecp.fr>
9 * This program is free software; you can redistribute it and/or modify
10 * it under the terms of the GNU General Public License as published by
11 * the Free Software Foundation; either version 2 of the License, or
12 * (at your option) any later version.
14 * This program is distributed in the hope that it will be useful,
15 * but WITHOUT ANY WARRANTY; without even the implied warranty of
16 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
17 * GNU General Public License for more details.
19 * You should have received a copy of the GNU General Public License
20 * along with this program; if not, write to the Free Software
21 * Foundation, Inc., 59 Temple Place - Suite 330, Boston, MA 02111, USA.
22 *****************************************************************************/
24 /*****************************************************************************
25 * InitIDCT : initialize data for SparseIDCT
26 *****************************************************************************/
27 static void InitIDCT ( void ** pp_idct_data )
32 *pp_idct_data = malloc( sizeof(dctelem_t) * 64 * 64 );
33 p_pre = (dctelem_t *) *pp_idct_data;
34 memset( p_pre, 0, 64 * 64 * sizeof(dctelem_t) );
36 for( i = 0 ; i < 64 ; i++ )
38 p_pre[i*64+i] = 1 << SPARSE_SCALE_FACTOR;
39 IDCT( &p_pre[i*64] ) ;
47 /*****************************************************************************
48 * SparseIDCT : IDCT function for sparse matrices
49 *****************************************************************************/
50 static inline void SparseIDCT( dctelem_t * p_block, void * p_idct_data,
59 dctelem_t * p_pre = (dctelem_t *) p_idct_data;
61 /* If DC Coefficient. */
62 if ( i_sparse_pos == 0 )
65 val=RIGHT_SHIFT((*p_block + 4), 3);
66 /* Compute int to assign. This speeds things up a bit */
67 v = ((val & 0xffff) | (val << 16));
68 dp[0] = v; dp[1] = v; dp[2] = v; dp[3] = v;
69 dp[4] = v; dp[5] = v; dp[6] = v; dp[7] = v;
70 dp[8] = v; dp[9] = v; dp[10] = v; dp[11] = v;
71 dp[12] = v; dp[13] = v; dp[14] = v; dp[15] = v;
72 dp[16] = v; dp[17] = v; dp[18] = v; dp[19] = v;
73 dp[20] = v; dp[21] = v; dp[22] = v; dp[23] = v;
74 dp[24] = v; dp[25] = v; dp[26] = v; dp[27] = v;
75 dp[28] = v; dp[29] = v; dp[30] = v; dp[31] = v;
78 /* Some other coefficient. */
79 p_dest = (s16*)p_block;
80 p_source = (s16*)&p_pre[i_sparse_pos * 64];
81 coeff = (int)p_dest[i_sparse_pos];
82 for( rr=0 ; rr < 4 ; rr++ )
84 p_dest[0] = (p_source[0] * coeff) >> SPARSE_SCALE_FACTOR;
85 p_dest[1] = (p_source[1] * coeff) >> SPARSE_SCALE_FACTOR;
86 p_dest[2] = (p_source[2] * coeff) >> SPARSE_SCALE_FACTOR;
87 p_dest[3] = (p_source[3] * coeff) >> SPARSE_SCALE_FACTOR;
88 p_dest[4] = (p_source[4] * coeff) >> SPARSE_SCALE_FACTOR;
89 p_dest[5] = (p_source[5] * coeff) >> SPARSE_SCALE_FACTOR;
90 p_dest[6] = (p_source[6] * coeff) >> SPARSE_SCALE_FACTOR;
91 p_dest[7] = (p_source[7] * coeff) >> SPARSE_SCALE_FACTOR;
92 p_dest[8] = (p_source[8] * coeff) >> SPARSE_SCALE_FACTOR;
93 p_dest[9] = (p_source[9] * coeff) >> SPARSE_SCALE_FACTOR;
94 p_dest[10] = (p_source[10] * coeff) >> SPARSE_SCALE_FACTOR;
95 p_dest[11] = (p_source[11] * coeff) >> SPARSE_SCALE_FACTOR;
96 p_dest[12] = (p_source[12] * coeff) >> SPARSE_SCALE_FACTOR;
97 p_dest[13] = (p_source[13] * coeff) >> SPARSE_SCALE_FACTOR;
98 p_dest[14] = (p_source[14] * coeff) >> SPARSE_SCALE_FACTOR;
99 p_dest[15] = (p_source[15] * coeff) >> SPARSE_SCALE_FACTOR;
105 /*****************************************************************************
107 *****************************************************************************/
108 static void SparseIDCTCopy( dctelem_t * p_block, yuv_data_t * p_dest,
109 int i_stride, void * p_idct_data, int i_sparse_pos )
111 SparseIDCT( p_block, p_idct_data, i_sparse_pos );
112 CopyBlock( p_block, p_dest, i_stride );
115 static void SparseIDCTAdd( dctelem_t * p_block, yuv_data_t * p_dest,
116 int i_stride, void * p_idct_data, int i_sparse_pos )
118 SparseIDCT( p_block, p_idct_data, i_sparse_pos );
119 AddBlock( p_block, p_dest, i_stride );