cnstream/easydk/3rdparty/libjpeg-turbo/simd/mips/jsimd_dspr2.S

/*
 * MIPS DSPr2 optimizations for libjpeg-turbo
 *
 * Copyright (C) 2013-2014, MIPS Technologies, Inc., California.
 *                          All Rights Reserved.
 * Authors:  Teodora Novkovic <teodora.novkovic@imgtec.com>
 *           Darko Laus       <darko.laus@imgtec.com>
 * Copyright (C) 2015, D. R. Commander.  All Rights Reserved.
 *
 * This software is provided 'as-is', without any express or implied
 * warranty.  In no event will the authors be held liable for any damages
 * arising from the use of this software.
 *
 * Permission is granted to anyone to use this software for any purpose,
 * including commercial applications, and to alter it and redistribute it
 * freely, subject to the following restrictions:
 *
 * 1. The origin of this software must not be misrepresented; you must not
 *    claim that you wrote the original software. If you use this software
 *    in a product, an acknowledgment in the product documentation would be
 *    appreciated but is not required.
 * 2. Altered source versions must be plainly marked as such, and must not be
 *    misrepresented as being the original software.
 * 3. This notice may not be removed or altered from any source distribution.
 */

#include "jsimd_dspr2_asm.h"


/*****************************************************************************/
LEAF_DSPR2(jsimd_c_null_convert_dspr2)
/*
 * a0     = cinfo->image_width
 * a1     = input_buf
 * a2     = output_buf
 * a3     = output_row
 * 16(sp) = num_rows
 * 20(sp) = cinfo->num_components
 *
 * Null conversion for compression
 */
    SAVE_REGS_ON_STACK 8, s0, s1

    lw          t9, 24(sp)      // t9 = num_rows
    lw          s0, 28(sp)      // s0 = cinfo->num_components
    andi        t0, a0, 3       // t0 = cinfo->image_width & 3
    beqz        t0, 4f          // no residual
     nop
0:
    addiu       t9, t9, -1
    bltz        t9, 7f
     li         t1, 0
1:
    sll         t3, t1, 2
    lwx         t5, t3(a2)      // t5 = outptr = output_buf[ci]
    lw          t2, 0(a1)       // t2 = inptr = *input_buf
    sll         t4, a3, 2
    lwx         t5, t4(t5)      // t5 = outptr = output_buf[ci][output_row]
    addu        t2, t2, t1
    addu        s1, t5, a0
    addu        t6, t5, t0
2:
    lbu         t3, 0(t2)
    addiu       t5, t5, 1
    sb          t3, -1(t5)
    bne         t6, t5, 2b
     addu       t2, t2, s0
3:
    lbu         t3, 0(t2)
    addu        t4, t2, s0
    addu        t7, t4, s0
    addu        t8, t7, s0
    addu        t2, t8, s0
    lbu         t4, 0(t4)
    lbu         t7, 0(t7)
    lbu         t8, 0(t8)
    addiu       t5, t5, 4
    sb          t3, -4(t5)
    sb          t4, -3(t5)
    sb          t7, -2(t5)
    bne         s1, t5, 3b
     sb         t8, -1(t5)
    addiu       t1, t1, 1
    bne         t1, s0, 1b
     nop
    addiu       a1, a1, 4
    bgez        t9, 0b
     addiu      a3, a3, 1
    b           7f
     nop
4:
    addiu       t9, t9, -1
    bltz        t9, 7f
     li         t1, 0
5:
    sll         t3, t1, 2
    lwx         t5, t3(a2)      // t5 = outptr = output_buf[ci]
    lw          t2, 0(a1)       // t2 = inptr = *input_buf
    sll         t4, a3, 2
    lwx         t5, t4(t5)      // t5 = outptr = output_buf[ci][output_row]
    addu        t2, t2, t1
    addu        s1, t5, a0
    addu        t6, t5, t0
6:
    lbu         t3, 0(t2)
    addu        t4, t2, s0
    addu        t7, t4, s0
    addu        t8, t7, s0
    addu        t2, t8, s0
    lbu         t4, 0(t4)
    lbu         t7, 0(t7)
    lbu         t8, 0(t8)
    addiu       t5, t5, 4
    sb          t3, -4(t5)
    sb          t4, -3(t5)
    sb          t7, -2(t5)
    bne         s1, t5, 6b
     sb         t8, -1(t5)
    addiu       t1, t1, 1
    bne         t1, s0, 5b
     nop
    addiu       a1, a1, 4
    bgez        t9, 4b
     addiu      a3, a3, 1
7:
    RESTORE_REGS_FROM_STACK 8, s0, s1

    j           ra
     nop

END(jsimd_c_null_convert_dspr2)


/*****************************************************************************/
/*
 * jsimd_extrgb_ycc_convert_dspr2
 * jsimd_extbgr_ycc_convert_dspr2
 * jsimd_extrgbx_ycc_convert_dspr2
 * jsimd_extbgrx_ycc_convert_dspr2
 * jsimd_extxbgr_ycc_convert_dspr2
 * jsimd_extxrgb_ycc_convert_dspr2
 *
 * Colorspace conversion RGB -> YCbCr
 */

.macro GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2  colorid, pixel_size, \
                                             r_offs, g_offs, b_offs

.macro DO_RGB_TO_YCC  r, g, b, inptr
    lbu         \r, \r_offs(\inptr)
    lbu         \g, \g_offs(\inptr)
    lbu         \b, \b_offs(\inptr)
    addiu       \inptr, \pixel_size
.endm

LEAF_DSPR2(jsimd_\colorid\()_ycc_convert_dspr2)
/*
 * a0     = cinfo->image_width
 * a1     = input_buf
 * a2     = output_buf
 * a3     = output_row
 * 16(sp) = num_rows
 */
    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    lw          t7, 48(sp)      // t7 = num_rows
    li          s0, 0x4c8b      // FIX(0.29900)
    li          s1, 0x9646      // FIX(0.58700)
    li          s2, 0x1d2f      // FIX(0.11400)
    li          s3, 0xffffd4cd  // -FIX(0.16874)
    li          s4, 0xffffab33  // -FIX(0.33126)
    li          s5, 0x8000      // FIX(0.50000)
    li          s6, 0xffff94d1  // -FIX(0.41869)
    li          s7, 0xffffeb2f  // -FIX(0.08131)
    li          t8, 0x807fff    // CBCR_OFFSET + ONE_HALF-1

0:
    addiu       t7, -1          // --num_rows
    lw          t6, 0(a1)       // t6 = input_buf[0]
    lw          t0, 0(a2)
    lw          t1, 4(a2)
    lw          t2, 8(a2)
    sll         t3, a3, 2
    lwx         t0, t3(t0)      // t0 = output_buf[0][output_row]
    lwx         t1, t3(t1)      // t1 = output_buf[1][output_row]
    lwx         t2, t3(t2)      // t2 = output_buf[2][output_row]

    addu        t9, t2, a0      // t9 = end address
    addiu       a3, 1

1:
    DO_RGB_TO_YCC t3, t4, t5, t6

    mtlo        s5, $ac0
    mtlo        t8, $ac1
    mtlo        t8, $ac2
    maddu       $ac0, s2, t5
    maddu       $ac1, s5, t5
    maddu       $ac2, s5, t3
    maddu       $ac0, s0, t3
    maddu       $ac1, s3, t3
    maddu       $ac2, s6, t4
    maddu       $ac0, s1, t4
    maddu       $ac1, s4, t4
    maddu       $ac2, s7, t5
    extr.w      t3, $ac0, 16
    extr.w      t4, $ac1, 16
    extr.w      t5, $ac2, 16
    sb          t3, 0(t0)
    sb          t4, 0(t1)
    sb          t5, 0(t2)
    addiu       t0, 1
    addiu       t2, 1
    bne         t2, t9, 1b
     addiu      t1, 1
    bgtz        t7, 0b
     addiu      a1, 4

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop
END(jsimd_\colorid\()_ycc_convert_dspr2)

.purgem DO_RGB_TO_YCC

.endm

/*-------------------------------------id -- pix R  G  B */
GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2 extrgb,  3, 0, 1, 2
GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2 extbgr,  3, 2, 1, 0
GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2 extrgbx, 4, 0, 1, 2
GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2 extbgrx, 4, 2, 1, 0
GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2 extxbgr, 4, 3, 2, 1
GENERATE_JSIMD_RGB_YCC_CONVERT_DSPR2 extxrgb, 4, 1, 2, 3


/*****************************************************************************/
/*
 * jsimd_ycc_extrgb_convert_dspr2
 * jsimd_ycc_extbgr_convert_dspr2
 * jsimd_ycc_extrgbx_convert_dspr2
 * jsimd_ycc_extbgrx_convert_dspr2
 * jsimd_ycc_extxbgr_convert_dspr2
 * jsimd_ycc_extxrgb_convert_dspr2
 *
 * Colorspace conversion YCbCr -> RGB
 */

.macro GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2  colorid, pixel_size, \
                                             r_offs, g_offs, b_offs, a_offs

.macro STORE_YCC_TO_RGB  scratch0 scratch1 scratch2 outptr
    sb          \scratch0, \r_offs(\outptr)
    sb          \scratch1, \g_offs(\outptr)
    sb          \scratch2, \b_offs(\outptr)
.if (\pixel_size == 4)
    li          t0, 0xFF
    sb          t0, \a_offs(\outptr)
.endif
    addiu       \outptr, \pixel_size
.endm

LEAF_DSPR2(jsimd_ycc_\colorid\()_convert_dspr2)
/*
 * a0     = cinfo->image_width
 * a1     = input_buf
 * a2     = input_row
 * a3     = output_buf
 * 16(sp) = num_rows
 */
    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    lw          s1, 48(sp)
    li          t3, 0x8000
    li          t4, 0x166e9     // FIX(1.40200)
    li          t5, 0x1c5a2     // FIX(1.77200)
    li          t6, 0xffff492e  // -FIX(0.71414)
    li          t7, 0xffffa7e6  // -FIX(0.34414)
    repl.ph     t8, 128

0:
    lw          s0, 0(a3)
    lw          t0, 0(a1)
    lw          t1, 4(a1)
    lw          t2, 8(a1)
    sll         s5, a2, 2
    addiu       s1, -1
    lwx         s2, s5(t0)
    lwx         s3, s5(t1)
    lwx         s4, s5(t2)
    addu        t9, s2, a0
    addiu       a2, 1

1:
    lbu         s7, 0(s4)       // cr
    lbu         s6, 0(s3)       // cb
    lbu         s5, 0(s2)       // y
    addiu       s2, 1
    addiu       s4, 1
    addiu       s7, -128
    addiu       s6, -128
    mul         t2, t7, s6
    mul         t0, t6, s7      // Crgtab[cr]
    sll         s7, 15
    mulq_rs.w   t1, t4, s7      // Crrtab[cr]
    sll         s6, 15
    addu        t2, t3          // Cbgtab[cb]
    addu        t2, t0

    mulq_rs.w   t0, t5, s6      // Cbbtab[cb]
    sra         t2, 16
    addu        t1, s5
    addu        t2, s5          // add y
    ins         t2, t1, 16, 16
    subu.ph     t2, t2, t8
    addu        t0, s5
    shll_s.ph   t2, t2, 8
    subu        t0, 128
    shra.ph     t2, t2, 8
    shll_s.w    t0, t0, 24
    addu.ph     t2, t2, t8      // clip & store
    sra         t0, t0, 24
    sra         t1, t2, 16
    addiu       t0, 128

    STORE_YCC_TO_RGB t1, t2, t0, s0

    bne         s2, t9, 1b
     addiu      s3, 1
    bgtz        s1, 0b
     addiu      a3, 4

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop
END(jsimd_ycc_\colorid\()_convert_dspr2)

.purgem STORE_YCC_TO_RGB

.endm

/*-------------------------------------id -- pix R  G  B  A */
GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2 extrgb,  3, 0, 1, 2, 3
GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2 extbgr,  3, 2, 1, 0, 3
GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2 extrgbx, 4, 0, 1, 2, 3
GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2 extbgrx, 4, 2, 1, 0, 3
GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2 extxbgr, 4, 3, 2, 1, 0
GENERATE_JSIMD_YCC_RGB_CONVERT_DSPR2 extxrgb, 4, 1, 2, 3, 0


/*****************************************************************************/
/*
 * jsimd_extrgb_gray_convert_dspr2
 * jsimd_extbgr_gray_convert_dspr2
 * jsimd_extrgbx_gray_convert_dspr2
 * jsimd_extbgrx_gray_convert_dspr2
 * jsimd_extxbgr_gray_convert_dspr2
 * jsimd_extxrgb_gray_convert_dspr2
 *
 * Colorspace conversion RGB -> GRAY
 */

.macro GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2  colorid, pixel_size, \
                                              r_offs, g_offs, b_offs

.macro DO_RGB_TO_GRAY  r, g, b, inptr
    lbu         \r, \r_offs(\inptr)
    lbu         \g, \g_offs(\inptr)
    lbu         \b, \b_offs(\inptr)
    addiu       \inptr, \pixel_size
.endm

LEAF_DSPR2(jsimd_\colorid\()_gray_convert_dspr2)
/*
 * a0     = cinfo->image_width
 * a1     = input_buf
 * a2     = output_buf
 * a3     = output_row
 * 16(sp) = num_rows
 */
    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    li          s0, 0x4c8b      // s0 = FIX(0.29900)
    li          s1, 0x9646      // s1 = FIX(0.58700)
    li          s2, 0x1d2f      // s2 = FIX(0.11400)
    li          s7, 0x8000      // s7 = FIX(0.50000)
    lw          s6, 48(sp)
    andi        t7, a0, 3

0:
    addiu       s6, -1          // s6 = num_rows
    lw          t0, 0(a1)
    lw          t1, 0(a2)
    sll         t3, a3, 2
    lwx         t1, t3(t1)
    addiu       a3, 1
    addu        t9, t1, a0
    subu        t8, t9, t7
    beq         t1, t8, 2f
     nop

1:
    DO_RGB_TO_GRAY t3, t4, t5, t0
    DO_RGB_TO_GRAY s3, s4, s5, t0

    mtlo        s7, $ac0
    maddu       $ac0, s2, t5
    maddu       $ac0, s1, t4
    maddu       $ac0, s0, t3
    mtlo        s7, $ac1
    maddu       $ac1, s2, s5
    maddu       $ac1, s1, s4
    maddu       $ac1, s0, s3
    extr.w      t6, $ac0, 16

    DO_RGB_TO_GRAY t3, t4, t5, t0
    DO_RGB_TO_GRAY s3, s4, s5, t0

    mtlo        s7, $ac0
    maddu       $ac0, s2, t5
    maddu       $ac0, s1, t4
    extr.w      t2, $ac1, 16
    maddu       $ac0, s0, t3
    mtlo        s7, $ac1
    maddu       $ac1, s2, s5
    maddu       $ac1, s1, s4
    maddu       $ac1, s0, s3
    extr.w      t5, $ac0, 16
    sb          t6, 0(t1)
    sb          t2, 1(t1)
    extr.w      t3, $ac1, 16
    addiu       t1, 4
    sb          t5, -2(t1)
    sb          t3, -1(t1)
    bne         t1, t8, 1b
     nop

2:
    beqz        t7, 4f
     nop

3:
    DO_RGB_TO_GRAY t3, t4, t5, t0

    mtlo        s7, $ac0
    maddu       $ac0, s2, t5
    maddu       $ac0, s1, t4
    maddu       $ac0, s0, t3
    extr.w      t6, $ac0, 16
    sb          t6, 0(t1)
    addiu       t1, 1
    bne         t1, t9, 3b
     nop

4:
    bgtz        s6, 0b
     addiu      a1, 4

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop
END(jsimd_\colorid\()_gray_convert_dspr2)

.purgem DO_RGB_TO_GRAY

.endm

/*-------------------------------------id --  pix R  G  B */
GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2 extrgb,  3, 0, 1, 2
GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2 extbgr,  3, 2, 1, 0
GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2 extrgbx, 4, 0, 1, 2
GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2 extbgrx, 4, 2, 1, 0
GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2 extxbgr, 4, 3, 2, 1
GENERATE_JSIMD_RGB_GRAY_CONVERT_DSPR2 extxrgb, 4, 1, 2, 3


/*****************************************************************************/
/*
 * jsimd_h2v2_merged_upsample_dspr2
 * jsimd_h2v2_extrgb_merged_upsample_dspr2
 * jsimd_h2v2_extrgbx_merged_upsample_dspr2
 * jsimd_h2v2_extbgr_merged_upsample_dspr2
 * jsimd_h2v2_extbgrx_merged_upsample_dspr2
 * jsimd_h2v2_extxbgr_merged_upsample_dspr2
 * jsimd_h2v2_extxrgb_merged_upsample_dspr2
 *
 * Merged h2v2 upsample routines
 */
.macro GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2  colorid, pixel_size, \
                                            r1_offs, g1_offs, \
                                            b1_offs, a1_offs, \
                                            r2_offs, g2_offs, \
                                            b2_offs, a2_offs

.macro STORE_H2V2_2_PIXELS  scratch0 scratch1 scratch2 scratch3 scratch4 \
                            scratch5 outptr
    sb          \scratch0, \r1_offs(\outptr)
    sb          \scratch1, \g1_offs(\outptr)
    sb          \scratch2, \b1_offs(\outptr)
    sb          \scratch3, \r2_offs(\outptr)
    sb          \scratch4, \g2_offs(\outptr)
    sb          \scratch5, \b2_offs(\outptr)
.if (\pixel_size == 8)
    li          \scratch0, 0xFF
    sb          \scratch0, \a1_offs(\outptr)
    sb          \scratch0, \a2_offs(\outptr)
.endif
    addiu       \outptr, \pixel_size
.endm

.macro STORE_H2V2_1_PIXEL  scratch0 scratch1 scratch2 outptr
    sb          \scratch0, \r1_offs(\outptr)
    sb          \scratch1, \g1_offs(\outptr)
    sb          \scratch2, \b1_offs(\outptr)

.if (\pixel_size == 8)
    li          t0, 0xFF
    sb          t0, \a1_offs(\outptr)
.endif
.endm

LEAF_DSPR2(jsimd_h2v2_\colorid\()_merged_upsample_dspr2)
/*
 * a0     = cinfo->output_width
 * a1     = input_buf
 * a2     = in_row_group_ctr
 * a3     = output_buf
 * 16(sp) = cinfo->sample_range_limit
 */
    SAVE_REGS_ON_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, ra

    lw          t9, 56(sp)      // cinfo->sample_range_limit
    lw          v0, 0(a1)
    lw          v1, 4(a1)
    lw          t0, 8(a1)
    sll         t1, a2, 3
    addiu       t2, t1, 4
    sll         t3, a2, 2
    lw          t4, 0(a3)       // t4 = output_buf[0]
    lwx         t1, t1(v0)      // t1 = input_buf[0][in_row_group_ctr*2]
    lwx         t2, t2(v0)      // t2 = input_buf[0][in_row_group_ctr*2 + 1]
    lwx         t5, t3(v1)      // t5 = input_buf[1][in_row_group_ctr]
    lwx         t6, t3(t0)      // t6 = input_buf[2][in_row_group_ctr]
    lw          t7, 4(a3)       // t7 = output_buf[1]
    li          s1, 0xe6ea
    addiu       t8, s1, 0x7fff    // t8 = 0x166e9 [FIX(1.40200)]
    addiu       s0, t8, 0x5eb9    // s0 = 0x1c5a2 [FIX(1.77200)]
    addiu       s1, zero, 0xa7e6  // s4 = 0xffffa7e6 [-FIX(0.34414)]
    xori        s2, s1, 0xeec8    // s3 = 0xffff492e [-FIX(0.71414)]
    srl         t3, a0, 1
    blez        t3, 2f
     addu       t0, t5, t3      // t0 = end address
 1:
    lbu         t3, 0(t5)
    lbu         s3, 0(t6)
    addiu       t5, t5, 1
    addiu       t3, t3, -128    // (cb - 128)
    addiu       s3, s3, -128    // (cr - 128)
    mult        $ac1, s1, t3
    madd        $ac1, s2, s3
    sll         s3, s3, 15
    sll         t3, t3, 15
    mulq_rs.w   s4, t8, s3      // s4 = (C1 * cr + ONE_HALF)>> SCALEBITS
    extr_r.w    s5, $ac1, 16
    mulq_rs.w   s6, s0, t3      // s6 = (C2 * cb + ONE_HALF)>> SCALEBITS
    lbu         v0, 0(t1)
    addiu       t6, t6, 1
    addiu       t1, t1, 2
    addu        t3, v0, s4      // y+cred
    addu        s3, v0, s5      // y+cgreen
    addu        v1, v0, s6      // y+cblue
    addu        t3, t9, t3      // y+cred
    addu        s3, t9, s3      // y+cgreen
    addu        v1, t9, v1      // y+cblue
    lbu         AT, 0(t3)
    lbu         s7, 0(s3)
    lbu         ra, 0(v1)
    lbu         v0, -1(t1)
    addu        t3, v0, s4      // y+cred
    addu        s3, v0, s5      // y+cgreen
    addu        v1, v0, s6      // y+cblue
    addu        t3, t9, t3      // y+cred
    addu        s3, t9, s3      // y+cgreen
    addu        v1, t9, v1      // y+cblue
    lbu         t3, 0(t3)
    lbu         s3, 0(s3)
    lbu         v1, 0(v1)
    lbu         v0, 0(t2)

    STORE_H2V2_2_PIXELS AT, s7, ra, t3, s3, v1, t4

    addu        t3, v0, s4      // y+cred
    addu        s3, v0, s5      // y+cgreen
    addu        v1, v0, s6      // y+cblue
    addu        t3, t9, t3      // y+cred
    addu        s3, t9, s3      // y+cgreen
    addu        v1, t9, v1      // y+cblue
    lbu         AT, 0(t3)
    lbu         s7, 0(s3)
    lbu         ra, 0(v1)
    lbu         v0, 1(t2)
    addiu       t2, t2, 2
    addu        t3, v0, s4      // y+cred
    addu        s3, v0, s5      // y+cgreen
    addu        v1, v0, s6      // y+cblue
    addu        t3, t9, t3      // y+cred
    addu        s3, t9, s3      // y+cgreen
    addu        v1, t9, v1      // y+cblue
    lbu         t3, 0(t3)
    lbu         s3, 0(s3)
    lbu         v1, 0(v1)

    STORE_H2V2_2_PIXELS AT, s7, ra, t3, s3, v1, t7

    bne         t0, t5, 1b
     nop
2:
    andi        t0, a0, 1
    beqz        t0, 4f
     lbu        t3, 0(t5)
    lbu         s3, 0(t6)
    addiu       t3, t3, -128    // (cb - 128)
    addiu       s3, s3, -128    // (cr - 128)
    mult        $ac1, s1, t3
    madd        $ac1, s2, s3
    sll         s3, s3, 15
    sll         t3, t3, 15
    lbu         v0, 0(t1)
    extr_r.w    s5, $ac1, 16
    mulq_rs.w   s4, t8, s3      // s4 = (C1 * cr + ONE_HALF)>> SCALEBITS
    mulq_rs.w   s6, s0, t3      // s6 = (C2 * cb + ONE_HALF)>> SCALEBITS
    addu        t3, v0, s4      // y+cred
    addu        s3, v0, s5      // y+cgreen
    addu        v1, v0, s6      // y+cblue
    addu        t3, t9, t3      // y+cred
    addu        s3, t9, s3      // y+cgreen
    addu        v1, t9, v1      // y+cblue
    lbu         t3, 0(t3)
    lbu         s3, 0(s3)
    lbu         v1, 0(v1)
    lbu         v0, 0(t2)

    STORE_H2V2_1_PIXEL t3, s3, v1, t4

    addu        t3, v0, s4      // y+cred
    addu        s3, v0, s5      // y+cgreen
    addu        v1, v0, s6      // y+cblue
    addu        t3, t9, t3      // y+cred
    addu        s3, t9, s3      // y+cgreen
    addu        v1, t9, v1      // y+cblue
    lbu         t3, 0(t3)
    lbu         s3, 0(s3)
    lbu         v1, 0(v1)

    STORE_H2V2_1_PIXEL t3, s3, v1, t7
4:
    RESTORE_REGS_FROM_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, ra

    j           ra
     nop

END(jsimd_h2v2_\colorid\()_merged_upsample_dspr2)

.purgem STORE_H2V2_1_PIXEL
.purgem STORE_H2V2_2_PIXELS
.endm

/*------------------------------------id -- pix R1 G1 B1 A1 R2 G2 B2 A2 */
GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2 extrgb,  6, 0, 1, 2, 6, 3, 4, 5, 6
GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2 extbgr,  6, 2, 1, 0, 3, 5, 4, 3, 6
GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2 extrgbx, 8, 0, 1, 2, 3, 4, 5, 6, 7
GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2 extbgrx, 8, 2, 1, 0, 3, 6, 5, 4, 7
GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2 extxbgr, 8, 3, 2, 1, 0, 7, 6, 5, 4
GENERATE_H2V2_MERGED_UPSAMPLE_DSPR2 extxrgb, 8, 1, 2, 3, 0, 5, 6, 7, 4


/*****************************************************************************/
/*
 * jsimd_h2v1_merged_upsample_dspr2
 * jsimd_h2v1_extrgb_merged_upsample_dspr2
 * jsimd_h2v1_extrgbx_merged_upsample_dspr2
 * jsimd_h2v1_extbgr_merged_upsample_dspr2
 * jsimd_h2v1_extbgrx_merged_upsample_dspr2
 * jsimd_h2v1_extxbgr_merged_upsample_dspr2
 * jsimd_h2v1_extxrgb_merged_upsample_dspr2
 *
 * Merged h2v1 upsample routines
 */

.macro GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2  colorid, pixel_size, \
                                            r1_offs, g1_offs, \
                                            b1_offs, a1_offs, \
                                            r2_offs, g2_offs, \
                                            b2_offs, a2_offs

.macro STORE_H2V1_2_PIXELS  scratch0 scratch1 scratch2 scratch3 scratch4 \
                            scratch5 outptr
    sb          \scratch0, \r1_offs(\outptr)
    sb          \scratch1, \g1_offs(\outptr)
    sb          \scratch2, \b1_offs(\outptr)
    sb          \scratch3, \r2_offs(\outptr)
    sb          \scratch4, \g2_offs(\outptr)
    sb          \scratch5, \b2_offs(\outptr)
.if (\pixel_size == 8)
    li          t0, 0xFF
    sb          t0, \a1_offs(\outptr)
    sb          t0, \a2_offs(\outptr)
.endif
    addiu       \outptr, \pixel_size
.endm

.macro STORE_H2V1_1_PIXEL  scratch0 scratch1 scratch2 outptr
    sb          \scratch0, \r1_offs(\outptr)
    sb          \scratch1, \g1_offs(\outptr)
    sb          \scratch2, \b1_offs(\outptr)
.if (\pixel_size == 8)
    li          t0, 0xFF
    sb          t0, \a1_offs(\outptr)
.endif
.endm

LEAF_DSPR2(jsimd_h2v1_\colorid\()_merged_upsample_dspr2)
/*
 * a0     = cinfo->output_width
 * a1     = input_buf
 * a2     = in_row_group_ctr
 * a3     = output_buf
 * 16(sp) = range_limit
 */
    SAVE_REGS_ON_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, ra

    li          t0, 0xe6ea
    lw          t1, 0(a1)         // t1 = input_buf[0]
    lw          t2, 4(a1)         // t2 = input_buf[1]
    lw          t3, 8(a1)         // t3 = input_buf[2]
    lw          t8, 56(sp)        // t8 = range_limit
    addiu       s1, t0, 0x7fff    // s1 = 0x166e9 [FIX(1.40200)]
    addiu       s2, s1, 0x5eb9    // s2 = 0x1c5a2 [FIX(1.77200)]
    addiu       s0, t0, 0x9916    // s0 = 0x8000
    addiu       s4, zero, 0xa7e6  // s4 = 0xffffa7e6 [-FIX(0.34414)]
    xori        s3, s4, 0xeec8    // s3 = 0xffff492e [-FIX(0.71414)]
    srl         t0, a0, 1
    sll         t4, a2, 2
    lwx         s5, t4(t1)      // s5 = inptr0
    lwx         s6, t4(t2)      // s6 = inptr1
    lwx         s7, t4(t3)      // s7 = inptr2
    lw          t7, 0(a3)       // t7 = outptr
    blez        t0, 2f
     addu       t9, s6, t0      // t9 = end address
1:
    lbu         t2, 0(s6)       // t2 = cb
    lbu         t0, 0(s7)       // t0 = cr
    lbu         t1, 0(s5)       // t1 = y
    addiu       t2, t2, -128    // t2 = cb - 128
    addiu       t0, t0, -128    // t0 = cr - 128
    mult        $ac1, s4, t2
    madd        $ac1, s3, t0
    sll         t0, t0, 15
    sll         t2, t2, 15
    mulq_rs.w   t0, s1, t0      // t0 = (C1*cr + ONE_HALF)>> SCALEBITS
    extr_r.w    t5, $ac1, 16
    mulq_rs.w   t6, s2, t2      // t6 = (C2*cb + ONE_HALF)>> SCALEBITS
    addiu       s7, s7, 1
    addiu       s6, s6, 1
    addu        t2, t1, t0      // t2 = y + cred
    addu        t3, t1, t5      // t3 = y + cgreen
    addu        t4, t1, t6      // t4 = y + cblue
    addu        t2, t8, t2
    addu        t3, t8, t3
    addu        t4, t8, t4
    lbu         t1, 1(s5)
    lbu         v0, 0(t2)
    lbu         v1, 0(t3)
    lbu         ra, 0(t4)
    addu        t2, t1, t0
    addu        t3, t1, t5
    addu        t4, t1, t6
    addu        t2, t8, t2
    addu        t3, t8, t3
    addu        t4, t8, t4
    lbu         t2, 0(t2)
    lbu         t3, 0(t3)
    lbu         t4, 0(t4)

    STORE_H2V1_2_PIXELS v0, v1, ra, t2, t3, t4, t7

    bne         t9, s6, 1b
     addiu      s5, s5, 2
2:
    andi        t0, a0, 1
    beqz        t0, 4f
     nop
3:
    lbu         t2, 0(s6)
    lbu         t0, 0(s7)
    lbu         t1, 0(s5)
    addiu       t2, t2, -128    // (cb - 128)
    addiu       t0, t0, -128    // (cr - 128)
    mul         t3, s4, t2
    mul         t4, s3, t0
    sll         t0, t0, 15
    sll         t2, t2, 15
    mulq_rs.w   t0, s1, t0      // (C1*cr + ONE_HALF)>> SCALEBITS
    mulq_rs.w   t6, s2, t2      // (C2*cb + ONE_HALF)>> SCALEBITS
    addu        t3, t3, s0
    addu        t3, t4, t3
    sra         t5, t3, 16      // (C4*cb + ONE_HALF + C3*cr)>> SCALEBITS
    addu        t2, t1, t0      // y + cred
    addu        t3, t1, t5      // y + cgreen
    addu        t4, t1, t6      // y + cblue
    addu        t2, t8, t2
    addu        t3, t8, t3
    addu        t4, t8, t4
    lbu         t2, 0(t2)
    lbu         t3, 0(t3)
    lbu         t4, 0(t4)

    STORE_H2V1_1_PIXEL t2, t3, t4, t7
4:
    RESTORE_REGS_FROM_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, ra

    j           ra
     nop

END(jsimd_h2v1_\colorid\()_merged_upsample_dspr2)

.purgem STORE_H2V1_1_PIXEL
.purgem STORE_H2V1_2_PIXELS
.endm

/*------------------------------------id -- pix R1 G1 B1 A1 R2 G2 B2 A2 */
GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2 extrgb,  6, 0, 1, 2, 6, 3, 4, 5, 6
GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2 extbgr,  6, 2, 1, 0, 3, 5, 4, 3, 6
GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2 extrgbx, 8, 0, 1, 2, 3, 4, 5, 6, 7
GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2 extbgrx, 8, 2, 1, 0, 3, 6, 5, 4, 7
GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2 extxbgr, 8, 3, 2, 1, 0, 7, 6, 5, 4
GENERATE_H2V1_MERGED_UPSAMPLE_DSPR2 extxrgb, 8, 1, 2, 3, 0, 5, 6, 7, 4


/*****************************************************************************/
/*
 * jsimd_h2v2_fancy_upsample_dspr2
 *
 * Fancy processing for the common case of 2:1 horizontal and 2:1 vertical.
 */
LEAF_DSPR2(jsimd_h2v2_fancy_upsample_dspr2)
/*
 * a0 = cinfo->max_v_samp_factor
 * a1 = downsampled_width
 * a2 = input_data
 * a3 = output_data_ptr
 */
    SAVE_REGS_ON_STACK 24, s0, s1, s2, s3, s4, s5

    li            s4, 0
    lw            s2, 0(a3)       // s2 = *output_data_ptr
0:
    li            t9, 2
    lw            s1, -4(a2)      // s1 = inptr1

1:
    lw            s0, 0(a2)       // s0 = inptr0
    lwx           s3, s4(s2)
    addiu         s5, a1, -2      // s5 = downsampled_width - 2
    srl           t4, s5, 1
    sll           t4, t4, 1
    lbu           t0, 0(s0)
    lbu           t1, 1(s0)
    lbu           t2, 0(s1)
    lbu           t3, 1(s1)
    addiu         s0, 2
    addiu         s1, 2
    addu          t8, s0, t4      // t8 = end address
    andi          s5, s5, 1       // s5 = residual
    sll           t4, t0, 1
    sll           t6, t1, 1
    addu          t0, t0, t4      // t0 = (*inptr0++) * 3
    addu          t1, t1, t6      // t1 = (*inptr0++) * 3
    addu          t7, t0, t2      // t7 = thiscolsum
    addu          t6, t1, t3      // t5 = nextcolsum
    sll           t0, t7, 2       // t0 = thiscolsum * 4
    subu          t1, t0, t7      // t1 = thiscolsum * 3
    shra_r.w      t0, t0, 4
    addiu         t1, 7
    addu          t1, t1, t6
    srl           t1, t1, 4
    sb            t0, 0(s3)
    sb            t1, 1(s3)
    beq           t8, s0, 22f     // skip to final iteration if width == 3
     addiu        s3, 2
2:
    lh            t0, 0(s0)       // t0 = A3|A2
    lh            t2, 0(s1)       // t2 = B3|B2
    addiu         s0, 2
    addiu         s1, 2
    preceu.ph.qbr t0, t0          // t0 = 0|A3|0|A2
    preceu.ph.qbr t2, t2          // t2 = 0|B3|0|B2
    shll.ph       t1, t0, 1
    sll           t3, t6, 1
    addu.ph       t0, t1, t0      // t0 = A3*3|A2*3
    addu          t3, t3, t6      // t3 = this * 3
    addu.ph       t0, t0, t2      // t0 = next2|next1
    addu          t1, t3, t7
    andi          t7, t0, 0xFFFF  // t7 = next1
    sll           t2, t7, 1
    addu          t2, t7, t2      // t2 = next1*3
    addu          t4, t2, t6
    srl           t6, t0, 16      // t6 = next2
    shra_r.w      t1, t1, 4       // t1 = (this*3 + last + 8) >> 4
    addu          t0, t3, t7
    addiu         t0, 7
    srl           t0, t0, 4       // t0 = (this*3 + next1 + 7) >> 4
    shra_r.w      t4, t4, 4       // t3 = (next1*3 + this + 8) >> 4
    addu          t2, t2, t6
    addiu         t2, 7
    srl           t2, t2, 4       // t2 = (next1*3 + next2 + 7) >> 4
    sb            t1, 0(s3)
    sb            t0, 1(s3)
    sb            t4, 2(s3)
    sb            t2, 3(s3)
    bne           t8, s0, 2b
     addiu        s3, 4
22:
    beqz          s5, 4f
     addu         t8, s0, s5
3:
    lbu           t0, 0(s0)
    lbu           t2, 0(s1)
    addiu         s0, 1
    addiu         s1, 1
    sll           t3, t6, 1
    sll           t1, t0, 1
    addu          t1, t0, t1      // t1 = inptr0 * 3
    addu          t3, t3, t6      // t3 = thiscolsum * 3
    addu          t5, t1, t2
    addu          t1, t3, t7
    shra_r.w      t1, t1, 4
    addu          t0, t3, t5
    addiu         t0, 7
    srl           t0, t0, 4
    sb            t1, 0(s3)
    sb            t0, 1(s3)
    addiu         s3, 2
    move          t7, t6
    bne           t8, s0, 3b
     move         t6, t5
4:
    sll           t0, t6, 2       // t0 = thiscolsum * 4
    subu          t1, t0, t6      // t1 = thiscolsum * 3
    addu          t1, t1, t7
    addiu         s4, 4
    shra_r.w      t1, t1, 4
    addiu         t0, 7
    srl           t0, t0, 4
    sb            t1, 0(s3)
    sb            t0, 1(s3)
    addiu         t9, -1
    addiu         s3, 2
    bnez          t9, 1b
     lw           s1, 4(a2)
    srl           t0, s4, 2
    subu          t0, a0, t0
    bgtz          t0, 0b
     addiu        a2, 4

    RESTORE_REGS_FROM_STACK 24, s0, s1, s2, s3, s4, s5

    j             ra
     nop
END(jsimd_h2v2_fancy_upsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_h2v1_fancy_upsample_dspr2)
/*
 * a0 = cinfo->max_v_samp_factor
 * a1 = downsampled_width
 * a2 = input_data
 * a3 = output_data_ptr
 */
    SAVE_REGS_ON_STACK 16, s0, s1, s2, s3

    .set at

    beqz          a0, 3f
     sll          t0, a0, 2
    lw            s1, 0(a3)
    li            s3, 0x10001
    addu          s0, s1, t0
0:
    addiu         t8, a1, -2
    srl           t9, t8, 2
    lw            t7, 0(a2)
    lw            s2, 0(s1)
    lbu           t0, 0(t7)
    lbu           t1, 1(t7)       // t1 = inptr[1]
    sll           t2, t0, 1
    addu          t2, t2, t0      // t2 = invalue*3
    addu          t2, t2, t1
    shra_r.w      t2, t2, 2
    sb            t0, 0(s2)
    sb            t2, 1(s2)
    beqz          t9, 11f
     addiu        s2, 2
1:
    ulw           t0, 0(t7)       // t0 = |P3|P2|P1|P0|
    ulw           t1, 1(t7)
    ulh           t2, 4(t7)       // t2 = |0|0|P5|P4|
    preceu.ph.qbl t3, t0          // t3 = |0|P3|0|P2|
    preceu.ph.qbr t0, t0          // t0 = |0|P1|0|P0|
    preceu.ph.qbr t2, t2          // t2 = |0|P5|0|P4|
    preceu.ph.qbl t4, t1          // t4 = |0|P4|0|P3|
    preceu.ph.qbr t1, t1          // t1 = |0|P2|0|P1|
    shll.ph       t5, t4, 1
    shll.ph       t6, t1, 1
    addu.ph       t5, t5, t4      // t5 = |P4*3|P3*3|
    addu.ph       t6, t6, t1      // t6 = |P2*3|P1*3|
    addu.ph       t4, t3, s3
    addu.ph       t0, t0, s3
    addu.ph       t4, t4, t5
    addu.ph       t0, t0, t6
    shrl.ph       t4, t4, 2       // t4 = |0|P3|0|P2|
    shrl.ph       t0, t0, 2       // t0 = |0|P1|0|P0|
    addu.ph       t2, t2, t5
    addu.ph       t3, t3, t6
    shra_r.ph     t2, t2, 2       // t2 = |0|P5|0|P4|
    shra_r.ph     t3, t3, 2       // t3 = |0|P3|0|P2|
    shll.ph       t2, t2, 8
    shll.ph       t3, t3, 8
    or            t2, t4, t2
    or            t3, t3, t0
    addiu         t9, -1
    usw           t3, 0(s2)
    usw           t2, 4(s2)
    addiu         s2, 8
    bgtz          t9, 1b
     addiu        t7, 4
11:
    andi          t8, 3
    beqz          t8, 22f
     addiu        t7, 1

2:
    lbu           t0, 0(t7)
    addiu         t7, 1
    sll           t1, t0, 1
    addu          t2, t0, t1      // t2 = invalue
    lbu           t3, -2(t7)
    lbu           t4, 0(t7)
    addiu         t3, 1
    addiu         t4, 2
    addu          t3, t3, t2
    addu          t4, t4, t2
    srl           t3, 2
    srl           t4, 2
    sb            t3, 0(s2)
    sb            t4, 1(s2)
    addiu         t8, -1
    bgtz          t8, 2b
     addiu        s2, 2

22:
    lbu           t0, 0(t7)
    lbu           t2, -1(t7)
    sll           t1, t0, 1
    addu          t1, t1, t0      // t1 = invalue * 3
    addu          t1, t1, t2
    addiu         t1, 1
    srl           t1, t1, 2
    sb            t1, 0(s2)
    sb            t0, 1(s2)
    addiu         s1, 4
    bne           s1, s0, 0b
     addiu        a2, 4
3:
    RESTORE_REGS_FROM_STACK 16, s0, s1, s2, s3

    j             ra
     nop
END(jsimd_h2v1_fancy_upsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_h2v1_downsample_dspr2)
/*
 * a0     = cinfo->image_width
 * a1     = cinfo->max_v_samp_factor
 * a2     = compptr->v_samp_factor
 * a3     = compptr->width_in_blocks
 * 16(sp) = input_data
 * 20(sp) = output_data
 */
    .set at

    SAVE_REGS_ON_STACK 24, s0, s1, s2, s3, s4

    beqz        a2, 7f
     lw         s1, 44(sp)      // s1 = output_data
    lw          s0, 40(sp)      // s0 = input_data
    srl         s2, a0, 2
    andi        t9, a0, 2
    srl         t7, t9, 1
    addu        s2, t7, s2
    sll         t0, a3, 3       // t0 = width_in_blocks*DCT
    srl         t7, t0, 1
    subu        s2, t7, s2
0:
    andi        t6, a0, 1       // t6 = temp_index
    addiu       t6, -1
    lw          t4, 0(s1)       // t4 = outptr
    lw          t5, 0(s0)       // t5 = inptr0
    li          s3, 0           // s3 = bias
    srl         t7, a0, 1       // t7 = image_width1
    srl         s4, t7, 2
    andi        t8, t7, 3
1:
    ulhu        t0, 0(t5)
    ulhu        t1, 2(t5)
    ulhu        t2, 4(t5)
    ulhu        t3, 6(t5)
    raddu.w.qb  t0, t0
    raddu.w.qb  t1, t1
    raddu.w.qb  t2, t2
    raddu.w.qb  t3, t3
    shra.ph     t0, t0, 1
    shra_r.ph   t1, t1, 1
    shra.ph     t2, t2, 1
    shra_r.ph   t3, t3, 1
    sb          t0, 0(t4)
    sb          t1, 1(t4)
    sb          t2, 2(t4)
    sb          t3, 3(t4)
    addiu       s4, -1
    addiu       t4, 4
    bgtz        s4, 1b
     addiu      t5, 8
    beqz        t8, 3f
     addu       s4, t4, t8
2:
    ulhu        t0, 0(t5)
    raddu.w.qb  t0, t0
    addqh.w     t0, t0, s3
    xori        s3, s3, 1
    sb          t0, 0(t4)
    addiu       t4, 1
    bne         t4, s4, 2b
     addiu      t5, 2
3:
    lbux        t1, t6(t5)
    sll         t1, 1
    addqh.w     t2, t1, s3      // t2 = pixval1
    xori        s3, s3, 1
    addqh.w     t3, t1, s3      // t3 = pixval2
    blez        s2, 5f
     append     t3, t2,  8
    addu        t5, t4, s2      // t5 = loop_end2
4:
    ush         t3, 0(t4)
    addiu       s2, -1
    bgtz        s2, 4b
     addiu      t4,  2
5:
    beqz        t9, 6f
     nop
    sb          t2, 0(t4)
6:
    addiu       s1, 4
    addiu       a2, -1
    bnez        a2, 0b
     addiu      s0, 4
7:
    RESTORE_REGS_FROM_STACK 24, s0, s1, s2, s3, s4

    j           ra
    nop
END(jsimd_h2v1_downsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_h2v2_downsample_dspr2)
/*
 * a0     = cinfo->image_width
 * a1     = cinfo->max_v_samp_factor
 * a2     = compptr->v_samp_factor
 * a3     = compptr->width_in_blocks
 * 16(sp) = input_data
 * 20(sp) = output_data
 */
    .set at

    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    beqz        a2, 8f
     lw         s1, 52(sp)      // s1 = output_data
    lw          s0, 48(sp)      // s0 = input_data

    andi        t6, a0, 1       // t6 = temp_index
    addiu       t6, -1
    srl         t7, a0, 1       // t7 = image_width1
    srl         s4, t7, 2
    andi        t8, t7, 3
    andi        t9, a0, 2
    srl         s2, a0, 2
    srl         t7, t9, 1
    addu        s2, t7, s2
    sll         t0, a3, 3       // s2 = width_in_blocks*DCT
    srl         t7, t0, 1
    subu        s2, t7, s2
0:
    lw          t4, 0(s1)       // t4 = outptr
    lw          t5, 0(s0)       // t5 = inptr0
    lw          s7, 4(s0)       // s7 = inptr1
    li          s6, 1           // s6 = bias
2:
    ulw         t0, 0(t5)       // t0 = |P3|P2|P1|P0|
    ulw         t1, 0(s7)       // t1 = |Q3|Q2|Q1|Q0|
    ulw         t2, 4(t5)
    ulw         t3, 4(s7)
    precrq.ph.w t7, t0, t1      // t2 = |P3|P2|Q3|Q2|
    ins         t0, t1, 16, 16  // t0 = |Q1|Q0|P1|P0|
    raddu.w.qb  t1, t7
    raddu.w.qb  t0, t0
    shra_r.w    t1, t1, 2
    addiu       t0, 1
    srl         t0, 2
    precrq.ph.w t7, t2, t3
    ins         t2, t3, 16, 16
    raddu.w.qb  t7, t7
    raddu.w.qb  t2, t2
    shra_r.w    t7, t7, 2
    addiu       t2, 1
    srl         t2, 2
    sb          t0, 0(t4)
    sb          t1, 1(t4)
    sb          t2, 2(t4)
    sb          t7, 3(t4)
    addiu       t4, 4
    addiu       t5, 8
    addiu       s4, s4, -1
    bgtz        s4, 2b
     addiu      s7, 8
    beqz        t8, 4f
     addu       t8, t4, t8
3:
    ulhu        t0, 0(t5)
    ulhu        t1, 0(s7)
    ins         t0, t1, 16, 16
    raddu.w.qb  t0, t0
    addu        t0, t0, s6
    srl         t0, 2
    xori        s6, s6, 3
    sb          t0, 0(t4)
    addiu       t5, 2
    addiu       t4, 1
    bne         t8, t4, 3b
     addiu      s7, 2
4:
    lbux        t1, t6(t5)
    sll         t1, 1
    lbux        t0, t6(s7)
    sll         t0, 1
    addu        t1, t1, t0
    addu        t3, t1, s6
    srl         t0, t3, 2       // t2 = pixval1
    xori        s6, s6, 3
    addu        t2, t1, s6
    srl         t1, t2, 2       // t3 = pixval2
    blez        s2, 6f
     append     t1, t0, 8
5:
    ush         t1, 0(t4)
    addiu       s2, -1
    bgtz        s2, 5b
     addiu      t4, 2
6:
    beqz        t9, 7f
     nop
    sb          t0, 0(t4)
7:
    addiu       s1, 4
    addiu       a2, -1
    bnez        a2, 0b
     addiu      s0, 8
8:
    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop
END(jsimd_h2v2_downsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_h2v2_smooth_downsample_dspr2)
/*
 * a0     = input_data
 * a1     = output_data
 * a2     = compptr->v_samp_factor
 * a3     = cinfo->max_v_samp_factor
 * 16(sp) = cinfo->smoothing_factor
 * 20(sp) = compptr->width_in_blocks
 * 24(sp) = cinfo->image_width
 */
    .set at

    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    lw          s7, 52(sp)      // compptr->width_in_blocks
    lw          s0, 56(sp)      // cinfo->image_width
    lw          s6, 48(sp)      // cinfo->smoothing_factor
    sll         s7, 3           // output_cols = width_in_blocks * DCTSIZE
    sll         v0, s7, 1
    subu        v0, v0, s0
    blez        v0, 2f
    move        v1, zero
    addiu       t0, a3, 2       // t0 = cinfo->max_v_samp_factor + 2
0:
    addiu       t1, a0, -4
    sll         t2, v1, 2
    lwx         t1, t2(t1)
    move        t3, v0
    addu        t1, t1, s0
    lbu         t2, -1(t1)
1:
    addiu       t3, t3, -1
    sb          t2, 0(t1)
    bgtz        t3, 1b
    addiu       t1, t1, 1
    addiu       v1, v1, 1
    bne         v1, t0, 0b
    nop
2:
    li          v0, 80
    mul         v0, s6, v0
    li          v1, 16384
    move        t4, zero
    move        t5, zero
    subu        t6, v1, v0      // t6 = 16384 - tmp_smoot_f * 80
    sll         t7, s6, 4       // t7 = tmp_smoot_f * 16
3:
/* Special case for first column: pretend column -1 is same as column 0 */
    sll         v0, t4, 2
    lwx         t8, v0(a1)      //  outptr = output_data[outrow]
    sll         v1, t5, 2
    addiu       t9, v1, 4
    addiu       s0, v1, -4
    addiu       s1, v1, 8
    lwx         s2, v1(a0)      // inptr0 = input_data[inrow]
    lwx         t9, t9(a0)      // inptr1 = input_data[inrow+1]
    lwx         s0, s0(a0)      // above_ptr = input_data[inrow-1]
    lwx         s1, s1(a0)      // below_ptr = input_data[inrow+2]
    lh          v0, 0(s2)
    lh          v1, 0(t9)
    lh          t0, 0(s0)
    lh          t1, 0(s1)
    ins         v0, v1, 16, 16
    ins         t0, t1, 16, 16
    raddu.w.qb  t2, v0
    raddu.w.qb  s3, t0
    lbu         v0, 0(s2)
    lbu         v1, 2(s2)
    lbu         t0, 0(t9)
    lbu         t1, 2(t9)
    addu        v0, v0, v1
    mult        $ac1, t2, t6
    addu        t0, t0, t1
    lbu         t2, 2(s0)
    addu        t0, t0, v0
    lbu         t3, 2(s1)
    addu        s3, t0, s3
    lbu         v0, 0(s0)
    lbu         t0, 0(s1)
    sll         s3, s3, 1
    addu        v0, v0, t2
    addu        t0, t0, t3
    addu        t0, t0, v0
    addu        s3, t0, s3
    madd        $ac1, s3, t7
    extr_r.w    v0, $ac1, 16
    addiu       t8, t8, 1
    addiu       s2, s2, 2
    addiu       t9, t9, 2
    addiu       s0, s0, 2
    addiu       s1, s1, 2
    sb          v0, -1(t8)
    addiu       s4, s7, -2
    and         s4, s4, 3
    addu        s5, s4, t8      // end address
4:
    lh          v0, 0(s2)
    lh          v1, 0(t9)
    lh          t0, 0(s0)
    lh          t1, 0(s1)
    ins         v0, v1, 16, 16
    ins         t0, t1, 16, 16
    raddu.w.qb  t2, v0
    raddu.w.qb  s3, t0
    lbu         v0, -1(s2)
    lbu         v1, 2(s2)
    lbu         t0, -1(t9)
    lbu         t1, 2(t9)
    addu        v0, v0, v1
    mult        $ac1, t2, t6
    addu        t0, t0, t1
    lbu         t2, 2(s0)
    addu        t0, t0, v0
    lbu         t3, 2(s1)
    addu        s3, t0, s3
    lbu         v0, -1(s0)
    lbu         t0, -1(s1)
    sll         s3, s3, 1
    addu        v0, v0, t2
    addu        t0, t0, t3
    addu        t0, t0, v0
    addu        s3, t0, s3
    madd        $ac1, s3, t7
    extr_r.w    t2, $ac1, 16
    addiu       t8, t8, 1
    addiu       s2, s2, 2
    addiu       t9, t9, 2
    addiu       s0, s0, 2
    sb          t2, -1(t8)
    bne         s5, t8, 4b
    addiu       s1, s1, 2
    addiu       s5, s7, -2
    subu        s5, s5, s4
    addu        s5, s5, t8      // end address
5:
    lh          v0, 0(s2)
    lh          v1, 0(t9)
    lh          t0, 0(s0)
    lh          t1, 0(s1)
    ins         v0, v1, 16, 16
    ins         t0, t1, 16, 16
    raddu.w.qb  t2, v0
    raddu.w.qb  s3, t0
    lbu         v0, -1(s2)
    lbu         v1, 2(s2)
    lbu         t0, -1(t9)
    lbu         t1, 2(t9)
    addu        v0, v0, v1
    mult        $ac1, t2, t6
    addu        t0, t0, t1
    lbu         t2, 2(s0)
    addu        t0, t0, v0
    lbu         t3, 2(s1)
    addu        s3, t0, s3
    lbu         v0, -1(s0)
    lbu         t0, -1(s1)
    sll         s3, s3, 1
    addu        v0, v0, t2
    addu        t0, t0, t3
    lh          v1, 2(t9)
    addu        t0, t0, v0
    lh          v0, 2(s2)
    addu        s3, t0, s3
    lh          t0, 2(s0)
    lh          t1, 2(s1)
    madd        $ac1, s3, t7
    extr_r.w    t2, $ac1, 16
    ins         t0, t1, 16, 16
    ins         v0, v1, 16, 16
    raddu.w.qb  s3, t0
    lbu         v1, 4(s2)
    lbu         t0, 1(t9)
    lbu         t1, 4(t9)
    sb          t2, 0(t8)
    raddu.w.qb  t3, v0
    lbu         v0, 1(s2)
    addu        t0, t0, t1
    mult        $ac1, t3, t6
    addu        v0, v0, v1
    lbu         t2, 4(s0)
    addu        t0, t0, v0
    lbu         v0, 1(s0)
    addu        s3, t0, s3
    lbu         t0, 1(s1)
    lbu         t3, 4(s1)
    addu        v0, v0, t2
    sll         s3, s3, 1
    addu        t0, t0, t3
    lh          v1, 4(t9)
    addu        t0, t0, v0
    lh          v0, 4(s2)
    addu        s3, t0, s3
    lh          t0, 4(s0)
    lh          t1, 4(s1)
    madd        $ac1, s3, t7
    extr_r.w    t2, $ac1, 16
    ins         t0, t1, 16, 16
    ins         v0, v1, 16, 16
    raddu.w.qb  s3, t0
    lbu         v1, 6(s2)
    lbu         t0, 3(t9)
    lbu         t1, 6(t9)
    sb          t2, 1(t8)
    raddu.w.qb  t3, v0
    lbu         v0, 3(s2)
    addu        t0, t0, t1
    mult        $ac1, t3, t6
    addu        v0, v0, v1
    lbu         t2, 6(s0)
    addu        t0, t0, v0
    lbu         v0, 3(s0)
    addu        s3, t0, s3
    lbu         t0, 3(s1)
    lbu         t3, 6(s1)
    addu        v0, v0, t2
    sll         s3, s3, 1
    addu        t0, t0, t3
    lh          v1, 6(t9)
    addu        t0, t0, v0
    lh          v0, 6(s2)
    addu        s3, t0, s3
    lh          t0, 6(s0)
    lh          t1, 6(s1)
    madd        $ac1, s3, t7
    extr_r.w    t3, $ac1, 16
    ins         t0, t1, 16, 16
    ins         v0, v1, 16, 16
    raddu.w.qb  s3, t0
    lbu         v1, 8(s2)
    lbu         t0, 5(t9)
    lbu         t1, 8(t9)
    sb          t3, 2(t8)
    raddu.w.qb  t2, v0
    lbu         v0, 5(s2)
    addu        t0, t0, t1
    mult        $ac1, t2, t6
    addu        v0, v0, v1
    lbu         t2, 8(s0)
    addu        t0, t0, v0
    lbu         v0, 5(s0)
    addu        s3, t0, s3
    lbu         t0, 5(s1)
    lbu         t3, 8(s1)
    addu        v0, v0, t2
    sll         s3, s3, 1
    addu        t0, t0, t3
    addiu       t8, t8, 4
    addu        t0, t0, v0
    addiu       s2, s2, 8
    addu        s3, t0, s3
    addiu       t9, t9, 8
    madd        $ac1, s3, t7
    extr_r.w    t1, $ac1, 16
    addiu       s0, s0, 8
    addiu       s1, s1, 8
    bne         s5, t8, 5b
    sb          t1, -1(t8)
/* Special case for last column */
    lh          v0, 0(s2)
    lh          v1, 0(t9)
    lh          t0, 0(s0)
    lh          t1, 0(s1)
    ins         v0, v1, 16, 16
    ins         t0, t1, 16, 16
    raddu.w.qb  t2, v0
    raddu.w.qb  s3, t0
    lbu         v0, -1(s2)
    lbu         v1, 1(s2)
    lbu         t0, -1(t9)
    lbu         t1, 1(t9)
    addu        v0, v0, v1
    mult        $ac1, t2, t6
    addu        t0, t0, t1
    lbu         t2, 1(s0)
    addu        t0, t0, v0
    lbu         t3, 1(s1)
    addu        s3, t0, s3
    lbu         v0, -1(s0)
    lbu         t0, -1(s1)
    sll         s3, s3, 1
    addu        v0, v0, t2
    addu        t0, t0, t3
    addu        t0, t0, v0
    addu        s3, t0, s3
    madd        $ac1, s3, t7
    extr_r.w    t0, $ac1, 16
    addiu       t5, t5, 2
    sb          t0, 0(t8)
    addiu       t4, t4, 1
    bne         t4, a2, 3b
    addiu       t5, t5, 2

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop

END(jsimd_h2v2_smooth_downsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_int_upsample_dspr2)
/*
 * a0     = upsample->h_expand[compptr->component_index]
 * a1     = upsample->v_expand[compptr->component_index]
 * a2     = input_data
 * a3     = output_data_ptr
 * 16(sp) = cinfo->output_width
 * 20(sp) = cinfo->max_v_samp_factor
 */
    .set at

    SAVE_REGS_ON_STACK 16, s0, s1, s2, s3

    lw          s0, 0(a3)       // s0 = output_data
    lw          s1, 32(sp)      // s1 = cinfo->output_width
    lw          s2, 36(sp)      // s2 = cinfo->max_v_samp_factor
    li          t6, 0           // t6 = inrow
    beqz        s2, 10f
     li         s3, 0           // s3 = outrow
0:
    addu        t0, a2, t6
    addu        t7, s0, s3
    lw          t3, 0(t0)       // t3 = inptr
    lw          t8, 0(t7)       // t8 = outptr
    beqz        s1, 4f
     addu       t5, t8, s1      // t5 = outend
1:
    lb          t2, 0(t3)       // t2 = invalue = *inptr++
    addiu       t3, 1
    beqz        a0, 3f
     move       t0, a0          // t0 = h_expand
2:
    sb          t2, 0(t8)
    addiu       t0, -1
    bgtz        t0, 2b
     addiu      t8, 1
3:
    bgt         t5, t8, 1b
     nop
4:
    addiu       t9, a1, -1      // t9 = v_expand - 1
    blez        t9, 9f
     nop
5:
    lw          t3, 0(s0)
    lw          t4, 4(s0)
    subu        t0, s1, 0xF
    blez        t0, 7f
     addu       t5, t3, s1      // t5 = end address
    andi        t7, s1, 0xF     // t7 = residual
    subu        t8, t5, t7
6:
    ulw         t0, 0(t3)
    ulw         t1, 4(t3)
    ulw         t2, 8(t3)
    usw         t0, 0(t4)
    ulw         t0, 12(t3)
    usw         t1, 4(t4)
    usw         t2, 8(t4)
    usw         t0, 12(t4)
    addiu       t3, 16
    bne         t3, t8, 6b
     addiu      t4, 16
    beqz        t7, 8f
     nop
7:
    lbu         t0, 0(t3)
    sb          t0, 0(t4)
    addiu       t3, 1
    bne         t3, t5, 7b
     addiu      t4, 1
8:
    addiu       t9, -1
    bgtz        t9, 5b
     addiu      s0, 8
9:
    addu        s3, s3, a1
    bne         s3, s2, 0b
     addiu      t6, 1
10:
    RESTORE_REGS_FROM_STACK 16, s0, s1, s2, s3

    j           ra
     nop
END(jsimd_int_upsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_h2v1_upsample_dspr2)
/*
 * a0 = cinfo->max_v_samp_factor
 * a1 = cinfo->output_width
 * a2 = input_data
 * a3 = output_data_ptr
 */
    lw          t7, 0(a3)       // t7 = output_data
    andi        t8, a1, 0xf     // t8 = residual
    sll         t0, a0, 2
    blez        a0, 4f
     addu       t9, t7, t0      // t9 = output_data end address
0:
    lw          t5, 0(t7)       // t5 = outptr
    lw          t6, 0(a2)       // t6 = inptr
    addu        t3, t5, a1      // t3 = outptr + output_width (end address)
    subu        t3, t8          // t3 = end address - residual
    beq         t5, t3, 2f
     move       t4, t8
1:
    ulw         t0, 0(t6)       // t0 = |P3|P2|P1|P0|
    ulw         t2, 4(t6)       // t2 = |P7|P6|P5|P4|
    srl         t1, t0, 16      // t1 = |X|X|P3|P2|
    ins         t0, t0, 16, 16  // t0 = |P1|P0|P1|P0|
    ins         t1, t1, 16, 16  // t1 = |P3|P2|P3|P2|
    ins         t0, t0, 8, 16   // t0 = |P1|P1|P0|P0|
    ins         t1, t1, 8, 16   // t1 = |P3|P3|P2|P2|
    usw         t0, 0(t5)
    usw         t1, 4(t5)
    srl         t0, t2, 16      // t0 = |X|X|P7|P6|
    ins         t2, t2, 16, 16  // t2 = |P5|P4|P5|P4|
    ins         t0, t0, 16, 16  // t0 = |P7|P6|P7|P6|
    ins         t2, t2, 8, 16   // t2 = |P5|P5|P4|P4|
    ins         t0, t0, 8, 16   // t0 = |P7|P7|P6|P6|
    usw         t2, 8(t5)
    usw         t0, 12(t5)
    addiu       t5, 16
    bne         t5, t3, 1b
     addiu      t6, 8
    beqz        t8, 3f
     move       t4, t8
2:
    lbu         t1, 0(t6)
    sb          t1, 0(t5)
    sb          t1, 1(t5)
    addiu       t4, -2
    addiu       t6, 1
    bgtz        t4, 2b
     addiu      t5, 2
3:
    addiu       t7, 4
    bne         t9, t7, 0b
     addiu      a2, 4
4:
    j           ra
     nop
END(jsimd_h2v1_upsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_h2v2_upsample_dspr2)
/*
 * a0 = cinfo->max_v_samp_factor
 * a1 = cinfo->output_width
 * a2 = input_data
 * a3 = output_data_ptr
 */
    lw          t7, 0(a3)
    blez        a0, 7f
     andi       t9, a1, 0xf     // t9 = residual
0:
    lw          t6, 0(a2)       // t6 = inptr
    lw          t5, 0(t7)       // t5 = outptr
    addu        t8, t5, a1      // t8 = outptr end address
    subu        t8, t9          // t8 = end address - residual
    beq         t5, t8, 2f
     move       t4, t9
1:
    ulw         t0, 0(t6)
    srl         t1, t0, 16
    ins         t0, t0, 16, 16
    ins         t0, t0, 8, 16
    ins         t1, t1, 16, 16
    ins         t1, t1, 8, 16
    ulw         t2, 4(t6)
    usw         t0, 0(t5)
    usw         t1, 4(t5)
    srl         t3, t2, 16
    ins         t2, t2, 16, 16
    ins         t2, t2, 8, 16
    ins         t3, t3, 16, 16
    ins         t3, t3, 8, 16
    usw         t2, 8(t5)
    usw         t3, 12(t5)
    addiu       t5, 16
    bne         t5, t8, 1b
     addiu      t6, 8
    beqz        t9, 3f
     move       t4, t9
2:
    lbu         t0, 0(t6)
    sb          t0, 0(t5)
    sb          t0, 1(t5)
    addiu       t4, -2
    addiu       t6, 1
    bgtz        t4, 2b
     addiu      t5, 2
3:
    lw          t6, 0(t7)       // t6 = outptr[0]
    lw          t5, 4(t7)       // t5 = outptr[1]
    addu        t4, t6, a1      // t4 = new end address
    beq         a1, t9, 5f
     subu       t8, t4, t9
4:
    ulw         t0, 0(t6)
    ulw         t1, 4(t6)
    ulw         t2, 8(t6)
    usw         t0, 0(t5)
    ulw         t0, 12(t6)
    usw         t1, 4(t5)
    usw         t2, 8(t5)
    usw         t0, 12(t5)
    addiu       t6, 16
    bne         t6, t8, 4b
     addiu      t5, 16
    beqz        t9, 6f
     nop
5:
    lbu         t0, 0(t6)
    sb          t0, 0(t5)
    addiu       t6, 1
    bne         t6, t4, 5b
     addiu      t5, 1
6:
    addiu       t7, 8
    addiu       a0, -2
    bgtz        a0, 0b
     addiu      a2, 4
7:
    j           ra
     nop
END(jsimd_h2v2_upsample_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_islow_dspr2)
/*
 * a0 = coef_block
 * a1 = compptr->dcttable
 * a2 = output
 * a3 = range_limit
 */
    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    addiu       sp, sp, -256
    move        v0, sp
    addiu       v1, zero, 8     // v1 = DCTSIZE = 8
1:
    lh          s4, 32(a0)      // s4 = inptr[16]
    lh          s5, 64(a0)      // s5 = inptr[32]
    lh          s6, 96(a0)      // s6 = inptr[48]
    lh          t1, 112(a0)     // t1 = inptr[56]
    lh          t7, 16(a0)      // t7 = inptr[8]
    lh          t5, 80(a0)      // t5 = inptr[40]
    lh          t3, 48(a0)      // t3 = inptr[24]
    or          s4, s4, t1
    or          s4, s4, t3
    or          s4, s4, t5
    or          s4, s4, t7
    or          s4, s4, s5
    or          s4, s4, s6
    bnez        s4, 2f
     addiu      v1, v1, -1
    lh          s5, 0(a1)       // quantptr[DCTSIZE*0]
    lh          s6, 0(a0)       // inptr[DCTSIZE*0]
    mul         s5, s5, s6      // DEQUANTIZE(inptr[0], quantptr[0])
    sll         s5, s5, 2
    sw          s5, 0(v0)
    sw          s5, 32(v0)
    sw          s5, 64(v0)
    sw          s5, 96(v0)
    sw          s5, 128(v0)
    sw          s5, 160(v0)
    sw          s5, 192(v0)
    b           3f
     sw         s5, 224(v0)
2:
    lh          t0, 112(a1)
    lh          t2, 48(a1)
    lh          t4, 80(a1)
    lh          t6, 16(a1)
    mul         t0, t0, t1      // DEQUANTIZE(inptr[DCTSIZE*7], quantptr[DCTSIZE*7])
    mul         t1, t2, t3      // DEQUANTIZE(inptr[DCTSIZE*3], quantptr[DCTSIZE*3])
    mul         t2, t4, t5      // DEQUANTIZE(inptr[DCTSIZE*5], quantptr[DCTSIZE*5])
    mul         t3, t6, t7      // DEQUANTIZE(inptr[DCTSIZE*1], quantptr[DCTSIZE*1])
    lh          t4, 32(a1)
    lh          t5, 32(a0)
    lh          t6, 96(a1)
    lh          t7, 96(a0)
    addu        s0, t0, t1       // z3 = tmp0 + tmp2
    addu        s1, t1, t2       // z2 = tmp1 + tmp2
    addu        s2, t2, t3       // z4 = tmp1 + tmp3
    addu        s3, s0, s2       // z3 + z4
    addiu       t9, zero, 9633   // FIX_1_175875602
    mul         s3, s3, t9       // z5 = MULTIPLY(z3 + z4, FIX_1_175875602)
    addu        t8, t0, t3       // z1 = tmp0 + tmp3
    addiu       t9, zero, 2446   // FIX_0_298631336
    mul         t0, t0, t9       // tmp0 = MULTIPLY(tmp0, FIX_0_298631336)
    addiu       t9, zero, 16819  // FIX_2_053119869
    mul         t2, t2, t9       // tmp1 = MULTIPLY(tmp1, FIX_2_053119869)
    addiu       t9, zero, 25172  // FIX_3_072711026
    mul         t1, t1, t9       // tmp2 = MULTIPLY(tmp2, FIX_3_072711026)
    addiu       t9, zero, 12299  // FIX_1_501321110
    mul         t3, t3, t9       // tmp3 = MULTIPLY(tmp3, FIX_1_501321110)
    addiu       t9, zero, 16069  // FIX_1_961570560
    mul         s0, s0, t9       // -z3 = MULTIPLY(z3, FIX_1_961570560)
    addiu       t9, zero, 3196   // FIX_0_390180644
    mul         s2, s2, t9       // -z4 = MULTIPLY(z4, FIX_0_390180644)
    addiu       t9, zero, 7373   // FIX_0_899976223
    mul         t8, t8, t9       // -z1 = MULTIPLY(z1, FIX_0_899976223)
    addiu       t9, zero, 20995  // FIX_2_562915447
    mul         s1, s1, t9       // -z2 = MULTIPLY(z2, FIX_2_562915447)
    subu        s0, s3, s0       // z3 += z5
    addu        t0, t0, s0       // tmp0 += z3
    addu        t1, t1, s0       // tmp2 += z3
    subu        s2, s3, s2       // z4 += z5
    addu        t2, t2, s2       // tmp1 += z4
    addu        t3, t3, s2       // tmp3 += z4
    subu        t0, t0, t8       // tmp0 += z1
    subu        t1, t1, s1       // tmp2 += z2
    subu        t2, t2, s1       // tmp1 += z2
    subu        t3, t3, t8       // tmp3 += z1
    mul         s0, t4, t5       // DEQUANTIZE(inptr[DCTSIZE*2], quantptr[DCTSIZE*2])
    addiu       t9, zero, 6270   // FIX_0_765366865
    mul         s1, t6, t7       // DEQUANTIZE(inptr[DCTSIZE*6], quantptr[DCTSIZE*6])
    lh          t4, 0(a1)
    lh          t5, 0(a0)
    lh          t6, 64(a1)
    lh          t7, 64(a0)
    mul         s2, t9, s0       // MULTIPLY(z2, FIX_0_765366865)
    mul         t5, t4, t5       // DEQUANTIZE(inptr[DCTSIZE*0], quantptr[DCTSIZE*0])
    mul         t6, t6, t7       // DEQUANTIZE(inptr[DCTSIZE*4], quantptr[DCTSIZE*4])
    addiu       t9, zero, 4433   // FIX_0_541196100
    addu        s3, s0, s1       // z2 + z3
    mul         s3, s3, t9       // z1 = MULTIPLY(z2 + z3, FIX_0_541196100)
    addiu       t9, zero, 15137  // FIX_1_847759065
    mul         t8, s1, t9       // MULTIPLY(z3, FIX_1_847759065)
    addu        t4, t5, t6
    subu        t5, t5, t6
    sll         t4, t4, 13      // tmp0 = (z2 + z3) << CONST_BITS
    sll         t5, t5, 13      // tmp1 = (z2 - z3) << CONST_BITS
    addu        t7, s3, s2      // tmp3 = z1 + MULTIPLY(z2, FIX_0_765366865)
    subu        t6, s3, t8      // tmp2 = z1 + MULTIPLY(z3, -FIX_1_847759065)
    addu        s0, t4, t7
    subu        s1, t4, t7
    addu        s2, t5, t6
    subu        s3, t5, t6
    addu        t4, s0, t3
    subu        s0, s0, t3
    addu        t3, s2, t1
    subu        s2, s2, t1
    addu        t1, s3, t2
    subu        s3, s3, t2
    addu        t2, s1, t0
    subu        s1, s1, t0
    shra_r.w    t4, t4, 11
    shra_r.w    t3, t3, 11
    shra_r.w    t1, t1, 11
    shra_r.w    t2, t2, 11
    shra_r.w    s1, s1, 11
    shra_r.w    s3, s3, 11
    shra_r.w    s2, s2, 11
    shra_r.w    s0, s0, 11
    sw          t4, 0(v0)
    sw          t3, 32(v0)
    sw          t1, 64(v0)
    sw          t2, 96(v0)
    sw          s1, 128(v0)
    sw          s3, 160(v0)
    sw          s2, 192(v0)
    sw          s0, 224(v0)
3:
    addiu       a1, a1, 2
    addiu       a0, a0, 2
    bgtz        v1, 1b
     addiu      v0, v0, 4
    move        v0, sp
    addiu       v1, zero, 8
4:
    lw          t0, 8(v0)       // z2 = (JLONG)wsptr[2]
    lw          t1, 24(v0)      // z3 = (JLONG)wsptr[6]
    lw          t2, 0(v0)       // (JLONG)wsptr[0]
    lw          t3, 16(v0)      // (JLONG)wsptr[4]
    lw          s4, 4(v0)       // (JLONG)wsptr[1]
    lw          s5, 12(v0)      // (JLONG)wsptr[3]
    lw          s6, 20(v0)      // (JLONG)wsptr[5]
    lw          s7, 28(v0)      // (JLONG)wsptr[7]
    or          s4, s4, t0
    or          s4, s4, t1
    or          s4, s4, t3
    or          s4, s4, s7
    or          s4, s4, s5
    or          s4, s4, s6
    bnez        s4, 5f
     addiu      v1, v1, -1
    shra_r.w    s5, t2, 5
    andi        s5, s5, 0x3ff
    lbux        s5, s5(a3)
    lw          s1, 0(a2)
    replv.qb    s5, s5
    usw         s5, 0(s1)
    usw         s5, 4(s1)
    b           6f
     nop
5:
    addu        t4, t0, t1       // z2 + z3
    addiu       t8, zero, 4433   // FIX_0_541196100
    mul         t5, t4, t8       // z1 = MULTIPLY(z2 + z3, FIX_0_541196100)
    addiu       t8, zero, 15137  // FIX_1_847759065
    mul         t1, t1, t8       // MULTIPLY(z3, FIX_1_847759065)
    addiu       t8, zero, 6270   // FIX_0_765366865
    mul         t0, t0, t8       // MULTIPLY(z2, FIX_0_765366865)
    addu        t4, t2, t3       // (JLONG)wsptr[0] + (JLONG)wsptr[4]
    subu        t2, t2, t3       // (JLONG)wsptr[0] - (JLONG)wsptr[4]
    sll         t4, t4, 13       // tmp0 = (wsptr[0] + wsptr[4]) << CONST_BITS
    sll         t2, t2, 13       // tmp1 = (wsptr[0] - wsptr[4]) << CONST_BITS
    subu        t1, t5, t1       // tmp2 = z1 + MULTIPLY(z3, -FIX_1_847759065)
    subu        t3, t2, t1       // tmp12 = tmp1 - tmp2
    addu        t2, t2, t1       // tmp11 = tmp1 + tmp2
    addu        t5, t5, t0       // tmp3 = z1 + MULTIPLY(z2, FIX_0_765366865)
    subu        t1, t4, t5       // tmp13 = tmp0 - tmp3
    addu        t0, t4, t5       // tmp10 = tmp0 + tmp3
    lw          t4, 28(v0)       // tmp0 = (JLONG)wsptr[7]
    lw          t6, 12(v0)       // tmp2 = (JLONG)wsptr[3]
    lw          t5, 20(v0)       // tmp1 = (JLONG)wsptr[5]
    lw          t7, 4(v0)        // tmp3 = (JLONG)wsptr[1]
    addu        s0, t4, t6       // z3 = tmp0 + tmp2
    addiu       t8, zero, 9633   // FIX_1_175875602
    addu        s1, t5, t7       // z4 = tmp1 + tmp3
    addu        s2, s0, s1       // z3 + z4
    mul         s2, s2, t8       // z5 = MULTIPLY(z3 + z4, FIX_1_175875602)
    addu        s3, t4, t7       // z1 = tmp0 + tmp3
    addu        t9, t5, t6       // z2 = tmp1 + tmp2
    addiu       t8, zero, 16069  // FIX_1_961570560
    mul         s0, s0, t8       // -z3 = MULTIPLY(z3, FIX_1_961570560)
    addiu       t8, zero, 3196   // FIX_0_390180644
    mul         s1, s1, t8       // -z4 = MULTIPLY(z4, FIX_0_390180644)
    addiu       t8, zero, 2446   // FIX_0_298631336
    mul         t4, t4, t8       // tmp0 = MULTIPLY(tmp0, FIX_0_298631336)
    addiu       t8, zero, 7373   // FIX_0_899976223
    mul         s3, s3, t8       // -z1 = MULTIPLY(z1, FIX_0_899976223)
    addiu       t8, zero, 16819  // FIX_2_053119869
    mul         t5, t5, t8       // tmp1 = MULTIPLY(tmp1, FIX_2_053119869)
    addiu       t8, zero, 20995  // FIX_2_562915447
    mul         t9, t9, t8       // -z2 = MULTIPLY(z2, FIX_2_562915447)
    addiu       t8, zero, 25172  // FIX_3_072711026
    mul         t6, t6, t8       // tmp2 = MULTIPLY(tmp2, FIX_3_072711026)
    addiu       t8, zero, 12299  // FIX_1_501321110
    mul         t7, t7, t8       // tmp3 = MULTIPLY(tmp3, FIX_1_501321110)
    subu        s0, s2, s0       // z3 += z5
    subu        s1, s2, s1       // z4 += z5
    addu        t4, t4, s0
    subu        t4, t4, s3      // tmp0
    addu        t5, t5, s1
    subu        t5, t5, t9      // tmp1
    addu        t6, t6, s0
    subu        t6, t6, t9      // tmp2
    addu        t7, t7, s1
    subu        t7, t7, s3      // tmp3
    addu        s0, t0, t7
    subu        t0, t0, t7
    addu        t7, t2, t6
    subu        t2, t2, t6
    addu        t6, t3, t5
    subu        t3, t3, t5
    addu        t5, t1, t4
    subu        t1, t1, t4
    shra_r.w    s0, s0, 18
    shra_r.w    t7, t7, 18
    shra_r.w    t6, t6, 18
    shra_r.w    t5, t5, 18
    shra_r.w    t1, t1, 18
    shra_r.w    t3, t3, 18
    shra_r.w    t2, t2, 18
    shra_r.w    t0, t0, 18
    andi        s0, s0, 0x3ff
    andi        t7, t7, 0x3ff
    andi        t6, t6, 0x3ff
    andi        t5, t5, 0x3ff
    andi        t1, t1, 0x3ff
    andi        t3, t3, 0x3ff
    andi        t2, t2, 0x3ff
    andi        t0, t0, 0x3ff
    lw          s1, 0(a2)
    lbux        s0, s0(a3)
    lbux        t7, t7(a3)
    lbux        t6, t6(a3)
    lbux        t5, t5(a3)
    lbux        t1, t1(a3)
    lbux        t3, t3(a3)
    lbux        t2, t2(a3)
    lbux        t0, t0(a3)
    sb          s0, 0(s1)
    sb          t7, 1(s1)
    sb          t6, 2(s1)
    sb          t5, 3(s1)
    sb          t1, 4(s1)
    sb          t3, 5(s1)
    sb          t2, 6(s1)
    sb          t0, 7(s1)
6:
    addiu       v0, v0, 32
    bgtz        v1, 4b
     addiu      a2, a2, 4
    addiu       sp, sp, 256

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop

END(jsimd_idct_islow_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_ifast_cols_dspr2)
/*
 * a0 = inptr
 * a1 = quantptr
 * a2 = wsptr
 * a3 = mips_idct_ifast_coefs
 */
    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    addiu         t9, a0, 16      // end address
    or            AT, a3, zero

0:
    lw            s0, 0(a1)       // quantptr[DCTSIZE*0]
    lw            t0, 0(a0)       // inptr[DCTSIZE*0]
    lw            t1, 16(a0)      // inptr[DCTSIZE*1]
    muleq_s.w.phl v0, t0, s0      // tmp0 ...
    lw            t2, 32(a0)      // inptr[DCTSIZE*2]
    lw            t3, 48(a0)      // inptr[DCTSIZE*3]
    lw            t4, 64(a0)      // inptr[DCTSIZE*4]
    lw            t5, 80(a0)      // inptr[DCTSIZE*5]
    muleq_s.w.phr t0, t0, s0      // ... tmp0 ...
    lw            t6, 96(a0)      // inptr[DCTSIZE*6]
    lw            t7, 112(a0)     // inptr[DCTSIZE*7]
    or            s4, t1, t2
    or            s5, t3, t4
    bnez          s4, 1f
     ins          t0, v0, 16, 16  // ... tmp0
    bnez          s5, 1f
     or           s6, t5, t6
    or            s6, s6, t7
    bnez          s6, 1f
     sw           t0, 0(a2)       // wsptr[DCTSIZE*0]
    sw            t0, 16(a2)      // wsptr[DCTSIZE*1]
    sw            t0, 32(a2)      // wsptr[DCTSIZE*2]
    sw            t0, 48(a2)      // wsptr[DCTSIZE*3]
    sw            t0, 64(a2)      // wsptr[DCTSIZE*4]
    sw            t0, 80(a2)      // wsptr[DCTSIZE*5]
    sw            t0, 96(a2)      // wsptr[DCTSIZE*6]
    sw            t0, 112(a2)     // wsptr[DCTSIZE*7]
    addiu         a0, a0, 4
    b             2f
     addiu        a1, a1, 4

1:
    lw            s1, 32(a1)      // quantptr[DCTSIZE*2]
    lw            s2, 64(a1)      // quantptr[DCTSIZE*4]
    muleq_s.w.phl v0, t2, s1      // tmp1 ...
    muleq_s.w.phr t2, t2, s1      // ... tmp1 ...
    lw            s0, 16(a1)      // quantptr[DCTSIZE*1]
    lw            s1, 48(a1)      // quantptr[DCTSIZE*3]
    lw            s3, 96(a1)      // quantptr[DCTSIZE*6]
    muleq_s.w.phl v1, t4, s2      // tmp2 ...
    muleq_s.w.phr t4, t4, s2      // ... tmp2 ...
    lw            s2, 80(a1)      // quantptr[DCTSIZE*5]
    lw            t8, 4(AT)       // FIX(1.414213562)
    ins           t2, v0, 16, 16  // ... tmp1
    muleq_s.w.phl v0, t6, s3      // tmp3 ...
    muleq_s.w.phr t6, t6, s3      // ... tmp3 ...
    ins           t4, v1, 16, 16  // ... tmp2
    addq.ph       s4, t0, t4      // tmp10
    subq.ph       s5, t0, t4      // tmp11
    ins           t6, v0, 16, 16  // ... tmp3
    subq.ph       s6, t2, t6      // tmp12 ...
    addq.ph       s7, t2, t6      // tmp13
    mulq_s.ph     s6, s6, t8      // ... tmp12 ...
    addq.ph       t0, s4, s7      // tmp0
    subq.ph       t6, s4, s7      // tmp3
    muleq_s.w.phl v0, t1, s0      // tmp4 ...
    muleq_s.w.phr t1, t1, s0      // ... tmp4 ...
    shll_s.ph     s6, s6, 1       // x2
    lw            s3, 112(a1)     // quantptr[DCTSIZE*7]
    subq.ph       s6, s6, s7      // ... tmp12
    muleq_s.w.phl v1, t7, s3      // tmp7 ...
    muleq_s.w.phr t7, t7, s3      // ... tmp7 ...
    ins           t1, v0, 16, 16  // ... tmp4
    addq.ph       t2, s5, s6      // tmp1
    subq.ph       t4, s5, s6      // tmp2
    muleq_s.w.phl v0, t5, s2      // tmp6 ...
    muleq_s.w.phr t5, t5, s2      // ... tmp6 ...
    ins           t7, v1, 16, 16  // ... tmp7
    addq.ph       s5, t1, t7      // z11
    subq.ph       s6, t1, t7      // z12
    muleq_s.w.phl v1, t3, s1      // tmp5 ...
    muleq_s.w.phr t3, t3, s1      // ... tmp5 ...
    ins           t5, v0, 16, 16  // ... tmp6
    ins           t3, v1, 16, 16  // ... tmp5
    addq.ph       s7, t5, t3      // z13
    subq.ph       v0, t5, t3      // z10
    addq.ph       t7, s5, s7      // tmp7
    subq.ph       s5, s5, s7      // tmp11 ...
    addq.ph       v1, v0, s6      // z5 ...
    mulq_s.ph     s5, s5, t8      // ... tmp11
    lw            t8, 8(AT)       // FIX(1.847759065)
    lw            s4, 0(AT)       // FIX(1.082392200)
    addq.ph       s0, t0, t7
    subq.ph       s1, t0, t7
    mulq_s.ph     v1, v1, t8      // ... z5
    shll_s.ph     s5, s5, 1       // x2
    lw            t8, 12(AT)      // FIX(-2.613125930)
    sw            s0, 0(a2)       // wsptr[DCTSIZE*0]
    shll_s.ph     v0, v0, 1       // x4
    mulq_s.ph     v0, v0, t8      // tmp12 ...
    mulq_s.ph     s4, s6, s4      // tmp10 ...
    shll_s.ph     v1, v1, 1       // x2
    addiu         a0, a0, 4
    addiu         a1, a1, 4
    sw            s1, 112(a2)     // wsptr[DCTSIZE*7]
    shll_s.ph     s6, v0, 1       // x4
    shll_s.ph     s4, s4, 1       // x2
    addq.ph       s6, s6, v1      // ... tmp12
    subq.ph       t5, s6, t7      // tmp6
    subq.ph       s4, s4, v1      // ... tmp10
    subq.ph       t3, s5, t5      // tmp5
    addq.ph       s2, t2, t5
    addq.ph       t1, s4, t3      // tmp4
    subq.ph       s3, t2, t5
    sw            s2, 16(a2)      // wsptr[DCTSIZE*1]
    sw            s3, 96(a2)      // wsptr[DCTSIZE*6]
    addq.ph       v0, t4, t3
    subq.ph       v1, t4, t3
    sw            v0, 32(a2)      // wsptr[DCTSIZE*2]
    sw            v1, 80(a2)      // wsptr[DCTSIZE*5]
    addq.ph       v0, t6, t1
    subq.ph       v1, t6, t1
    sw            v0, 64(a2)      // wsptr[DCTSIZE*4]
    sw            v1, 48(a2)      // wsptr[DCTSIZE*3]

2:
    bne           a0, t9, 0b
     addiu        a2, a2, 4

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j             ra
     nop

END(jsimd_idct_ifast_cols_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_ifast_rows_dspr2)
/*
 * a0 = wsptr
 * a1 = output_buf
 * a2 = output_col
 * a3 = mips_idct_ifast_coefs
 */
    SAVE_REGS_ON_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, s8, a3

    addiu         t9, a0, 128     // end address
    lui           s8, 0x8080
    ori           s8, s8, 0x8080

0:
    lw            AT, 36(sp)      // restore $a3 (mips_idct_ifast_coefs)
    lw            t0, 0(a0)       // wsptr[DCTSIZE*0+0/1]  b a
    lw            s0, 16(a0)      // wsptr[DCTSIZE*1+0/1]  B A
    lw            t2, 4(a0)       // wsptr[DCTSIZE*0+2/3]  d c
    lw            s2, 20(a0)      // wsptr[DCTSIZE*1+2/3]  D C
    lw            t4, 8(a0)       // wsptr[DCTSIZE*0+4/5]  f e
    lw            s4, 24(a0)      // wsptr[DCTSIZE*1+4/5]  F E
    lw            t6, 12(a0)      // wsptr[DCTSIZE*0+6/7]  h g
    lw            s6, 28(a0)      // wsptr[DCTSIZE*1+6/7]  H G
    precrq.ph.w   t1, s0, t0      // B b
    ins           t0, s0, 16, 16  // A a
    bnez          t1, 1f
     or           s0, t2, s2
    bnez          s0, 1f
     or           s0, t4, s4
    bnez          s0, 1f
     or           s0, t6, s6
    bnez          s0, 1f
     shll_s.ph    s0, t0, 2       // A a
    lw            a3, 0(a1)
    lw            AT, 4(a1)
    precrq.ph.w   t0, s0, s0      // A A
    ins           s0, s0, 16, 16  // a a
    addu          a3, a3, a2
    addu          AT, AT, a2
    precrq.qb.ph  t0, t0, t0      // A A A A
    precrq.qb.ph  s0, s0, s0      // a a a a
    addu.qb       s0, s0, s8
    addu.qb       t0, t0, s8
    sw            s0, 0(a3)
    sw            s0, 4(a3)
    sw            t0, 0(AT)
    sw            t0, 4(AT)
    addiu         a0, a0, 32
    bne           a0, t9, 0b
     addiu        a1, a1, 8
    b             2f
     nop

1:
    precrq.ph.w   t3, s2, t2
    ins           t2, s2, 16, 16
    precrq.ph.w   t5, s4, t4
    ins           t4, s4, 16, 16
    precrq.ph.w   t7, s6, t6
    ins           t6, s6, 16, 16
    lw            t8, 4(AT)       // FIX(1.414213562)
    addq.ph       s4, t0, t4      // tmp10
    subq.ph       s5, t0, t4      // tmp11
    subq.ph       s6, t2, t6      // tmp12 ...
    addq.ph       s7, t2, t6      // tmp13
    mulq_s.ph     s6, s6, t8      // ... tmp12 ...
    addq.ph       t0, s4, s7      // tmp0
    subq.ph       t6, s4, s7      // tmp3
    shll_s.ph     s6, s6, 1       // x2
    subq.ph       s6, s6, s7      // ... tmp12
    addq.ph       t2, s5, s6      // tmp1
    subq.ph       t4, s5, s6      // tmp2
    addq.ph       s5, t1, t7      // z11
    subq.ph       s6, t1, t7      // z12
    addq.ph       s7, t5, t3      // z13
    subq.ph       v0, t5, t3      // z10
    addq.ph       t7, s5, s7      // tmp7
    subq.ph       s5, s5, s7      // tmp11 ...
    addq.ph       v1, v0, s6      // z5 ...
    mulq_s.ph     s5, s5, t8      // ... tmp11
    lw            t8, 8(AT)       // FIX(1.847759065)
    lw            s4, 0(AT)       // FIX(1.082392200)
    addq.ph       s0, t0, t7      // tmp0 + tmp7
    subq.ph       s7, t0, t7      // tmp0 - tmp7
    mulq_s.ph     v1, v1, t8      // ... z5
    lw            a3, 0(a1)
    lw            t8, 12(AT)      // FIX(-2.613125930)
    shll_s.ph     s5, s5, 1       // x2
    addu          a3, a3, a2
    shll_s.ph     v0, v0, 1       // x4
    mulq_s.ph     v0, v0, t8      // tmp12 ...
    mulq_s.ph     s4, s6, s4      // tmp10 ...
    shll_s.ph     v1, v1, 1       // x2
    addiu         a0, a0, 32
    addiu         a1, a1, 8
    shll_s.ph     s6, v0, 1       // x4
    shll_s.ph     s4, s4, 1       // x2
    addq.ph       s6, s6, v1      // ... tmp12
    shll_s.ph     s0, s0, 2
    subq.ph       t5, s6, t7      // tmp6
    subq.ph       s4, s4, v1      // ... tmp10
    subq.ph       t3, s5, t5      // tmp5
    shll_s.ph     s7, s7, 2
    addq.ph       t1, s4, t3      // tmp4
    addq.ph       s1, t2, t5      // tmp1 + tmp6
    subq.ph       s6, t2, t5      // tmp1 - tmp6
    addq.ph       s2, t4, t3      // tmp2 + tmp5
    subq.ph       s5, t4, t3      // tmp2 - tmp5
    addq.ph       s4, t6, t1      // tmp3 + tmp4
    subq.ph       s3, t6, t1      // tmp3 - tmp4
    shll_s.ph     s1, s1, 2
    shll_s.ph     s2, s2, 2
    shll_s.ph     s3, s3, 2
    shll_s.ph     s4, s4, 2
    shll_s.ph     s5, s5, 2
    shll_s.ph     s6, s6, 2
    precrq.ph.w   t0, s1, s0      // B A
    ins           s0, s1, 16, 16  // b a
    precrq.ph.w   t2, s3, s2      // D C
    ins           s2, s3, 16, 16  // d c
    precrq.ph.w   t4, s5, s4      // F E
    ins           s4, s5, 16, 16  // f e
    precrq.ph.w   t6, s7, s6      // H G
    ins           s6, s7, 16, 16  // h g
    precrq.qb.ph  t0, t2, t0      // D C B A
    precrq.qb.ph  s0, s2, s0      // d c b a
    precrq.qb.ph  t4, t6, t4      // H G F E
    precrq.qb.ph  s4, s6, s4      // h g f e
    addu.qb       s0, s0, s8
    addu.qb       s4, s4, s8
    sw            s0, 0(a3)       // outptr[0/1/2/3]       d c b a
    sw            s4, 4(a3)       // outptr[4/5/6/7]       h g f e
    lw            a3, -4(a1)
    addu.qb       t0, t0, s8
    addu          a3, a3, a2
    addu.qb       t4, t4, s8
    sw            t0, 0(a3)       // outptr[0/1/2/3]       D C B A
    bne           a0, t9, 0b
     sw           t4, 4(a3)       // outptr[4/5/6/7]       H G F E

2:

    RESTORE_REGS_FROM_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, s8, a3

    j             ra
     nop

END(jsimd_idct_ifast_rows_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_fdct_islow_dspr2)
/*
 * a0 = data
 */
    SAVE_REGS_ON_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, s8

    lui         t0, 6437
    ori         t0, 2260
    lui         t1, 9633
    ori         t1, 11363
    lui         t2, 0xd39e
    ori         t2, 0xe6dc
    lui         t3, 0xf72d
    ori         t3, 9633
    lui         t4, 2261
    ori         t4, 9633
    lui         t5, 0xd39e
    ori         t5, 6437
    lui         t6, 9633
    ori         t6, 0xd39d
    lui         t7, 0xe6dc
    ori         t7, 2260
    lui         t8, 4433
    ori         t8, 10703
    lui         t9, 0xd630
    ori         t9, 4433
    li          s8, 8
    move        a1, a0
1:
    lw          s0, 0(a1)       // tmp0 = 1|0
    lw          s1, 4(a1)       // tmp1 = 3|2
    lw          s2, 8(a1)       // tmp2 = 5|4
    lw          s3, 12(a1)      // tmp3 = 7|6
    packrl.ph   s1, s1, s1      // tmp1 = 2|3
    packrl.ph   s3, s3, s3      // tmp3 = 6|7
    subq.ph     s7, s1, s2      // tmp7 = 2-5|3-4 = t5|t4
    subq.ph     s5, s0, s3      // tmp5 = 1-6|0-7 = t6|t7
    mult        $0, $0          // ac0  = 0
    dpa.w.ph    $ac0, s7, t0    // ac0 += t5*  6437 + t4*  2260
    dpa.w.ph    $ac0, s5, t1    // ac0 += t6*  9633 + t7* 11363
    mult        $ac1, $0, $0    // ac1  = 0
    dpa.w.ph    $ac1, s7, t2    // ac1 += t5*-11362 + t4* -6436
    dpa.w.ph    $ac1, s5, t3    // ac1 += t6* -2259 + t7*  9633
    mult        $ac2, $0, $0    // ac2  = 0
    dpa.w.ph    $ac2, s7, t4    // ac2 += t5*  2261 + t4*  9633
    dpa.w.ph    $ac2, s5, t5    // ac2 += t6*-11362 + t7*  6437
    mult        $ac3, $0, $0    // ac3  = 0
    dpa.w.ph    $ac3, s7, t6    // ac3 += t5*  9633 + t4*-11363
    dpa.w.ph    $ac3, s5, t7    // ac3 += t6* -6436 + t7*  2260
    addq.ph     s6, s1, s2      // tmp6 = 2+5|3+4 = t2|t3
    addq.ph     s4, s0, s3      // tmp4 = 1+6|0+7 = t1|t0
    extr_r.w    s0, $ac0, 11    // tmp0 = (ac0 + 1024) >> 11
    extr_r.w    s1, $ac1, 11    // tmp1 = (ac1 + 1024) >> 11
    extr_r.w    s2, $ac2, 11    // tmp2 = (ac2 + 1024) >> 11
    extr_r.w    s3, $ac3, 11    // tmp3 = (ac3 + 1024) >> 11
    addq.ph     s5, s4, s6      // tmp5 = t1+t2|t0+t3 = t11|t10
    subq.ph     s7, s4, s6      // tmp7 = t1-t2|t0-t3 = t12|t13
    sh          s0, 2(a1)
    sh          s1, 6(a1)
    sh          s2, 10(a1)
    sh          s3, 14(a1)
    mult        $0, $0          // ac0  = 0
    dpa.w.ph    $ac0, s7, t8    // ac0 += t12*  4433 + t13* 10703
    mult        $ac1, $0, $0    // ac1  = 0
    dpa.w.ph    $ac1, s7, t9    // ac1 += t12*-10704 + t13*  4433
    sra         s4, s5, 16      // tmp4 = t11
    addiu       a1, a1, 16
    addiu       s8, s8, -1
    extr_r.w    s0, $ac0, 11    // tmp0 = (ac0 + 1024) >> 11
    extr_r.w    s1, $ac1, 11    // tmp1 = (ac1 + 1024) >> 11
    addu        s2, s5, s4      // tmp2 = t10 + t11
    subu        s3, s5, s4      // tmp3 = t10 - t11
    sll         s2, s2, 2       // tmp2 = (t10 + t11) << 2
    sll         s3, s3, 2       // tmp3 = (t10 - t11) << 2
    sh          s2, -16(a1)
    sh          s3, -8(a1)
    sh          s0, -12(a1)
    bgtz        s8, 1b
     sh         s1, -4(a1)
    li          t0, 2260
    li          t1, 11363
    li          t2, 9633
    li          t3, 6436
    li          t4, 6437
    li          t5, 2261
    li          t6, 11362
    li          t7, 2259
    li          t8, 4433
    li          t9, 10703
    li          a1, 10704
    li          s8, 8

2:
    lh          a2, 0(a0)       // 0
    lh          a3, 16(a0)      // 8
    lh          v0, 32(a0)      // 16
    lh          v1, 48(a0)      // 24
    lh          s4, 64(a0)      // 32
    lh          s5, 80(a0)      // 40
    lh          s6, 96(a0)      // 48
    lh          s7, 112(a0)     // 56
    addu        s2, v0, s5      // tmp2 = 16 + 40
    subu        s5, v0, s5      // tmp5 = 16 - 40
    addu        s3, v1, s4      // tmp3 = 24 + 32
    subu        s4, v1, s4      // tmp4 = 24 - 32
    addu        s0, a2, s7      // tmp0 =  0 + 56
    subu        s7, a2, s7      // tmp7 =  0 - 56
    addu        s1, a3, s6      // tmp1 =  8 + 48
    subu        s6, a3, s6      // tmp6 =  8 - 48
    addu        a2, s0, s3      // tmp10 = tmp0 + tmp3
    subu        v1, s0, s3      // tmp13 = tmp0 - tmp3
    addu        a3, s1, s2      // tmp11 = tmp1 + tmp2
    subu        v0, s1, s2      // tmp12 = tmp1 - tmp2
    mult        s7, t1          // ac0  = tmp7 * c1
    madd        s4, t0          // ac0 += tmp4 * c0
    madd        s5, t4          // ac0 += tmp5 * c4
    madd        s6, t2          // ac0 += tmp6 * c2
    mult        $ac1, s7, t2    // ac1  = tmp7 * c2
    msub        $ac1, s4, t3    // ac1 -= tmp4 * c3
    msub        $ac1, s5, t6    // ac1 -= tmp5 * c6
    msub        $ac1, s6, t7    // ac1 -= tmp6 * c7
    mult        $ac2, s7, t4    // ac2  = tmp7 * c4
    madd        $ac2, s4, t2    // ac2 += tmp4 * c2
    madd        $ac2, s5, t5    // ac2 += tmp5 * c5
    msub        $ac2, s6, t6    // ac2 -= tmp6 * c6
    mult        $ac3, s7, t0    // ac3  = tmp7 * c0
    msub        $ac3, s4, t1    // ac3 -= tmp4 * c1
    madd        $ac3, s5, t2    // ac3 += tmp5 * c2
    msub        $ac3, s6, t3    // ac3 -= tmp6 * c3
    extr_r.w    s0, $ac0, 15    // tmp0 = (ac0 + 16384) >> 15
    extr_r.w    s1, $ac1, 15    // tmp1 = (ac1 + 16384) >> 15
    extr_r.w    s2, $ac2, 15    // tmp2 = (ac2 + 16384) >> 15
    extr_r.w    s3, $ac3, 15    // tmp3 = (ac3 + 16384) >> 15
    addiu       s8, s8, -1
    addu        s4, a2, a3      // tmp4 = tmp10 + tmp11
    subu        s5, a2, a3      // tmp5 = tmp10 - tmp11
    sh          s0, 16(a0)
    sh          s1, 48(a0)
    sh          s2, 80(a0)
    sh          s3, 112(a0)
    mult        v0, t8          // ac0  = tmp12 * c8
    madd        v1, t9          // ac0 += tmp13 * c9
    mult        $ac1, v1, t8    // ac1  = tmp13 * c8
    msub        $ac1, v0, a1    // ac1 -= tmp12 * c10
    addiu       a0, a0, 2
    extr_r.w    s6, $ac0, 15    // tmp6 = (ac0 + 16384) >> 15
    extr_r.w    s7, $ac1, 15    // tmp7 = (ac1 + 16384) >> 15
    shra_r.w    s4, s4, 2       // tmp4 = (tmp4 + 2) >> 2
    shra_r.w    s5, s5, 2       // tmp5 = (tmp5 + 2) >> 2
    sh          s4, -2(a0)
    sh          s5, 62(a0)
    sh          s6, 30(a0)
    bgtz        s8, 2b
     sh         s7, 94(a0)

    RESTORE_REGS_FROM_STACK 40, s0, s1, s2, s3, s4, s5, s6, s7, s8

    jr          ra
     nop

END(jsimd_fdct_islow_dspr2)


/**************************************************************************/
LEAF_DSPR2(jsimd_fdct_ifast_dspr2)
/*
 * a0 = data
 */
    .set at

    SAVE_REGS_ON_STACK 8, s0, s1

    li          a1, 0x014e014e  // FIX_1_306562965 (334 << 16)|(334 & 0xffff)
    li          a2, 0x008b008b  // FIX_0_541196100 (139 << 16)|(139 & 0xffff)
    li          a3, 0x00620062  // FIX_0_382683433 (98 << 16) |(98 & 0xffff)
    li          s1, 0x00b500b5  // FIX_0_707106781 (181 << 16)|(181 & 0xffff)

    move        v0, a0
    addiu       v1, v0, 128     // end address

0:
    lw          t0, 0(v0)       // tmp0 = 1|0
    lw          t1, 4(v0)       // tmp1 = 3|2
    lw          t2, 8(v0)       // tmp2 = 5|4
    lw          t3, 12(v0)      // tmp3 = 7|6
    packrl.ph   t1, t1, t1      // tmp1 = 2|3
    packrl.ph   t3, t3, t3      // tmp3 = 6|7
    subq.ph     t7, t1, t2      // tmp7 = 2-5|3-4 = t5|t4
    subq.ph     t5, t0, t3      // tmp5 = 1-6|0-7 = t6|t7
    addq.ph     t6, t1, t2      // tmp6 = 2+5|3+4 = t2|t3
    addq.ph     t4, t0, t3      // tmp4 = 1+6|0+7 = t1|t0
    addq.ph     t8, t4, t6      // tmp5 = t1+t2|t0+t3 = t11|t10
    subq.ph     t9, t4, t6      // tmp7 = t1-t2|t0-t3 = t12|t13
    sra         t4, t8, 16      // tmp4 = t11
    mult        $0, $0          // ac0  = 0
    dpa.w.ph    $ac0, t9, s1
    mult        $ac1, $0, $0    // ac1  = 0
    dpa.w.ph    $ac1, t7, a3    // ac1 += t4*98 + t5*98
    dpsx.w.ph   $ac1, t5, a3    // ac1 += t6*98 + t7*98
    mult        $ac2, $0, $0    // ac2  = 0
    dpa.w.ph    $ac2, t7, a2    // ac2 += t4*139 + t5*139
    mult        $ac3, $0, $0    // ac3  = 0
    dpa.w.ph    $ac3, t5, a1    // ac3 += t6*334 + t7*334
    precrq.ph.w t0, t5, t7      // t0 = t5|t6
    addq.ph     t2, t8, t4      // tmp2 = t10 + t11
    subq.ph     t3, t8, t4      // tmp3 = t10 - t11
    extr.w      t4, $ac0, 8
    mult        $0, $0          // ac0  = 0
    dpa.w.ph    $ac0, t0, s1    // ac0 += t5*181 + t6*181
    extr.w      t0, $ac1, 8     // t0 = z5
    extr.w      t1, $ac2, 8     // t1 = MULTIPLY(tmp10, 139)
    extr.w      t7, $ac3, 8     // t2 = MULTIPLY(tmp12, 334)
    extr.w      t8, $ac0, 8     // t8 = z3 = MULTIPLY(tmp11, 181)
    add         t6, t1, t0      // t6 = z2
    add         t7, t7, t0      // t7 = z4
    subq.ph     t0, t5, t8      // t0 = z13 = tmp7 - z3
    addq.ph     t8, t5, t8      // t9 = z11 = tmp7 + z3
    addq.ph     t1, t0, t6      // t1 = z13 + z2
    subq.ph     t6, t0, t6      // t6 = z13 - z2
    addq.ph     t0, t8, t7      // t0 = z11 + z4
    subq.ph     t7, t8, t7      // t7 = z11 - z4
    addq.ph     t5, t4, t9
    subq.ph     t4, t9, t4
    sh          t2, 0(v0)
    sh          t5, 4(v0)
    sh          t3, 8(v0)
    sh          t4, 12(v0)
    sh          t1, 10(v0)
    sh          t6, 6(v0)
    sh          t0, 2(v0)
    sh          t7, 14(v0)
    addiu       v0, 16
    bne         v1, v0, 0b
     nop
    move        v0, a0
    addiu       v1, v0, 16

1:
    lh          t0, 0(v0)       // 0
    lh          t1, 16(v0)      // 8
    lh          t2, 32(v0)      // 16
    lh          t3, 48(v0)      // 24
    lh          t4, 64(v0)      // 32
    lh          t5, 80(v0)      // 40
    lh          t6, 96(v0)      // 48
    lh          t7, 112(v0)     // 56
    add         t8, t0, t7      // t8 = tmp0
    sub         t7, t0, t7      // t7 = tmp7
    add         t0, t1, t6      // t0 = tmp1
    sub         t1, t1, t6      // t1 = tmp6
    add         t6, t2, t5      // t6 = tmp2
    sub         t5, t2, t5      // t5 = tmp5
    add         t2, t3, t4      // t2 = tmp3
    sub         t3, t3, t4      // t3 = tmp4
    add         t4, t8, t2      // t4 = tmp10 = tmp0 + tmp3
    sub         t8, t8, t2      // t8 = tmp13 = tmp0 - tmp3
    sub         s0, t0, t6      // s0 = tmp12 = tmp1 - tmp2
    ins         t8, s0, 16, 16  // t8 = tmp12|tmp13
    add         t2, t0, t6      // t2 = tmp11 = tmp1 + tmp2
    mult        $0, $0          // ac0  = 0
    dpa.w.ph    $ac0, t8, s1    // ac0 += t12*181 + t13*181
    add         s0, t4, t2      // t8 = tmp10+tmp11
    sub         t4, t4, t2      // t4 = tmp10-tmp11
    sh          s0, 0(v0)
    sh          t4, 64(v0)
    extr.w      t2, $ac0, 8     // z1 = MULTIPLY(tmp12+tmp13, FIX_0_707106781)
    addq.ph     t4, t8, t2      // t9 = tmp13 + z1
    subq.ph     t8, t8, t2      // t2 = tmp13 - z1
    sh          t4, 32(v0)
    sh          t8, 96(v0)
    add         t3, t3, t5      // t3 = tmp10 = tmp4 + tmp5
    add         t0, t5, t1      // t0 = tmp11 = tmp5 + tmp6
    add         t1, t1, t7      // t1 = tmp12 = tmp6 + tmp7
    andi        t4, a1, 0xffff
    mul         s0, t1, t4
    sra         s0, s0, 8       // s0 = z4 = MULTIPLY(tmp12, FIX_1_306562965)
    ins         t1, t3, 16, 16  // t1 = tmp10|tmp12
    mult        $0, $0          // ac0  = 0
    mulsa.w.ph  $ac0, t1, a3    // ac0 += t10*98 - t12*98
    extr.w      t8, $ac0, 8     // z5 = MULTIPLY(tmp10-tmp12, FIX_0_382683433)
    add         t2, t7, t8      // t2 = tmp7 + z5
    sub         t7, t7, t8      // t7 = tmp7 - z5
    andi        t4, a2, 0xffff
    mul         t8, t3, t4
    sra         t8, t8, 8       // t8 = z2 = MULTIPLY(tmp10, FIX_0_541196100)
    andi        t4, s1, 0xffff
    mul         t6, t0, t4
    sra         t6, t6, 8       // t6 = z3 = MULTIPLY(tmp11, FIX_0_707106781)
    add         t0, t6, t8      // t0 = z3 + z2
    sub         t1, t6, t8      // t1 = z3 - z2
    add         t3, t6, s0      // t3 = z3 + z4
    sub         t4, t6, s0      // t4 = z3 - z4
    sub         t5, t2, t1      // t5 = dataptr[5]
    sub         t6, t7, t0      // t6 = dataptr[3]
    add         t3, t2, t3      // t3 = dataptr[1]
    add         t4, t7, t4      // t4 = dataptr[7]
    sh          t5, 80(v0)
    sh          t6, 48(v0)
    sh          t3, 16(v0)
    sh          t4, 112(v0)
    addiu       v0, 2
    bne         v0, v1, 1b
     nop

    RESTORE_REGS_FROM_STACK 8, s0, s1

    j           ra
     nop
END(jsimd_fdct_ifast_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_quantize_dspr2)
/*
 * a0 = coef_block
 * a1 = divisors
 * a2 = workspace
 */
    .set at

    SAVE_REGS_ON_STACK 16, s0, s1, s2

    addiu       v0, a2, 124     // v0 = workspace_end
    lh          t0, 0(a2)
    lh          t1, 0(a1)
    lh          t2, 128(a1)
    sra         t3, t0, 15
    sll         t3, t3, 1
    addiu       t3, t3, 1
    mul         t0, t0, t3
    lh          t4, 384(a1)
    lh          t5, 130(a1)
    lh          t6, 2(a2)
    lh          t7, 2(a1)
    lh          t8, 386(a1)

1:
    andi        t1, 0xffff
    add         t9, t0, t2
    andi        t9, 0xffff
    mul         v1, t9, t1
    sra         s0, t6, 15
    sll         s0, s0, 1
    addiu       s0, s0, 1
    addiu       t9, t4, 16
    srav        v1, v1, t9
    mul         v1, v1, t3
    mul         t6, t6, s0
    andi        t7, 0xffff
    addiu       a2, a2, 4
    addiu       a1, a1, 4
    add         s1, t6, t5
    andi        s1, 0xffff
    sh          v1, 0(a0)

    mul         s2, s1, t7
    addiu       s1, t8, 16
    srav        s2, s2, s1
    mul         s2, s2, s0
    lh          t0, 0(a2)
    lh          t1, 0(a1)
    sra         t3, t0, 15
    sll         t3, t3, 1
    addiu       t3, t3, 1
    mul         t0, t0, t3
    lh          t2, 128(a1)
    lh          t4, 384(a1)
    lh          t5, 130(a1)
    lh          t8, 386(a1)
    lh          t6, 2(a2)
    lh          t7, 2(a1)
    sh          s2, 2(a0)
    lh          t0, 0(a2)
    sra         t3, t0, 15
    sll         t3, t3, 1
    addiu       t3, t3, 1
    mul         t0, t0, t3
    bne         a2, v0, 1b
     addiu      a0, a0, 4

    andi        t1, 0xffff
    add         t9, t0, t2
    andi        t9, 0xffff
    mul         v1, t9, t1
    sra         s0, t6, 15
    sll         s0, s0, 1
    addiu       s0, s0, 1
    addiu       t9, t4, 16
    srav        v1, v1, t9
    mul         v1, v1, t3
    mul         t6, t6, s0
    andi        t7, 0xffff
    sh          v1, 0(a0)
    add         s1, t6, t5
    andi        s1, 0xffff
    mul         s2, s1, t7
    addiu       s1, t8, 16
    addiu       a2, a2, 4
    addiu       a1, a1, 4
    srav        s2, s2, s1
    mul         s2, s2, s0
    sh          s2, 2(a0)

    RESTORE_REGS_FROM_STACK 16, s0, s1, s2

    j           ra
     nop

END(jsimd_quantize_dspr2)


#ifndef __mips_soft_float

/*****************************************************************************/
LEAF_DSPR2(jsimd_quantize_float_dspr2)
/*
 * a0 = coef_block
 * a1 = divisors
 * a2 = workspace
 */
    .set at

    li          t1, 0x46800100  // integer representation 16384.5
    mtc1        t1, f0
    li          t0, 63
0:
    lwc1        f2, 0(a2)
    lwc1        f10, 0(a1)
    lwc1        f4, 4(a2)
    lwc1        f12, 4(a1)
    lwc1        f6, 8(a2)
    lwc1        f14, 8(a1)
    lwc1        f8, 12(a2)
    lwc1        f16, 12(a1)
    madd.s      f2, f0, f2, f10
    madd.s      f4, f0, f4, f12
    madd.s      f6, f0, f6, f14
    madd.s      f8, f0, f8, f16
    lwc1        f10, 16(a1)
    lwc1        f12, 20(a1)
    trunc.w.s   f2, f2
    trunc.w.s   f4, f4
    trunc.w.s   f6, f6
    trunc.w.s   f8, f8
    lwc1        f14, 24(a1)
    lwc1        f16, 28(a1)
    mfc1        t1, f2
    mfc1        t2, f4
    mfc1        t3, f6
    mfc1        t4, f8
    lwc1        f2, 16(a2)
    lwc1        f4, 20(a2)
    lwc1        f6, 24(a2)
    lwc1        f8, 28(a2)
    madd.s      f2, f0, f2, f10
    madd.s      f4, f0, f4, f12
    madd.s      f6, f0, f6, f14
    madd.s      f8, f0, f8, f16
    addiu       t1, t1, -16384
    addiu       t2, t2, -16384
    addiu       t3, t3, -16384
    addiu       t4, t4, -16384
    trunc.w.s   f2, f2
    trunc.w.s   f4, f4
    trunc.w.s   f6, f6
    trunc.w.s   f8, f8
    sh          t1, 0(a0)
    sh          t2, 2(a0)
    sh          t3, 4(a0)
    sh          t4, 6(a0)
    mfc1        t1, f2
    mfc1        t2, f4
    mfc1        t3, f6
    mfc1        t4, f8
    addiu       t0, t0, -8
    addiu       a2, a2, 32
    addiu       a1, a1, 32
    addiu       t1, t1, -16384
    addiu       t2, t2, -16384
    addiu       t3, t3, -16384
    addiu       t4, t4, -16384
    sh          t1, 8(a0)
    sh          t2, 10(a0)
    sh          t3, 12(a0)
    sh          t4, 14(a0)
    bgez        t0, 0b
     addiu      a0, a0, 16

    j           ra
     nop

END(jsimd_quantize_float_dspr2)

#endif


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_2x2_dspr2)
/*
 * a0 = compptr->dct_table
 * a1 = coef_block
 * a2 = output_buf
 * a3 = output_col
 */
    .set at

    SAVE_REGS_ON_STACK 24, s0, s1, s2, s3, s4, s5

    addiu       sp, sp, -40
    move        v0, sp
    addiu       s2, zero, 29692
    addiu       s3, zero, -10426
    addiu       s4, zero, 6967
    addiu       s5, zero, -5906
    lh          t0, 0(a1)       // t0 = inptr[DCTSIZE*0]
    lh          t5, 0(a0)       // t5 = quantptr[DCTSIZE*0]
    lh          t1, 48(a1)      // t1 = inptr[DCTSIZE*3]
    lh          t6, 48(a0)      // t6 = quantptr[DCTSIZE*3]
    mul         t4, t5, t0
    lh          t0, 16(a1)      // t0 = inptr[DCTSIZE*1]
    lh          t5, 16(a0)      // t5 = quantptr[DCTSIZE*1]
    mul         t6, t6, t1
    mul         t5, t5, t0
    lh          t2, 80(a1)      // t2 = inptr[DCTSIZE*5]
    lh          t7, 80(a0)      // t7 = quantptr[DCTSIZE*5]
    lh          t3, 112(a1)     // t3 = inptr[DCTSIZE*7]
    lh          t8, 112(a0)     // t8 = quantptr[DCTSIZE*7]
    mul         t7, t7, t2
    mult        zero, zero
    mul         t8, t8, t3
    li          s0, 0x73FCD746  // s0 = (29692 << 16) | (-10426 & 0xffff)
    li          s1, 0x1B37E8EE  // s1 = (6967 << 16) | (-5906 & 0xffff)
    ins         t6, t5, 16, 16  // t6 = t5|t6
    sll         t4, t4, 15
    dpa.w.ph    $ac0, t6, s0
    lh          t1, 2(a1)
    lh          t6, 2(a0)
    ins         t8, t7, 16, 16  // t8 = t7|t8
    dpa.w.ph    $ac0, t8, s1
    mflo        t0, $ac0
    mul         t5, t6, t1
    lh          t1, 18(a1)
    lh          t6, 18(a0)
    lh          t2, 50(a1)
    lh          t7, 50(a0)
    mul         t6, t6, t1
    subu        t8, t4, t0
    mul         t7, t7, t2
    addu        t0, t4, t0
    shra_r.w    t0, t0, 13
    lh          t1, 82(a1)
    lh          t2, 82(a0)
    lh          t3, 114(a1)
    lh          t4, 114(a0)
    shra_r.w    t8, t8, 13
    mul         t1, t1, t2
    mul         t3, t3, t4
    sw          t0, 0(v0)
    sw          t8, 20(v0)
    sll         t4, t5, 15
    ins         t7, t6, 16, 16
    mult        zero, zero
    dpa.w.ph    $ac0, t7, s0
    ins         t3, t1, 16, 16
    lh          t1, 6(a1)
    lh          t6, 6(a0)
    dpa.w.ph    $ac0, t3, s1
    mflo        t0, $ac0
    mul         t5, t6, t1
    lh          t1, 22(a1)
    lh          t6, 22(a0)
    lh          t2, 54(a1)
    lh          t7, 54(a0)
    mul         t6, t6, t1
    subu        t8, t4, t0
    mul         t7, t7, t2
    addu        t0, t4, t0
    shra_r.w    t0, t0, 13
    lh          t1, 86(a1)
    lh          t2, 86(a0)
    lh          t3, 118(a1)
    lh          t4, 118(a0)
    shra_r.w    t8, t8, 13
    mul         t1, t1, t2
    mul         t3, t3, t4
    sw          t0, 4(v0)
    sw          t8, 24(v0)
    sll         t4, t5, 15
    ins         t7, t6, 16, 16
    mult        zero, zero
    dpa.w.ph    $ac0, t7, s0
    ins         t3, t1, 16, 16
    lh          t1, 10(a1)
    lh          t6, 10(a0)
    dpa.w.ph    $ac0, t3, s1
    mflo        t0, $ac0
    mul         t5, t6, t1
    lh          t1, 26(a1)
    lh          t6, 26(a0)
    lh          t2, 58(a1)
    lh          t7, 58(a0)
    mul         t6, t6, t1
    subu        t8, t4, t0
    mul         t7, t7, t2
    addu        t0, t4, t0
    shra_r.w    t0, t0, 13
    lh          t1, 90(a1)
    lh          t2, 90(a0)
    lh          t3, 122(a1)
    lh          t4, 122(a0)
    shra_r.w    t8, t8, 13
    mul         t1, t1, t2
    mul         t3, t3, t4
    sw          t0, 8(v0)
    sw          t8, 28(v0)
    sll         t4, t5, 15
    ins         t7, t6, 16, 16
    mult        zero, zero
    dpa.w.ph    $ac0, t7, s0
    ins         t3, t1, 16, 16
    lh          t1, 14(a1)
    lh          t6, 14(a0)
    dpa.w.ph    $ac0, t3, s1
    mflo        t0, $ac0
    mul         t5, t6, t1
    lh          t1, 30(a1)
    lh          t6, 30(a0)
    lh          t2, 62(a1)
    lh          t7, 62(a0)
    mul         t6, t6, t1
    subu        t8, t4, t0
    mul         t7, t7, t2
    addu        t0, t4, t0
    shra_r.w    t0, t0, 13
    lh          t1, 94(a1)
    lh          t2, 94(a0)
    lh          t3, 126(a1)
    lh          t4, 126(a0)
    shra_r.w    t8, t8, 13
    mul         t1, t1, t2
    mul         t3, t3, t4
    sw          t0, 12(v0)
    sw          t8, 32(v0)
    sll         t4, t5, 15
    ins         t7, t6, 16, 16
    mult        zero, zero
    dpa.w.ph    $ac0, t7, s0
    ins         t3, t1, 16, 16
    dpa.w.ph    $ac0, t3, s1
    mflo        t0, $ac0
    lw          t9, 0(a2)
    lw          t3, 0(v0)
    lw          t7, 4(v0)
    lw          t1, 8(v0)
    addu        t9, t9, a3
    sll         t3, t3, 15
    subu        t8, t4, t0
    addu        t0, t4, t0
    shra_r.w    t0, t0, 13
    shra_r.w    t8, t8, 13
    sw          t0, 16(v0)
    sw          t8, 36(v0)
    lw          t5, 12(v0)
    lw          t6, 16(v0)
    mult        t7, s2
    madd        t1, s3
    madd        t5, s4
    madd        t6, s5
    lw          t5, 24(v0)
    lw          t7, 28(v0)
    mflo        t0, $ac0
    lw          t8, 32(v0)
    lw          t2, 36(v0)
    mult        $ac1, t5, s2
    madd        $ac1, t7, s3
    madd        $ac1, t8, s4
    madd        $ac1, t2, s5
    addu        t1, t3, t0
    subu        t6, t3, t0
    shra_r.w    t1, t1, 20
    shra_r.w    t6, t6, 20
    mflo        t4, $ac1
    shll_s.w    t1, t1, 24
    shll_s.w    t6, t6, 24
    sra         t1, t1, 24
    sra         t6, t6, 24
    addiu       t1, t1, 128
    addiu       t6, t6, 128
    lw          t0, 20(v0)
    sb          t1, 0(t9)
    sb          t6, 1(t9)
    sll         t0, t0, 15
    lw          t9, 4(a2)
    addu        t1, t0, t4
    subu        t6, t0, t4
    addu        t9, t9, a3
    shra_r.w    t1, t1, 20
    shra_r.w    t6, t6, 20
    shll_s.w    t1, t1, 24
    shll_s.w    t6, t6, 24
    sra         t1, t1, 24
    sra         t6, t6, 24
    addiu       t1, t1, 128
    addiu       t6, t6, 128
    sb          t1, 0(t9)
    sb          t6, 1(t9)
    addiu       sp, sp, 40

    RESTORE_REGS_FROM_STACK 24, s0, s1, s2, s3, s4, s5

    j           ra
     nop

END(jsimd_idct_2x2_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_4x4_dspr2)
/*
 * a0     = compptr->dct_table
 * a1     = coef_block
 * a2     = output_buf
 * a3     = output_col
 * 16(sp) = workspace[DCTSIZE*4];  // buffers data between passes
 */
    .set at

    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    lw          v1, 48(sp)
    move        t0, a1
    move        t1, v1
    li          t9, 4
    li          s0, 0x2e75f93e
    li          s1, 0x21f9ba79
    li          s2, 0xecc2efb0
    li          s3, 0x52031ccd

0:
    lh          s6, 32(t0)      // inptr[DCTSIZE*2]
    lh          t6, 32(a0)      // quantptr[DCTSIZE*2]
    lh          s7, 96(t0)      // inptr[DCTSIZE*6]
    lh          t7, 96(a0)      // quantptr[DCTSIZE*6]
    mul         t6, s6, t6      // z2 = (inptr[DCTSIZE*2] * quantptr[DCTSIZE*2])
    lh          s4, 0(t0)       // inptr[DCTSIZE*0]
    mul         t7, s7, t7      // z3 = (inptr[DCTSIZE*6] * quantptr[DCTSIZE*6])
    lh          s5, 0(a0)       // quantptr[0]
    li          s6, 15137
    li          s7, 6270
    mul         t2, s4, s5      // tmp0 = (inptr[0] * quantptr[0])
    mul         t6, s6, t6      // z2 = (inptr[DCTSIZE*2] * quantptr[DCTSIZE*2])
    lh          t5, 112(t0)     // inptr[DCTSIZE*7]
    mul         t7, s7, t7      // z3 = (inptr[DCTSIZE*6] * quantptr[DCTSIZE*6])
    lh          s4, 112(a0)     // quantptr[DCTSIZE*7]
    lh          v0, 80(t0)      // inptr[DCTSIZE*5]
    lh          s5, 80(a0)      // quantptr[DCTSIZE*5]
    lh          s6, 48(a0)      // quantptr[DCTSIZE*3]
    sll         t2, t2, 14      // tmp0 <<= (CONST_BITS+1)
    lh          s7, 16(a0)      // quantptr[DCTSIZE*1]
    lh          t8, 16(t0)      // inptr[DCTSIZE*1]
    subu        t6, t6, t7      // tmp2 = MULTIPLY(z2, t5) - MULTIPLY(z3, t6)
    lh          t7, 48(t0)      // inptr[DCTSIZE*3]
    mul         t5, s4, t5      // z1 = (inptr[DCTSIZE*7] * quantptr[DCTSIZE*7])
    mul         v0, s5, v0      // z2 = (inptr[DCTSIZE*5] * quantptr[DCTSIZE*5])
    mul         t7, s6, t7      // z3 = (inptr[DCTSIZE*3] * quantptr[DCTSIZE*3])
    mul         t8, s7, t8      // z4 = (inptr[DCTSIZE*1] * quantptr[DCTSIZE*1])
    addu        t3, t2, t6      // tmp10 = tmp0 + z2
    subu        t4, t2, t6      // tmp10 = tmp0 - z2
    mult        $ac0, zero, zero
    mult        $ac1, zero, zero
    ins         t5, v0, 16, 16
    ins         t7, t8, 16, 16
    addiu       t9, t9, -1
    dpa.w.ph    $ac0, t5, s0
    dpa.w.ph    $ac0, t7, s1
    dpa.w.ph    $ac1, t5, s2
    dpa.w.ph    $ac1, t7, s3
    mflo        s4, $ac0
    mflo        s5, $ac1
    addiu       a0, a0, 2
    addiu       t1, t1, 4
    addiu       t0, t0, 2
    addu        t6, t4, s4
    subu        t5, t4, s4
    addu        s6, t3, s5
    subu        s7, t3, s5
    shra_r.w    t6, t6, 12      // DESCALE(tmp12 + temp1, 12)
    shra_r.w    t5, t5, 12      // DESCALE(tmp12 - temp1, 12)
    shra_r.w    s6, s6, 12      // DESCALE(tmp10 + temp2, 12)
    shra_r.w    s7, s7, 12      // DESCALE(tmp10 - temp2, 12)
    sw          t6, 28(t1)
    sw          t5, 60(t1)
    sw          s6, -4(t1)
    bgtz        t9, 0b
     sw         s7, 92(t1)
    // second loop three pass
    li          t9, 3
1:
    lh          s6, 34(t0)      // inptr[DCTSIZE*2]
    lh          t6, 34(a0)      // quantptr[DCTSIZE*2]
    lh          s7, 98(t0)      // inptr[DCTSIZE*6]
    lh          t7, 98(a0)      // quantptr[DCTSIZE*6]
    mul         t6, s6, t6      // z2 = (inptr[DCTSIZE*2] * quantptr[DCTSIZE*2])
    lh          s4, 2(t0)       // inptr[DCTSIZE*0]
    mul         t7, s7, t7      // z3 = (inptr[DCTSIZE*6] * quantptr[DCTSIZE*6])
    lh          s5, 2(a0)       // quantptr[DCTSIZE*0]
    li          s6, 15137
    li          s7, 6270
    mul         t2, s4, s5      // tmp0 = (inptr[0] * quantptr[0])
    mul         v0, s6, t6      // z2 = (inptr[DCTSIZE*2] * quantptr[DCTSIZE*2])
    lh          t5, 114(t0)     // inptr[DCTSIZE*7]
    mul         t7, s7, t7      // z3 = (inptr[DCTSIZE*6] * quantptr[DCTSIZE*6])
    lh          s4, 114(a0)     // quantptr[DCTSIZE*7]
    lh          s5, 82(a0)      // quantptr[DCTSIZE*5]
    lh          t6, 82(t0)      // inptr[DCTSIZE*5]
    sll         t2, t2, 14      // tmp0 <<= (CONST_BITS+1)
    lh          s6, 50(a0)      // quantptr[DCTSIZE*3]
    lh          t8, 18(t0)      // inptr[DCTSIZE*1]
    subu        v0, v0, t7      // tmp2 = MULTIPLY(z2, t5) - MULTIPLY(z3, t6)
    lh          t7, 50(t0)      // inptr[DCTSIZE*3]
    lh          s7, 18(a0)      // quantptr[DCTSIZE*1]
    mul         t5, s4, t5      // z1 = (inptr[DCTSIZE*7] * quantptr[DCTSIZE*7])
    mul         t6, s5, t6      // z2 = (inptr[DCTSIZE*5] * quantptr[DCTSIZE*5])
    mul         t7, s6, t7      // z3 = (inptr[DCTSIZE*3] * quantptr[DCTSIZE*3])
    mul         t8, s7, t8      // z4 = (inptr[DCTSIZE*1] * quantptr[DCTSIZE*1])
    addu        t3, t2, v0      // tmp10 = tmp0 + z2
    subu        t4, t2, v0      // tmp10 = tmp0 - z2
    mult        $ac0, zero, zero
    mult        $ac1, zero, zero
    ins         t5, t6, 16, 16
    ins         t7, t8, 16, 16
    dpa.w.ph    $ac0, t5, s0
    dpa.w.ph    $ac0, t7, s1
    dpa.w.ph    $ac1, t5, s2
    dpa.w.ph    $ac1, t7, s3
    mflo        t5, $ac0
    mflo        t6, $ac1
    addiu       t9, t9, -1
    addiu       t0, t0, 2
    addiu       a0, a0, 2
    addiu       t1, t1, 4
    addu        s5, t4, t5
    subu        s4, t4, t5
    addu        s6, t3, t6
    subu        s7, t3, t6
    shra_r.w    s5, s5, 12      // DESCALE(tmp12 + temp1, 12)
    shra_r.w    s4, s4, 12      // DESCALE(tmp12 - temp1, 12)
    shra_r.w    s6, s6, 12      // DESCALE(tmp10 + temp2, 12)
    shra_r.w    s7, s7, 12      // DESCALE(tmp10 - temp2, 12)
    sw          s5, 32(t1)
    sw          s4, 64(t1)
    sw          s6, 0(t1)
    bgtz        t9, 1b
     sw         s7, 96(t1)
    move        t1, v1
    li          s4, 15137
    lw          s6, 8(t1)       // wsptr[2]
    li          s5, 6270
    lw          s7, 24(t1)      // wsptr[6]
    mul         s4, s4, s6      // MULTIPLY((JLONG)wsptr[2], FIX_1_847759065)
    lw          t2, 0(t1)       // wsptr[0]
    mul         s5, s5, s7      // MULTIPLY((JLONG)wsptr[6], -FIX_0_765366865)
    lh          t5, 28(t1)      // wsptr[7]
    lh          t6, 20(t1)      // wsptr[5]
    lh          t7, 12(t1)      // wsptr[3]
    lh          t8, 4(t1)       // wsptr[1]
    ins         t5, t6, 16, 16
    ins         t7, t8, 16, 16
    mult        $ac0, zero, zero
    dpa.w.ph    $ac0, t5, s0
    dpa.w.ph    $ac0, t7, s1
    mult        $ac1, zero, zero
    dpa.w.ph    $ac1, t5, s2
    dpa.w.ph    $ac1, t7, s3
    sll         t2, t2, 14      // tmp0 = ((JLONG)wsptr[0]) << (CONST_BITS+1)
    mflo        s6, $ac0
    // MULTIPLY(wsptr[2], FIX_1_847759065 + MULTIPLY(wsptr[6], -FIX_0_765366865)
    subu        s4, s4, s5
    addu        t3, t2, s4      // tmp10 = tmp0 + z2
    mflo        s7, $ac1
    subu        t4, t2, s4      // tmp10 = tmp0 - z2
    addu        t7, t4, s6
    subu        t8, t4, s6
    addu        t5, t3, s7
    subu        t6, t3, s7
    shra_r.w    t5, t5, 19      // DESCALE(tmp10 + temp2, 19)
    shra_r.w    t6, t6, 19      // DESCALE(tmp10 - temp2, 19)
    shra_r.w    t7, t7, 19      // DESCALE(tmp12 + temp1, 19)
    shra_r.w    t8, t8, 19      // DESCALE(tmp12 - temp1, 19)
    sll         s4, t9, 2
    lw          v0, 0(a2)       // output_buf[ctr]
    shll_s.w    t5, t5, 24
    shll_s.w    t6, t6, 24
    shll_s.w    t7, t7, 24
    shll_s.w    t8, t8, 24
    sra         t5, t5, 24
    sra         t6, t6, 24
    sra         t7, t7, 24
    sra         t8, t8, 24
    addu        v0, v0, a3      // outptr = output_buf[ctr] + output_col
    addiu       t5, t5, 128
    addiu       t6, t6, 128
    addiu       t7, t7, 128
    addiu       t8, t8, 128
    sb          t5, 0(v0)
    sb          t7, 1(v0)
    sb          t8, 2(v0)
    sb          t6, 3(v0)
    // 2
    li          s4, 15137
    lw          s6, 40(t1)      // wsptr[2]
    li          s5, 6270
    lw          s7, 56(t1)      // wsptr[6]
    mul         s4, s4, s6      // MULTIPLY((JLONG)wsptr[2], FIX_1_847759065)
    lw          t2, 32(t1)      // wsptr[0]
    mul         s5, s5, s7      // MULTIPLY((JLONG)wsptr[6], -FIX_0_765366865)
    lh          t5, 60(t1)      // wsptr[7]
    lh          t6, 52(t1)      // wsptr[5]
    lh          t7, 44(t1)      // wsptr[3]
    lh          t8, 36(t1)      // wsptr[1]
    ins         t5, t6, 16, 16
    ins         t7, t8, 16, 16
    mult        $ac0, zero, zero
    dpa.w.ph    $ac0, t5, s0
    dpa.w.ph    $ac0, t7, s1
    mult        $ac1, zero, zero
    dpa.w.ph    $ac1, t5, s2
    dpa.w.ph    $ac1, t7, s3
    sll         t2, t2, 14      // tmp0 = ((JLONG)wsptr[0]) << (CONST_BITS+1)
    mflo        s6, $ac0
    // MULTIPLY(wsptr[2], FIX_1_847759065 + MULTIPLY(wsptr[6], -FIX_0_765366865)
    subu        s4, s4, s5
    addu        t3, t2, s4      // tmp10 = tmp0 + z2
    mflo        s7, $ac1
    subu        t4, t2, s4      // tmp10 = tmp0 - z2
    addu        t7, t4, s6
    subu        t8, t4, s6
    addu        t5, t3, s7
    subu        t6, t3, s7
    shra_r.w    t5, t5, 19      // DESCALE(tmp10 + temp2, CONST_BITS-PASS1_BITS+1)
    shra_r.w    t6, t6, 19      // DESCALE(tmp10 - temp2, CONST_BITS-PASS1_BITS+1)
    shra_r.w    t7, t7, 19      // DESCALE(tmp12 + temp1, CONST_BITS-PASS1_BITS+1)
    shra_r.w    t8, t8, 19      // DESCALE(tmp12 - temp1, CONST_BITS-PASS1_BITS+1)
    sll         s4, t9, 2
    lw          v0, 4(a2)       // output_buf[ctr]
    shll_s.w    t5, t5, 24
    shll_s.w    t6, t6, 24
    shll_s.w    t7, t7, 24
    shll_s.w    t8, t8, 24
    sra         t5, t5, 24
    sra         t6, t6, 24
    sra         t7, t7, 24
    sra         t8, t8, 24
    addu        v0, v0, a3      // outptr = output_buf[ctr] + output_col
    addiu       t5, t5, 128
    addiu       t6, t6, 128
    addiu       t7, t7, 128
    addiu       t8, t8, 128
    sb          t5, 0(v0)
    sb          t7, 1(v0)
    sb          t8, 2(v0)
    sb          t6, 3(v0)
    // 3
    li          s4, 15137
    lw          s6, 72(t1)      // wsptr[2]
    li          s5, 6270
    lw          s7, 88(t1)      // wsptr[6]
    mul         s4, s4, s6      // MULTIPLY((JLONG)wsptr[2], FIX_1_847759065)
    lw          t2, 64(t1)      // wsptr[0]
    mul         s5, s5, s7      // MULTIPLY((JLONG)wsptr[6], -FIX_0_765366865)
    lh          t5, 92(t1)      // wsptr[7]
    lh          t6, 84(t1)      // wsptr[5]
    lh          t7, 76(t1)      // wsptr[3]
    lh          t8, 68(t1)      // wsptr[1]
    ins         t5, t6, 16, 16
    ins         t7, t8, 16, 16
    mult        $ac0, zero, zero
    dpa.w.ph    $ac0, t5, s0
    dpa.w.ph    $ac0, t7, s1
    mult        $ac1, zero, zero
    dpa.w.ph    $ac1, t5, s2
    dpa.w.ph    $ac1, t7, s3
    sll         t2, t2, 14      // tmp0 = ((JLONG)wsptr[0]) << (CONST_BITS+1)
    mflo        s6, $ac0
    // MULTIPLY(wsptr[2], FIX_1_847759065 + MULTIPLY(wsptr[6], -FIX_0_765366865)
    subu        s4, s4, s5
    addu        t3, t2, s4      // tmp10 = tmp0 + z2
    mflo        s7, $ac1
    subu        t4, t2, s4      // tmp10 = tmp0 - z2
    addu        t7, t4, s6
    subu        t8, t4, s6
    addu        t5, t3, s7
    subu        t6, t3, s7
    shra_r.w    t5, t5, 19      // DESCALE(tmp10 + temp2, 19)
    shra_r.w    t6, t6, 19      // DESCALE(tmp10 - temp2, 19)
    shra_r.w    t7, t7, 19      // DESCALE(tmp12 + temp1, 19)
    shra_r.w    t8, t8, 19      // DESCALE(tmp12 - temp1, 19)
    sll         s4, t9, 2
    lw          v0, 8(a2)       // output_buf[ctr]
    shll_s.w    t5, t5, 24
    shll_s.w    t6, t6, 24
    shll_s.w    t7, t7, 24
    shll_s.w    t8, t8, 24
    sra         t5, t5, 24
    sra         t6, t6, 24
    sra         t7, t7, 24
    sra         t8, t8, 24
    addu        v0, v0, a3      // outptr = output_buf[ctr] + output_col
    addiu       t5, t5, 128
    addiu       t6, t6, 128
    addiu       t7, t7, 128
    addiu       t8, t8, 128
    sb          t5, 0(v0)
    sb          t7, 1(v0)
    sb          t8, 2(v0)
    sb          t6, 3(v0)
    li          s4, 15137
    lw          s6, 104(t1)     // wsptr[2]
    li          s5, 6270
    lw          s7, 120(t1)     // wsptr[6]
    mul         s4, s4, s6      // MULTIPLY((JLONG)wsptr[2], FIX_1_847759065)
    lw          t2, 96(t1)      // wsptr[0]
    mul         s5, s5, s7      // MULTIPLY((JLONG)wsptr[6], -FIX_0_765366865)
    lh          t5, 124(t1)     // wsptr[7]
    lh          t6, 116(t1)     // wsptr[5]
    lh          t7, 108(t1)     // wsptr[3]
    lh          t8, 100(t1)     // wsptr[1]
    ins         t5, t6, 16, 16
    ins         t7, t8, 16, 16
    mult        $ac0, zero, zero
    dpa.w.ph    $ac0, t5, s0
    dpa.w.ph    $ac0, t7, s1
    mult        $ac1, zero, zero
    dpa.w.ph    $ac1, t5, s2
    dpa.w.ph    $ac1, t7, s3
    sll         t2, t2, 14      // tmp0 = ((JLONG)wsptr[0]) << (CONST_BITS+1)
    mflo        s6, $ac0
    // MULTIPLY(wsptr[2], FIX_1_847759065 + MULTIPLY(wsptr[6], -FIX_0_765366865)
    subu        s4, s4, s5
    addu        t3, t2, s4      // tmp10 = tmp0 + z2;
    mflo        s7, $ac1
    subu        t4, t2, s4      // tmp10 = tmp0 - z2;
    addu        t7, t4, s6
    subu        t8, t4, s6
    addu        t5, t3, s7
    subu        t6, t3, s7
    shra_r.w    t5, t5, 19      // DESCALE(tmp10 + temp2, 19)
    shra_r.w    t6, t6, 19      // DESCALE(tmp10 - temp2, 19)
    shra_r.w    t7, t7, 19      // DESCALE(tmp12 + temp1, 19)
    shra_r.w    t8, t8, 19      // DESCALE(tmp12 - temp1, 19)
    sll         s4, t9, 2
    lw          v0, 12(a2)      // output_buf[ctr]
    shll_s.w    t5, t5, 24
    shll_s.w    t6, t6, 24
    shll_s.w    t7, t7, 24
    shll_s.w    t8, t8, 24
    sra         t5, t5, 24
    sra         t6, t6, 24
    sra         t7, t7, 24
    sra         t8, t8, 24
    addu        v0, v0, a3      // outptr = output_buf[ctr] + output_col
    addiu       t5, t5, 128
    addiu       t6, t6, 128
    addiu       t7, t7, 128
    addiu       t8, t8, 128
    sb          t5, 0(v0)
    sb          t7, 1(v0)
    sb          t8, 2(v0)
    sb          t6, 3(v0)

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop
END(jsimd_idct_4x4_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_6x6_dspr2)
/*
 * a0 = compptr->dct_table
 * a1 = coef_block
 * a2 = output_buf
 * a3 = output_col
 */
    .set at

    SAVE_REGS_ON_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    addiu       sp, sp, -144
    move        v0, sp
    addiu       v1, v0, 24
    addiu       t9, zero, 5793
    addiu       s0, zero, 10033
    addiu       s1, zero, 2998

1:
    lh          s2, 0(a0)       // q0 = quantptr[ 0]
    lh          s3, 32(a0)      // q1 = quantptr[16]
    lh          s4, 64(a0)      // q2 = quantptr[32]
    lh          t2, 64(a1)      // tmp2 = inptr[32]
    lh          t1, 32(a1)      // tmp1 = inptr[16]
    lh          t0, 0(a1)       // tmp0 = inptr[ 0]
    mul         t2, t2, s4      // tmp2 = tmp2 * q2
    mul         t1, t1, s3      // tmp1 = tmp1 * q1
    mul         t0, t0, s2      // tmp0 = tmp0 * q0
    lh          t6, 16(a1)      // z1 = inptr[ 8]
    lh          t8, 80(a1)      // z3 = inptr[40]
    lh          t7, 48(a1)      // z2 = inptr[24]
    lh          s2, 16(a0)      // q0 = quantptr[ 8]
    lh          s4, 80(a0)      // q2 = quantptr[40]
    lh          s3, 48(a0)      // q1 = quantptr[24]
    mul         t2, t2, t9      // tmp2 = tmp2 * 5793
    mul         t1, t1, s0      // tmp1 = tmp1 * 10033
    sll         t0, t0, 13      // tmp0 = tmp0 << 13
    mul         t6, t6, s2      // z1 = z1 * q0
    mul         t8, t8, s4      // z3 = z3 * q2
    mul         t7, t7, s3      // z2 = z2 * q1
    addu        t3, t0, t2      // tmp10 = tmp0 + tmp2
    sll         t2, t2, 1       // tmp2 = tmp2 << 2
    subu        t4, t0, t2      // tmp11 = tmp0 - tmp2;
    subu        t5, t3, t1      // tmp12 = tmp10 - tmp1
    addu        t3, t3, t1      // tmp10 = tmp10 + tmp1
    addu        t1, t6, t8      // tmp1 = z1 + z3
    mul         t1, t1, s1      // tmp1 = tmp1 * 2998
    shra_r.w    t4, t4, 11      // tmp11 = (tmp11 + 1024) >> 11
    subu        t2, t6, t8      // tmp2 = z1 - z3
    subu        t2, t2, t7      // tmp2 = tmp2 - z2
    sll         t2, t2, 2       // tmp2 = tmp2 << 2
    addu        t0, t6, t7      // tmp0 = z1 + z2
    sll         t0, t0, 13      // tmp0 = tmp0 << 13
    subu        s2, t8, t7      // q0 = z3 - z2
    sll         s2, s2, 13      // q0 = q0 << 13
    addu        t0, t0, t1      // tmp0 = tmp0 + tmp1
    addu        t1, s2, t1      // tmp1 = q0 + tmp1
    addu        s2, t4, t2      // q0 = tmp11 + tmp2
    subu        s3, t4, t2      // q1 = tmp11 - tmp2
    addu        t6, t3, t0      // z1 = tmp10 + tmp0
    subu        t7, t3, t0      // z2 = tmp10 - tmp0
    addu        t4, t5, t1      // tmp11 = tmp12 + tmp1
    subu        t5, t5, t1      // tmp12 = tmp12 - tmp1
    shra_r.w    t6, t6, 11      // z1 = (z1 + 1024) >> 11
    shra_r.w    t7, t7, 11      // z2 = (z2 + 1024) >> 11
    shra_r.w    t4, t4, 11      // tmp11 = (tmp11 + 1024) >> 11
    shra_r.w    t5, t5, 11      // tmp12 = (tmp12 + 1024) >> 11
    sw          s2, 24(v0)
    sw          s3, 96(v0)
    sw          t6, 0(v0)
    sw          t7, 120(v0)
    sw          t4, 48(v0)
    sw          t5, 72(v0)
    addiu       v0, v0, 4
    addiu       a1, a1, 2
    bne         v0, v1, 1b
     addiu      a0, a0, 2

    /* Pass 2: process 6 rows from work array, store into output array. */
    move        v0, sp
    addiu       v1, v0, 144

2:
    lw          t0, 0(v0)
    lw          t2, 16(v0)
    lw          s5, 0(a2)
    addiu       t0, t0, 16
    sll         t0, t0, 13
    mul         t3, t2, t9
    lw          t6, 4(v0)
    lw          t8, 20(v0)
    lw          t7, 12(v0)
    addu        s5, s5, a3
    addu        s6, t6, t8
    mul         s6, s6, s1
    addu        t1, t0, t3
    subu        t4, t0, t3
    subu        t4, t4, t3
    lw          t3, 8(v0)
    mul         t0, t3, s0
    addu        s7, t6, t7
    sll         s7, s7, 13
    addu        s7, s6, s7
    subu        t2, t8, t7
    sll         t2, t2, 13
    addu        t2, s6, t2
    subu        s6, t6, t7
    subu        s6, s6, t8
    sll         s6, s6, 13
    addu        t3, t1, t0
    subu        t5, t1, t0
    addu        t6, t3, s7
    subu        t3, t3, s7
    addu        t7, t4, s6
    subu        t4, t4, s6
    addu        t8, t5, t2
    subu        t5, t5, t2
    shll_s.w    t6, t6, 6
    shll_s.w    t3, t3, 6
    shll_s.w    t7, t7, 6
    shll_s.w    t4, t4, 6
    shll_s.w    t8, t8, 6
    shll_s.w    t5, t5, 6
    sra         t6, t6, 24
    addiu       t6, t6, 128
    sra         t3, t3, 24
    addiu       t3, t3, 128
    sb          t6, 0(s5)
    sra         t7, t7, 24
    addiu       t7, t7, 128
    sb          t3, 5(s5)
    sra         t4, t4, 24
    addiu       t4, t4, 128
    sb          t7, 1(s5)
    sra         t8, t8, 24
    addiu       t8, t8, 128
    sb          t4, 4(s5)
    addiu       v0, v0, 24
    sra         t5, t5, 24
    addiu       t5, t5, 128
    sb          t8, 2(s5)
    addiu       a2, a2,  4
    bne         v0, v1, 2b
     sb         t5, 3(s5)

    addiu       sp, sp, 144

    RESTORE_REGS_FROM_STACK 32, s0, s1, s2, s3, s4, s5, s6, s7

    j           ra
     nop

END(jsimd_idct_6x6_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_12x12_pass1_dspr2)
/*
 * a0 = compptr->dct_table
 * a1 = coef_block
 * a2 = workspace
 */
    SAVE_REGS_ON_STACK 16, s0, s1, s2, s3

    li          a3, 8

1:
    // odd part
    lh          t0, 48(a1)
    lh          t1, 48(a0)
    lh          t2, 16(a1)
    lh          t3, 16(a0)
    lh          t4, 80(a1)
    lh          t5, 80(a0)
    lh          t6, 112(a1)
    lh          t7, 112(a0)
    mul         t0, t0, t1      // z2
    mul         t1, t2, t3      // z1
    mul         t2, t4, t5      // z3
    mul         t3, t6, t7      // z4
    li          t4, 10703       // FIX(1.306562965)
    li          t5, 4433        // FIX_0_541196100
    li          t6, 7053        // FIX(0.860918669)
    mul         t4, t0, t4      // tmp11
    mul         t5, t0, t5      // -tmp14
    addu        t7, t1, t2      // tmp10
    addu        t8, t7, t3      // tmp10 + z4
    mul         t6, t6, t8      // tmp15
    li          t8, 2139        // FIX(0.261052384)
    mul         t8, t7, t8      // MULTIPLY(tmp10, FIX(0.261052384))
    li          t7, 2295        // FIX(0.280143716)
    mul         t7, t1, t7      // MULTIPLY(z1, FIX(0.280143716))
    addu        t9, t2, t3      // z3 + z4
    li          s0, 8565        // FIX(1.045510580)
    mul         t9, t9, s0      // -tmp13
    li          s0, 12112       // FIX(1.478575242)
    mul         s0, t2, s0      // MULTIPLY(z3, FIX(1.478575242)
    li          s1, 12998       // FIX(1.586706681)
    mul         s1, t3, s1      // MULTIPLY(z4, FIX(1.586706681))
    li          s2, 5540        // FIX(0.676326758)
    mul         s2, t1, s2      // MULTIPLY(z1, FIX(0.676326758))
    li          s3, 16244       // FIX(1.982889723)
    mul         s3, t3, s3      // MULTIPLY(z4, FIX(1.982889723))
    subu        t1, t1, t3      // z1-=z4
    subu        t0, t0, t2      // z2-=z3
    addu        t2, t0, t1      // z1+z2
    li          t3, 4433        // FIX_0_541196100
    mul         t2, t2, t3      // z3
    li          t3, 6270        // FIX_0_765366865
    mul         t1, t1, t3      // MULTIPLY(z1, FIX_0_765366865)
    li          t3, 15137       // FIX_0_765366865
    mul         t0, t0, t3      // MULTIPLY(z2, FIX_1_847759065)
    addu        t8, t6, t8      // tmp12
    addu        t3, t8, t4      // tmp12 + tmp11
    addu        t3, t3, t7      // tmp10
    subu        t8, t8, t9      // tmp12 + tmp13
    addu        s0, t5, s0
    subu        t8, t8, s0      // tmp12
    subu        t9, t6, t9
    subu        s1, s1, t4
    addu        t9, t9, s1      // tmp13
    subu        t6, t6, t5
    subu        t6, t6, s2
    subu        t6, t6, s3      // tmp15
    // even part start
    lh          t4, 64(a1)
    lh          t5, 64(a0)
    lh          t7, 32(a1)
    lh          s0, 32(a0)
    lh          s1, 0(a1)
    lh          s2, 0(a0)
    lh          s3, 96(a1)
    lh          v0, 96(a0)
    mul         t4, t4, t5      // DEQUANTIZE(inptr[DCTSIZE*4], quantptr[DCTSIZE*4])
    mul         t5, t7, s0      // DEQUANTIZE(inptr[DCTSIZE*2], quantptr[DCTSIZE*2])
    mul         t7, s1, s2      // DEQUANTIZE(inptr[DCTSIZE*0], quantptr[DCTSIZE*0])
    mul         s0, s3, v0      // DEQUANTIZE(inptr[DCTSIZE*6], quantptr[DCTSIZE*6])
    // odd part end
    addu        t1, t2, t1      // tmp11
    subu        t0, t2, t0      // tmp14
    // update counter and pointers
    addiu       a3, a3, -1
    addiu       a0, a0, 2
    addiu       a1, a1, 2
    // even part rest
    li          s1, 10033
    li          s2, 11190
    mul         t4, t4, s1      // z4
    mul         s1, t5, s2      // z4
    sll         t5, t5, 13      // z1
    sll         t7, t7, 13
    addiu       t7, t7, 1024    // z3
    sll         s0, s0, 13      // z2
    addu        s2, t7, t4      // tmp10
    subu        t4, t7, t4      // tmp11
    subu        s3, t5, s0      // tmp12
    addu        t2, t7, s3      // tmp21
    subu        s3, t7, s3      // tmp24
    addu        t7, s1, s0      // tmp12
    addu        v0, s2, t7      // tmp20
    subu        s2, s2, t7      // tmp25
    subu        s1, s1, t5      // z4 - z1
    subu        s1, s1, s0      // tmp12
    addu        s0, t4, s1      // tmp22
    subu        t4, t4, s1      // tmp23
    // final output stage
    addu        t5, v0, t3
    subu        v0, v0, t3
    addu        t3, t2, t1
    subu        t2, t2, t1
    addu        t1, s0, t8
    subu        s0, s0, t8
    addu        t8, t4, t9
    subu        t4, t4, t9
    addu        t9, s3, t0
    subu        s3, s3, t0
    addu        t0, s2, t6
    subu        s2, s2, t6
    sra         t5, t5, 11
    sra         t3, t3, 11
    sra         t1, t1, 11
    sra         t8, t8, 11
    sra         t9, t9, 11
    sra         t0, t0, 11
    sra         s2, s2, 11
    sra         s3, s3, 11
    sra         t4, t4, 11
    sra         s0, s0, 11
    sra         t2, t2, 11
    sra         v0, v0, 11
    sw          t5, 0(a2)
    sw          t3, 32(a2)
    sw          t1, 64(a2)
    sw          t8, 96(a2)
    sw          t9, 128(a2)
    sw          t0, 160(a2)
    sw          s2, 192(a2)
    sw          s3, 224(a2)
    sw          t4, 256(a2)
    sw          s0, 288(a2)
    sw          t2, 320(a2)
    sw          v0, 352(a2)
    bgtz        a3, 1b
     addiu      a2, a2, 4

    RESTORE_REGS_FROM_STACK 16, s0, s1, s2, s3

    j           ra
     nop

END(jsimd_idct_12x12_pass1_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_idct_12x12_pass2_dspr2)
/*
 * a0 = workspace
 * a1 = output
 */
    SAVE_REGS_ON_STACK 16, s0, s1, s2, s3

    li          a3, 12

1:
    // Odd part
    lw          t0, 12(a0)
    lw          t1, 4(a0)
    lw          t2, 20(a0)
    lw          t3, 28(a0)
    li          t4, 10703       // FIX(1.306562965)
    li          t5, 4433        // FIX_0_541196100
    mul         t4, t0, t4      // tmp11
    mul         t5, t0, t5      // -tmp14
    addu        t6, t1, t2      // tmp10
    li          t7, 2139        // FIX(0.261052384)
    mul         t7, t6, t7      // MULTIPLY(tmp10, FIX(0.261052384))
    addu        t6, t6, t3      // tmp10 + z4
    li          t8, 7053        // FIX(0.860918669)
    mul         t6, t6, t8      // tmp15
    li          t8, 2295        // FIX(0.280143716)
    mul         t8, t1, t8      // MULTIPLY(z1, FIX(0.280143716))
    addu        t9, t2, t3      // z3 + z4
    li          s0, 8565        // FIX(1.045510580)
    mul         t9, t9, s0      // -tmp13
    li          s0, 12112       // FIX(1.478575242)
    mul         s0, t2, s0      // MULTIPLY(z3, FIX(1.478575242))
    li          s1, 12998       // FIX(1.586706681)
    mul         s1, t3, s1      // MULTIPLY(z4, FIX(1.586706681))
    li          s2, 5540        // FIX(0.676326758)
    mul         s2, t1, s2      // MULTIPLY(z1, FIX(0.676326758))
    li          s3, 16244       // FIX(1.982889723)
    mul         s3, t3, s3      // MULTIPLY(z4, FIX(1.982889723))
    subu        t1, t1, t3      // z1 -= z4
    subu        t0, t0, t2      // z2 -= z3
    addu        t2, t1, t0      // z1 + z2
    li          t3, 4433        // FIX_0_541196100
    mul         t2, t2, t3      // z3
    li          t3, 6270        // FIX_0_765366865
    mul         t1, t1, t3      // MULTIPLY(z1, FIX_0_765366865)
    li          t3, 15137       // FIX_1_847759065
    mul         t0, t0, t3      // MULTIPLY(z2, FIX_1_847759065)
    addu        t3, t6, t7      // tmp12
    addu        t7, t3, t4
    addu        t7, t7, t8      // tmp10
    subu        t3, t3, t9
    subu        t3, t3, t5
    subu        t3, t3, s0      // tmp12
    subu        t9, t6, t9
    subu        t9, t9, t4
    addu        t9, t9, s1      // tmp13
    subu        t6, t6, t5
    subu        t6, t6, s2
    subu        t6, t6, s3      // tmp15
    addu        t1, t2, t1      // tmp11
    subu        t0, t2, t0      // tmp14
    // even part
    lw          t2, 16(a0)      // z4
    lw          t4, 8(a0)       // z1
    lw          t5, 0(a0)       // z3
    lw          t8, 24(a0)      // z2
    li          s0, 10033       // FIX(1.224744871)
    li          s1, 11190       // FIX(1.366025404)
    mul         t2, t2, s0      // z4
    mul         s0, t4, s1      // z4
    addiu       t5, t5, 0x10
    sll         t5, t5, 13      // z3
    sll         t4, t4, 13      // z1
    sll         t8, t8, 13      // z2
    subu        s1, t4, t8      // tmp12
    addu        s2, t5, t2      // tmp10
    subu        t2, t5, t2      // tmp11
    addu        s3, t5, s1      // tmp21
    subu        s1, t5, s1      // tmp24
    addu        t5, s0, t8      // tmp12
    addu        v0, s2, t5      // tmp20
    subu        t5, s2, t5      // tmp25
    subu        t4, s0, t4
    subu        t4, t4, t8      // tmp12
    addu        t8, t2, t4      // tmp22
    subu        t2, t2, t4      // tmp23
    // increment counter and pointers
    addiu       a3, a3, -1
    addiu       a0, a0, 32
    // Final stage
    addu        t4, v0, t7
    subu        v0, v0, t7
    addu        t7, s3, t1
    subu        s3, s3, t1
    addu        t1, t8, t3
    subu        t8, t8, t3
    addu        t3, t2, t9
    subu        t2, t2, t9
    addu        t9, s1, t0
    subu        s1, s1, t0
    addu        t0, t5, t6
    subu        t5, t5, t6
    sll         t4, t4, 4
    sll         t7, t7, 4
    sll         t1, t1, 4
    sll         t3, t3, 4
    sll         t9, t9, 4
    sll         t0, t0, 4
    sll         t5, t5, 4
    sll         s1, s1, 4
    sll         t2, t2, 4
    sll         t8, t8, 4
    sll         s3, s3, 4
    sll         v0, v0, 4
    shll_s.w    t4, t4, 2
    shll_s.w    t7, t7, 2
    shll_s.w    t1, t1, 2
    shll_s.w    t3, t3, 2
    shll_s.w    t9, t9, 2
    shll_s.w    t0, t0, 2
    shll_s.w    t5, t5, 2
    shll_s.w    s1, s1, 2
    shll_s.w    t2, t2, 2
    shll_s.w    t8, t8, 2
    shll_s.w    s3, s3, 2
    shll_s.w    v0, v0, 2
    srl         t4, t4, 24
    srl         t7, t7, 24
    srl         t1, t1, 24
    srl         t3, t3, 24
    srl         t9, t9, 24
    srl         t0, t0, 24
    srl         t5, t5, 24
    srl         s1, s1, 24
    srl         t2, t2, 24
    srl         t8, t8, 24
    srl         s3, s3, 24
    srl         v0, v0, 24
    lw          t6, 0(a1)
    addiu       t4, t4, 0x80
    addiu       t7, t7, 0x80
    addiu       t1, t1, 0x80
    addiu       t3, t3, 0x80
    addiu       t9, t9, 0x80
    addiu       t0, t0, 0x80
    addiu       t5, t5, 0x80
    addiu       s1, s1, 0x80
    addiu       t2, t2, 0x80
    addiu       t8, t8, 0x80
    addiu       s3, s3, 0x80
    addiu       v0, v0, 0x80
    sb          t4, 0(t6)
    sb          t7, 1(t6)
    sb          t1, 2(t6)
    sb          t3, 3(t6)
    sb          t9, 4(t6)
    sb          t0, 5(t6)
    sb          t5, 6(t6)
    sb          s1, 7(t6)
    sb          t2, 8(t6)
    sb          t8, 9(t6)
    sb          s3, 10(t6)
    sb          v0, 11(t6)
    bgtz        a3, 1b
     addiu      a1, a1, 4

    RESTORE_REGS_FROM_STACK 16, s0, s1, s2, s3

    jr          ra
     nop

END(jsimd_idct_12x12_pass2_dspr2)


/*****************************************************************************/
LEAF_DSPR2(jsimd_convsamp_dspr2)
/*
 * a0 = sample_data
 * a1 = start_col
 * a2 = workspace
 */
    lw            t0, 0(a0)
    li            t7, 0xff80ff80
    addu          t0, t0, a1
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    lw            t0, 4(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 0(a2)
    usw           t4, 4(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 8(a2)
    usw           t6, 12(a2)

    lw            t0, 8(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 16(a2)
    usw           t4, 20(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 24(a2)
    usw           t6, 28(a2)

    lw            t0, 12(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 32(a2)
    usw           t4, 36(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 40(a2)
    usw           t6, 44(a2)

    lw            t0, 16(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 48(a2)
    usw           t4, 52(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 56(a2)
    usw           t6, 60(a2)

    lw            t0, 20(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 64(a2)
    usw           t4, 68(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 72(a2)
    usw           t6, 76(a2)

    lw            t0, 24(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 80(a2)
    usw           t4, 84(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 88(a2)
    usw           t6, 92(a2)

    lw            t0, 28(a0)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu          t0, t0, a1
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    ulw           t1, 0(t0)
    ulw           t2, 4(t0)
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 96(a2)
    usw           t4, 100(a2)
    preceu.ph.qbr t3, t1
    preceu.ph.qbl t4, t1
    usw           t5, 104(a2)
    usw           t6, 108(a2)
    preceu.ph.qbr t5, t2
    preceu.ph.qbl t6, t2
    addu.ph       t3, t3, t7
    addu.ph       t4, t4, t7
    addu.ph       t5, t5, t7
    addu.ph       t6, t6, t7
    usw           t3, 112(a2)
    usw           t4, 116(a2)
    usw           t5, 120(a2)
    usw           t6, 124(a2)

    j             ra
     nop

END(jsimd_convsamp_dspr2)


#ifndef __mips_soft_float

/*****************************************************************************/
LEAF_DSPR2(jsimd_convsamp_float_dspr2)
/*
 * a0 = sample_data
 * a1 = start_col
 * a2 = workspace
 */
    .set at

    lw          t0, 0(a0)
    addu        t0, t0, a1
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 4(a0)
    swc1        f2, 0(a2)
    swc1        f4, 4(a2)
    swc1        f6, 8(a2)
    addu        t0, t0, a1
    swc1        f8, 12(a2)
    swc1        f10, 16(a2)
    swc1        f12, 20(a2)
    swc1        f14, 24(a2)
    swc1        f16, 28(a2)
    // elemr 1
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 8(a0)
    swc1        f2, 32(a2)
    swc1        f4, 36(a2)
    swc1        f6, 40(a2)
    addu        t0, t0, a1
    swc1        f8, 44(a2)
    swc1        f10, 48(a2)
    swc1        f12, 52(a2)
    swc1        f14, 56(a2)
    swc1        f16, 60(a2)
    // elemr 2
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 12(a0)
    swc1        f2, 64(a2)
    swc1        f4, 68(a2)
    swc1        f6, 72(a2)
    addu        t0, t0, a1
    swc1        f8, 76(a2)
    swc1        f10, 80(a2)
    swc1        f12, 84(a2)
    swc1        f14, 88(a2)
    swc1        f16, 92(a2)
    //  elemr 3
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 16(a0)
    swc1        f2, 96(a2)
    swc1        f4, 100(a2)
    swc1        f6, 104(a2)
    addu        t0, t0, a1
    swc1        f8, 108(a2)
    swc1        f10, 112(a2)
    swc1        f12, 116(a2)
    swc1        f14, 120(a2)
    swc1        f16, 124(a2)
    // elemr 4
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 20(a0)
    swc1        f2, 128(a2)
    swc1        f4, 132(a2)
    swc1        f6, 136(a2)
    addu        t0, t0, a1
    swc1        f8, 140(a2)
    swc1        f10, 144(a2)
    swc1        f12, 148(a2)
    swc1        f14, 152(a2)
    swc1        f16, 156(a2)
    // elemr 5
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 24(a0)
    swc1        f2, 160(a2)
    swc1        f4, 164(a2)
    swc1        f6, 168(a2)
    addu        t0, t0, a1
    swc1        f8, 172(a2)
    swc1        f10, 176(a2)
    swc1        f12, 180(a2)
    swc1        f14, 184(a2)
    swc1        f16, 188(a2)
    // elemr 6
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    lw          t0, 28(a0)
    swc1        f2, 192(a2)
    swc1        f4, 196(a2)
    swc1        f6, 200(a2)
    addu        t0, t0, a1
    swc1        f8, 204(a2)
    swc1        f10, 208(a2)
    swc1        f12, 212(a2)
    swc1        f14, 216(a2)
    swc1        f16, 220(a2)
    // elemr 7
    lbu         t1, 0(t0)
    lbu         t2, 1(t0)
    lbu         t3, 2(t0)
    lbu         t4, 3(t0)
    lbu         t5, 4(t0)
    lbu         t6, 5(t0)
    lbu         t7, 6(t0)
    lbu         t8, 7(t0)
    addiu       t1, t1, -128
    addiu       t2, t2, -128
    addiu       t3, t3, -128
    addiu       t4, t4, -128
    addiu       t5, t5, -128
    addiu       t6, t6, -128
    addiu       t7, t7, -128
    addiu       t8, t8, -128
    mtc1        t1, f2
    mtc1        t2, f4
    mtc1        t3, f6
    mtc1        t4, f8
    mtc1        t5, f10
    mtc1        t6, f12
    mtc1        t7, f14
    mtc1        t8, f16
    cvt.s.w     f2, f2
    cvt.s.w     f4, f4
    cvt.s.w     f6, f6
    cvt.s.w     f8, f8
    cvt.s.w     f10, f10
    cvt.s.w     f12, f12
    cvt.s.w     f14, f14
    cvt.s.w     f16, f16
    swc1        f2, 224(a2)
    swc1        f4, 228(a2)
    swc1        f6, 232(a2)
    swc1        f8, 236(a2)
    swc1        f10, 240(a2)
    swc1        f12, 244(a2)
    swc1        f14, 248(a2)
    swc1        f16, 252(a2)

    j           ra
     nop

END(jsimd_convsamp_float_dspr2)

#endif

/*****************************************************************************/