blob: a75483a3dd730460ee1cc3a5accd33319ffa645a [file] [log] [blame]
/* ------------------------------------------------------------------
* Copyright (C) 1998-2009 PacketVideo
*
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either
* express or implied.
* See the License for the specific language governing permissions
* and limitations under the License.
* -------------------------------------------------------------------
*/
/*
------------------------------------------------------------------------------
INPUT AND OUTPUT DEFINITIONS
Inputs:
[input_variable_name] = [description of the input to module, its type
definition, and length (when applicable)]
Local Stores/Buffers/Pointers Needed:
[local_store_name] = [description of the local store, its type
definition, and length (when applicable)]
[local_buffer_name] = [description of the local buffer, its type
definition, and length (when applicable)]
[local_ptr_name] = [description of the local pointer, its type
definition, and length (when applicable)]
Global Stores/Buffers/Pointers Needed:
[global_store_name] = [description of the global store, its type
definition, and length (when applicable)]
[global_buffer_name] = [description of the global buffer, its type
definition, and length (when applicable)]
[global_ptr_name] = [description of the global pointer, its type
definition, and length (when applicable)]
Outputs:
[return_variable_name] = [description of data/pointer returned
by module, its type definition, and length
(when applicable)]
Pointers and Buffers Modified:
[variable_bfr_ptr] points to the [describe where the
variable_bfr_ptr points to, its type definition, and length
(when applicable)]
[variable_bfr] contents are [describe the new contents of
variable_bfr]
Local Stores Modified:
[local_store_name] = [describe new contents, its type
definition, and length (when applicable)]
Global Stores Modified:
[global_store_name] = [describe new contents, its type
definition, and length (when applicable)]
------------------------------------------------------------------------------
FUNCTION DESCRIPTION
------------------------------------------------------------------------------
REQUIREMENTS
------------------------------------------------------------------------------
REFERENCES
------------------------------------------------------------------------------
PSEUDO-CODE
------------------------------------------------------------------------------
RESOURCES USED
When the code is written for a specific target processor the
the resources used should be documented below.
STACK USAGE: [stack count for this module] + [variable to represent
stack usage for each subroutine called]
where: [stack usage variable] = stack usage for [subroutine
name] (see [filename].ext)
DATA MEMORY USED: x words
PROGRAM MEMORY USED: x words
CLOCK CYCLES: [cycle count equation for this module] + [variable
used to represent cycle count for each subroutine
called]
where: [cycle count variable] = cycle count for [subroutine
name] (see [filename].ext)
------------------------------------------------------------------------------
*/
/*----------------------------------------------------------------------------
; INCLUDES
----------------------------------------------------------------------------*/
#include "mp4dec_lib.h"
#include "idct.h"
#include "motion_comp.h"
#define OSCL_DISABLE_WARNING_CONV_POSSIBLE_LOSS_OF_DATA
/*----------------------------------------------------------------------------
; MACROS
; Define module specific macros here
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; DEFINES
; Include all pre-processor statements here. Include conditional
; compile variables also.
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; LOCAL FUNCTION DEFINITIONS
; Function Prototype declaration
----------------------------------------------------------------------------*/
/* private prototypes */
static void idctrow(int16 *blk, uint8 *pred, uint8 *dst, int width);
static void idctrow_intra(int16 *blk, PIXEL *, int width);
static void idctcol(int16 *blk);
#ifdef FAST_IDCT
// mapping from nz_coefs to functions to be used
// ARM4 does not allow global data when they are not constant hence
// an array of function pointers cannot be considered as array of constants
// (actual addresses are only known when the dll is loaded).
// So instead of arrays of function pointers, we'll store here
// arrays of rows or columns and then call the idct function
// corresponding to such the row/column number:
static void (*const idctcolVCA[10][4])(int16*) =
{
{&idctcol1, &idctcol0, &idctcol0, &idctcol0},
{&idctcol1, &idctcol1, &idctcol0, &idctcol0},
{&idctcol2, &idctcol1, &idctcol0, &idctcol0},
{&idctcol3, &idctcol1, &idctcol0, &idctcol0},
{&idctcol3, &idctcol2, &idctcol0, &idctcol0},
{&idctcol3, &idctcol2, &idctcol1, &idctcol0},
{&idctcol3, &idctcol2, &idctcol1, &idctcol1},
{&idctcol3, &idctcol2, &idctcol2, &idctcol1},
{&idctcol3, &idctcol3, &idctcol2, &idctcol1},
{&idctcol4, &idctcol3, &idctcol2, &idctcol1}
};
static void (*const idctrowVCA[10])(int16*, uint8*, uint8*, int) =
{
&idctrow1,
&idctrow2,
&idctrow2,
&idctrow2,
&idctrow2,
&idctrow3,
&idctrow4,
&idctrow4,
&idctrow4,
&idctrow4
};
static void (*const idctcolVCA2[16])(int16*) =
{
&idctcol0, &idctcol4, &idctcol3, &idctcol4,
&idctcol2, &idctcol4, &idctcol3, &idctcol4,
&idctcol1, &idctcol4, &idctcol3, &idctcol4,
&idctcol2, &idctcol4, &idctcol3, &idctcol4
};
static void (*const idctrowVCA2[8])(int16*, uint8*, uint8*, int) =
{
&idctrow1, &idctrow4, &idctrow3, &idctrow4,
&idctrow2, &idctrow4, &idctrow3, &idctrow4
};
static void (*const idctrowVCA_intra[10])(int16*, PIXEL *, int) =
{
&idctrow1_intra,
&idctrow2_intra,
&idctrow2_intra,
&idctrow2_intra,
&idctrow2_intra,
&idctrow3_intra,
&idctrow4_intra,
&idctrow4_intra,
&idctrow4_intra,
&idctrow4_intra
};
static void (*const idctrowVCA2_intra[8])(int16*, PIXEL *, int) =
{
&idctrow1_intra, &idctrow4_intra, &idctrow3_intra, &idctrow4_intra,
&idctrow2_intra, &idctrow4_intra, &idctrow3_intra, &idctrow4_intra
};
#endif
/*----------------------------------------------------------------------------
; LOCAL STORE/BUFFER/POINTER DEFINITIONS
; Variable declaration - defined here and used outside this module
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; EXTERNAL FUNCTION REFERENCES
; Declare functions defined elsewhere and referenced in this module
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; EXTERNAL GLOBAL STORE/BUFFER/POINTER REFERENCES
; Declare variables used in this module but defined elsewhere
----------------------------------------------------------------------------*/
/*----------------------------------------------------------------------------
; FUNCTION CODE
----------------------------------------------------------------------------*/
void MBlockIDCT(VideoDecData *video)
{
Vop *currVop = video->currVop;
MacroBlock *mblock = video->mblock;
PIXEL *c_comp;
PIXEL *cu_comp;
PIXEL *cv_comp;
int x_pos = video->mbnum_col;
int y_pos = video->mbnum_row;
int width, width_uv;
int32 offset;
width = video->width;
width_uv = width >> 1;
offset = (int32)(y_pos << 4) * width + (x_pos << 4);
c_comp = currVop->yChan + offset;
cu_comp = currVop->uChan + (offset >> 2) + (x_pos << 2);
cv_comp = currVop->vChan + (offset >> 2) + (x_pos << 2);
BlockIDCT_intra(mblock, c_comp, 0, width);
BlockIDCT_intra(mblock, c_comp + 8, 1, width);
BlockIDCT_intra(mblock, c_comp + (width << 3), 2, width);
BlockIDCT_intra(mblock, c_comp + (width << 3) + 8, 3, width);
BlockIDCT_intra(mblock, cu_comp, 4, width_uv);
BlockIDCT_intra(mblock, cv_comp, 5, width_uv);
}
void BlockIDCT_intra(
MacroBlock *mblock, PIXEL *c_comp, int comp, int width)
{
/*----------------------------------------------------------------------------
; Define all local variables
----------------------------------------------------------------------------*/
int16 *coeff_in = mblock->block[comp];
#ifdef INTEGER_IDCT
#ifdef FAST_IDCT /* VCA IDCT using nzcoefs and bitmaps*/
int i, bmapr;
int nz_coefs = mblock->no_coeff[comp];
uint8 *bitmapcol = mblock->bitmapcol[comp];
uint8 bitmaprow = mblock->bitmaprow[comp];
/*----------------------------------------------------------------------------
; Function body here
----------------------------------------------------------------------------*/
if (nz_coefs <= 10)
{
bmapr = (nz_coefs - 1);
(*(idctcolVCA[bmapr]))(coeff_in);
(*(idctcolVCA[bmapr][1]))(coeff_in + 1);
(*(idctcolVCA[bmapr][2]))(coeff_in + 2);
(*(idctcolVCA[bmapr][3]))(coeff_in + 3);
(*idctrowVCA_intra[nz_coefs-1])(coeff_in, c_comp, width);
}
else
{
i = 8;
while (i--)
{
bmapr = (int)bitmapcol[i];
if (bmapr)
{
if ((bmapr&0xf) == 0) /* 07/18/01 */
{
(*(idctcolVCA2[bmapr>>4]))(coeff_in + i);
}
else
{
idctcol(coeff_in + i);
}
}
}
if ((bitmapcol[4] | bitmapcol[5] | bitmapcol[6] | bitmapcol[7]) == 0)
{
bitmaprow >>= 4;
(*(idctrowVCA2_intra[(int)bitmaprow]))(coeff_in, c_comp, width);
}
else
{
idctrow_intra(coeff_in, c_comp, width);
}
}
#else
void idct_intra(int *block, uint8 *comp, int width);
idct_intra(coeff_in, c_comp, width);
#endif
#else
void idctref_intra(int *block, uint8 *comp, int width);
idctref_intra(coeff_in, c_comp, width);
#endif
/*----------------------------------------------------------------------------
; Return nothing or data or data pointer
----------------------------------------------------------------------------*/
return;
}
/* 08/04/05, no residue, just copy from pred to output */
void Copy_Blk_to_Vop(uint8 *dst, uint8 *pred, int width)
{
/* copy 4 bytes at a time */
width -= 4;
*((uint32*)dst) = *((uint32*)pred);
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
*((uint32*)(dst += width)) = *((uint32*)(pred += 12));
*((uint32*)(dst += 4)) = *((uint32*)(pred += 4));
return ;
}
/* 08/04/05 compute IDCT and add prediction at the end */
void BlockIDCT(
uint8 *dst, /* destination */
uint8 *pred, /* prediction block, pitch 16 */
int16 *coeff_in, /* DCT data, size 64 */
int width, /* width of dst */
int nz_coefs,
uint8 *bitmapcol,
uint8 bitmaprow
)
{
#ifdef INTEGER_IDCT
#ifdef FAST_IDCT /* VCA IDCT using nzcoefs and bitmaps*/
int i, bmapr;
/*----------------------------------------------------------------------------
; Function body here
----------------------------------------------------------------------------*/
if (nz_coefs <= 10)
{
bmapr = (nz_coefs - 1);
(*(idctcolVCA[bmapr]))(coeff_in);
(*(idctcolVCA[bmapr][1]))(coeff_in + 1);
(*(idctcolVCA[bmapr][2]))(coeff_in + 2);
(*(idctcolVCA[bmapr][3]))(coeff_in + 3);
(*idctrowVCA[nz_coefs-1])(coeff_in, pred, dst, width);
return ;
}
else
{
i = 8;
while (i--)
{
bmapr = (int)bitmapcol[i];
if (bmapr)
{
if ((bmapr&0xf) == 0) /* 07/18/01 */
{
(*(idctcolVCA2[bmapr>>4]))(coeff_in + i);
}
else
{
idctcol(coeff_in + i);
}
}
}
if ((bitmapcol[4] | bitmapcol[5] | bitmapcol[6] | bitmapcol[7]) == 0)
{
(*(idctrowVCA2[bitmaprow>>4]))(coeff_in, pred, dst, width);
}
else
{
idctrow(coeff_in, pred, dst, width);
}
return ;
}
#else // FAST_IDCT
void idct(int *block, uint8 *pred, uint8 *dst, int width);
idct(coeff_in, pred, dst, width);
return;
#endif // FAST_IDCT
#else // INTEGER_IDCT
void idctref(int *block, uint8 *pred, uint8 *dst, int width);
idctref(coeff_in, pred, dst, width);
return;
#endif // INTEGER_IDCT
}
/*----------------------------------------------------------------------------
; End Function: block_idct
----------------------------------------------------------------------------*/
/****************************************************************************/
/*
------------------------------------------------------------------------------
FUNCTION NAME: idctrow
------------------------------------------------------------------------------
INPUT AND OUTPUT DEFINITIONS FOR idctrow
Inputs:
[input_variable_name] = [description of the input to module, its type
definition, and length (when applicable)]
Local Stores/Buffers/Pointers Needed:
[local_store_name] = [description of the local store, its type
definition, and length (when applicable)]
[local_buffer_name] = [description of the local buffer, its type
definition, and length (when applicable)]
[local_ptr_name] = [description of the local pointer, its type
definition, and length (when applicable)]
Global Stores/Buffers/Pointers Needed:
[global_store_name] = [description of the global store, its type
definition, and length (when applicable)]
[global_buffer_name] = [description of the global buffer, its type
definition, and length (when applicable)]
[global_ptr_name] = [description of the global pointer, its type
definition, and length (when applicable)]
Outputs:
[return_variable_name] = [description of data/pointer returned
by module, its type definition, and length
(when applicable)]
Pointers and Buffers Modified:
[variable_bfr_ptr] points to the [describe where the
variable_bfr_ptr points to, its type definition, and length
(when applicable)]
[variable_bfr] contents are [describe the new contents of
variable_bfr]
Local Stores Modified:
[local_store_name] = [describe new contents, its type
definition, and length (when applicable)]
Global Stores Modified:
[global_store_name] = [describe new contents, its type
definition, and length (when applicable)]
------------------------------------------------------------------------------
FUNCTION DESCRIPTION FOR idctrow
------------------------------------------------------------------------------
REQUIREMENTS FOR idctrow
------------------------------------------------------------------------------
REFERENCES FOR idctrow
------------------------------------------------------------------------------
PSEUDO-CODE FOR idctrow
------------------------------------------------------------------------------
RESOURCES USED FOR idctrow
When the code is written for a specific target processor the
the resources used should be documented below.
STACK USAGE: [stack count for this module] + [variable to represent
stack usage for each subroutine called]
where: [stack usage variable] = stack usage for [subroutine
name] (see [filename].ext)
DATA MEMORY USED: x words
PROGRAM MEMORY USED: x words
CLOCK CYCLES: [cycle count equation for this module] + [variable
used to represent cycle count for each subroutine
called]
where: [cycle count variable] = cycle count for [subroutine
name] (see [filename].ext)
------------------------------------------------------------------------------
*/
/*----------------------------------------------------------------------------
; Function Code FOR idctrow
----------------------------------------------------------------------------*/
void idctrow(
int16 *blk, uint8 *pred, uint8 *dst, int width
)
{
/*----------------------------------------------------------------------------
; Define all local variables
----------------------------------------------------------------------------*/
int32 x0, x1, x2, x3, x4, x5, x6, x7, x8;
int i = 8;
uint32 pred_word, dst_word;
int res, res2;
/*----------------------------------------------------------------------------
; Function body here
----------------------------------------------------------------------------*/
/* row (horizontal) IDCT
*
* 7 pi 1 dst[k] = sum c[l] * src[l] * cos( -- *
* ( k + - ) * l ) l=0 8 2
*
* where: c[0] = 128 c[1..7] = 128*sqrt(2) */
/* preset the offset, such that we can take advantage pre-offset addressing mode */
width -= 4;
dst -= width;
pred -= 12;
blk -= 8;
while (i--)
{
x1 = (int32)blk[12] << 8;
blk[12] = 0;
x2 = blk[14];
blk[14] = 0;
x3 = blk[10];
blk[10] = 0;
x4 = blk[9];
blk[9] = 0;
x5 = blk[15];
blk[15] = 0;
x6 = blk[13];
blk[13] = 0;
x7 = blk[11];
blk[11] = 0;
x0 = ((*(blk += 8)) << 8) + 8192;
blk[0] = 0; /* for proper rounding in the fourth stage */
/* first stage */
x8 = W7 * (x4 + x5) + 4;
x4 = (x8 + (W1 - W7) * x4) >> 3;
x5 = (x8 - (W1 + W7) * x5) >> 3;
x8 = W3 * (x6 + x7) + 4;
x6 = (x8 - (W3 - W5) * x6) >> 3;
x7 = (x8 - (W3 + W5) * x7) >> 3;
/* second stage */
x8 = x0 + x1;
x0 -= x1;
x1 = W6 * (x3 + x2) + 4;
x2 = (x1 - (W2 + W6) * x2) >> 3;
x3 = (x1 + (W2 - W6) * x3) >> 3;
x1 = x4 + x6;
x4 -= x6;
x6 = x5 + x7;
x5 -= x7;
/* third stage */
x7 = x8 + x3;
x8 -= x3;
x3 = x0 + x2;
x0 -= x2;
x2 = (181 * (x4 + x5) + 128) >> 8;
x4 = (181 * (x4 - x5) + 128) >> 8;
/* fourth stage */
pred_word = *((uint32*)(pred += 12)); /* read 4 bytes from pred */
res = (x7 + x1) >> 14;
ADD_AND_CLIP1(res);
res2 = (x3 + x2) >> 14;
ADD_AND_CLIP2(res2);
dst_word = (res2 << 8) | res;
res = (x0 + x4) >> 14;
ADD_AND_CLIP3(res);
dst_word |= (res << 16);
res = (x8 + x6) >> 14;
ADD_AND_CLIP4(res);
dst_word |= (res << 24);
*((uint32*)(dst += width)) = dst_word; /* save 4 bytes to dst */
pred_word = *((uint32*)(pred += 4)); /* read 4 bytes from pred */
res = (x8 - x6) >> 14;
ADD_AND_CLIP1(res);
res2 = (x0 - x4) >> 14;
ADD_AND_CLIP2(res2);
dst_word = (res2 << 8) | res;
res = (x3 - x2) >> 14;
ADD_AND_CLIP3(res);
dst_word |= (res << 16);
res = (x7 - x1) >> 14;
ADD_AND_CLIP4(res);
dst_word |= (res << 24);
*((uint32*)(dst += 4)) = dst_word; /* save 4 bytes to dst */
}
/*----------------------------------------------------------------------------
; Return nothing or data or data pointer
----------------------------------------------------------------------------*/
return;
}
void idctrow_intra(
int16 *blk, PIXEL *comp, int width
)
{
/*----------------------------------------------------------------------------
; Define all local variables
----------------------------------------------------------------------------*/
int32 x0, x1, x2, x3, x4, x5, x6, x7, x8, temp;
int i = 8;
int offset = width;
int32 word;
/*----------------------------------------------------------------------------
; Function body here
----------------------------------------------------------------------------*/
/* row (horizontal) IDCT
*
* 7 pi 1 dst[k] = sum c[l] * src[l] * cos( -- *
* ( k + - ) * l ) l=0 8 2
*
* where: c[0] = 128 c[1..7] = 128*sqrt(2) */
while (i--)
{
x1 = (int32)blk[4] << 8;
blk[4] = 0;
x2 = blk[6];
blk[6] = 0;
x3 = blk[2];
blk[2] = 0;
x4 = blk[1];
blk[1] = 0;
x5 = blk[7];
blk[7] = 0;
x6 = blk[5];
blk[5] = 0;
x7 = blk[3];
blk[3] = 0;
#ifndef FAST_IDCT
/* shortcut */ /* covered by idctrow1 01/9/2001 */
if (!(x1 | x2 | x3 | x4 | x5 | x6 | x7))
{
blk[0] = blk[1] = blk[2] = blk[3] = blk[4] = blk[5] = blk[6] = blk[7] = (blk[0] + 32) >> 6;
return;
}
#endif
x0 = ((int32)blk[0] << 8) + 8192;
blk[0] = 0; /* for proper rounding in the fourth stage */
/* first stage */
x8 = W7 * (x4 + x5) + 4;
x4 = (x8 + (W1 - W7) * x4) >> 3;
x5 = (x8 - (W1 + W7) * x5) >> 3;
x8 = W3 * (x6 + x7) + 4;
x6 = (x8 - (W3 - W5) * x6) >> 3;
x7 = (x8 - (W3 + W5) * x7) >> 3;
/* second stage */
x8 = x0 + x1;
x0 -= x1;
x1 = W6 * (x3 + x2) + 4;
x2 = (x1 - (W2 + W6) * x2) >> 3;
x3 = (x1 + (W2 - W6) * x3) >> 3;
x1 = x4 + x6;
x4 -= x6;
x6 = x5 + x7;
x5 -= x7;
/* third stage */
x7 = x8 + x3;
x8 -= x3;
x3 = x0 + x2;
x0 -= x2;
x2 = (181 * (x4 + x5) + 128) >> 8;
x4 = (181 * (x4 - x5) + 128) >> 8;
/* fourth stage */
word = ((x7 + x1) >> 14);
CLIP_RESULT(word)
temp = ((x3 + x2) >> 14);
CLIP_RESULT(temp)
word = word | (temp << 8);
temp = ((x0 + x4) >> 14);
CLIP_RESULT(temp)
word = word | (temp << 16);
temp = ((x8 + x6) >> 14);
CLIP_RESULT(temp)
word = word | (temp << 24);
*((int32*)(comp)) = word;
word = ((x8 - x6) >> 14);
CLIP_RESULT(word)
temp = ((x0 - x4) >> 14);
CLIP_RESULT(temp)
word = word | (temp << 8);
temp = ((x3 - x2) >> 14);
CLIP_RESULT(temp)
word = word | (temp << 16);
temp = ((x7 - x1) >> 14);
CLIP_RESULT(temp)
word = word | (temp << 24);
*((int32*)(comp + 4)) = word;
comp += offset;
blk += B_SIZE;
}
/*----------------------------------------------------------------------------
; Return nothing or data or data pointer
----------------------------------------------------------------------------*/
return;
}
/*----------------------------------------------------------------------------
; End Function: idctrow
----------------------------------------------------------------------------*/
/****************************************************************************/
/*
------------------------------------------------------------------------------
FUNCTION NAME: idctcol
------------------------------------------------------------------------------
INPUT AND OUTPUT DEFINITIONS FOR idctcol
Inputs:
[input_variable_name] = [description of the input to module, its type
definition, and length (when applicable)]
Local Stores/Buffers/Pointers Needed:
[local_store_name] = [description of the local store, its type
definition, and length (when applicable)]
[local_buffer_name] = [description of the local buffer, its type
definition, and length (when applicable)]
[local_ptr_name] = [description of the local pointer, its type
definition, and length (when applicable)]
Global Stores/Buffers/Pointers Needed:
[global_store_name] = [description of the global store, its type
definition, and length (when applicable)]
[global_buffer_name] = [description of the global buffer, its type
definition, and length (when applicable)]
[global_ptr_name] = [description of the global pointer, its type
definition, and length (when applicable)]
Outputs:
[return_variable_name] = [description of data/pointer returned
by module, its type definition, and length
(when applicable)]
Pointers and Buffers Modified:
[variable_bfr_ptr] points to the [describe where the
variable_bfr_ptr points to, its type definition, and length
(when applicable)]
[variable_bfr] contents are [describe the new contents of
variable_bfr]
Local Stores Modified:
[local_store_name] = [describe new contents, its type
definition, and length (when applicable)]
Global Stores Modified:
[global_store_name] = [describe new contents, its type
definition, and length (when applicable)]
------------------------------------------------------------------------------
FUNCTION DESCRIPTION FOR idctcol
------------------------------------------------------------------------------
REQUIREMENTS FOR idctcol
------------------------------------------------------------------------------
REFERENCES FOR idctcol
------------------------------------------------------------------------------
PSEUDO-CODE FOR idctcol
------------------------------------------------------------------------------
RESOURCES USED FOR idctcol
When the code is written for a specific target processor the
the resources used should be documented below.
STACK USAGE: [stack count for this module] + [variable to represent
stack usage for each subroutine called]
where: [stack usage variable] = stack usage for [subroutine
name] (see [filename].ext)
DATA MEMORY USED: x words
PROGRAM MEMORY USED: x words
CLOCK CYCLES: [cycle count equation for this module] + [variable
used to represent cycle count for each subroutine
called]
where: [cycle count variable] = cycle count for [subroutine
name] (see [filename].ext)
------------------------------------------------------------------------------
*/
/*----------------------------------------------------------------------------
; Function Code FOR idctcol
----------------------------------------------------------------------------*/
void idctcol(
int16 *blk
)
{
/*----------------------------------------------------------------------------
; Define all local variables
----------------------------------------------------------------------------*/
int32 x0, x1, x2, x3, x4, x5, x6, x7, x8;
/*----------------------------------------------------------------------------
; Function body here
----------------------------------------------------------------------------*/
/* column (vertical) IDCT
*
* 7 pi 1 dst[8*k] = sum c[l] * src[8*l] *
* cos( -- * ( k + - ) * l ) l=0 8 2
*
* where: c[0] = 1/1024 c[1..7] = (1/1024)*sqrt(2) */
x1 = (int32)blk[32] << 11;
x2 = blk[48];
x3 = blk[16];
x4 = blk[8];
x5 = blk[56];
x6 = blk[40];
x7 = blk[24];
#ifndef FAST_IDCT
/* shortcut */ /* covered by idctcolumn1 01/9/2001 */
if (!(x1 | x2 | x3 | x4 | x5 | x6 | x7))
{
blk[0] = blk[8] = blk[16] = blk[24] = blk[32] = blk[40] = blk[48] = blk[56]
= blk[0] << 3;
return;
}
#endif
x0 = ((int32)blk[0] << 11) + 128;
/* first stage */
x8 = W7 * (x4 + x5);
x4 = x8 + (W1 - W7) * x4;
x5 = x8 - (W1 + W7) * x5;
x8 = W3 * (x6 + x7);
x6 = x8 - (W3 - W5) * x6;
x7 = x8 - (W3 + W5) * x7;
/* second stage */
x8 = x0 + x1;
x0 -= x1;
x1 = W6 * (x3 + x2);
x2 = x1 - (W2 + W6) * x2;
x3 = x1 + (W2 - W6) * x3;
x1 = x4 + x6;
x4 -= x6;
x6 = x5 + x7;
x5 -= x7;
/* third stage */
x7 = x8 + x3;
x8 -= x3;
x3 = x0 + x2;
x0 -= x2;
x2 = (181 * (x4 + x5) + 128) >> 8;
x4 = (181 * (x4 - x5) + 128) >> 8;
/* fourth stage */
blk[0] = (x7 + x1) >> 8;
blk[8] = (x3 + x2) >> 8;
blk[16] = (x0 + x4) >> 8;
blk[24] = (x8 + x6) >> 8;
blk[32] = (x8 - x6) >> 8;
blk[40] = (x0 - x4) >> 8;
blk[48] = (x3 - x2) >> 8;
blk[56] = (x7 - x1) >> 8;
/*----------------------------------------------------------------------------
; Return nothing or data or data pointer
----------------------------------------------------------------------------*/
return;
}
/*----------------------------------------------------------------------------
; End Function: idctcol
----------------------------------------------------------------------------*/