ref: 5d3327e891da9b8163f0150d3e52cc9f89016353
dir: /vpx_dsp/mips/itrans32_dspr2.c/
/* * Copyright (c) 2013 The WebM project authors. All Rights Reserved. * * Use of this source code is governed by a BSD-style license * that can be found in the LICENSE file in the root of the source * tree. An additional intellectual property rights grant can be found * in the file PATENTS. All contributing project authors may * be found in the AUTHORS file in the root of the source tree. */ #include <assert.h> #include <stdio.h> #include "./vpx_config.h" #include "vpx_dsp/mips/inv_txfm_dspr2.h" #include "vpx_dsp/txfm_common.h" #if HAVE_DSPR2 static void idct32_rows_dspr2(const int16_t *input, int16_t *output, uint32_t no_rows) { int16_t step1_0, step1_1, step1_2, step1_3, step1_4, step1_5, step1_6; int16_t step1_7, step1_8, step1_9, step1_10, step1_11, step1_12, step1_13; int16_t step1_14, step1_15, step1_16, step1_17, step1_18, step1_19, step1_20; int16_t step1_21, step1_22, step1_23, step1_24, step1_25, step1_26, step1_27; int16_t step1_28, step1_29, step1_30, step1_31; int16_t step2_0, step2_1, step2_2, step2_3, step2_4, step2_5, step2_6; int16_t step2_7, step2_8, step2_9, step2_10, step2_11, step2_12, step2_13; int16_t step2_14, step2_15, step2_16, step2_17, step2_18, step2_19, step2_20; int16_t step2_21, step2_22, step2_23, step2_24, step2_25, step2_26, step2_27; int16_t step2_28, step2_29, step2_30, step2_31; int16_t step3_8, step3_9, step3_10, step3_11, step3_12, step3_13, step3_14; int16_t step3_15, step3_16, step3_17, step3_18, step3_19, step3_20, step3_21; int16_t step3_22, step3_23, step3_24, step3_25, step3_26, step3_27, step3_28; int16_t step3_29, step3_30, step3_31; int temp0, temp1, temp2, temp3; int load1, load2, load3, load4; int result1, result2; int temp21; int i; const int const_2_power_13 = 8192; const int32_t *input_int; for (i = no_rows; i--; ) { input_int = (const int32_t *)input; if (!(input_int[0] | input_int[1] | input_int[2] | input_int[3] | input_int[4] | input_int[5] | input_int[6] | input_int[7] | input_int[8] | input_int[9] | input_int[10] | input_int[11] | input_int[12] | input_int[13] | input_int[14] | input_int[15])) { input += 32; __asm__ __volatile__ ( "sh $zero, 0(%[output]) \n\t" "sh $zero, 64(%[output]) \n\t" "sh $zero, 128(%[output]) \n\t" "sh $zero, 192(%[output]) \n\t" "sh $zero, 256(%[output]) \n\t" "sh $zero, 320(%[output]) \n\t" "sh $zero, 384(%[output]) \n\t" "sh $zero, 448(%[output]) \n\t" "sh $zero, 512(%[output]) \n\t" "sh $zero, 576(%[output]) \n\t" "sh $zero, 640(%[output]) \n\t" "sh $zero, 704(%[output]) \n\t" "sh $zero, 768(%[output]) \n\t" "sh $zero, 832(%[output]) \n\t" "sh $zero, 896(%[output]) \n\t" "sh $zero, 960(%[output]) \n\t" "sh $zero, 1024(%[output]) \n\t" "sh $zero, 1088(%[output]) \n\t" "sh $zero, 1152(%[output]) \n\t" "sh $zero, 1216(%[output]) \n\t" "sh $zero, 1280(%[output]) \n\t" "sh $zero, 1344(%[output]) \n\t" "sh $zero, 1408(%[output]) \n\t" "sh $zero, 1472(%[output]) \n\t" "sh $zero, 1536(%[output]) \n\t" "sh $zero, 1600(%[output]) \n\t" "sh $zero, 1664(%[output]) \n\t" "sh $zero, 1728(%[output]) \n\t" "sh $zero, 1792(%[output]) \n\t" "sh $zero, 1856(%[output]) \n\t" "sh $zero, 1920(%[output]) \n\t" "sh $zero, 1984(%[output]) \n\t" : : [output] "r" (output) ); output += 1; continue; } /* prefetch row */ prefetch_load((const uint8_t *)(input + 32)); prefetch_load((const uint8_t *)(input + 48)); __asm__ __volatile__ ( "lh %[load1], 2(%[input]) \n\t" "lh %[load2], 62(%[input]) \n\t" "lh %[load3], 34(%[input]) \n\t" "lh %[load4], 30(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_31_64] \n\t" "msub $ac1, %[load2], %[cospi_1_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_1_64] \n\t" "madd $ac3, %[load2], %[cospi_31_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_15_64] \n\t" "msub $ac2, %[load4], %[cospi_17_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_17_64] \n\t" "madd $ac1, %[load4], %[cospi_15_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp3], %[temp2] \n\t" "sub %[load2], %[temp0], %[temp1] \n\t" "madd $ac1, %[load1], %[cospi_28_64] \n\t" "msub $ac1, %[load2], %[cospi_4_64] \n\t" "madd $ac3, %[load1], %[cospi_4_64] \n\t" "madd $ac3, %[load2], %[cospi_28_64] \n\t" "extp %[step1_17], $ac1, 31 \n\t" "extp %[step1_30], $ac3, 31 \n\t" "add %[step1_16], %[temp0], %[temp1] \n\t" "add %[step1_31], %[temp2], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step1_16] "=r" (step1_16), [step1_17] "=r" (step1_17), [step1_30] "=r" (step1_30), [step1_31] "=r" (step1_31) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_31_64] "r" (cospi_31_64), [cospi_1_64] "r" (cospi_1_64), [cospi_4_64] "r" (cospi_4_64), [cospi_17_64] "r" (cospi_17_64), [cospi_15_64] "r" (cospi_15_64), [cospi_28_64] "r" (cospi_28_64) ); __asm__ __volatile__ ( "lh %[load1], 18(%[input]) \n\t" "lh %[load2], 46(%[input]) \n\t" "lh %[load3], 50(%[input]) \n\t" "lh %[load4], 14(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_23_64] \n\t" "msub $ac1, %[load2], %[cospi_9_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_9_64] \n\t" "madd $ac3, %[load2], %[cospi_23_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_7_64] \n\t" "msub $ac2, %[load4], %[cospi_25_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_25_64] \n\t" "madd $ac1, %[load4], %[cospi_7_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp1], %[temp0] \n\t" "sub %[load2], %[temp2], %[temp3] \n\t" "msub $ac1, %[load1], %[cospi_28_64] \n\t" "msub $ac1, %[load2], %[cospi_4_64] \n\t" "msub $ac3, %[load1], %[cospi_4_64] \n\t" "madd $ac3, %[load2], %[cospi_28_64] \n\t" "extp %[step1_18], $ac1, 31 \n\t" "extp %[step1_29], $ac3, 31 \n\t" "add %[step1_19], %[temp0], %[temp1] \n\t" "add %[step1_28], %[temp2], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step1_18] "=r" (step1_18), [step1_19] "=r" (step1_19), [step1_28] "=r" (step1_28), [step1_29] "=r" (step1_29) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_23_64] "r" (cospi_23_64), [cospi_9_64] "r" (cospi_9_64), [cospi_4_64] "r" (cospi_4_64), [cospi_7_64] "r" (cospi_7_64), [cospi_25_64] "r" (cospi_25_64), [cospi_28_64] "r" (cospi_28_64) ); __asm__ __volatile__ ( "lh %[load1], 10(%[input]) \n\t" "lh %[load2], 54(%[input]) \n\t" "lh %[load3], 42(%[input]) \n\t" "lh %[load4], 22(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_27_64] \n\t" "msub $ac1, %[load2], %[cospi_5_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_5_64] \n\t" "madd $ac3, %[load2], %[cospi_27_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_11_64] \n\t" "msub $ac2, %[load4], %[cospi_21_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_21_64] \n\t" "madd $ac1, %[load4], %[cospi_11_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp0], %[temp1] \n\t" "sub %[load2], %[temp3], %[temp2] \n\t" "madd $ac1, %[load2], %[cospi_12_64] \n\t" "msub $ac1, %[load1], %[cospi_20_64] \n\t" "madd $ac3, %[load1], %[cospi_12_64] \n\t" "madd $ac3, %[load2], %[cospi_20_64] \n\t" "extp %[step1_21], $ac1, 31 \n\t" "extp %[step1_26], $ac3, 31 \n\t" "add %[step1_20], %[temp0], %[temp1] \n\t" "add %[step1_27], %[temp2], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step1_20] "=r" (step1_20), [step1_21] "=r" (step1_21), [step1_26] "=r" (step1_26), [step1_27] "=r" (step1_27) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_27_64] "r" (cospi_27_64), [cospi_5_64] "r" (cospi_5_64), [cospi_11_64] "r" (cospi_11_64), [cospi_21_64] "r" (cospi_21_64), [cospi_12_64] "r" (cospi_12_64), [cospi_20_64] "r" (cospi_20_64) ); __asm__ __volatile__ ( "lh %[load1], 26(%[input]) \n\t" "lh %[load2], 38(%[input]) \n\t" "lh %[load3], 58(%[input]) \n\t" "lh %[load4], 6(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_19_64] \n\t" "msub $ac1, %[load2], %[cospi_13_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_13_64] \n\t" "madd $ac3, %[load2], %[cospi_19_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_3_64] \n\t" "msub $ac2, %[load4], %[cospi_29_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_29_64] \n\t" "madd $ac1, %[load4], %[cospi_3_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp1], %[temp0] \n\t" "sub %[load2], %[temp2], %[temp3] \n\t" "msub $ac1, %[load1], %[cospi_12_64] \n\t" "msub $ac1, %[load2], %[cospi_20_64] \n\t" "msub $ac3, %[load1], %[cospi_20_64] \n\t" "madd $ac3, %[load2], %[cospi_12_64] \n\t" "extp %[step1_22], $ac1, 31 \n\t" "extp %[step1_25], $ac3, 31 \n\t" "add %[step1_23], %[temp0], %[temp1] \n\t" "add %[step1_24], %[temp2], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step1_22] "=r" (step1_22), [step1_23] "=r" (step1_23), [step1_24] "=r" (step1_24), [step1_25] "=r" (step1_25) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_19_64] "r" (cospi_19_64), [cospi_13_64] "r" (cospi_13_64), [cospi_3_64] "r" (cospi_3_64), [cospi_29_64] "r" (cospi_29_64), [cospi_12_64] "r" (cospi_12_64), [cospi_20_64] "r" (cospi_20_64) ); __asm__ __volatile__ ( "lh %[load1], 4(%[input]) \n\t" "lh %[load2], 60(%[input]) \n\t" "lh %[load3], 36(%[input]) \n\t" "lh %[load4], 28(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_30_64] \n\t" "msub $ac1, %[load2], %[cospi_2_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_2_64] \n\t" "madd $ac3, %[load2], %[cospi_30_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_14_64] \n\t" "msub $ac2, %[load4], %[cospi_18_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_18_64] \n\t" "madd $ac1, %[load4], %[cospi_14_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp0], %[temp1] \n\t" "sub %[load2], %[temp3], %[temp2] \n\t" "msub $ac1, %[load1], %[cospi_8_64] \n\t" "madd $ac1, %[load2], %[cospi_24_64] \n\t" "madd $ac3, %[load1], %[cospi_24_64] \n\t" "madd $ac3, %[load2], %[cospi_8_64] \n\t" "extp %[step2_9], $ac1, 31 \n\t" "extp %[step2_14], $ac3, 31 \n\t" "add %[step2_8], %[temp0], %[temp1] \n\t" "add %[step2_15], %[temp2], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step2_8] "=r" (step2_8), [step2_9] "=r" (step2_9), [step2_14] "=r" (step2_14), [step2_15] "=r" (step2_15) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_30_64] "r" (cospi_30_64), [cospi_2_64] "r" (cospi_2_64), [cospi_14_64] "r" (cospi_14_64), [cospi_18_64] "r" (cospi_18_64), [cospi_8_64] "r" (cospi_8_64), [cospi_24_64] "r" (cospi_24_64) ); __asm__ __volatile__ ( "lh %[load1], 20(%[input]) \n\t" "lh %[load2], 44(%[input]) \n\t" "lh %[load3], 52(%[input]) \n\t" "lh %[load4], 12(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_22_64] \n\t" "msub $ac1, %[load2], %[cospi_10_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_10_64] \n\t" "madd $ac3, %[load2], %[cospi_22_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_6_64] \n\t" "msub $ac2, %[load4], %[cospi_26_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_26_64] \n\t" "madd $ac1, %[load4], %[cospi_6_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp1], %[temp0] \n\t" "sub %[load2], %[temp2], %[temp3] \n\t" "msub $ac1, %[load1], %[cospi_24_64] \n\t" "msub $ac1, %[load2], %[cospi_8_64] \n\t" "madd $ac3, %[load2], %[cospi_24_64] \n\t" "msub $ac3, %[load1], %[cospi_8_64] \n\t" "extp %[step2_10], $ac1, 31 \n\t" "extp %[step2_13], $ac3, 31 \n\t" "add %[step2_11], %[temp0], %[temp1] \n\t" "add %[step2_12], %[temp2], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step2_10] "=r" (step2_10), [step2_11] "=r" (step2_11), [step2_12] "=r" (step2_12), [step2_13] "=r" (step2_13) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_22_64] "r" (cospi_22_64), [cospi_10_64] "r" (cospi_10_64), [cospi_6_64] "r" (cospi_6_64), [cospi_26_64] "r" (cospi_26_64), [cospi_8_64] "r" (cospi_8_64), [cospi_24_64] "r" (cospi_24_64) ); __asm__ __volatile__ ( "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "sub %[temp0], %[step2_14], %[step2_13] \n\t" "sub %[temp0], %[temp0], %[step2_9] \n\t" "add %[temp0], %[temp0], %[step2_10] \n\t" "madd $ac0, %[temp0], %[cospi_16_64] \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "sub %[temp1], %[step2_14], %[step2_13] \n\t" "add %[temp1], %[temp1], %[step2_9] \n\t" "sub %[temp1], %[temp1], %[step2_10] \n\t" "madd $ac1, %[temp1], %[cospi_16_64] \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "sub %[temp0], %[step2_15], %[step2_12] \n\t" "sub %[temp0], %[temp0], %[step2_8] \n\t" "add %[temp0], %[temp0], %[step2_11] \n\t" "madd $ac2, %[temp0], %[cospi_16_64] \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[temp1], %[step2_15], %[step2_12] \n\t" "add %[temp1], %[temp1], %[step2_8] \n\t" "sub %[temp1], %[temp1], %[step2_11] \n\t" "madd $ac3, %[temp1], %[cospi_16_64] \n\t" "add %[step3_8], %[step2_8], %[step2_11] \n\t" "add %[step3_9], %[step2_9], %[step2_10] \n\t" "add %[step3_14], %[step2_13], %[step2_14] \n\t" "add %[step3_15], %[step2_12], %[step2_15] \n\t" "extp %[step3_10], $ac0, 31 \n\t" "extp %[step3_13], $ac1, 31 \n\t" "extp %[step3_11], $ac2, 31 \n\t" "extp %[step3_12], $ac3, 31 \n\t" : [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [step3_8] "=r" (step3_8), [step3_9] "=r" (step3_9), [step3_10] "=r" (step3_10), [step3_11] "=r" (step3_11), [step3_12] "=r" (step3_12), [step3_13] "=r" (step3_13), [step3_14] "=r" (step3_14), [step3_15] "=r" (step3_15) : [const_2_power_13] "r" (const_2_power_13), [step2_8] "r" (step2_8), [step2_9] "r" (step2_9), [step2_10] "r" (step2_10), [step2_11] "r" (step2_11), [step2_12] "r" (step2_12), [step2_13] "r" (step2_13), [step2_14] "r" (step2_14), [step2_15] "r" (step2_15), [cospi_16_64] "r" (cospi_16_64) ); step2_18 = step1_17 - step1_18; step2_29 = step1_30 - step1_29; __asm__ __volatile__ ( "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "msub $ac0, %[step2_18], %[cospi_8_64] \n\t" "madd $ac0, %[step2_29], %[cospi_24_64] \n\t" "extp %[step3_18], $ac0, 31 \n\t" : [step3_18] "=r" (step3_18) : [const_2_power_13] "r" (const_2_power_13), [step2_18] "r" (step2_18), [step2_29] "r" (step2_29), [cospi_24_64] "r" (cospi_24_64), [cospi_8_64] "r" (cospi_8_64) ); temp21 = step2_18 * cospi_24_64 + step2_29 * cospi_8_64; step3_29 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; step2_19 = step1_16 - step1_19; step2_28 = step1_31 - step1_28; __asm__ __volatile__ ( "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "msub $ac0, %[step2_19], %[cospi_8_64] \n\t" "madd $ac0, %[step2_28], %[cospi_24_64] \n\t" "extp %[step3_19], $ac0, 31 \n\t" : [step3_19] "=r" (step3_19) : [const_2_power_13] "r" (const_2_power_13), [step2_19] "r" (step2_19), [step2_28] "r" (step2_28), [cospi_24_64] "r" (cospi_24_64), [cospi_8_64] "r" (cospi_8_64) ); temp21 = step2_19 * cospi_24_64 + step2_28 * cospi_8_64; step3_28 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; step3_16 = step1_16 + step1_19; step3_17 = step1_17 + step1_18; step3_30 = step1_29 + step1_30; step3_31 = step1_28 + step1_31; step2_20 = step1_23 - step1_20; step2_27 = step1_24 - step1_27; __asm__ __volatile__ ( "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "msub $ac0, %[step2_20], %[cospi_24_64] \n\t" "msub $ac0, %[step2_27], %[cospi_8_64] \n\t" "extp %[step3_20], $ac0, 31 \n\t" : [step3_20] "=r" (step3_20) : [const_2_power_13] "r" (const_2_power_13), [step2_20] "r" (step2_20), [step2_27] "r" (step2_27), [cospi_24_64] "r" (cospi_24_64), [cospi_8_64] "r" (cospi_8_64) ); temp21 = -step2_20 * cospi_8_64 + step2_27 * cospi_24_64; step3_27 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; step2_21 = step1_22 - step1_21; step2_26 = step1_25 - step1_26; __asm__ __volatile__ ( "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "msub $ac1, %[step2_21], %[cospi_24_64] \n\t" "msub $ac1, %[step2_26], %[cospi_8_64] \n\t" "extp %[step3_21], $ac1, 31 \n\t" : [step3_21] "=r" (step3_21) : [const_2_power_13] "r" (const_2_power_13), [step2_21] "r" (step2_21), [step2_26] "r" (step2_26), [cospi_24_64] "r" (cospi_24_64), [cospi_8_64] "r" (cospi_8_64) ); temp21 = -step2_21 * cospi_8_64 + step2_26 * cospi_24_64; step3_26 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; step3_22 = step1_21 + step1_22; step3_23 = step1_20 + step1_23; step3_24 = step1_24 + step1_27; step3_25 = step1_25 + step1_26; step2_16 = step3_16 + step3_23; step2_17 = step3_17 + step3_22; step2_18 = step3_18 + step3_21; step2_19 = step3_19 + step3_20; step2_20 = step3_19 - step3_20; step2_21 = step3_18 - step3_21; step2_22 = step3_17 - step3_22; step2_23 = step3_16 - step3_23; step2_24 = step3_31 - step3_24; step2_25 = step3_30 - step3_25; step2_26 = step3_29 - step3_26; step2_27 = step3_28 - step3_27; step2_28 = step3_28 + step3_27; step2_29 = step3_29 + step3_26; step2_30 = step3_30 + step3_25; step2_31 = step3_31 + step3_24; __asm__ __volatile__ ( "lh %[load1], 0(%[input]) \n\t" "lh %[load2], 32(%[input]) \n\t" "lh %[load3], 16(%[input]) \n\t" "lh %[load4], 48(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "add %[result1], %[load1], %[load2] \n\t" "sub %[result2], %[load1], %[load2] \n\t" "madd $ac1, %[result1], %[cospi_16_64] \n\t" "madd $ac2, %[result2], %[cospi_16_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "extp %[temp1], $ac2, 31 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac3, %[load3], %[cospi_24_64] \n\t" "msub $ac3, %[load4], %[cospi_8_64] \n\t" "extp %[temp2], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "madd $ac1, %[load3], %[cospi_8_64] \n\t" "madd $ac1, %[load4], %[cospi_24_64] \n\t" "extp %[temp3], $ac1, 31 \n\t" "add %[step1_0], %[temp0], %[temp3] \n\t" "add %[step1_1], %[temp1], %[temp2] \n\t" "sub %[step1_2], %[temp1], %[temp2] \n\t" "sub %[step1_3], %[temp0], %[temp3] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [result1] "=&r" (result1), [result2] "=&r" (result2), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step1_0] "=r" (step1_0), [step1_1] "=r" (step1_1), [step1_2] "=r" (step1_2), [step1_3] "=r" (step1_3) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_16_64] "r" (cospi_16_64), [cospi_24_64] "r" (cospi_24_64), [cospi_8_64] "r" (cospi_8_64) ); __asm__ __volatile__ ( "lh %[load1], 8(%[input]) \n\t" "lh %[load2], 56(%[input]) \n\t" "lh %[load3], 40(%[input]) \n\t" "lh %[load4], 24(%[input]) \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "madd $ac1, %[load1], %[cospi_28_64] \n\t" "msub $ac1, %[load2], %[cospi_4_64] \n\t" "extp %[temp0], $ac1, 31 \n\t" "madd $ac3, %[load1], %[cospi_4_64] \n\t" "madd $ac3, %[load2], %[cospi_28_64] \n\t" "extp %[temp3], $ac3, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac2 \n\t" "mthi $zero, $ac2 \n\t" "madd $ac2, %[load3], %[cospi_12_64] \n\t" "msub $ac2, %[load4], %[cospi_20_64] \n\t" "extp %[temp1], $ac2, 31 \n\t" "madd $ac1, %[load3], %[cospi_20_64] \n\t" "madd $ac1, %[load4], %[cospi_12_64] \n\t" "extp %[temp2], $ac1, 31 \n\t" "mtlo %[const_2_power_13], $ac1 \n\t" "mthi $zero, $ac1 \n\t" "mtlo %[const_2_power_13], $ac3 \n\t" "mthi $zero, $ac3 \n\t" "sub %[load1], %[temp3], %[temp2] \n\t" "sub %[load1], %[load1], %[temp0] \n\t" "add %[load1], %[load1], %[temp1] \n\t" "sub %[load2], %[temp0], %[temp1] \n\t" "sub %[load2], %[load2], %[temp2] \n\t" "add %[load2], %[load2], %[temp3] \n\t" "madd $ac1, %[load1], %[cospi_16_64] \n\t" "madd $ac3, %[load2], %[cospi_16_64] \n\t" "extp %[step1_5], $ac1, 31 \n\t" "extp %[step1_6], $ac3, 31 \n\t" "add %[step1_4], %[temp0], %[temp1] \n\t" "add %[step1_7], %[temp3], %[temp2] \n\t" : [load1] "=&r" (load1), [load2] "=&r" (load2), [load3] "=&r" (load3), [load4] "=&r" (load4), [temp0] "=&r" (temp0), [temp1] "=&r" (temp1), [temp2] "=&r" (temp2), [temp3] "=&r" (temp3), [step1_4] "=r" (step1_4), [step1_5] "=r" (step1_5), [step1_6] "=r" (step1_6), [step1_7] "=r" (step1_7) : [const_2_power_13] "r" (const_2_power_13), [input] "r" (input), [cospi_20_64] "r" (cospi_20_64), [cospi_12_64] "r" (cospi_12_64), [cospi_4_64] "r" (cospi_4_64), [cospi_28_64] "r" (cospi_28_64), [cospi_16_64] "r" (cospi_16_64) ); step2_0 = step1_0 + step1_7; step2_1 = step1_1 + step1_6; step2_2 = step1_2 + step1_5; step2_3 = step1_3 + step1_4; step2_4 = step1_3 - step1_4; step2_5 = step1_2 - step1_5; step2_6 = step1_1 - step1_6; step2_7 = step1_0 - step1_7; step1_0 = step2_0 + step3_15; step1_1 = step2_1 + step3_14; step1_2 = step2_2 + step3_13; step1_3 = step2_3 + step3_12; step1_4 = step2_4 + step3_11; step1_5 = step2_5 + step3_10; step1_6 = step2_6 + step3_9; step1_7 = step2_7 + step3_8; step1_8 = step2_7 - step3_8; step1_9 = step2_6 - step3_9; step1_10 = step2_5 - step3_10; step1_11 = step2_4 - step3_11; step1_12 = step2_3 - step3_12; step1_13 = step2_2 - step3_13; step1_14 = step2_1 - step3_14; step1_15 = step2_0 - step3_15; __asm__ __volatile__ ( "sub %[temp0], %[step2_27], %[step2_20] \n\t" "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "madd $ac0, %[temp0], %[cospi_16_64] \n\t" "extp %[step1_20], $ac0, 31 \n\t" : [temp0] "=&r" (temp0), [step1_20] "=r" (step1_20) : [const_2_power_13] "r" (const_2_power_13), [step2_20] "r" (step2_20), [step2_27] "r" (step2_27), [cospi_16_64] "r" (cospi_16_64) ); temp21 = (step2_20 + step2_27) * cospi_16_64; step1_27 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; __asm__ __volatile__ ( "sub %[temp0], %[step2_26], %[step2_21] \n\t" "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "madd $ac0, %[temp0], %[cospi_16_64] \n\t" "extp %[step1_21], $ac0, 31 \n\t" : [temp0] "=&r" (temp0), [step1_21] "=r" (step1_21) : [const_2_power_13] "r" (const_2_power_13), [step2_26] "r" (step2_26), [step2_21] "r" (step2_21), [cospi_16_64] "r" (cospi_16_64) ); temp21 = (step2_21 + step2_26) * cospi_16_64; step1_26 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; __asm__ __volatile__ ( "sub %[temp0], %[step2_25], %[step2_22] \n\t" "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "madd $ac0, %[temp0], %[cospi_16_64] \n\t" "extp %[step1_22], $ac0, 31 \n\t" : [temp0] "=&r" (temp0), [step1_22] "=r" (step1_22) : [const_2_power_13] "r" (const_2_power_13), [step2_25] "r" (step2_25), [step2_22] "r" (step2_22), [cospi_16_64] "r" (cospi_16_64) ); temp21 = (step2_22 + step2_25) * cospi_16_64; step1_25 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; __asm__ __volatile__ ( "sub %[temp0], %[step2_24], %[step2_23] \n\t" "mtlo %[const_2_power_13], $ac0 \n\t" "mthi $zero, $ac0 \n\t" "madd $ac0, %[temp0], %[cospi_16_64] \n\t" "extp %[step1_23], $ac0, 31 \n\t" : [temp0] "=&r" (temp0), [step1_23] "=r" (step1_23) : [const_2_power_13] "r" (const_2_power_13), [step2_24] "r" (step2_24), [step2_23] "r" (step2_23), [cospi_16_64] "r" (cospi_16_64) ); temp21 = (step2_23 + step2_24) * cospi_16_64; step1_24 = (temp21 + DCT_CONST_ROUNDING) >> DCT_CONST_BITS; // final stage output[0 * 32] = step1_0 + step2_31; output[1 * 32] = step1_1 + step2_30; output[2 * 32] = step1_2 + step2_29; output[3 * 32] = step1_3 + step2_28; output[4 * 32] = step1_4 + step1_27; output[5 * 32] = step1_5 + step1_26; output[6 * 32] = step1_6 + step1_25; output[7 * 32] = step1_7 + step1_24; output[8 * 32] = step1_8 + step1_23; output[9 * 32] = step1_9 + step1_22; output[10 * 32] = step1_10 + step1_21; output[11 * 32] = step1_11 + step1_20; output[12 * 32] = step1_12 + step2_19; output[13 * 32] = step1_13 + step2_18; output[14 * 32] = step1_14 + step2_17; output[15 * 32] = step1_15 + step2_16; output[16 * 32] = step1_15 - step2_16; output[17 * 32] = step1_14 - step2_17; output[18 * 32] = step1_13 - step2_18; output[19 * 32] = step1_12 - step2_19; output[20 * 32] = step1_11 - step1_20; output[21 * 32] = step1_10 - step1_21; output[22 * 32] = step1_9 - step1_22; output[23 * 32] = step1_8 - step1_23; output[24 * 32] = step1_7 - step1_24; output[25 * 32] = step1_6 - step1_25; output[26 * 32] = step1_5 - step1_26; output[27 * 32] = step1_4 - step1_27; output[28 * 32] = step1_3 - step2_28; output[29 * 32] = step1_2 - step2_29; output[30 * 32] = step1_1 - step2_30; output[31 * 32] = step1_0 - step2_31; input += 32; output += 1; } } void vpx_idct32x32_1024_add_dspr2(const int16_t *input, uint8_t *dest, int dest_stride) { DECLARE_ALIGNED(32, int16_t, out[32 * 32]); 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 idct32_rows_dspr2(input, outptr, 32); // Columns vpx_idct32_cols_add_blk_dspr2(out, dest, dest_stride); } void vpx_idct32x32_34_add_dspr2(const int16_t *input, uint8_t *dest, int stride) { DECLARE_ALIGNED(32, int16_t, out[32 * 32]); int16_t *outptr = out; uint32_t i; uint32_t pos = 45; /* bit positon for extract from acc */ __asm__ __volatile__ ( "wrdsp %[pos], 1 \n\t" : : [pos] "r" (pos) ); // Rows idct32_rows_dspr2(input, outptr, 8); outptr += 8; __asm__ __volatile__ ( "sw $zero, 0(%[outptr]) \n\t" "sw $zero, 4(%[outptr]) \n\t" "sw $zero, 8(%[outptr]) \n\t" "sw $zero, 12(%[outptr]) \n\t" "sw $zero, 16(%[outptr]) \n\t" "sw $zero, 20(%[outptr]) \n\t" "sw $zero, 24(%[outptr]) \n\t" "sw $zero, 28(%[outptr]) \n\t" "sw $zero, 32(%[outptr]) \n\t" "sw $zero, 36(%[outptr]) \n\t" "sw $zero, 40(%[outptr]) \n\t" "sw $zero, 44(%[outptr]) \n\t" : : [outptr] "r" (outptr) ); for (i = 0; i < 31; ++i) { outptr += 32; __asm__ __volatile__ ( "sw $zero, 0(%[outptr]) \n\t" "sw $zero, 4(%[outptr]) \n\t" "sw $zero, 8(%[outptr]) \n\t" "sw $zero, 12(%[outptr]) \n\t" "sw $zero, 16(%[outptr]) \n\t" "sw $zero, 20(%[outptr]) \n\t" "sw $zero, 24(%[outptr]) \n\t" "sw $zero, 28(%[outptr]) \n\t" "sw $zero, 32(%[outptr]) \n\t" "sw $zero, 36(%[outptr]) \n\t" "sw $zero, 40(%[outptr]) \n\t" "sw $zero, 44(%[outptr]) \n\t" : : [outptr] "r" (outptr) ); } // Columns vpx_idct32_cols_add_blk_dspr2(out, dest, stride); } void vpx_idct32x32_1_add_dspr2(const int16_t *input, uint8_t *dest, int stride) { int r, out; int32_t a1, absa1; int32_t vector_a1; int32_t t1, t2, t3, t4; int32_t vector_1, vector_2, vector_3, vector_4; uint32_t pos = 45; /* 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[0]); __asm__ __volatile__ ( "addi %[out], %[out], 32 \n\t" "sra %[a1], %[out], 6 \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 = 32; r--;) { __asm__ __volatile__ ( "lw %[t1], 0(%[dest]) \n\t" "lw %[t2], 4(%[dest]) \n\t" "lw %[t3], 8(%[dest]) \n\t" "lw %[t4], 12(%[dest]) \n\t" "subu_s.qb %[vector_1], %[t1], %[vector_a1] \n\t" "subu_s.qb %[vector_2], %[t2], %[vector_a1] \n\t" "subu_s.qb %[vector_3], %[t3], %[vector_a1] \n\t" "subu_s.qb %[vector_4], %[t4], %[vector_a1] \n\t" "sw %[vector_1], 0(%[dest]) \n\t" "sw %[vector_2], 4(%[dest]) \n\t" "sw %[vector_3], 8(%[dest]) \n\t" "sw %[vector_4], 12(%[dest]) \n\t" "lw %[t1], 16(%[dest]) \n\t" "lw %[t2], 20(%[dest]) \n\t" "lw %[t3], 24(%[dest]) \n\t" "lw %[t4], 28(%[dest]) \n\t" "subu_s.qb %[vector_1], %[t1], %[vector_a1] \n\t" "subu_s.qb %[vector_2], %[t2], %[vector_a1] \n\t" "subu_s.qb %[vector_3], %[t3], %[vector_a1] \n\t" "subu_s.qb %[vector_4], %[t4], %[vector_a1] \n\t" "sw %[vector_1], 16(%[dest]) \n\t" "sw %[vector_2], 20(%[dest]) \n\t" "sw %[vector_3], 24(%[dest]) \n\t" "sw %[vector_4], 28(%[dest]) \n\t" "add %[dest], %[dest], %[stride] \n\t" : [t1] "=&r" (t1), [t2] "=&r" (t2), [t3] "=&r" (t3), [t4] "=&r" (t4), [vector_1] "=&r" (vector_1), [vector_2] "=&r" (vector_2), [vector_3] "=&r" (vector_3), [vector_4] "=&r" (vector_4), [dest] "+&r" (dest) : [stride] "r" (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 = 32; r--;) { __asm__ __volatile__ ( "lw %[t1], 0(%[dest]) \n\t" "lw %[t2], 4(%[dest]) \n\t" "lw %[t3], 8(%[dest]) \n\t" "lw %[t4], 12(%[dest]) \n\t" "addu_s.qb %[vector_1], %[t1], %[vector_a1] \n\t" "addu_s.qb %[vector_2], %[t2], %[vector_a1] \n\t" "addu_s.qb %[vector_3], %[t3], %[vector_a1] \n\t" "addu_s.qb %[vector_4], %[t4], %[vector_a1] \n\t" "sw %[vector_1], 0(%[dest]) \n\t" "sw %[vector_2], 4(%[dest]) \n\t" "sw %[vector_3], 8(%[dest]) \n\t" "sw %[vector_4], 12(%[dest]) \n\t" "lw %[t1], 16(%[dest]) \n\t" "lw %[t2], 20(%[dest]) \n\t" "lw %[t3], 24(%[dest]) \n\t" "lw %[t4], 28(%[dest]) \n\t" "addu_s.qb %[vector_1], %[t1], %[vector_a1] \n\t" "addu_s.qb %[vector_2], %[t2], %[vector_a1] \n\t" "addu_s.qb %[vector_3], %[t3], %[vector_a1] \n\t" "addu_s.qb %[vector_4], %[t4], %[vector_a1] \n\t" "sw %[vector_1], 16(%[dest]) \n\t" "sw %[vector_2], 20(%[dest]) \n\t" "sw %[vector_3], 24(%[dest]) \n\t" "sw %[vector_4], 28(%[dest]) \n\t" "add %[dest], %[dest], %[stride] \n\t" : [t1] "=&r" (t1), [t2] "=&r" (t2), [t3] "=&r" (t3), [t4] "=&r" (t4), [vector_1] "=&r" (vector_1), [vector_2] "=&r" (vector_2), [vector_3] "=&r" (vector_3), [vector_4] "=&r" (vector_4), [dest] "+&r" (dest) : [stride] "r" (stride), [vector_a1] "r" (vector_a1) ); } } } #endif // #if HAVE_DSPR2