(Patch by Teodora Novkovic <teodora.petrovic@gmail.com>, originally reviewed at https://codereview.chromium.org/14929006/).

Added MIPS DSPr2 optimization for BGRAConvolve2D routine. 

The following routines are optimized: 
- ConvolveVertically 
- ConvolveHorizontally 

Performance gain measured on Malta 74Kc board: 
- our standalone test/bench application ~45% 
- chromium unit_test ~20% 

R=hubbe@chromium.org, senorblanco@chromium.org

Review URL: https://chromiumcodereview.appspot.com/15742005

git-svn-id: svn://svn.chromium.org/chrome/trunk/src@202197 0039d316-1c4b-4281-b951-d872f2087c98

8 files changed, 522 insertions, 5 deletions

diff --git a/AUTHORS b/AUTHORS
index 6d67acc..5518a71 100644
--- a/AUTHORS
+++ b/AUTHORS
@@ -260,3 +260,4 @@ Dongseong Hwang <dongseong.hwang@intel.com>
 Umar Hansa <umar.hansa@gmail.com>
 Radu Stavila <stavila@adobe.com>
 Rijubrata Bhaumik <rijubrata.bhaumik@intel.com>
+Teodora Novkovic <teodora.petrovic@gmail.com>
diff --git a/skia/ext/convolver.cc b/skia/ext/convolver.cc
index a7824aa..4b40ffd 100644
--- a/skia/ext/convolver.cc
+++ b/skia/ext/convolver.cc
@@ -7,6 +7,7 @@
 #include "base/logging.h"
 #include "skia/ext/convolver.h"
 #include "skia/ext/convolver_SSE2.h"
+#include "skia/ext/convolver_mips_dspr2.h"
 #include "third_party/skia/include/core/SkSize.h"
 #include "third_party/skia/include/core/SkTypes.h"
 
@@ -347,7 +348,8 @@ typedef void (*Convolve4RowsHorizontally_pointer)(
 typedef void (*ConvolveHorizontally_pointer)(
     const unsigned char* src_data,
     const ConvolutionFilter1D& filter,
-    unsigned char* out_row);
+    unsigned char* out_row,
+    bool has_alpha);
 
 struct ConvolveProcs {
   // This is how many extra pixels may be read by the
@@ -367,6 +369,10 @@ void SetupSIMD(ConvolveProcs *procs) {
     procs->convolve_4rows_horizontally = &Convolve4RowsHorizontally_SSE2;
     procs->convolve_horizontally = &ConvolveHorizontally_SSE2;
   }
+#elif defined SIMD_MIPS_DSPR2
+  procs->extra_horizontal_reads = 3;
+  procs->convolve_vertically = &ConvolveVertically_mips_dspr2;
+  procs->convolve_horizontally = &ConvolveHorizontally_mips_dspr2;
 #endif
 }
 
@@ -464,7 +470,7 @@ void BGRAConvolve2D(const unsigned char* source_data,
             avoid_simd_rows) {
           simd.convolve_horizontally(
               &source_data[next_x_row * source_byte_row_stride],
-              filter_x, row_buffer.AdvanceRow());
+              filter_x, row_buffer.AdvanceRow(), source_has_alpha);
         } else {
           if (source_has_alpha) {
             ConvolveHorizontally<true>(
diff --git a/skia/ext/convolver.h b/skia/ext/convolver.h
index 4248bc6..dd99a72 100644
--- a/skia/ext/convolver.h
+++ b/skia/ext/convolver.h
@@ -20,6 +20,10 @@
 #define SIMD_PADDING 8  // 8 * int16
 #endif
 
+#if defined (ARCH_CPU_MIPS_FAMILY) && \
+    defined(__mips_dsp) && (__mips_dsp_rev >= 2)
+#define SIMD_MIPS_DSPR2 1
+#endif
 // avoid confusion with Mac OS X's math library (Carbon)
 #if defined(__APPLE__)
 #undef FloatToFixed
diff --git a/skia/ext/convolver_SSE2.cc b/skia/ext/convolver_SSE2.cc
index a823edc..a77a1f4 100644
--- a/skia/ext/convolver_SSE2.cc
+++ b/skia/ext/convolver_SSE2.cc
@@ -16,7 +16,8 @@ namespace skia {
 // |src_data| and continues for the num_values() of the filter.
 void ConvolveHorizontally_SSE2(const unsigned char* src_data,
                                const ConvolutionFilter1D& filter,
-                               unsigned char* out_row) {
+                               unsigned char* out_row,
+                               bool /*has_alpha*/) {
   int num_values = filter.num_values();
 
   int filter_offset, filter_length;
diff --git a/skia/ext/convolver_SSE2.h b/skia/ext/convolver_SSE2.h
index 6b79dae..cf60406 100644
--- a/skia/ext/convolver_SSE2.h
+++ b/skia/ext/convolver_SSE2.h
@@ -20,7 +20,8 @@ void Convolve4RowsHorizontally_SSE2(const unsigned char* src_data[4],
                                     unsigned char* out_row[4]);
 void ConvolveHorizontally_SSE2(const unsigned char* src_data,
                                const ConvolutionFilter1D& filter,
-                               unsigned char* out_row);
+                               unsigned char* out_row,
+                               bool has_alpha);
 }  // namespace skia
 
 #endif  // SKIA_EXT_CONVOLVER_SSE2_H_
diff --git a/skia/ext/convolver_mips_dspr2.cc b/skia/ext/convolver_mips_dspr2.cc
new file mode 100644
index 0000000..955abef
--- /dev/null
+++ b/skia/ext/convolver_mips_dspr2.cc
@@ -0,0 +1,478 @@
+// Copyright (c) 2013 The Chromium Authors. All rights reserved.
+// Use of this source code is governed by a BSD-style license that can be
+// found in the LICENSE file.
+
+#include <algorithm>
+#include "skia/ext/convolver.h"
+#include "skia/ext/convolver_mips_dspr2.h"
+#include "third_party/skia/include/core/SkTypes.h"
+
+namespace skia {
+// Convolves horizontally along a single row. The row data is given in
+// |src_data| and continues for the num_values() of the filter.
+void ConvolveHorizontally_mips_dspr2(const unsigned char* src_data,
+                                     const ConvolutionFilter1D& filter,
+                                     unsigned char* out_row,
+                                     bool has_alpha) {
+#if SIMD_MIPS_DSPR2
+  int row_to_filter = 0;
+  int num_values = filter.num_values();
+  if (has_alpha) {
+    for (int out_x = 0; out_x < num_values; out_x++) {
+      // Get the filter that determines the current output pixel.
+      int filter_offset, filter_length;
+      const ConvolutionFilter1D::Fixed* filter_values =
+        filter.FilterForValue(out_x, &filter_offset, &filter_length);
+      int filter_x = 0;
+
+      __asm__ __volatile__ (
+        ".set push                                  \n"
+        ".set noreorder                             \n"
+
+        "beqz            %[filter_len], 3f          \n"
+        " sll            $t0, %[filter_offset], 2   \n"
+        "addu            %[rtf], %[src_data], $t0   \n"
+        "mtlo            $0, $ac0                   \n"
+        "mtlo            $0, $ac1                   \n"
+        "mtlo            $0, $ac2                   \n"
+        "mtlo            $0, $ac3                   \n"
+        "srl             $t7, %[filter_len], 2      \n"
+        "beqz            $t7, 2f                    \n"
+        " li             %[fx], 0                   \n"
+
+        "11:                                        \n"
+        "addu            $t4, %[filter_val], %[fx]  \n"
+        "sll             $t5, %[fx], 1              \n"
+        "ulw             $t6, 0($t4)                \n" // t6 = |cur[1]|cur[0]|
+        "ulw             $t8, 4($t4)                \n" // t8 = |cur[3]|cur[2]|
+        "addu            $t0, %[rtf], $t5           \n"
+        "lw              $t1, 0($t0)                \n" // t1 = |a0|b0|g0|r0|
+        "lw              $t2, 4($t0)                \n" // t2 = |a1|b1|g1|r1|
+        "lw              $t3, 8($t0)                \n" // t3 = |a2|b2|g2|r2|
+        "lw              $t4, 12($t0)               \n" // t4 = |a3|b3|g3|r3|
+        "precrq.qb.ph    $t0, $t2, $t1              \n" // t0 = |a1|g1|a0|g0|
+        "precr.qb.ph     $t5, $t2, $t1              \n" // t5 = |b1|r1|b0|r0|
+        "preceu.ph.qbla  $t1, $t0                   \n" // t1 = |0|a1|0|a0|
+        "preceu.ph.qbra  $t2, $t0                   \n" // t2 = |0|g1|0|g0|
+        "preceu.ph.qbla  $t0, $t5                   \n" // t0 = |0|b1|0|b0|
+        "preceu.ph.qbra  $t5, $t5                   \n" // t5 = |0|r1|0|r0|
+        "dpa.w.ph        $ac0, $t1, $t6             \n" // ac0+(cur*a1)+(cur*a0)
+        "dpa.w.ph        $ac1, $t0, $t6             \n" // ac1+(cur*b1)+(cur*b0)
+        "dpa.w.ph        $ac2, $t2, $t6             \n" // ac2+(cur*g1)+(cur*g0)
+        "dpa.w.ph        $ac3, $t5, $t6             \n" // ac3+(cur*r1)+(cur*r0)
+        "precrq.qb.ph    $t0, $t4, $t3              \n" // t0 = |a3|g3|a2|g2|
+        "precr.qb.ph     $t5, $t4, $t3              \n" // t5 = |b3|r3|b2|r2|
+        "preceu.ph.qbla  $t1, $t0                   \n" // t1 = |0|a3|0|a2|
+        "preceu.ph.qbra  $t2, $t0                   \n" // t2 = |0|g3|0|g2|
+        "preceu.ph.qbla  $t0, $t5                   \n" // t0 = |0|b3|0|b2|
+        "preceu.ph.qbra  $t5, $t5                   \n" // t5 = |0|r3|0|r2|
+        "dpa.w.ph        $ac0, $t1, $t8             \n" // ac0+(cur*a3)+(cur*a2)
+        "dpa.w.ph        $ac1, $t0, $t8             \n" // ac1+(cur*b3)+(cur*b2)
+        "dpa.w.ph        $ac2, $t2, $t8             \n" // ac2+(cur*g3)+(cur*g2)
+        "dpa.w.ph        $ac3, $t5, $t8             \n" // ac3+(cur*r3)+(cur*r2)
+        "addiu           $t7, $t7, -1               \n"
+        "bgtz            $t7, 11b                   \n"
+        " addiu          %[fx], %[fx], 8            \n"
+
+        "2:                                         \n"
+        "andi            $t7, %[filter_len], 0x3    \n" // residual
+        "beqz            $t7, 3f                    \n"
+        " nop                                       \n"
+
+        "21:                                        \n"
+        "sll             $t1, %[fx], 1              \n"
+        "addu            $t2, %[filter_val], %[fx]  \n"
+        "addu            $t0, %[rtf], $t1           \n"
+        "lh              $t6, 0($t2)                \n" // t6 = filter_val[fx]
+        "lbu             $t1, 0($t0)                \n" // t1 = row[fx * 4 + 0]
+        "lbu             $t2, 1($t0)                \n" // t2 = row[fx * 4 + 1]
+        "lbu             $t3, 2($t0)                \n" // t3 = row[fx * 4 + 2]
+        "lbu             $t4, 3($t0)                \n" // t4 = row[fx * 4 + 2]
+        "maddu           $ac3, $t6, $t1             \n"
+        "maddu           $ac2, $t6, $t2             \n"
+        "maddu           $ac1, $t6, $t3             \n"
+        "maddu           $ac0, $t6, $t4             \n"
+        "addiu           $t7, $t7, -1               \n"
+        "bgtz            $t7, 21b                   \n"
+        " addiu          %[fx], %[fx], 2            \n"
+
+        "3:                                         \n"
+        "extrv.w         $t0, $ac0, %[kShiftBits]   \n" // a >> kShiftBits
+        "extrv.w         $t1, $ac1, %[kShiftBits]   \n" // b >> kShiftBits
+        "extrv.w         $t2, $ac2, %[kShiftBits]   \n" // g >> kShiftBits
+        "extrv.w         $t3, $ac3, %[kShiftBits]   \n" // r >> kShiftBits
+        "sll             $t5, %[out_x], 2           \n"
+        "repl.ph         $t6, 128                   \n" // t6 = | 128 | 128 |
+        "addu            $t5, %[out_row], $t5       \n"
+        "append          $t2, $t3, 16               \n"
+        "append          $t0, $t1, 16               \n"
+        "subu.ph         $t1, $t0, $t6              \n"
+        "shll_s.ph       $t1, $t1, 8                \n"
+        "shra.ph         $t1, $t1, 8                \n"
+        "addu.ph         $t1, $t1, $t6              \n"
+        "subu.ph         $t3, $t2, $t6              \n"
+        "shll_s.ph       $t3, $t3, 8                \n"
+        "shra.ph         $t3, $t3, 8                \n"
+        "addu.ph         $t3, $t3, $t6              \n"
+        "precr.qb.ph     $t0, $t1, $t3              \n"
+        "usw             $t0, 0($t5)                \n"
+
+        ".set pop                                   \n"
+      : [fx] "+r" (filter_x), [out_x] "+r" (out_x), [out_row] "+r" (out_row),
+        [rtf] "+r" (row_to_filter)
+      : [filter_val] "r" (filter_values), [filter_len] "r" (filter_length),
+        [kShiftBits] "r" (ConvolutionFilter1D::kShiftBits),
+        [filter_offset] "r" (filter_offset), [src_data] "r" (src_data)
+      : "lo", "hi", "$ac1lo", "$ac1hi", "$ac2lo", "$ac2hi", "$ac3lo", "$ac3hi",
+        "t0", "t1", "t2", "t3", "t4", "t5", "t6", "t7", "t8"
+      );
+    }
+  } else {
+    for (int out_x = 0; out_x < num_values; out_x++) {
+      // Get the filter that determines the current output pixel.
+      int filter_offset, filter_length;
+      const ConvolutionFilter1D::Fixed* filter_values =
+        filter.FilterForValue(out_x, &filter_offset, &filter_length);
+      int filter_x = 0;
+      __asm__ __volatile__ (
+        ".set push                                  \n"
+        ".set noreorder                             \n"
+
+        "beqz            %[filter_len], 3f          \n"
+        " sll            $t0, %[filter_offset], 2   \n"
+        "addu            %[rtf], %[src_data], $t0   \n"
+        "mtlo            $0, $ac1                   \n"
+        "mtlo            $0, $ac2                   \n"
+        "mtlo            $0, $ac3                   \n"
+        "srl             $t7, %[filter_len], 2      \n"
+        "beqz            $t7, 2f                    \n"
+        " li             %[fx], 0                   \n"
+
+        "11:                                        \n"
+        "addu            $t4, %[filter_val], %[fx]  \n"
+        "sll             $t5, %[fx], 1              \n"
+        "ulw             $t6, 0($t4)                \n" // t6 = |cur[1]|cur[0]|
+        "ulw             $t8, 4($t4)                \n" // t8 = |cur[3]|cur[2]|
+        "addu            $t0, %[rtf], $t5           \n"
+        "lw              $t1, 0($t0)                \n" // t1 = |a0|b0|g0|r0|
+        "lw              $t2, 4($t0)                \n" // t2 = |a1|b1|g1|r1|
+        "lw              $t3, 8($t0)                \n" // t3 = |a2|b2|g2|r2|
+        "lw              $t4, 12($t0)               \n" // t4 = |a3|b3|g3|r3|
+        "precrq.qb.ph    $t0, $t2, $t1              \n" // t0 = |a1|g1|a0|g0|
+        "precr.qb.ph     $t5, $t2, $t1              \n" // t5 = |b1|r1|b0|r0|
+        "preceu.ph.qbra  $t2, $t0                   \n" // t2 = |0|g1|0|g0|
+        "preceu.ph.qbla  $t0, $t5                   \n" // t0 = |0|b1|0|b0|
+        "preceu.ph.qbra  $t5, $t5                   \n" // t5 = |0|r1|0|r0|
+        "dpa.w.ph        $ac1, $t0, $t6             \n" // ac1+(cur*b1)+(cur*b0)
+        "dpa.w.ph        $ac2, $t2, $t6             \n" // ac2+(cur*g1)+(cur*g0)
+        "dpa.w.ph        $ac3, $t5, $t6             \n" // ac3+(cur*r1)+(cur*r0)
+        "precrq.qb.ph    $t0, $t4, $t3              \n" // t0 = |a3|g3|a2|g2|
+        "precr.qb.ph     $t5, $t4, $t3              \n" // t5 = |b3|r3|b2|r2|
+        "preceu.ph.qbra  $t2, $t0                   \n" // t2 = |0|g3|0|g2|
+        "preceu.ph.qbla  $t0, $t5                   \n" // t0 = |0|b3|0|b2|
+        "preceu.ph.qbra  $t5, $t5                   \n" // t5 = |0|r3|0|r2|
+        "dpa.w.ph        $ac1, $t0, $t8             \n" // ac1+(cur*b3)+(cur*b2)
+        "dpa.w.ph        $ac2, $t2, $t8             \n" // ac2+(cur*g3)+(cur*g2)
+        "dpa.w.ph        $ac3, $t5, $t8             \n" // ac3+(cur*r3)+(cur*r2)
+        "addiu           $t7, $t7, -1               \n"
+        "bgtz            $t7, 11b                   \n"
+        " addiu          %[fx], %[fx], 8            \n"
+
+        "2:                                         \n"
+        "andi            $t7, %[filter_len], 0x3    \n" // residual
+        "beqz            $t7, 3f                    \n"
+        " nop                                       \n"
+
+        "21:                                        \n"
+        "sll             $t1, %[fx], 1              \n"
+        "addu            $t2, %[filter_val], %[fx]  \n"
+        "addu            $t0, %[rtf], $t1           \n"
+        "lh              $t6, 0($t2)                \n" // t6 = filter_val[fx]
+        "lbu             $t1, 0($t0)                \n" // t1 = row[fx * 4 + 0]
+        "lbu             $t2, 1($t0)                \n" // t2 = row[fx * 4 + 1]
+        "lbu             $t3, 2($t0)                \n" // t3 = row[fx * 4 + 2]
+        "maddu           $ac3, $t6, $t1             \n"
+        "maddu           $ac2, $t6, $t2             \n"
+        "maddu           $ac1, $t6, $t3             \n"
+        "addiu           $t7, $t7, -1               \n"
+        "bgtz            $t7, 21b                   \n"
+        " addiu          %[fx], %[fx], 2            \n"
+
+        "3:                                         \n"
+        "extrv.w         $t1, $ac1, %[kShiftBits]   \n" // b >> kShiftBits
+        "extrv.w         $t2, $ac2, %[kShiftBits]   \n" // g >> kShiftBits
+        "extrv.w         $t3, $ac3, %[kShiftBits]   \n" // r >> kShiftBits
+        "repl.ph         $t6, 128                   \n" // t6 = | 128 | 128 |
+        "sll             $t8, %[out_x], 2           \n"
+        "addu            $t8, %[out_row], $t8       \n"
+        "append          $t2, $t3, 16               \n"
+        "andi            $t1, 0xFFFF                \n"
+        "subu.ph         $t5, $t1, $t6              \n"
+        "shll_s.ph       $t5, $t5, 8                \n"
+        "shra.ph         $t5, $t5, 8                \n"
+        "addu.ph         $t5, $t5, $t6              \n"
+        "subu.ph         $t4, $t2, $t6              \n"
+        "shll_s.ph       $t4, $t4, 8                \n"
+        "shra.ph         $t4, $t4, 8                \n"
+        "addu.ph         $t4, $t4, $t6              \n"
+        "precr.qb.ph     $t0, $t5, $t4              \n"
+        "usw             $t0, 0($t8)                \n"
+
+        ".set pop                                   \n"
+      : [fx] "+r" (filter_x), [out_x] "+r" (out_x), [out_row] "+r" (out_row),
+        [rtf] "+r" (row_to_filter)
+      : [filter_val] "r" (filter_values), [filter_len] "r" (filter_length),
+        [kShiftBits] "r" (ConvolutionFilter1D::kShiftBits),
+        [filter_offset] "r" (filter_offset), [src_data] "r" (src_data)
+      : "lo", "hi", "$ac1lo", "$ac1hi", "$ac2lo", "$ac2hi", "$ac3lo", "$ac3hi",
+        "t0", "t1", "t2", "t3", "t4", "t5", "t6", "t7", "t8"
+      );
+    }
+  }
+#endif
+}
+void ConvolveVertically_mips_dspr2(const ConvolutionFilter1D::Fixed* filter_val,
+                                   int filter_length,
+                                   unsigned char* const* source_data_rows,
+                                   int pixel_width,
+                                   unsigned char* out_row,
+                                   bool has_alpha) {
+#if SIMD_MIPS_DSPR2
+  // We go through each column in the output and do a vertical convolution,
+  // generating one output pixel each time.
+  int byte_offset;
+  int cnt;
+  int filter_y;
+  if (has_alpha) {
+    for (int out_x = 0; out_x < pixel_width; out_x++) {
+      __asm__ __volatile__ (
+        ".set push                                   \n"
+        ".set noreorder                              \n"
+
+        "beqz            %[filter_len], 3f           \n"
+        " sll            %[offset], %[out_x], 2      \n"
+        "mtlo            $0, $ac0                    \n"
+        "mtlo            $0, $ac1                    \n"
+        "mtlo            $0, $ac2                    \n"
+        "mtlo            $0, $ac3                    \n"
+        "srl             %[cnt], %[filter_len], 2    \n"
+        "beqz            %[cnt], 2f                  \n"
+        " li             %[fy], 0                    \n"
+
+        "11:                                         \n"
+        "sll             $t1, %[fy], 1               \n"
+        "addu            $t0, %[src_data_rows], $t1  \n"
+        "lw              $t1, 0($t0)                 \n"
+        "lw              $t2, 4($t0)                 \n"
+        "lw              $t3, 8($t0)                 \n"
+        "lw              $t4, 12($t0)                \n"
+        "addu            $t1, $t1, %[offset]         \n"
+        "addu            $t2, $t2, %[offset]         \n"
+        "addu            $t3, $t3, %[offset]         \n"
+        "addu            $t4, $t4, %[offset]         \n"
+        "lw              $t1, 0($t1)                 \n" // t1 = |a0|b0|g0|r0|
+        "lw              $t2, 0($t2)                 \n" // t2 = |a1|b1|g1|r1|
+        "lw              $t3, 0($t3)                 \n" // t3 = |a0|b0|g0|r0|
+        "lw              $t4, 0($t4)                 \n" // t4 = |a1|b1|g1|r1|
+        "precrq.qb.ph    $t5, $t2, $t1               \n" // t5 = |a1|g1|a0|g0|
+        "precr.qb.ph     $t6, $t2, $t1               \n" // t6 = |b1|r1|b0|r0|
+        "preceu.ph.qbla  $t0, $t5                    \n" // t0 = |0|a1|0|a0|
+        "preceu.ph.qbra  $t1, $t5                    \n" // t1 = |0|g1|0|g0|
+        "preceu.ph.qbla  $t2, $t6                    \n" // t2 = |0|b1|0|b0|
+        "preceu.ph.qbra  $t5, $t6                    \n" // t5 = |0|r1|0|r0|
+        "addu            $t6, %[filter_val], %[fy]   \n"
+        "ulw             $t7, 0($t6)                 \n" // t7 = |cur_1|cur_0|
+        "ulw             $t6, 4($t6)                 \n" // t6 = |cur_3|cur_2|
+        "dpa.w.ph        $ac0, $t5, $t7              \n" // (cur*r1)+(cur*r0)
+        "dpa.w.ph        $ac1, $t1, $t7              \n" // (cur*g1)+(cur*g0)
+        "dpa.w.ph        $ac2, $t2, $t7              \n" // (cur*b1)+(cur*b0)
+        "dpa.w.ph        $ac3, $t0, $t7              \n" // (cur*a1)+(cur*a0)
+        "precrq.qb.ph    $t5, $t4, $t3               \n" // t5 = |a3|g3|a2|g2|
+        "precr.qb.ph     $t7, $t4, $t3               \n" // t7 = |b3|r3|b2|r2|
+        "preceu.ph.qbla  $t0, $t5                    \n" // t0 = |0|a3|0|a2|
+        "preceu.ph.qbra  $t1, $t5                    \n" // t1 = |0|g3|0|g2|
+        "preceu.ph.qbla  $t2, $t7                    \n" // t2 = |0|b3|0|b2|
+        "preceu.ph.qbra  $t5, $t7                    \n" // t5 = |0|r3|0|r2|
+        "dpa.w.ph        $ac0, $t5, $t6              \n" // (cur*r3)+(cur*r2)
+        "dpa.w.ph        $ac1, $t1, $t6              \n" // (cur*g3)+(cur*g2)
+        "dpa.w.ph        $ac2, $t2, $t6              \n" // (cur*b3)+(cur*b2)
+        "dpa.w.ph        $ac3, $t0, $t6              \n" // (cur*a3)+(cur*a2)
+        "addiu           %[cnt], %[cnt], -1          \n"
+        "bgtz            %[cnt], 11b                 \n"
+        " addiu          %[fy], %[fy], 8             \n"
+
+        "2:                                          \n"
+        "andi            %[cnt], %[filter_len], 0x3  \n" // residual
+        "beqz            %[cnt], 3f                  \n"
+        " nop                                        \n"
+
+        "21:                                         \n"
+        "addu            $t0, %[filter_val], %[fy]   \n"
+        "lh              $t4, 0($t0)                 \n" // t4=filter_val[fx]
+        "sll             $t1, %[fy], 1               \n"
+        "addu            $t0, %[src_data_rows], $t1  \n"
+        "lw              $t1, 0($t0)                 \n"
+        "addu            $t0, $t1, %[offset]         \n"
+        "lbu             $t1, 0($t0)                 \n" // t1 = row[fx*4 + 0]
+        "lbu             $t2, 1($t0)                 \n" // t2 = row[fx*4 + 1]
+        "lbu             $t3, 2($t0)                 \n" // t3 = row[fx*4 + 2]
+        "lbu             $t0, 3($t0)                 \n" // t4 = row[fx*4 + 2]
+        "maddu           $ac0, $t4, $t1              \n"
+        "maddu           $ac1, $t4, $t2              \n"
+        "maddu           $ac2, $t4, $t3              \n"
+        "maddu           $ac3, $t4, $t0              \n"
+        "addiu           %[cnt], %[cnt], -1          \n"
+        "bgtz            %[cnt], 21b                 \n"
+        " addiu          %[fy], %[fy], 2             \n"
+
+        "3:                                          \n"
+        "extrv.w         $t3, $ac0, %[kShiftBits]    \n" // a >> kShiftBits
+        "extrv.w         $t2, $ac1, %[kShiftBits]    \n" // b >> kShiftBits
+        "extrv.w         $t1, $ac2, %[kShiftBits]    \n" // g >> kShiftBits
+        "extrv.w         $t0, $ac3, %[kShiftBits]    \n" // r >> kShiftBits
+        "repl.ph         $t4, 128                    \n" // t4 = | 128 | 128 |
+        "addu            $t5, %[out_row], %[offset]  \n"
+        "append          $t2, $t3, 16                \n" // t2 = |0|g|0|r|
+        "append          $t0, $t1, 16                \n" // t0 = |0|a|0|b|
+        "subu.ph         $t1, $t0, $t4               \n"
+        "shll_s.ph       $t1, $t1, 8                 \n"
+        "shra.ph         $t1, $t1, 8                 \n"
+        "addu.ph         $t1, $t1, $t4               \n" // Clamp(a)|Clamp(b)
+        "subu.ph         $t2, $t2, $t4               \n"
+        "shll_s.ph       $t2, $t2, 8                 \n"
+        "shra.ph         $t2, $t2, 8                 \n"
+        "addu.ph         $t2, $t2, $t4               \n" // Clamp(g)|Clamp(r)
+        "andi            $t3, $t1, 0xFF              \n" // t3 = ClampTo8(b)
+        "cmp.lt.ph       $t3, $t2                    \n" // cmp b, g, r
+        "pick.ph         $t0, $t2, $t3               \n"
+        "andi            $t3, $t0, 0xFF              \n"
+        "srl             $t4, $t0, 16                \n"
+        "cmp.lt.ph       $t3, $t4                    \n"
+        "pick.ph         $t0, $t4, $t3               \n" // t0 = max_color_ch
+        "srl             $t3, $t1, 16                \n" // t1 = ClampTo8(a)
+        "cmp.lt.ph       $t3, $t0                    \n"
+        "pick.ph         $t0, $t0, $t3               \n"
+        "ins             $t1, $t0, 16, 8             \n"
+        "precr.qb.ph     $t0, $t1, $t2               \n" // t0 = |a|b|g|r|
+        "usw             $t0, 0($t5)                 \n"
+
+        ".set pop                                    \n"
+      : [filter_val] "+r" (filter_val), [filter_len] "+r" (filter_length),
+        [offset] "+r" (byte_offset), [fy] "+r" (filter_y), [cnt] "+r" (cnt),
+        [out_x] "+r" (out_x), [pixel_width] "+r" (pixel_width)
+      : [src_data_rows] "r" (source_data_rows), [out_row] "r" (out_row),
+        [kShiftBits] "r" (ConvolutionFilter1D::kShiftBits)
+      : "lo", "hi", "$ac1lo", "$ac1hi", "$ac2lo", "$ac2hi", "$ac3lo", "$ac3hi",
+        "t0", "t1", "t2", "t3", "t4", "t5", "t6","t7", "memory"
+      );
+    }
+  } else {
+    for (int out_x = 0; out_x < pixel_width; out_x++) {
+      __asm__ __volatile__ (
+        ".set push                                   \n"
+        ".set noreorder                              \n"
+
+        "beqz            %[filter_len], 3f           \n"
+        " sll            %[offset], %[out_x], 2      \n"
+        "mtlo            $0, $ac0                    \n"
+        "mtlo            $0, $ac1                    \n"
+        "mtlo            $0, $ac2                    \n"
+        "srl             %[cnt], %[filter_len], 2    \n"
+        "beqz            %[cnt], 2f                  \n"
+        " li             %[fy], 0                    \n"
+
+        "11:                                         \n"
+        "sll             $t1, %[fy], 1               \n"
+        "addu            $t0, %[src_data_rows], $t1  \n"
+        "lw              $t1, 0($t0)                 \n"
+        "lw              $t2, 4($t0)                 \n"
+        "lw              $t3, 8($t0)                 \n"
+        "lw              $t4, 12($t0)                \n"
+        "addu            $t1, $t1, %[offset]         \n"
+        "addu            $t2, $t2, %[offset]         \n"
+        "addu            $t3, $t3, %[offset]         \n"
+        "addu            $t4, $t4, %[offset]         \n"
+        "lw              $t1, 0($t1)                 \n" // t1 = |a0|b0|g0|r0|
+        "lw              $t2, 0($t2)                 \n" // t2 = |a1|b1|g1|r1|
+        "lw              $t3, 0($t3)                 \n" // t3 = |a0|b0|g0|r0|
+        "lw              $t4, 0($t4)                 \n" // t4 = |a1|b1|g1|r1|
+        "precrq.qb.ph    $t5, $t2, $t1               \n" // t5 = |a1|g1|a0|g0|
+        "precr.qb.ph     $t6, $t2, $t1               \n" // t6 = |b1|r1|b0|r0|
+        "preceu.ph.qbra  $t1, $t5                    \n" // t1 = |0|g1|0|g0|
+        "preceu.ph.qbla  $t2, $t6                    \n" // t2 = |0|b1|0|b0|
+        "preceu.ph.qbra  $t5, $t6                    \n" // t5 = |0|r1|0|r0|
+        "addu            $t6, %[filter_val], %[fy]   \n"
+        "ulw             $t0, 0($t6)                 \n" // t0 = |cur_1|cur_0|
+        "ulw             $t6, 4($t6)                 \n" // t6 = |cur_1|cur_0|
+        "dpa.w.ph        $ac0, $t5, $t0              \n" // (cur*r1)+(cur*r0)
+        "dpa.w.ph        $ac1, $t1, $t0              \n" // (cur*g1)+(cur*g0)
+        "dpa.w.ph        $ac2, $t2, $t0              \n" // (cur*b1)+(cur*b0)
+        "precrq.qb.ph    $t5, $t4, $t3               \n" // t5 = |a3|g3|a2|g2|
+        "precr.qb.ph     $t0, $t4, $t3               \n" // t0 = |b3|r3|b2|r2|
+        "preceu.ph.qbra  $t1, $t5                    \n" // t1 = |0|g3|0|g2|
+        "preceu.ph.qbla  $t2, $t0                    \n" // t2 = |0|b3|0|b2|
+        "preceu.ph.qbra  $t5, $t0                    \n" // t5 = |0|r3|0|r2|
+        "dpa.w.ph        $ac0, $t5, $t6              \n" // (cur*r1)+(cur*r0)
+        "dpa.w.ph        $ac1, $t1, $t6              \n" // (cur*g1)+(cur*g0)
+        "dpa.w.ph        $ac2, $t2, $t6              \n" // (cur*b1)+(cur*b0)
+        "addiu           %[cnt], %[cnt], -1          \n"
+        "bgtz            %[cnt], 11b                 \n"
+        " addiu          %[fy], %[fy], 8             \n"
+
+        "2:                                          \n"
+        "andi            %[cnt], %[filter_len], 0x3  \n" // residual
+        "beqz            %[cnt], 3f                  \n"
+        " nop                                        \n"
+
+        "21:                                         \n"
+        "addu            $t0, %[filter_val], %[fy]   \n"
+        "lh              $t4, 0($t0)                 \n" // filter_val[fx]
+        "sll             $t1, %[fy], 1               \n"
+        "addu            $t0, %[src_data_rows], $t1  \n"
+        "lw              $t1, 0($t0)                 \n"
+        "addu            $t0, $t1, %[offset]         \n"
+        "lbu             $t1, 0($t0)                 \n" // t1 = row[fx*4 + 0]
+        "lbu             $t2, 1($t0)                 \n" // t2 = row[fx*4 + 1]
+        "lbu             $t3, 2($t0)                 \n" // t3 = row[fx*4 + 2]
+        "maddu           $ac0, $t4, $t1              \n"
+        "maddu           $ac1, $t4, $t2              \n"
+        "maddu           $ac2, $t4, $t3              \n"
+        "addiu           %[cnt], %[cnt], -1          \n"
+        "bgtz            %[cnt], 21b                 \n"
+        " addiu          %[fy], %[fy], 2             \n"
+
+        "3:                                          \n"
+        "extrv.w         $t3, $ac0, %[kShiftBits]    \n" // r >> kShiftBits
+        "extrv.w         $t2, $ac1, %[kShiftBits]    \n" // g >> kShiftBits
+        "extrv.w         $t1, $ac2, %[kShiftBits]    \n" // b >> kShiftBits
+        "repl.ph         $t6, 128                    \n" // t6 = | 128 | 128 |
+        "addu            $t5, %[out_row], %[offset]  \n"
+        "append          $t2, $t3, 16                \n" // t2 = |0|g|0|r|
+        "andi            $t1, $t1, 0xFFFF            \n"
+        "subu.ph         $t1, $t1, $t6               \n"
+        "shll_s.ph       $t1, $t1, 8                 \n"
+        "shra.ph         $t1, $t1, 8                 \n"
+        "addu.ph         $t1, $t1, $t6               \n" // Clamp(a)|Clamp(b)
+        "subu.ph         $t2, $t2, $t6               \n"
+        "shll_s.ph       $t2, $t2, 8                 \n"
+        "shra.ph         $t2, $t2, 8                 \n"
+        "addu.ph         $t2, $t2, $t6               \n" // Clamp(g)|Clamp(r)
+        "li              $t0, 0xFF                   \n"
+        "ins             $t1, $t0, 16, 8             \n"
+        "precr.qb.ph     $t0, $t1, $t2               \n" // t0 = |a|b|g|r|
+        "usw             $t0, 0($t5)                 \n"
+
+        ".set pop                                    \n"
+      : [filter_val] "+r" (filter_val), [filter_len] "+r" (filter_length),
+        [offset] "+r" (byte_offset), [fy] "+r" (filter_y), [cnt] "+r" (cnt),
+        [out_x] "+r" (out_x), [pixel_width] "+r" (pixel_width)
+      : [src_data_rows] "r" (source_data_rows), [out_row] "r" (out_row),
+        [kShiftBits] "r" (ConvolutionFilter1D::kShiftBits)
+      : "lo", "hi", "$ac1lo", "$ac1hi", "$ac2lo", "$ac2hi", "$ac3lo", "$ac3hi",
+        "t0", "t1", "t2", "t3", "t4", "t5", "t6", "memory"
+      );
+    }
+  }
+#endif
+}
+} // namespace skia
diff --git a/skia/ext/convolver_mips_dspr2.h b/skia/ext/convolver_mips_dspr2.h
new file mode 100644
index 0000000..a5e59f9
--- /dev/null
+++ b/skia/ext/convolver_mips_dspr2.h
@@ -0,0 +1,25 @@
+// Copyright (c) 2013 The Chromium Authors. All rights reserved.
+// Use of this source code is governed by a BSD-style license that can be
+// found in the LICENSE file.
+
+#ifndef SKIA_EXT_CONVOLVER_MIPS_DSPR2_H_
+#define SKIA_EXT_CONVOLVER_MIPS_DSPR2_H_
+
+#include "skia/ext/convolver.h"
+
+namespace skia {
+
+void ConvolveVertically_mips_dspr2(const ConvolutionFilter1D::Fixed* filter_val,
+                                   int filter_length,
+                                   unsigned char* const* source_data_rows,
+                                   int pixel_width,
+                                   unsigned char* out_row,
+                                   bool has_alpha);
+
+void ConvolveHorizontally_mips_dspr2(const unsigned char* src_data,
+                                     const ConvolutionFilter1D& filter,
+                                     unsigned char* out_row,
+                                     bool has_alpha);
+}  // namespace skia
+
+#endif  // SKIA_EXT_CONVOLVER_MIPS_DSPR2_H_
diff --git a/skia/skia.gyp b/skia/skia.gyp
index 79773ec4..4b79090 100644
--- a/skia/skia.gyp
+++ b/skia/skia.gyp
@@ -747,6 +747,7 @@
             '../third_party/skia/src/opts/SkBitmapProcState_opts_none.cpp',
             '../third_party/skia/src/opts/SkBlitRow_opts_none.cpp',
             '../third_party/skia/src/opts/SkUtils_opts_none.cpp',
+            'ext/convolver_mips_dspr2.cc',
           ],
         }],
       ],
@@ -790,7 +791,7 @@
             ],
           },
         }],
-        [ 'target_arch != "arm"', {
+        [ 'target_arch != "arm" and target_arch != "mipsel"', {
           'sources': [
             '../third_party/skia/src/opts/SkBitmapProcState_opts_SSSE3.cpp',
           ],