move more geometry to simd

src/utils/SkPatchUtils.cpp

 /*
  * Copyright 2014 Google Inc.
  *
  * Use of this source code is governed by a BSD-style license that can be
  * found in the LICENSE file.
  */
 
 #include "SkPatchUtils.h"
 
 #include "SkColorPriv.h"
@@ skipping 12 matching lines @@
  * For the cubic case the first difference gives as a result a quadratic polynomial to which we can
  * apply again forward differences and get linear function to which we can apply again forward
  * differences to get a constant difference. This is why we keep an array of size 4, the 0th
  * position keeps the sampled value while the next ones keep the quadratic, linear and constant
  * difference values.
  */
 
 class FwDCubicEvaluator {
     
 public:
-    FwDCubicEvaluator()
-    : fMax(0)
-    , fCurrent(0)
-    , fDivisions(0) {
-        memset(fPoints, 0, 4 * sizeof(SkPoint));
-        memset(fPoints, 0, 4 * sizeof(SkPoint));
-        memset(fPoints, 0, 4 * sizeof(SkPoint));
-    }
     
     /**
      * Receives the 4 control points of the cubic bezier.
      */
-    FwDCubicEvaluator(SkPoint a, SkPoint b, SkPoint c, SkPoint d) {
-        fPoints[0] = a;
-        fPoints[1] = b;
-        fPoints[2] = c;
-        fPoints[3] = d;
-        
-        SkCubicToCoeff(fPoints, fCoefs);
+    
+    explicit FwDCubicEvaluator(const SkPoint points[4])
+            : fCoefs(points) {
+        memcpy(fPoints, points, 4 * sizeof(SkPoint));
         
         this->restart(1);
     }
     
-    explicit FwDCubicEvaluator(const SkPoint points[4])  {
-        memcpy(fPoints, points, 4 * sizeof(SkPoint));
-        
-        SkCubicToCoeff(fPoints, fCoefs);
-        
-        this->restart(1);
-    }
-    
     /**
      * Restarts the forward differences evaluator to the first value of t = 0.
      */
     void restart(int divisions)  {
         fDivisions = divisions;
-        SkScalar h  = 1.f / fDivisions;
         fCurrent    = 0;
         fMax        = fDivisions + 1;
-        fFwDiff[0]  = fCoefs[3];
-        SkScalar h2 = h * h;
-        SkScalar h3 = h2 * h;
-        
-        fFwDiff[3].set(6.f * fCoefs[0].x() * h3, 6.f * fCoefs[0].y() * h3); //6ah^3
-        fFwDiff[2].set(fFwDiff[3].x() + 2.f * fCoefs[1].x() * h2, //6ah^3 + 2bh^2
-                       fFwDiff[3].y() + 2.f * fCoefs[1].y() * h2);
-        fFwDiff[1].set(fCoefs[0].x() * h3 + fCoefs[1].x() * h2 + fCoefs[2].x() * h,//ah^3 + bh^2 +ch
-                       fCoefs[0].y() * h3 + fCoefs[1].y() * h2 + fCoefs[2].y() * h);
+        Sk2s h  = Sk2s(1.f / fDivisions);
+        Sk2s h2 = h * h;
+        Sk2s h3 = h2 * h;
+        Sk2s fwDiff3 = Sk2s(6) * fCoefs.fA * h3;
+        fFwDiff[3] = to_point(fwDiff3);
+        fFwDiff[2] = to_point(fwDiff3 + times_2(fCoefs.fB) * h2);
+        fFwDiff[1] = to_point(fCoefs.fA * h3 + fCoefs.fB * h2 + fCoefs.fC * h);
+        fFwDiff[0] = to_point(fCoefs.fD);
     }
     
     /**
      * Check if the evaluator is still within the range of 0<=t<=1
      */
     bool done() const {
         return fCurrent > fMax;
     }
     
     /**
      * Call next to obtain the SkPoint sampled and move to the next one.
      */
     SkPoint next() {
         SkPoint point = fFwDiff[0];
         fFwDiff[0]    += fFwDiff[1];
         fFwDiff[1]    += fFwDiff[2];
         fFwDiff[2]    += fFwDiff[3];
         fCurrent++;
         return point;
     }
     
     const SkPoint* getCtrlPoints() const {
         return fPoints;
     }
     
 private:
+    SkCubicCoeff fCoefs;
     int fMax, fCurrent, fDivisions;
-    SkPoint fFwDiff[4], fCoefs[4], fPoints[4];
+    SkPoint fFwDiff[4], fPoints[4];
 };
 
 ////////////////////////////////////////////////////////////////////////////////
 
 // size in pixels of each partition per axis, adjust this knob
 static const int kPartitionSize = 10;
 
 /**
  * Calculate the approximate arc length given a bezier curve's control points.
  */
@@ skipping 203 matching lines @@
                 data->fIndices[i + 4] = data->fIndices[i + 2];
                 data->fIndices[i + 5] = (x + 1) * stride + y;
             }
             v = SkScalarClampMax(v + 1.f / lodY, 1);
         }
         u = SkScalarClampMax(u + 1.f / lodX, 1);
     }
     return true;
 
 }