// Copyright 2014 Google Inc. All Rights Reserved.

//

// Use of this source code is governed by a BSD-style license

// that can be found in the COPYING file in the root of the source

// tree. An additional intellectual property rights grant can be found

// in the file PATENTS. All contributing project authors may

// be found in the AUTHORS file in the root of the source tree.

// -----------------------------------------------------------------------------

//

// MIPS version of dsp functions

//

// Author(s):  Djordje Pesut    (djordje.pesut@imgtec.com)

//             Jovan Zelincevic (jovan.zelincevic@imgtec.com)

#include "src/dsp/dsp.h"

#if defined(WEBP_USE_MIPS32)

#include "src/dsp/mips_macro.h"

static const int kC1 = 20091 + (1 << 16);

static const int kC2 = 35468;

static WEBP_INLINE int abs_mips32(int x) {

  const int sign = x >> 31;

  return (x ^ sign) - sign;

}

// 4 pixels in, 2 pixels out

static WEBP_INLINE void do_filter2(uint8_t* p, int step) {

  const int p1 = p[-2 * step], p0 = p[-step], q0 = p[0], q1 = p[step];

  const int a = 3 * (q0 - p0) + VP8ksclip1[p1 - q1];

  const int a1 = VP8ksclip2[(a + 4) >> 3];

  const int a2 = VP8ksclip2[(a + 3) >> 3];

  p[-step] = VP8kclip1[p0 + a2];

  p[    0] = VP8kclip1[q0 - a1];

}

// 4 pixels in, 4 pixels out

static WEBP_INLINE void do_filter4(uint8_t* p, int step) {

  const int p1 = p[-2 * step], p0 = p[-step], q0 = p[0], q1 = p[step];

  const int a = 3 * (q0 - p0);

  const int a1 = VP8ksclip2[(a + 4) >> 3];

  const int a2 = VP8ksclip2[(a + 3) >> 3];

  const int a3 = (a1 + 1) >> 1;

  p[-2 * step] = VP8kclip1[p1 + a3];

  p[-    step] = VP8kclip1[p0 + a2];

  p[        0] = VP8kclip1[q0 - a1];

  p[     step] = VP8kclip1[q1 - a3];

}

// 6 pixels in, 6 pixels out

static WEBP_INLINE void do_filter6(uint8_t* p, int step) {

  const int p2 = p[-3 * step], p1 = p[-2 * step], p0 = p[-step];

  const int q0 = p[0], q1 = p[step], q2 = p[2 * step];

  const int a = VP8ksclip1[3 * (q0 - p0) + VP8ksclip1[p1 - q1]];

  // a is in [-128,127], a1 in [-27,27], a2 in [-18,18] and a3 in [-9,9]

  const int a1 = (27 * a + 63) >> 7;  // eq. to ((3 * a + 7) * 9) >> 7

  const int a2 = (18 * a + 63) >> 7;  // eq. to ((2 * a + 7) * 9) >> 7

  const int a3 = (9  * a + 63) >> 7;  // eq. to ((1 * a + 7) * 9) >> 7

  p[-3 * step] = VP8kclip1[p2 + a3];

  p[-2 * step] = VP8kclip1[p1 + a2];

  p[-    step] = VP8kclip1[p0 + a1];

  p[        0] = VP8kclip1[q0 - a1];

  p[     step] = VP8kclip1[q1 - a2];

  p[ 2 * step] = VP8kclip1[q2 - a3];

}

static WEBP_INLINE int hev(const uint8_t* p, int step, int thresh) {

  const int p1 = p[-2 * step], p0 = p[-step], q0 = p[0], q1 = p[step];

  return (abs_mips32(p1 - p0) > thresh) || (abs_mips32(q1 - q0) > thresh);

}

static WEBP_INLINE int needs_filter(const uint8_t* p, int step, int t) {

  const int p1 = p[-2 * step], p0 = p[-step], q0 = p[0], q1 = p[step];

  return ((4 * abs_mips32(p0 - q0) + abs_mips32(p1 - q1)) <= t);

}

static WEBP_INLINE int needs_filter2(const uint8_t* p,

                                     int step, int t, int it) {

  const int p3 = p[-4 * step], p2 = p[-3 * step];

  const int p1 = p[-2 * step], p0 = p[-step];

  const int q0 = p[0], q1 = p[step], q2 = p[2 * step], q3 = p[3 * step];

  if ((4 * abs_mips32(p0 - q0) + abs_mips32(p1 - q1)) > t) {

    return 0;

  }

  return abs_mips32(p3 - p2) <= it && abs_mips32(p2 - p1) <= it &&

         abs_mips32(p1 - p0) <= it && abs_mips32(q3 - q2) <= it &&

         abs_mips32(q2 - q1) <= it && abs_mips32(q1 - q0) <= it;

}

static WEBP_INLINE void FilterLoop26(uint8_t* p,

                                     int hstride, int vstride, int size,

                                     int thresh, int ithresh, int hev_thresh) {

  const int thresh2 = 2 * thresh + 1;

  while (size-- > 0) {

    if (needs_filter2(p, hstride, thresh2, ithresh)) {

      if (hev(p, hstride, hev_thresh)) {

        do_filter2(p, hstride);

      } else {

        do_filter6(p, hstride);

      }

    }

    p += vstride;

  }

}

static WEBP_INLINE void FilterLoop24(uint8_t* p,

                                     int hstride, int vstride, int size,

                                     int thresh, int ithresh, int hev_thresh) {

  const int thresh2 = 2 * thresh + 1;

  while (size-- > 0) {

    if (needs_filter2(p, hstride, thresh2, ithresh)) {

      if (hev(p, hstride, hev_thresh)) {

        do_filter2(p, hstride);

      } else {

        do_filter4(p, hstride);

      }

    }

    p += vstride;

  }

}

// on macroblock edges

static void VFilter16(uint8_t* p, int stride,

                      int thresh, int ithresh, int hev_thresh) {

  FilterLoop26(p, stride, 1, 16, thresh, ithresh, hev_thresh);

}

static void HFilter16(uint8_t* p, int stride,

                      int thresh, int ithresh, int hev_thresh) {

  FilterLoop26(p, 1, stride, 16, thresh, ithresh, hev_thresh);

}

// 8-pixels wide variant, for chroma filtering

static void VFilter8(uint8_t* u, uint8_t* v, int stride,

                     int thresh, int ithresh, int hev_thresh) {

  FilterLoop26(u, stride, 1, 8, thresh, ithresh, hev_thresh);

  FilterLoop26(v, stride, 1, 8, thresh, ithresh, hev_thresh);

}

static void HFilter8(uint8_t* u, uint8_t* v, int stride,

                     int thresh, int ithresh, int hev_thresh) {

  FilterLoop26(u, 1, stride, 8, thresh, ithresh, hev_thresh);

  FilterLoop26(v, 1, stride, 8, thresh, ithresh, hev_thresh);

}

static void VFilter8i(uint8_t* u, uint8_t* v, int stride,

                      int thresh, int ithresh, int hev_thresh) {

  FilterLoop24(u + 4 * stride, stride, 1, 8, thresh, ithresh, hev_thresh);

  FilterLoop24(v + 4 * stride, stride, 1, 8, thresh, ithresh, hev_thresh);

}

static void HFilter8i(uint8_t* u, uint8_t* v, int stride,

                      int thresh, int ithresh, int hev_thresh) {

  FilterLoop24(u + 4, 1, stride, 8, thresh, ithresh, hev_thresh);

  FilterLoop24(v + 4, 1, stride, 8, thresh, ithresh, hev_thresh);

}

// on three inner edges

static void VFilter16i(uint8_t* p, int stride,

                       int thresh, int ithresh, int hev_thresh) {

  int k;

  for (k = 3; k > 0; --k) {

    p += 4 * stride;

    FilterLoop24(p, stride, 1, 16, thresh, ithresh, hev_thresh);

  }

}

static void HFilter16i(uint8_t* p, int stride,

                       int thresh, int ithresh, int hev_thresh) {

  int k;

  for (k = 3; k > 0; --k) {

    p += 4;

    FilterLoop24(p, 1, stride, 16, thresh, ithresh, hev_thresh);

  }

}

//------------------------------------------------------------------------------

// Simple In-loop filtering (Paragraph 15.2)

static void SimpleVFilter16(uint8_t* p, int stride, int thresh) {

  int i;

  const int thresh2 = 2 * thresh + 1;

  for (i = 0; i < 16; ++i) {

    if (needs_filter(p + i, stride, thresh2)) {

      do_filter2(p + i, stride);

    }

  }

}

static void SimpleHFilter16(uint8_t* p, int stride, int thresh) {

  int i;

  const int thresh2 = 2 * thresh + 1;

  for (i = 0; i < 16; ++i) {

    if (needs_filter(p + i * stride, 1, thresh2)) {

      do_filter2(p + i * stride, 1);

    }

  }

}

static void SimpleVFilter16i(uint8_t* p, int stride, int thresh) {

  int k;

  for (k = 3; k > 0; --k) {

    p += 4 * stride;

    SimpleVFilter16(p, stride, thresh);

  }

}

static void SimpleHFilter16i(uint8_t* p, int stride, int thresh) {

  int k;

  for (k = 3; k > 0; --k) {

    p += 4;

    SimpleHFilter16(p, stride, thresh);

  }

}

static void TransformOne(const int16_t* in, uint8_t* dst) {

  int temp0, temp1, temp2, temp3, temp4;

  int temp5, temp6, temp7, temp8, temp9;

  int temp10, temp11, temp12, temp13, temp14;

  int temp15, temp16, temp17, temp18;

  int16_t* p_in = (int16_t*)in;

  // loops unrolled and merged to avoid usage of tmp buffer

  // and to reduce number of stalls. MUL macro is written

  // in assembler and inlined

  __asm__ volatile(

    "lh       %[temp0],  0(%[in])                      \n\t"

    "lh       %[temp8],  16(%[in])                     \n\t"

    "lh       %[temp4],  8(%[in])                      \n\t"

    "lh       %[temp12], 24(%[in])                     \n\t"

    "addu     %[temp16], %[temp0],  %[temp8]           \n\t"

    "subu     %[temp0],  %[temp0],  %[temp8]           \n\t"

    "mul      %[temp8],  %[temp4],  %[kC2]             \n\t"

    "mul      %[temp17], %[temp12], %[kC1]             \n\t"

    "mul      %[temp4],  %[temp4],  %[kC1]             \n\t"

    "mul      %[temp12], %[temp12], %[kC2]             \n\t"

    "lh       %[temp1],  2(%[in])                      \n\t"

    "lh       %[temp5],  10(%[in])                     \n\t"

    "lh       %[temp9],  18(%[in])                     \n\t"

    "lh       %[temp13], 26(%[in])                     \n\t"

    "sra      %[temp8],  %[temp8],  16                 \n\t"

    "sra      %[temp17], %[temp17], 16                 \n\t"

    "sra      %[temp4],  %[temp4],  16                 \n\t"

    "sra      %[temp12], %[temp12], 16                 \n\t"

    "lh       %[temp2],  4(%[in])                      \n\t"

    "lh       %[temp6],  12(%[in])                     \n\t"

    "lh       %[temp10], 20(%[in])                     \n\t"

    "lh       %[temp14], 28(%[in])                     \n\t"

    "subu     %[temp17], %[temp8],  %[temp17]          \n\t"

    "addu     %[temp4],  %[temp4],  %[temp12]          \n\t"

    "addu     %[temp8],  %[temp16], %[temp4]           \n\t"

    "subu     %[temp4],  %[temp16], %[temp4]           \n\t"

    "addu     %[temp16], %[temp1],  %[temp9]           \n\t"

    "subu     %[temp1],  %[temp1],  %[temp9]           \n\t"

    "lh       %[temp3],  6(%[in])                      \n\t"

    "lh       %[temp7],  14(%[in])                     \n\t"

    "lh       %[temp11], 22(%[in])                     \n\t"

    "lh       %[temp15], 30(%[in])                     \n\t"

    "addu     %[temp12], %[temp0],  %[temp17]          \n\t"

    "subu     %[temp0],  %[temp0],  %[temp17]          \n\t"

    "mul      %[temp9],  %[temp5],  %[kC2]             \n\t"

    "mul      %[temp17], %[temp13], %[kC1]             \n\t"

    "mul      %[temp5],  %[temp5],  %[kC1]             \n\t"

    "mul      %[temp13], %[temp13], %[kC2]             \n\t"

    "sra      %[temp9],  %[temp9],  16                 \n\t"

    "sra      %[temp17], %[temp17], 16                 \n\t"

    "subu     %[temp17], %[temp9],  %[temp17]          \n\t"

    "sra      %[temp5],  %[temp5],  16                 \n\t"

    "sra      %[temp13], %[temp13], 16                 \n\t"

    "addu     %[temp5],  %[temp5],  %[temp13]          \n\t"

    "addu     %[temp13], %[temp1],  %[temp17]          \n\t"

    "subu     %[temp1],  %[temp1],  %[temp17]          \n\t"

    "mul      %[temp17], %[temp14], %[kC1]             \n\t"

    "mul      %[temp14], %[temp14], %[kC2]             \n\t"

    "addu     %[temp9],  %[temp16], %[temp5]           \n\t"

    "subu     %[temp5],  %[temp16], %[temp5]           \n\t"

    "addu     %[temp16], %[temp2],  %[temp10]          \n\t"

    "subu     %[temp2],  %[temp2],  %[temp10]          \n\t"

    "mul      %[temp10], %[temp6],  %[kC2]             \n\t"

    "mul      %[temp6],  %[temp6],  %[kC1]             \n\t"

    "sra      %[temp17], %[temp17], 16                 \n\t"

    "sra      %[temp14], %[temp14], 16                 \n\t"

    "sra      %[temp10], %[temp10], 16                 \n\t"

    "sra      %[temp6],  %[temp6],  16                 \n\t"

    "subu     %[temp17], %[temp10], %[temp17]          \n\t"

    "addu     %[temp6],  %[temp6],  %[temp14]          \n\t"

    "addu     %[temp10], %[temp16], %[temp6]           \n\t"

    "subu     %[temp6],  %[temp16], %[temp6]           \n\t"

    "addu     %[temp14], %[temp2],  %[temp17]          \n\t"

    "subu     %[temp2],  %[temp2],  %[temp17]          \n\t"

    "mul      %[temp17], %[temp15], %[kC1]             \n\t"

    "mul      %[temp15], %[temp15], %[kC2]             \n\t"

    "addu     %[temp16], %[temp3],  %[temp11]          \n\t"

    "subu     %[temp3],  %[temp3],  %[temp11]          \n\t"

    "mul      %[temp11], %[temp7],  %[kC2]             \n\t"

    "mul      %[temp7],  %[temp7],  %[kC1]             \n\t"

    "addiu    %[temp8],  %[temp8],  4                  \n\t"

    "addiu    %[temp12], %[temp12], 4                  \n\t"

    "addiu    %[temp0],  %[temp0],  4                  \n\t"

    "addiu    %[temp4],  %[temp4],  4                  \n\t"

    "sra      %[temp17], %[temp17], 16                 \n\t"

    "sra      %[temp15], %[temp15], 16                 \n\t"

    "sra      %[temp11], %[temp11], 16                 \n\t"

    "sra      %[temp7],  %[temp7],  16                 \n\t"

    "subu     %[temp17], %[temp11], %[temp17]          \n\t"

    "addu     %[temp7],  %[temp7],  %[temp15]          \n\t"

    "addu     %[temp15], %[temp3],  %[temp17]          \n\t"

    "subu     %[temp3],  %[temp3],  %[temp17]          \n\t"

    "addu     %[temp11], %[temp16], %[temp7]           \n\t"

    "subu     %[temp7],  %[temp16], %[temp7]           \n\t"

    "addu     %[temp16], %[temp8],  %[temp10]          \n\t"

    "subu     %[temp8],  %[temp8],  %[temp10]          \n\t"

    "mul      %[temp10], %[temp9],  %[kC2]             \n\t"

    "mul      %[temp17], %[temp11], %[kC1]             \n\t"

    "mul      %[temp9],  %[temp9],  %[kC1]             \n\t"

    "mul      %[temp11], %[temp11], %[kC2]             \n\t"

    "sra      %[temp10], %[temp10], 16                 \n\t"

    "sra      %[temp17], %[temp17], 16                 \n\t"

    "sra      %[temp9],  %[temp9],  16                 \n\t"

    "sra      %[temp11], %[temp11], 16                 \n\t"

    "subu     %[temp17], %[temp10], %[temp17]          \n\t"

    "addu     %[temp11], %[temp9],  %[temp11]          \n\t"

    "addu     %[temp10], %[temp12], %[temp14]          \n\t"

    "subu     %[temp12], %[temp12], %[temp14]          \n\t"

    "mul      %[temp14], %[temp13], %[kC2]             \n\t"

    "mul      %[temp9],  %[temp15], %[kC1]             \n\t"

    "mul      %[temp13], %[temp13], %[kC1]             \n\t"

    "mul      %[temp15], %[temp15], %[kC2]             \n\t"

    "sra      %[temp14], %[temp14], 16                 \n\t"

    "sra      %[temp9],  %[temp9],  16                 \n\t"

    "sra      %[temp13], %[temp13], 16                 \n\t"

    "sra      %[temp15], %[temp15], 16                 \n\t"

    "subu     %[temp9],  %[temp14], %[temp9]           \n\t"

    "addu     %[temp15], %[temp13], %[temp15]          \n\t"

    "addu     %[temp14], %[temp0],  %[temp2]           \n\t"

    "subu     %[temp0],  %[temp0],  %[temp2]           \n\t"

    "mul      %[temp2],  %[temp1],  %[kC2]             \n\t"

    "mul      %[temp13], %[temp3],  %[kC1]             \n\t"

    "mul      %[temp1],  %[temp1],  %[kC1]             \n\t"

    "mul      %[temp3],  %[temp3],  %[kC2]             \n\t"

    "sra      %[temp2],  %[temp2],  16                 \n\t"

    "sra      %[temp13], %[temp13], 16                 \n\t"

    "sra      %[temp1],  %[temp1],  16                 \n\t"

    "sra      %[temp3],  %[temp3],  16                 \n\t"

    "subu     %[temp13], %[temp2],  %[temp13]          \n\t"

    "addu     %[temp3],  %[temp1],  %[temp3]           \n\t"

    "addu     %[temp2],  %[temp4],  %[temp6]           \n\t"

    "subu     %[temp4],  %[temp4],  %[temp6]           \n\t"

    "mul      %[temp6],  %[temp5],  %[kC2]             \n\t"

    "mul      %[temp1],  %[temp7],  %[kC1]             \n\t"

    "mul      %[temp5],  %[temp5],  %[kC1]             \n\t"

    "mul      %[temp7],  %[temp7],  %[kC2]             \n\t"

    "sra      %[temp6],  %[temp6],  16                 \n\t"

    "sra      %[temp1],  %[temp1],  16                 \n\t"

    "sra      %[temp5],  %[temp5],  16                 \n\t"

    "sra      %[temp7],  %[temp7],  16                 \n\t"

    "subu     %[temp1],  %[temp6],  %[temp1]           \n\t"

    "addu     %[temp7],  %[temp5],  %[temp7]           \n\t"

    "addu     %[temp5],  %[temp16], %[temp11]          \n\t"

    "subu     %[temp16], %[temp16], %[temp11]          \n\t"

    "addu     %[temp11], %[temp8],  %[temp17]          \n\t"

    "subu     %[temp8],  %[temp8],  %[temp17]          \n\t"

    "sra      %[temp5],  %[temp5],  3                  \n\t"

    "sra      %[temp16], %[temp16], 3                  \n\t"

    "sra      %[temp11], %[temp11], 3                  \n\t"

    "sra      %[temp8],  %[temp8],  3                  \n\t"

    "addu     %[temp17], %[temp10], %[temp15]          \n\t"

    "subu     %[temp10], %[temp10], %[temp15]          \n\t"

    "addu     %[temp15], %[temp12], %[temp9]           \n\t"

    "subu     %[temp12], %[temp12], %[temp9]           \n\t"

    "sra      %[temp17], %[temp17], 3                  \n\t"

    "sra      %[temp10], %[temp10], 3                  \n\t"

    "sra      %[temp15], %[temp15], 3                  \n\t"

    "sra      %[temp12], %[temp12], 3                  \n\t"

    "addu     %[temp9],  %[temp14], %[temp3]           \n\t"

    "subu     %[temp14], %[temp14], %[temp3]           \n\t"

    "addu     %[temp3],  %[temp0],  %[temp13]          \n\t"

    "subu     %[temp0],  %[temp0],  %[temp13]          \n\t"

    "sra      %[temp9],  %[temp9],  3                  \n\t"

    "sra      %[temp14], %[temp14], 3                  \n\t"

    "sra      %[temp3],  %[temp3],  3                  \n\t"

    "sra      %[temp0],  %[temp0],  3                  \n\t"

    "addu     %[temp13], %[temp2],  %[temp7]           \n\t"

    "subu     %[temp2],  %[temp2],  %[temp7]           \n\t"

    "addu     %[temp7],  %[temp4],  %[temp1]           \n\t"

    "subu     %[temp4],  %[temp4],  %[temp1]           \n\t"

    "sra      %[temp13], %[temp13], 3                  \n\t"

    "sra      %[temp2],  %[temp2],  3                  \n\t"

    "sra      %[temp7],  %[temp7],  3                  \n\t"

    "sra      %[temp4],  %[temp4],  3                  \n\t"

    "addiu    %[temp6],  $zero,     255                \n\t"

    "lbu      %[temp1],  0+0*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp1],  %[temp1],  %[temp5]           \n\t"

    "sra      %[temp5],  %[temp1],  8                  \n\t"

    "sra      %[temp18], %[temp1],  31                 \n\t"

    "beqz     %[temp5],  1f                            \n\t"

    "xor      %[temp1],  %[temp1],  %[temp1]           \n\t"

    "movz     %[temp1],  %[temp6],  %[temp18]          \n\t"

  "1:                                                  \n\t"

    "lbu      %[temp18], 1+0*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp1],  0+0*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp18], %[temp18], %[temp11]          \n\t"

    "sra      %[temp11], %[temp18], 8                  \n\t"

    "sra      %[temp1],  %[temp18], 31                 \n\t"

    "beqz     %[temp11], 2f                            \n\t"

    "xor      %[temp18], %[temp18], %[temp18]          \n\t"

    "movz     %[temp18], %[temp6],  %[temp1]           \n\t"

  "2:                                                  \n\t"

    "lbu      %[temp1],  2+0*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp18], 1+0*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp1],  %[temp1],  %[temp8]           \n\t"

    "sra      %[temp8],  %[temp1],  8                  \n\t"

    "sra      %[temp18], %[temp1],  31                 \n\t"

    "beqz     %[temp8],  3f                            \n\t"

    "xor      %[temp1],  %[temp1],  %[temp1]           \n\t"

    "movz     %[temp1],  %[temp6],  %[temp18]          \n\t"

  "3:                                                  \n\t"

    "lbu      %[temp18], 3+0*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp1],  2+0*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp18], %[temp18], %[temp16]          \n\t"

    "sra      %[temp16], %[temp18], 8                  \n\t"

    "sra      %[temp1],  %[temp18], 31                 \n\t"

    "beqz     %[temp16], 4f                            \n\t"

    "xor      %[temp18], %[temp18], %[temp18]          \n\t"

    "movz     %[temp18], %[temp6],  %[temp1]           \n\t"

  "4:                                                  \n\t"

    "sb       %[temp18], 3+0*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp5],  0+1*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp8],  1+1*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp11], 2+1*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp16], 3+1*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp5],  %[temp5],  %[temp17]          \n\t"

    "addu     %[temp8],  %[temp8],  %[temp15]          \n\t"

    "addu     %[temp11], %[temp11], %[temp12]          \n\t"

    "addu     %[temp16], %[temp16], %[temp10]          \n\t"

    "sra      %[temp18], %[temp5],  8                  \n\t"

    "sra      %[temp1],  %[temp5],  31                 \n\t"

    "beqz     %[temp18], 5f                            \n\t"

    "xor      %[temp5],  %[temp5],  %[temp5]           \n\t"

    "movz     %[temp5],  %[temp6],  %[temp1]           \n\t"

  "5:                                                  \n\t"

    "sra      %[temp18], %[temp8],  8                  \n\t"

    "sra      %[temp1],  %[temp8],  31                 \n\t"

    "beqz     %[temp18], 6f                            \n\t"

    "xor      %[temp8],  %[temp8],  %[temp8]           \n\t"

    "movz     %[temp8],  %[temp6],  %[temp1]           \n\t"

  "6:                                                  \n\t"

    "sra      %[temp18], %[temp11], 8                  \n\t"

    "sra      %[temp1],  %[temp11], 31                 \n\t"

    "sra      %[temp17], %[temp16], 8                  \n\t"

    "sra      %[temp15], %[temp16], 31                 \n\t"

    "beqz     %[temp18], 7f                            \n\t"

    "xor      %[temp11], %[temp11], %[temp11]          \n\t"

    "movz     %[temp11], %[temp6],  %[temp1]           \n\t"

  "7:                                                  \n\t"

    "beqz     %[temp17], 8f                            \n\t"

    "xor      %[temp16], %[temp16], %[temp16]          \n\t"

    "movz     %[temp16], %[temp6],  %[temp15]          \n\t"

  "8:                                                  \n\t"

    "sb       %[temp5],  0+1*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp8],  1+1*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp11], 2+1*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp16], 3+1*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp5],  0+2*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp8],  1+2*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp11], 2+2*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp16], 3+2*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp5],  %[temp5],  %[temp9]           \n\t"

    "addu     %[temp8],  %[temp8],  %[temp3]           \n\t"

    "addu     %[temp11], %[temp11], %[temp0]           \n\t"

    "addu     %[temp16], %[temp16], %[temp14]          \n\t"

    "sra      %[temp18], %[temp5],  8                  \n\t"

    "sra      %[temp1],  %[temp5],  31                 \n\t"

    "sra      %[temp17], %[temp8],  8                  \n\t"

    "sra      %[temp15], %[temp8],  31                 \n\t"

    "sra      %[temp12], %[temp11], 8                  \n\t"

    "sra      %[temp10], %[temp11], 31                 \n\t"

    "sra      %[temp9],  %[temp16], 8                  \n\t"

    "sra      %[temp3],  %[temp16], 31                 \n\t"

    "beqz     %[temp18], 9f                            \n\t"

    "xor      %[temp5],  %[temp5],  %[temp5]           \n\t"

    "movz     %[temp5],  %[temp6],  %[temp1]           \n\t"

  "9:                                                  \n\t"

    "beqz     %[temp17], 10f                           \n\t"

    "xor      %[temp8],  %[temp8],  %[temp8]           \n\t"

    "movz     %[temp8],  %[temp6],  %[temp15]          \n\t"

  "10:                                                 \n\t"

    "beqz     %[temp12], 11f                           \n\t"

    "xor      %[temp11], %[temp11], %[temp11]          \n\t"

    "movz     %[temp11], %[temp6],  %[temp10]          \n\t"

  "11:                                                 \n\t"

    "beqz     %[temp9],  12f                           \n\t"

    "xor      %[temp16], %[temp16], %[temp16]          \n\t"

    "movz     %[temp16], %[temp6],  %[temp3]           \n\t"

  "12:                                                 \n\t"

    "sb       %[temp5],  0+2*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp8],  1+2*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp11], 2+2*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp16], 3+2*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp5],  0+3*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp8],  1+3*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp11], 2+3*" XSTR(BPS) "(%[dst])     \n\t"

    "lbu      %[temp16], 3+3*" XSTR(BPS) "(%[dst])     \n\t"

    "addu     %[temp5],  %[temp5],  %[temp13]          \n\t"

    "addu     %[temp8],  %[temp8],  %[temp7]           \n\t"

    "addu     %[temp11], %[temp11], %[temp4]           \n\t"

    "addu     %[temp16], %[temp16], %[temp2]           \n\t"

    "sra      %[temp18], %[temp5],  8                  \n\t"

    "sra      %[temp1],  %[temp5],  31                 \n\t"

    "sra      %[temp17], %[temp8],  8                  \n\t"

    "sra      %[temp15], %[temp8],  31                 \n\t"

    "sra      %[temp12], %[temp11], 8                  \n\t"

    "sra      %[temp10], %[temp11], 31                 \n\t"

    "sra      %[temp9],  %[temp16], 8                  \n\t"

    "sra      %[temp3],  %[temp16], 31                 \n\t"

    "beqz     %[temp18], 13f                           \n\t"

    "xor      %[temp5],  %[temp5],  %[temp5]           \n\t"

    "movz     %[temp5],  %[temp6],  %[temp1]           \n\t"

  "13:                                                 \n\t"

    "beqz     %[temp17], 14f                           \n\t"

    "xor      %[temp8],  %[temp8],  %[temp8]           \n\t"

    "movz     %[temp8],  %[temp6],  %[temp15]          \n\t"

  "14:                                                 \n\t"

    "beqz     %[temp12], 15f                           \n\t"

    "xor      %[temp11], %[temp11], %[temp11]          \n\t"

    "movz     %[temp11], %[temp6],  %[temp10]          \n\t"

  "15:                                                 \n\t"

    "beqz     %[temp9],  16f                           \n\t"

    "xor      %[temp16], %[temp16], %[temp16]          \n\t"

    "movz     %[temp16], %[temp6],  %[temp3]           \n\t"

  "16:                                                 \n\t"

    "sb       %[temp5],  0+3*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp8],  1+3*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp11], 2+3*" XSTR(BPS) "(%[dst])     \n\t"

    "sb       %[temp16], 3+3*" XSTR(BPS) "(%[dst])     \n\t"

    : [temp0]"=&r"(temp0), [temp1]"=&r"(temp1), [temp2]"=&r"(temp2),

      [temp3]"=&r"(temp3), [temp4]"=&r"(temp4), [temp5]"=&r"(temp5),

      [temp6]"=&r"(temp6), [temp7]"=&r"(temp7), [temp8]"=&r"(temp8),

      [temp9]"=&r"(temp9), [temp10]"=&r"(temp10), [temp11]"=&r"(temp11),

      [temp12]"=&r"(temp12), [temp13]"=&r"(temp13), [temp14]"=&r"(temp14),

      [temp15]"=&r"(temp15), [temp16]"=&r"(temp16), [temp17]"=&r"(temp17),

      [temp18]"=&r"(temp18)

    : [in]"r"(p_in), [kC1]"r"(kC1), [kC2]"r"(kC2), [dst]"r"(dst)

    : "memory", "hi", "lo"

  );

}

static void TransformTwo(const int16_t* in, uint8_t* dst, int do_two) {

  TransformOne(in, dst);

  if (do_two) {

    TransformOne(in + 16, dst + 4);

  }

}

//------------------------------------------------------------------------------

// Entry point

extern void VP8DspInitMIPS32(void);

WEBP_TSAN_IGNORE_FUNCTION void VP8DspInitMIPS32(void) {

  VP8InitClipTables();

  VP8Transform = TransformTwo;

  VP8VFilter16 = VFilter16;

  VP8HFilter16 = HFilter16;

  VP8VFilter8 = VFilter8;

  VP8HFilter8 = HFilter8;

  VP8VFilter16i = VFilter16i;

  VP8HFilter16i = HFilter16i;

  VP8VFilter8i = VFilter8i;

  VP8HFilter8i = HFilter8i;

  VP8SimpleVFilter16 = SimpleVFilter16;

  VP8SimpleHFilter16 = SimpleHFilter16;

  VP8SimpleVFilter16i = SimpleVFilter16i;

  VP8SimpleHFilter16i = SimpleHFilter16i;

}

#else  // !WEBP_USE_MIPS32

WEBP_DSP_INIT_STUB(VP8DspInitMIPS32)

#endif  // WEBP_USE_MIPS32