| OLD | NEW |
| 1 /* | 1 /* |
| 2 * Copyright 2008 The Android Open Source Project | 2 * Copyright 2008 The Android Open Source Project |
| 3 * | 3 * |
| 4 * Use of this source code is governed by a BSD-style license that can be | 4 * Use of this source code is governed by a BSD-style license that can be |
| 5 * found in the LICENSE file. | 5 * found in the LICENSE file. |
| 6 */ | 6 */ |
| 7 | 7 |
| 8 #include "SkMathPriv.h" | 8 #include "SkMathPriv.h" |
| 9 #include "SkCordic.h" | 9 #include "SkCordic.h" |
| 10 #include "SkFloatBits.h" | 10 #include "SkFloatBits.h" |
| (...skipping 87 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 98 #ifdef SK_DEBUGx | 98 #ifdef SK_DEBUGx |
| 99 int32_t tmp = hi >> shift; | 99 int32_t tmp = hi >> shift; |
| 100 SkASSERT(tmp == 0 || tmp == -1); | 100 SkASSERT(tmp == 0 || tmp == -1); |
| 101 #endif | 101 #endif |
| 102 // we want (hi << (32 - shift)) | (lo >> shift) but rounded | 102 // we want (hi << (32 - shift)) | (lo >> shift) but rounded |
| 103 int roundBit = (lo >> (shift - 1)) & 1; | 103 int roundBit = (lo >> (shift - 1)) & 1; |
| 104 return ((hi << (32 - shift)) | (lo >> shift)) + roundBit; | 104 return ((hi << (32 - shift)) | (lo >> shift)) + roundBit; |
| 105 } | 105 } |
| 106 } | 106 } |
| 107 | 107 |
| 108 SkFixed SkFixedMul_portable(SkFixed a, SkFixed b) { | |
| 109 #if 0 | |
| 110 Sk64 tmp; | |
| 111 | |
| 112 tmp.setMul(a, b); | |
| 113 tmp.shiftRight(16); | |
| 114 return tmp.fLo; | |
| 115 #elif defined(SkLONGLONG) | |
| 116 return static_cast<SkFixed>((SkLONGLONG)a * b >> 16); | |
| 117 #else | |
| 118 int sa = SkExtractSign(a); | |
| 119 int sb = SkExtractSign(b); | |
| 120 // now make them positive | |
| 121 a = SkApplySign(a, sa); | |
| 122 b = SkApplySign(b, sb); | |
| 123 | |
| 124 uint32_t ah = a >> 16; | |
| 125 uint32_t al = a & 0xFFFF; | |
| 126 uint32_t bh = b >> 16; | |
| 127 uint32_t bl = b & 0xFFFF; | |
| 128 | |
| 129 uint32_t R = ah * b + al * bh + (al * bl >> 16); | |
| 130 | |
| 131 return SkApplySign(R, sa ^ sb); | |
| 132 #endif | |
| 133 } | |
| 134 | |
| 135 SkFract SkFractMul_portable(SkFract a, SkFract b) { | |
| 136 #if 0 | |
| 137 Sk64 tmp; | |
| 138 tmp.setMul(a, b); | |
| 139 return tmp.getFract(); | |
| 140 #elif defined(SkLONGLONG) | |
| 141 return static_cast<SkFract>((SkLONGLONG)a * b >> 30); | |
| 142 #else | |
| 143 int sa = SkExtractSign(a); | |
| 144 int sb = SkExtractSign(b); | |
| 145 // now make them positive | |
| 146 a = SkApplySign(a, sa); | |
| 147 b = SkApplySign(b, sb); | |
| 148 | |
| 149 uint32_t ah = a >> 16; | |
| 150 uint32_t al = a & 0xFFFF; | |
| 151 uint32_t bh = b >> 16; | |
| 152 uint32_t bl = b & 0xFFFF; | |
| 153 | |
| 154 uint32_t A = ah * bh; | |
| 155 uint32_t B = ah * bl + al * bh; | |
| 156 uint32_t C = al * bl; | |
| 157 | |
| 158 /* [ A ] | |
| 159 [ B ] | |
| 160 [ C ] | |
| 161 */ | |
| 162 uint32_t Lo = C + (B << 16); | |
| 163 uint32_t Hi = A + (B >>16) + (Lo < C); | |
| 164 | |
| 165 SkASSERT((Hi >> 29) == 0); // else overflow | |
| 166 | |
| 167 int32_t R = (Hi << 2) + (Lo >> 30); | |
| 168 | |
| 169 return SkApplySign(R, sa ^ sb); | |
| 170 #endif | |
| 171 } | |
| 172 | |
| 173 int SkFixedMulCommon(SkFixed a, int b, int bias) { | 108 int SkFixedMulCommon(SkFixed a, int b, int bias) { |
| 174 // this function only works if b is 16bits | 109 // this function only works if b is 16bits |
| 175 SkASSERT(b == (int16_t)b); | 110 SkASSERT(b == (int16_t)b); |
| 176 SkASSERT(b >= 0); | 111 SkASSERT(b >= 0); |
| 177 | 112 |
| 178 int sa = SkExtractSign(a); | 113 int sa = SkExtractSign(a); |
| 179 a = SkApplySign(a, sa); | 114 a = SkApplySign(a, sa); |
| 180 uint32_t ah = a >> 16; | 115 uint32_t ah = a >> 16; |
| 181 uint32_t al = a & 0xFFFF; | 116 uint32_t al = a & 0xFFFF; |
| 182 uint32_t R = ah * b + ((al * b + bias) >> 16); | 117 uint32_t R = ah * b + ((al * b + bias) >> 16); |
| (...skipping 325 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 508 } | 443 } |
| 509 | 444 |
| 510 /////////////////////////////////////////////////////////////////////////////// | 445 /////////////////////////////////////////////////////////////////////////////// |
| 511 | 446 |
| 512 SkFixed SkFixedTan(SkFixed radians) { return SkCordicTan(radians); } | 447 SkFixed SkFixedTan(SkFixed radians) { return SkCordicTan(radians); } |
| 513 SkFixed SkFixedASin(SkFixed x) { return SkCordicASin(x); } | 448 SkFixed SkFixedASin(SkFixed x) { return SkCordicASin(x); } |
| 514 SkFixed SkFixedACos(SkFixed x) { return SkCordicACos(x); } | 449 SkFixed SkFixedACos(SkFixed x) { return SkCordicACos(x); } |
| 515 SkFixed SkFixedATan2(SkFixed y, SkFixed x) { return SkCordicATan2(y, x); } | 450 SkFixed SkFixedATan2(SkFixed y, SkFixed x) { return SkCordicATan2(y, x); } |
| 516 SkFixed SkFixedExp(SkFixed x) { return SkCordicExp(x); } | 451 SkFixed SkFixedExp(SkFixed x) { return SkCordicExp(x); } |
| 517 SkFixed SkFixedLog(SkFixed x) { return SkCordicLog(x); } | 452 SkFixed SkFixedLog(SkFixed x) { return SkCordicLog(x); } |
| OLD | NEW |