aom_dsp/mips/itrans4_dspr2.c - aom - Git at Google

 /*
  * Copyright (c) 2016, Alliance for Open Media. All rights reserved
  *
  * This source code is subject to the terms of the BSD 2 Clause License and
  * the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
  * was not distributed with this source code in the LICENSE file, you can
  * obtain it at www.aomedia.org/license/software. If the Alliance for Open
  * Media Patent License 1.0 was not distributed with this source code in the
  * PATENTS file, you can obtain it at www.aomedia.org/license/patent.
  */

 #include "./aom_config.h"
 #include "./aom_dsp_rtcd.h"
 #include "aom_dsp/mips/inv_txfm_dspr2.h"
 #include "aom_dsp/txfm_common.h"

 #if HAVE_DSPR2
 void aom_idct4_rows_dspr2(const int16_t *input, int16_t *output) {
   int16_t step_0, step_1, step_2, step_3;
   int Temp0, Temp1, Temp2, Temp3;
   const int const_2_power_13 = 8192;
   int i;

   for (i = 4; i--;) {
     __asm__ __volatile__(
         /*
           temp_1 = (input[0] + input[2]) * cospi_16_64;
           step_0 = dct_const_round_shift(temp_1);

           temp_2 = (input[0] - input[2]) * cospi_16_64;
           step_1 = dct_const_round_shift(temp_2);
         */
         "lh       %[Temp0],             0(%[input])                     \n\t"
         "lh       %[Temp1],             4(%[input])                     \n\t"
         "mtlo     %[const_2_power_13],  $ac0                            \n\t"
         "mthi     $zero,                $ac0                            \n\t"
         "mtlo     %[const_2_power_13],  $ac1                            \n\t"
         "mthi     $zero,                $ac1                            \n\t"
         "add      %[Temp2],             %[Temp0],       %[Temp1]        \n\t"
         "sub      %[Temp3],             %[Temp0],       %[Temp1]        \n\t"
         "madd     $ac0,                 %[Temp2],       %[cospi_16_64]  \n\t"
         "lh       %[Temp0],             2(%[input])                     \n\t"
         "lh       %[Temp1],             6(%[input])                     \n\t"
         "extp     %[step_0],            $ac0,           31              \n\t"
         "mtlo     %[const_2_power_13],  $ac0                            \n\t"
         "mthi     $zero,                $ac0                            \n\t"

         "madd     $ac1,                 %[Temp3],       %[cospi_16_64]  \n\t"
         "extp     %[step_1],            $ac1,           31              \n\t"
         "mtlo     %[const_2_power_13],  $ac1                            \n\t"
         "mthi     $zero,                $ac1                            \n\t"

         /*
           temp1 = input[1] * cospi_24_64 - input[3] * cospi_8_64;
           step_2 = dct_const_round_shift(temp1);
         */
         "madd     $ac0,                 %[Temp0],       %[cospi_24_64]  \n\t"
         "msub     $ac0,                 %[Temp1],       %[cospi_8_64]   \n\t"
         "extp     %[step_2],            $ac0,           31              \n\t"

         /*
           temp2 = input[1] * cospi_8_64 + input[3] * cospi_24_64;
           step_3 = dct_const_round_shift(temp2);
         */
         "madd     $ac1,                 %[Temp0],       %[cospi_8_64]   \n\t"
         "madd     $ac1,                 %[Temp1],       %[cospi_24_64]  \n\t"
         "extp     %[step_3],            $ac1,           31              \n\t"

         /*
           output[0]  = step_0 + step_3;
           output[4]  = step_1 + step_2;
           output[8]  = step_1 - step_2;
           output[12] = step_0 - step_3;
         */
         "add      %[Temp0],             %[step_0],      %[step_3]       \n\t"
         "sh       %[Temp0],             0(%[output])                    \n\t"

         "add      %[Temp1],             %[step_1],      %[step_2]       \n\t"
         "sh       %[Temp1],             8(%[output])                    \n\t"

         "sub      %[Temp2],             %[step_1],      %[step_2]       \n\t"
         "sh       %[Temp2],             16(%[output])                   \n\t"

         "sub      %[Temp3],             %[step_0],      %[step_3]       \n\t"
         "sh       %[Temp3],             24(%[output])                   \n\t"

         : [Temp0] "=&r"(Temp0), [Temp1] "=&r"(Temp1), [Temp2] "=&r"(Temp2),
           [Temp3] "=&r"(Temp3), [step_0] "=&r"(step_0), [step_1] "=&r"(step_1),
           [step_2] "=&r"(step_2), [step_3] "=&r"(step_3), [output] "+r"(output)
         : [const_2_power_13] "r"(const_2_power_13),
           [cospi_8_64] "r"(cospi_8_64), [cospi_16_64] "r"(cospi_16_64),
           [cospi_24_64] "r"(cospi_24_64), [input] "r"(input));

     input += 4;
     output += 1;
   }
 }

 void aom_idct4_columns_add_blk_dspr2(int16_t *input, uint8_t *dest,
                                      int dest_stride) {
   int16_t step_0, step_1, step_2, step_3;
   int Temp0, Temp1, Temp2, Temp3;
   const int const_2_power_13 = 8192;
   int i;
   uint8_t *dest_pix;
   uint8_t *cm = aom_ff_cropTbl;

   /* prefetch aom_ff_cropTbl */
   prefetch_load(aom_ff_cropTbl);
   prefetch_load(aom_ff_cropTbl + 32);
   prefetch_load(aom_ff_cropTbl + 64);
   prefetch_load(aom_ff_cropTbl + 96);
   prefetch_load(aom_ff_cropTbl + 128);
   prefetch_load(aom_ff_cropTbl + 160);
   prefetch_load(aom_ff_cropTbl + 192);
   prefetch_load(aom_ff_cropTbl + 224);

   for (i = 0; i < 4; ++i) {
     dest_pix = (dest + i);

     __asm__ __volatile__(
         /*
           temp_1 = (input[0] + input[2]) * cospi_16_64;
           step_0 = dct_const_round_shift(temp_1);

           temp_2 = (input[0] - input[2]) * cospi_16_64;
           step_1 = dct_const_round_shift(temp_2);
         */
         "lh       %[Temp0],             0(%[input])                     \n\t"
         "lh       %[Temp1],             4(%[input])                     \n\t"
         "mtlo     %[const_2_power_13],  $ac0                            \n\t"
         "mthi     $zero,                $ac0                            \n\t"
         "mtlo     %[const_2_power_13],  $ac1                            \n\t"
         "mthi     $zero,                $ac1                            \n\t"
         "add      %[Temp2],             %[Temp0],       %[Temp1]        \n\t"
         "sub      %[Temp3],             %[Temp0],       %[Temp1]        \n\t"
         "madd     $ac0,                 %[Temp2],       %[cospi_16_64]  \n\t"
         "lh       %[Temp0],             2(%[input])                     \n\t"
         "lh       %[Temp1],             6(%[input])                     \n\t"
         "extp     %[step_0],            $ac0,           31              \n\t"
         "mtlo     %[const_2_power_13],  $ac0                            \n\t"
         "mthi     $zero,                $ac0                            \n\t"

         "madd     $ac1,                 %[Temp3],       %[cospi_16_64]  \n\t"
         "extp     %[step_1],            $ac1,           31              \n\t"
         "mtlo     %[const_2_power_13],  $ac1                            \n\t"
         "mthi     $zero,                $ac1                            \n\t"

         /*
           temp1 = input[1] * cospi_24_64 - input[3] * cospi_8_64;
           step_2 = dct_const_round_shift(temp1);
         */
         "madd     $ac0,                 %[Temp0],       %[cospi_24_64]  \n\t"
         "msub     $ac0,                 %[Temp1],       %[cospi_8_64]   \n\t"
         "extp     %[step_2],            $ac0,           31              \n\t"

         /*
           temp2 = input[1] * cospi_8_64 + input[3] * cospi_24_64;
           step_3 = dct_const_round_shift(temp2);
         */
         "madd     $ac1,                 %[Temp0],       %[cospi_8_64]   \n\t"
         "madd     $ac1,                 %[Temp1],       %[cospi_24_64]  \n\t"
         "extp     %[step_3],            $ac1,           31              \n\t"

         /*
           output[0]  = step_0 + step_3;
           output[4]  = step_1 + step_2;
           output[8]  = step_1 - step_2;
           output[12] = step_0 - step_3;
         */
         "add      %[Temp0],             %[step_0],      %[step_3]       \n\t"
         "addi     %[Temp0],             %[Temp0],       8               \n\t"
         "sra      %[Temp0],             %[Temp0],       4               \n\t"
         "lbu      %[Temp1],             0(%[dest_pix])                  \n\t"
         "add      %[Temp1],             %[Temp1],       %[Temp0]        \n\t"
         "add      %[Temp0],             %[step_1],      %[step_2]       \n\t"
         "lbux     %[Temp2],             %[Temp1](%[cm])                 \n\t"
         "sb       %[Temp2],             0(%[dest_pix])                  \n\t"
         "addu     %[dest_pix],          %[dest_pix],    %[dest_stride]  \n\t"

         "addi     %[Temp0],             %[Temp0],       8               \n\t"
         "sra      %[Temp0],             %[Temp0],       4               \n\t"
         "lbu      %[Temp1],             0(%[dest_pix])                  \n\t"
         "add      %[Temp1],             %[Temp1],       %[Temp0]        \n\t"
         "sub      %[Temp0],             %[step_1],      %[step_2]       \n\t"
         "lbux     %[Temp2],             %[Temp1](%[cm])                 \n\t"
         "sb       %[Temp2],             0(%[dest_pix])                  \n\t"
         "addu     %[dest_pix],          %[dest_pix],    %[dest_stride]  \n\t"

         "addi     %[Temp0],             %[Temp0],       8               \n\t"
         "sra      %[Temp0],             %[Temp0],       4               \n\t"
         "lbu      %[Temp1],             0(%[dest_pix])                  \n\t"
         "add      %[Temp1],             %[Temp1],       %[Temp0]        \n\t"
         "sub      %[Temp0],             %[step_0],      %[step_3]       \n\t"
         "lbux     %[Temp2],             %[Temp1](%[cm])                 \n\t"
         "sb       %[Temp2],             0(%[dest_pix])                  \n\t"
         "addu     %[dest_pix],          %[dest_pix],    %[dest_stride]  \n\t"

         "addi     %[Temp0],             %[Temp0],       8               \n\t"
         "sra      %[Temp0],             %[Temp0],       4               \n\t"
         "lbu      %[Temp1],             0(%[dest_pix])                  \n\t"
         "add      %[Temp1],             %[Temp1],       %[Temp0]        \n\t"
         "lbux     %[Temp2],             %[Temp1](%[cm])                 \n\t"
         "sb       %[Temp2],             0(%[dest_pix])                  \n\t"

         : [Temp0] "=&r"(Temp0), [Temp1] "=&r"(Temp1), [Temp2] "=&r"(Temp2),
           [Temp3] "=&r"(Temp3), [step_0] "=&r"(step_0), [step_1] "=&r"(step_1),
           [step_2] "=&r"(step_2), [step_3] "=&r"(step_3),
           [dest_pix] "+r"(dest_pix)
         : [const_2_power_13] "r"(const_2_power_13),
           [cospi_8_64] "r"(cospi_8_64), [cospi_16_64] "r"(cospi_16_64),
           [cospi_24_64] "r"(cospi_24_64), [input] "r"(input), [cm] "r"(cm),
           [dest_stride] "r"(dest_stride));

     input += 4;
   }
 }

 void aom_idct4x4_16_add_dspr2(const int16_t *input, uint8_t *dest,
                               int dest_stride) {
   DECLARE_ALIGNED(32, int16_t, out[4 * 4]);
   int16_t *outptr = out;
   uint32_t pos = 45;

   /* bit positon for extract from acc */
   __asm__ __volatile__("wrdsp      %[pos],     1           \n\t"
                        :
                        : [pos] "r"(pos));

   // Rows
   aom_idct4_rows_dspr2(input, outptr);

   // Columns
   aom_idct4_columns_add_blk_dspr2(&out[0], dest, dest_stride);
 }

 void aom_idct4x4_1_add_dspr2(const int16_t *input, uint8_t *dest,
                              int dest_stride) {
   int a1, absa1;
   int r;
   int32_t out;
   int t2, vector_a1, vector_a;
   uint32_t pos = 45;
   int16_t input_dc = input[0];

   /* bit positon for extract from acc */
   __asm__ __volatile__("wrdsp      %[pos],     1           \n\t"

                        :
                        : [pos] "r"(pos));

   out = DCT_CONST_ROUND_SHIFT_TWICE_COSPI_16_64(input_dc);
   __asm__ __volatile__(
       "addi     %[out],     %[out],    8       \n\t"
       "sra      %[a1],      %[out],    4       \n\t"

       : [out] "+r"(out), [a1] "=r"(a1)
       :);

   if (a1 < 0) {
     /* use quad-byte
      * input and output memory are four byte aligned */
     __asm__ __volatile__(
         "abs        %[absa1],     %[a1]         \n\t"
         "replv.qb   %[vector_a1], %[absa1]      \n\t"

         : [absa1] "=r"(absa1), [vector_a1] "=r"(vector_a1)
         : [a1] "r"(a1));

     for (r = 4; r--;) {
       __asm__ __volatile__(
           "lw             %[t2],          0(%[dest])                      \n\t"
           "subu_s.qb      %[vector_a],    %[t2],          %[vector_a1]    \n\t"
           "sw             %[vector_a],    0(%[dest])                      \n\t"
           "add            %[dest],        %[dest],        %[dest_stride]  \n\t"

           : [t2] "=&r"(t2), [vector_a] "=&r"(vector_a), [dest] "+&r"(dest)
           : [dest_stride] "r"(dest_stride), [vector_a1] "r"(vector_a1));
     }
   } else {
     /* use quad-byte
      * input and output memory are four byte aligned */
     __asm__ __volatile__("replv.qb       %[vector_a1],   %[a1]     \n\t"
                          : [vector_a1] "=r"(vector_a1)
                          : [a1] "r"(a1));

     for (r = 4; r--;) {
       __asm__ __volatile__(
           "lw           %[t2],          0(%[dest])                        \n\t"
           "addu_s.qb    %[vector_a],    %[t2],            %[vector_a1]    \n\t"
           "sw           %[vector_a],    0(%[dest])                        \n\t"
           "add          %[dest],        %[dest],          %[dest_stride]  \n\t"

           : [t2] "=&r"(t2), [vector_a] "=&r"(vector_a), [dest] "+&r"(dest)
           : [dest_stride] "r"(dest_stride), [vector_a1] "r"(vector_a1));
     }
   }
 }

 void iadst4_dspr2(const int16_t *input, int16_t *output) {
   int s0, s1, s2, s3, s4, s5, s6, s7;
   int x0, x1, x2, x3;

   x0 = input[0];
   x1 = input[1];
   x2 = input[2];
   x3 = input[3];

   if (!(x0 | x1 | x2 | x3)) {
     output[0] = output[1] = output[2] = output[3] = 0;
     return;
   }

   s0 = sinpi_1_9 * x0;
   s1 = sinpi_2_9 * x0;
   s2 = sinpi_3_9 * x1;
   s3 = sinpi_4_9 * x2;
   s4 = sinpi_1_9 * x2;
   s5 = sinpi_2_9 * x3;
   s6 = sinpi_4_9 * x3;
   s7 = x0 - x2 + x3;

   x0 = s0 + s3 + s5;
   x1 = s1 - s4 - s6;
   x2 = sinpi_3_9 * s7;
   x3 = s2;

   s0 = x0 + x3;
   s1 = x1 + x3;
   s2 = x2;
   s3 = x0 + x1 - x3;

   // 1-D transform scaling factor is sqrt(2).
   // The overall dynamic range is 14b (input) + 14b (multiplication scaling)
   // + 1b (addition) = 29b.
   // Hence the output bit depth is 15b.
   output[0] = dct_const_round_shift(s0);
   output[1] = dct_const_round_shift(s1);
   output[2] = dct_const_round_shift(s2);
   output[3] = dct_const_round_shift(s3);
 }
 #endif  // #if HAVE_DSPR2
	/*
	* Copyright (c) 2016, Alliance for Open Media. All rights reserved
	*
	* This source code is subject to the terms of the BSD 2 Clause License and
	* the Alliance for Open Media Patent License 1.0. If the BSD 2 Clause License
	* was not distributed with this source code in the LICENSE file, you can
	* obtain it at www.aomedia.org/license/software. If the Alliance for Open
	* Media Patent License 1.0 was not distributed with this source code in the
	* PATENTS file, you can obtain it at www.aomedia.org/license/patent.
	*/

	#include "./aom_config.h"
	#include "./aom_dsp_rtcd.h"
	#include "aom_dsp/mips/inv_txfm_dspr2.h"
	#include "aom_dsp/txfm_common.h"

	#if HAVE_DSPR2
	void aom_idct4_rows_dspr2(const int16_t input, int16_t output) {
	int16_t step_0, step_1, step_2, step_3;
	int Temp0, Temp1, Temp2, Temp3;
	const int const_2_power_13 = 8192;
	int i;

	for (i = 4; i--;) {
	__asm__ __volatile__(
	/*
	temp_1 = (input[0] + input[2]) * cospi_16_64;
	step_0 = dct_const_round_shift(temp_1);

	temp_2 = (input[0] - input[2]) * cospi_16_64;
	step_1 = dct_const_round_shift(temp_2);
	*/
	"lh %[Temp0], 0(%[input]) \n\t"
	"lh %[Temp1], 4(%[input]) \n\t"
	"mtlo %[const_2_power_13], $ac0 \n\t"
	"mthi $zero, $ac0 \n\t"
	"mtlo %[const_2_power_13], $ac1 \n\t"
	"mthi $zero, $ac1 \n\t"
	"add %[Temp2], %[Temp0], %[Temp1] \n\t"
	"sub %[Temp3], %[Temp0], %[Temp1] \n\t"
	"madd $ac0, %[Temp2], %[cospi_16_64] \n\t"
	"lh %[Temp0], 2(%[input]) \n\t"
	"lh %[Temp1], 6(%[input]) \n\t"
	"extp %[step_0], $ac0, 31 \n\t"
	"mtlo %[const_2_power_13], $ac0 \n\t"
	"mthi $zero, $ac0 \n\t"

	"madd $ac1, %[Temp3], %[cospi_16_64] \n\t"
	"extp %[step_1], $ac1, 31 \n\t"
	"mtlo %[const_2_power_13], $ac1 \n\t"
	"mthi $zero, $ac1 \n\t"

	/*
	temp1 = input[1] * cospi_24_64 - input[3] * cospi_8_64;
	step_2 = dct_const_round_shift(temp1);
	*/
	"madd $ac0, %[Temp0], %[cospi_24_64] \n\t"
	"msub $ac0, %[Temp1], %[cospi_8_64] \n\t"
	"extp %[step_2], $ac0, 31 \n\t"

	/*
	temp2 = input[1] * cospi_8_64 + input[3] * cospi_24_64;
	step_3 = dct_const_round_shift(temp2);
	*/
	"madd $ac1, %[Temp0], %[cospi_8_64] \n\t"
	"madd $ac1, %[Temp1], %[cospi_24_64] \n\t"
	"extp %[step_3], $ac1, 31 \n\t"

	/*
	output[0] = step_0 + step_3;
	output[4] = step_1 + step_2;
	output[8] = step_1 - step_2;
	output[12] = step_0 - step_3;
	*/
	"add %[Temp0], %[step_0], %[step_3] \n\t"
	"sh %[Temp0], 0(%[output]) \n\t"

	"add %[Temp1], %[step_1], %[step_2] \n\t"
	"sh %[Temp1], 8(%[output]) \n\t"

	"sub %[Temp2], %[step_1], %[step_2] \n\t"
	"sh %[Temp2], 16(%[output]) \n\t"

	"sub %[Temp3], %[step_0], %[step_3] \n\t"
	"sh %[Temp3], 24(%[output]) \n\t"

	: [Temp0] "=&r"(Temp0), [Temp1] "=&r"(Temp1), [Temp2] "=&r"(Temp2),
	[Temp3] "=&r"(Temp3), [step_0] "=&r"(step_0), [step_1] "=&r"(step_1),
	[step_2] "=&r"(step_2), [step_3] "=&r"(step_3), [output] "+r"(output)
	: [const_2_power_13] "r"(const_2_power_13),
	[cospi_8_64] "r"(cospi_8_64), [cospi_16_64] "r"(cospi_16_64),
	[cospi_24_64] "r"(cospi_24_64), [input] "r"(input));

	input += 4;
	output += 1;
	}
	}

	void aom_idct4_columns_add_blk_dspr2(int16_t input, uint8_t dest,
	int dest_stride) {
	int16_t step_0, step_1, step_2, step_3;
	int Temp0, Temp1, Temp2, Temp3;
	const int const_2_power_13 = 8192;
	int i;
	uint8_t *dest_pix;
	uint8_t *cm = aom_ff_cropTbl;

	/* prefetch aom_ff_cropTbl */
	prefetch_load(aom_ff_cropTbl);
	prefetch_load(aom_ff_cropTbl + 32);
	prefetch_load(aom_ff_cropTbl + 64);
	prefetch_load(aom_ff_cropTbl + 96);
	prefetch_load(aom_ff_cropTbl + 128);
	prefetch_load(aom_ff_cropTbl + 160);
	prefetch_load(aom_ff_cropTbl + 192);
	prefetch_load(aom_ff_cropTbl + 224);

	for (i = 0; i < 4; ++i) {
	dest_pix = (dest + i);

	__asm__ __volatile__(
	/*
	temp_1 = (input[0] + input[2]) * cospi_16_64;
	step_0 = dct_const_round_shift(temp_1);

	temp_2 = (input[0] - input[2]) * cospi_16_64;
	step_1 = dct_const_round_shift(temp_2);
	*/
	"lh %[Temp0], 0(%[input]) \n\t"
	"lh %[Temp1], 4(%[input]) \n\t"
	"mtlo %[const_2_power_13], $ac0 \n\t"
	"mthi $zero, $ac0 \n\t"
	"mtlo %[const_2_power_13], $ac1 \n\t"
	"mthi $zero, $ac1 \n\t"
	"add %[Temp2], %[Temp0], %[Temp1] \n\t"
	"sub %[Temp3], %[Temp0], %[Temp1] \n\t"
	"madd $ac0, %[Temp2], %[cospi_16_64] \n\t"
	"lh %[Temp0], 2(%[input]) \n\t"
	"lh %[Temp1], 6(%[input]) \n\t"
	"extp %[step_0], $ac0, 31 \n\t"
	"mtlo %[const_2_power_13], $ac0 \n\t"
	"mthi $zero, $ac0 \n\t"

	"madd $ac1, %[Temp3], %[cospi_16_64] \n\t"
	"extp %[step_1], $ac1, 31 \n\t"
	"mtlo %[const_2_power_13], $ac1 \n\t"
	"mthi $zero, $ac1 \n\t"

	/*
	temp1 = input[1] * cospi_24_64 - input[3] * cospi_8_64;
	step_2 = dct_const_round_shift(temp1);
	*/
	"madd $ac0, %[Temp0], %[cospi_24_64] \n\t"
	"msub $ac0, %[Temp1], %[cospi_8_64] \n\t"
	"extp %[step_2], $ac0, 31 \n\t"

	/*
	temp2 = input[1] * cospi_8_64 + input[3] * cospi_24_64;
	step_3 = dct_const_round_shift(temp2);
	*/
	"madd $ac1, %[Temp0], %[cospi_8_64] \n\t"
	"madd $ac1, %[Temp1], %[cospi_24_64] \n\t"
	"extp %[step_3], $ac1, 31 \n\t"

	/*
	output[0] = step_0 + step_3;
	output[4] = step_1 + step_2;
	output[8] = step_1 - step_2;
	output[12] = step_0 - step_3;
	*/
	"add %[Temp0], %[step_0], %[step_3] \n\t"
	"addi %[Temp0], %[Temp0], 8 \n\t"
	"sra %[Temp0], %[Temp0], 4 \n\t"
	"lbu %[Temp1], 0(%[dest_pix]) \n\t"
	"add %[Temp1], %[Temp1], %[Temp0] \n\t"
	"add %[Temp0], %[step_1], %[step_2] \n\t"
	"lbux %[Temp2], %[Temp1](%[cm]) \n\t"
	"sb %[Temp2], 0(%[dest_pix]) \n\t"
	"addu %[dest_pix], %[dest_pix], %[dest_stride] \n\t"

	"addi %[Temp0], %[Temp0], 8 \n\t"
	"sra %[Temp0], %[Temp0], 4 \n\t"
	"lbu %[Temp1], 0(%[dest_pix]) \n\t"
	"add %[Temp1], %[Temp1], %[Temp0] \n\t"
	"sub %[Temp0], %[step_1], %[step_2] \n\t"
	"lbux %[Temp2], %[Temp1](%[cm]) \n\t"
	"sb %[Temp2], 0(%[dest_pix]) \n\t"
	"addu %[dest_pix], %[dest_pix], %[dest_stride] \n\t"

	"addi %[Temp0], %[Temp0], 8 \n\t"
	"sra %[Temp0], %[Temp0], 4 \n\t"
	"lbu %[Temp1], 0(%[dest_pix]) \n\t"
	"add %[Temp1], %[Temp1], %[Temp0] \n\t"
	"sub %[Temp0], %[step_0], %[step_3] \n\t"
	"lbux %[Temp2], %[Temp1](%[cm]) \n\t"
	"sb %[Temp2], 0(%[dest_pix]) \n\t"
	"addu %[dest_pix], %[dest_pix], %[dest_stride] \n\t"

	"addi %[Temp0], %[Temp0], 8 \n\t"
	"sra %[Temp0], %[Temp0], 4 \n\t"
	"lbu %[Temp1], 0(%[dest_pix]) \n\t"
	"add %[Temp1], %[Temp1], %[Temp0] \n\t"
	"lbux %[Temp2], %[Temp1](%[cm]) \n\t"
	"sb %[Temp2], 0(%[dest_pix]) \n\t"

	: [Temp0] "=&r"(Temp0), [Temp1] "=&r"(Temp1), [Temp2] "=&r"(Temp2),
	[Temp3] "=&r"(Temp3), [step_0] "=&r"(step_0), [step_1] "=&r"(step_1),
	[step_2] "=&r"(step_2), [step_3] "=&r"(step_3),
	[dest_pix] "+r"(dest_pix)
	: [const_2_power_13] "r"(const_2_power_13),
	[cospi_8_64] "r"(cospi_8_64), [cospi_16_64] "r"(cospi_16_64),
	[cospi_24_64] "r"(cospi_24_64), [input] "r"(input), [cm] "r"(cm),
	[dest_stride] "r"(dest_stride));

	input += 4;
	}
	}

	void aom_idct4x4_16_add_dspr2(const int16_t input, uint8_t dest,
	int dest_stride) {
	DECLARE_ALIGNED(32, int16_t, out[4 * 4]);
	int16_t *outptr = out;
	uint32_t pos = 45;

	/* bit positon for extract from acc */
	__asm__ __volatile__("wrdsp %[pos], 1 \n\t"
	:
	: [pos] "r"(pos));

	// Rows
	aom_idct4_rows_dspr2(input, outptr);

	// Columns
	aom_idct4_columns_add_blk_dspr2(&out[0], dest, dest_stride);
	}

	void aom_idct4x4_1_add_dspr2(const int16_t input, uint8_t dest,
	int dest_stride) {
	int a1, absa1;
	int r;
	int32_t out;
	int t2, vector_a1, vector_a;
	uint32_t pos = 45;
	int16_t input_dc = input[0];

	/* bit positon for extract from acc */
	__asm__ __volatile__("wrdsp %[pos], 1 \n\t"

	:
	: [pos] "r"(pos));

	out = DCT_CONST_ROUND_SHIFT_TWICE_COSPI_16_64(input_dc);
	__asm__ __volatile__(
	"addi %[out], %[out], 8 \n\t"
	"sra %[a1], %[out], 4 \n\t"

	: [out] "+r"(out), [a1] "=r"(a1)
	:);

	if (a1 < 0) {
	/* use quad-byte
	* input and output memory are four byte aligned */
	__asm__ __volatile__(
	"abs %[absa1], %[a1] \n\t"
	"replv.qb %[vector_a1], %[absa1] \n\t"

	: [absa1] "=r"(absa1), [vector_a1] "=r"(vector_a1)
	: [a1] "r"(a1));

	for (r = 4; r--;) {
	__asm__ __volatile__(
	"lw %[t2], 0(%[dest]) \n\t"
	"subu_s.qb %[vector_a], %[t2], %[vector_a1] \n\t"
	"sw %[vector_a], 0(%[dest]) \n\t"
	"add %[dest], %[dest], %[dest_stride] \n\t"

	: [t2] "=&r"(t2), [vector_a] "=&r"(vector_a), [dest] "+&r"(dest)
	: [dest_stride] "r"(dest_stride), [vector_a1] "r"(vector_a1));
	}
	} else {
	/* use quad-byte
	* input and output memory are four byte aligned */
	__asm__ __volatile__("replv.qb %[vector_a1], %[a1] \n\t"
	: [vector_a1] "=r"(vector_a1)
	: [a1] "r"(a1));

	for (r = 4; r--;) {
	__asm__ __volatile__(
	"lw %[t2], 0(%[dest]) \n\t"
	"addu_s.qb %[vector_a], %[t2], %[vector_a1] \n\t"
	"sw %[vector_a], 0(%[dest]) \n\t"
	"add %[dest], %[dest], %[dest_stride] \n\t"

	: [t2] "=&r"(t2), [vector_a] "=&r"(vector_a), [dest] "+&r"(dest)
	: [dest_stride] "r"(dest_stride), [vector_a1] "r"(vector_a1));
	}
	}
	}

	void iadst4_dspr2(const int16_t input, int16_t output) {
	int s0, s1, s2, s3, s4, s5, s6, s7;
	int x0, x1, x2, x3;

	x0 = input[0];
	x1 = input[1];
	x2 = input[2];
	x3 = input[3];

	if (!(x0 \| x1 \| x2 \| x3)) {
	output[0] = output[1] = output[2] = output[3] = 0;
	return;
	}

	s0 = sinpi_1_9 * x0;
	s1 = sinpi_2_9 * x0;
	s2 = sinpi_3_9 * x1;
	s3 = sinpi_4_9 * x2;
	s4 = sinpi_1_9 * x2;
	s5 = sinpi_2_9 * x3;
	s6 = sinpi_4_9 * x3;
	s7 = x0 - x2 + x3;

	x0 = s0 + s3 + s5;
	x1 = s1 - s4 - s6;
	x2 = sinpi_3_9 * s7;
	x3 = s2;

	s0 = x0 + x3;
	s1 = x1 + x3;
	s2 = x2;
	s3 = x0 + x1 - x3;

	// 1-D transform scaling factor is sqrt(2).
	// The overall dynamic range is 14b (input) + 14b (multiplication scaling)
	// + 1b (addition) = 29b.
	// Hence the output bit depth is 15b.
	output[0] = dct_const_round_shift(s0);
	output[1] = dct_const_round_shift(s1);
	output[2] = dct_const_round_shift(s2);
	output[3] = dct_const_round_shift(s3);
	}
	#endif // #if HAVE_DSPR2