OLD | NEW |
(Empty) | |
| 1 /* |
| 2 * Copyright 2015 Google Inc. |
| 3 * |
| 4 * Use of this source code is governed by a BSD-style license that can be |
| 5 * found in the LICENSE file. |
| 6 */ |
| 7 |
| 8 // It is important _not_ to put header guards here. |
| 9 // This file will be intentionally included three times. |
| 10 |
| 11 #include "SkTypes.h" // Keep this before any #ifdef for skbug.com/3362 |
| 12 |
| 13 #if defined(SK2X_PREAMBLE) |
| 14 #include <immintrin.h> |
| 15 template <typename T> struct SkScalarToSIMD; |
| 16 template <> struct SkScalarToSIMD< float> { typedef __m128 Type; }; |
| 17 template <> struct SkScalarToSIMD<double> { typedef __m128d Type; }; |
| 18 |
| 19 |
| 20 #elif defined(SK2X_PRIVATE) |
| 21 typename SkScalarToSIMD<T>::Type fVec; |
| 22 /*implicit*/ Sk2x(const typename SkScalarToSIMD<T>::Type vec) { fVec = vec;
} |
| 23 |
| 24 #else |
| 25 |
| 26 #define M(...) template <> inline __VA_ARGS__ Sk2x<float>:: |
| 27 |
| 28 M() Sk2x() {} |
| 29 M() Sk2x(float val) { fVec = _mm_set1_ps(val); } |
| 30 M() Sk2x(float a, float b) { fVec = _mm_set_ps(b,a,b,a); } |
| 31 M(Sk2f&) operator=(const Sk2f& o) { fVec = o.fVec; return *this; } |
| 32 |
| 33 M(Sk2f) Load(const float vals[2]) { |
| 34 return _mm_castsi128_ps(_mm_loadl_epi64((const __m128i*)vals)); |
| 35 } |
| 36 M(void) store(float vals[2]) const { _mm_storel_pi((__m64*)vals, fVec); } |
| 37 |
| 38 M(Sk2f) add(const Sk2f& o) const { return _mm_add_ps(fVec, o.fVec); } |
| 39 M(Sk2f) subtract(const Sk2f& o) const { return _mm_sub_ps(fVec, o.fVec); } |
| 40 M(Sk2f) multiply(const Sk2f& o) const { return _mm_mul_ps(fVec, o.fVec); } |
| 41 |
| 42 M(Sk2f) Min(const Sk2f& a, const Sk2f& b) { return _mm_min_ps(a.fVec, b.fVec); } |
| 43 M(Sk2f) Max(const Sk2f& a, const Sk2f& b) { return _mm_max_ps(a.fVec, b.fVec); } |
| 44 |
| 45 M(Sk2f) rsqrt() const { return _mm_rsqrt_ps(fVec); } |
| 46 M(Sk2f) sqrt() const { return _mm_sqrt_ps (fVec); } |
| 47 |
| 48 #undef M |
| 49 |
| 50 #define M(...) template <> inline __VA_ARGS__ Sk2x<double>:: |
| 51 |
| 52 M() Sk2x() {} |
| 53 M() Sk2x(double val) { fVec = _mm_set1_pd(val); } |
| 54 M() Sk2x(double a, double b) { fVec = _mm_set_pd(b, a); } |
| 55 M(Sk2d&) operator=(const Sk2d& o) { fVec = o.fVec; return *this; } |
| 56 |
| 57 M(Sk2d) Load(const double vals[2]) { return _mm_loadu_pd(vals); } |
| 58 M(void) store(double vals[2]) const { _mm_storeu_pd(vals, fVec); } |
| 59 |
| 60 M(Sk2d) add(const Sk2d& o) const { return _mm_add_pd(fVec, o.fVec); } |
| 61 M(Sk2d) subtract(const Sk2d& o) const { return _mm_sub_pd(fVec, o.fVec); } |
| 62 M(Sk2d) multiply(const Sk2d& o) const { return _mm_mul_pd(fVec, o.fVec); } |
| 63 |
| 64 M(Sk2d) Min(const Sk2d& a, const Sk2d& b) { return _mm_min_pd(a.fVec, b.fVec); } |
| 65 M(Sk2d) Max(const Sk2d& a, const Sk2d& b) { return _mm_max_pd(a.fVec, b.fVec); } |
| 66 |
| 67 // There is no _mm_rsqrt_pd, so we do Sk2d::rsqrt() in floats. |
| 68 M(Sk2d) rsqrt() const { return _mm_cvtps_pd(_mm_rsqrt_ps(_mm_cvtpd_ps(fVec))); } |
| 69 M(Sk2d) sqrt() const { return _mm_sqrt_pd(fVec); } |
| 70 |
| 71 #undef M |
| 72 |
| 73 #endif |
OLD | NEW |