| OLD | NEW |
| 1 /* | 1 /* |
| 2 * Copyright 2012 Google Inc. | 2 * Copyright 2012 Google Inc. |
| 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 #include "SkFloatBits.h" | 7 #include "SkFloatBits.h" |
| 8 #include "SkOpCoincidence.h" | 8 #include "SkOpCoincidence.h" |
| 9 #include "SkPathOpsTypes.h" | 9 #include "SkPathOpsTypes.h" |
| 10 | 10 |
| 11 static bool arguments_denormalized(float a, float b, int epsilon) { | 11 static bool arguments_denormalized(float a, float b, int epsilon) { |
| 12 float denormalizedCheck = FLT_EPSILON * epsilon / 2; | 12 float denormalizedCheck = FLT_EPSILON * epsilon / 2; |
| 13 return fabsf(a) <= denormalizedCheck && fabsf(b) <= denormalizedCheck; | 13 return fabsf(a) <= denormalizedCheck && fabsf(b) <= denormalizedCheck; |
| 14 } | 14 } |
| 15 | 15 |
| 16 // from http://randomascii.wordpress.com/2012/02/25/comparing-floating-point-num
bers-2012-edition/ | 16 // from http://randomascii.wordpress.com/2012/02/25/comparing-floating-point-num
bers-2012-edition/ |
| 17 // FIXME: move to SkFloatBits.h | 17 // FIXME: move to SkFloatBits.h |
| 18 static bool equal_ulps(float a, float b, int epsilon, int depsilon) { | 18 static bool equal_ulps(float a, float b, int epsilon, int depsilon) { |
| 19 if (arguments_denormalized(a, b, depsilon)) { | 19 if (arguments_denormalized(a, b, depsilon)) { |
| 20 return true; | 20 return true; |
| 21 } | 21 } |
| 22 int aBits = SkFloatAs2sCompliment(a); | 22 int aBits = SkFloatAs2sCompliment(a); |
| 23 int bBits = SkFloatAs2sCompliment(b); | 23 int bBits = SkFloatAs2sCompliment(b); |
| 24 // Find the difference in ULPs. | 24 // Find the difference in ULPs. |
| 25 return aBits < bBits + epsilon && bBits < aBits + epsilon; | 25 return aBits < bBits + epsilon && bBits < aBits + epsilon; |
| 26 } | 26 } |
| 27 | 27 |
| 28 static bool equal_ulps_no_normal_check(float a, float b, int epsilon, int depsil
on) { |
| 29 int aBits = SkFloatAs2sCompliment(a); |
| 30 int bBits = SkFloatAs2sCompliment(b); |
| 31 // Find the difference in ULPs. |
| 32 return aBits < bBits + epsilon && bBits < aBits + epsilon; |
| 33 } |
| 34 |
| 28 static bool equal_ulps_pin(float a, float b, int epsilon, int depsilon) { | 35 static bool equal_ulps_pin(float a, float b, int epsilon, int depsilon) { |
| 29 if (!SkScalarIsFinite(a) || !SkScalarIsFinite(b)) { | 36 if (!SkScalarIsFinite(a) || !SkScalarIsFinite(b)) { |
| 30 return false; | 37 return false; |
| 31 } | 38 } |
| 32 if (arguments_denormalized(a, b, depsilon)) { | 39 if (arguments_denormalized(a, b, depsilon)) { |
| 33 return true; | 40 return true; |
| 34 } | 41 } |
| 35 int aBits = SkFloatAs2sCompliment(a); | 42 int aBits = SkFloatAs2sCompliment(a); |
| 36 int bBits = SkFloatAs2sCompliment(b); | 43 int bBits = SkFloatAs2sCompliment(b); |
| 37 // Find the difference in ULPs. | 44 // Find the difference in ULPs. |
| (...skipping 75 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 113 | 120 |
| 114 bool AlmostDequalUlps(double a, double b) { | 121 bool AlmostDequalUlps(double a, double b) { |
| 115 return AlmostDequalUlps(SkDoubleToScalar(a), SkDoubleToScalar(b)); | 122 return AlmostDequalUlps(SkDoubleToScalar(a), SkDoubleToScalar(b)); |
| 116 } | 123 } |
| 117 | 124 |
| 118 bool AlmostEqualUlps(float a, float b) { | 125 bool AlmostEqualUlps(float a, float b) { |
| 119 const int UlpsEpsilon = 16; | 126 const int UlpsEpsilon = 16; |
| 120 return equal_ulps(a, b, UlpsEpsilon, UlpsEpsilon); | 127 return equal_ulps(a, b, UlpsEpsilon, UlpsEpsilon); |
| 121 } | 128 } |
| 122 | 129 |
| 130 bool AlmostEqualUlpsNoNormalCheck(float a, float b) { |
| 131 const int UlpsEpsilon = 16; |
| 132 return equal_ulps_no_normal_check(a, b, UlpsEpsilon, UlpsEpsilon); |
| 133 } |
| 134 |
| 123 bool AlmostEqualUlps_Pin(float a, float b) { | 135 bool AlmostEqualUlps_Pin(float a, float b) { |
| 124 const int UlpsEpsilon = 16; | 136 const int UlpsEpsilon = 16; |
| 125 return equal_ulps_pin(a, b, UlpsEpsilon, UlpsEpsilon); | 137 return equal_ulps_pin(a, b, UlpsEpsilon, UlpsEpsilon); |
| 126 } | 138 } |
| 127 | 139 |
| 128 bool NotAlmostEqualUlps(float a, float b) { | 140 bool NotAlmostEqualUlps(float a, float b) { |
| 129 const int UlpsEpsilon = 16; | 141 const int UlpsEpsilon = 16; |
| 130 return not_equal_ulps(a, b, UlpsEpsilon); | 142 return not_equal_ulps(a, b, UlpsEpsilon); |
| 131 } | 143 } |
| 132 | 144 |
| (...skipping 72 matching lines...) Expand 10 before | Expand all | Expand 10 after Loading... |
| 205 if (approximately_zero_cubed(x)) { | 217 if (approximately_zero_cubed(x)) { |
| 206 return 0; | 218 return 0; |
| 207 } | 219 } |
| 208 double result = halley_cbrt3d(fabs(x)); | 220 double result = halley_cbrt3d(fabs(x)); |
| 209 if (x < 0) { | 221 if (x < 0) { |
| 210 result = -result; | 222 result = -result; |
| 211 } | 223 } |
| 212 return result; | 224 return result; |
| 213 } | 225 } |
| 214 | 226 |
| 215 SkOpGlobalState::SkOpGlobalState(SkOpCoincidence* coincidence, SkOpContourHead*
head | 227 SkOpGlobalState::SkOpGlobalState(SkOpContourHead* head, |
| 228 SkChunkAlloc* allocator |
| 216 SkDEBUGPARAMS(bool debugSkipAssert) | 229 SkDEBUGPARAMS(bool debugSkipAssert) |
| 217 SkDEBUGPARAMS(const char* testName)) | 230 SkDEBUGPARAMS(const char* testName)) |
| 218 : fCoincidence(coincidence) | 231 : fAllocator(allocator) |
| 232 , fCoincidence(nullptr) |
| 219 , fContourHead(head) | 233 , fContourHead(head) |
| 220 , fNested(0) | 234 , fNested(0) |
| 221 , fWindingFailed(false) | 235 , fWindingFailed(false) |
| 222 , fAngleCoincidence(false) | 236 , fAngleCoincidence(false) |
| 223 , fPhase(kIntersecting) | 237 , fPhase(kIntersecting) |
| 224 SkDEBUGPARAMS(fDebugTestName(testName)) | 238 SkDEBUGPARAMS(fDebugTestName(testName)) |
| 225 SkDEBUGPARAMS(fAngleID(0)) | 239 SkDEBUGPARAMS(fAngleID(0)) |
| 226 SkDEBUGPARAMS(fCoinID(0)) | 240 SkDEBUGPARAMS(fCoinID(0)) |
| 227 SkDEBUGPARAMS(fContourID(0)) | 241 SkDEBUGPARAMS(fContourID(0)) |
| 228 SkDEBUGPARAMS(fPtTID(0)) | 242 SkDEBUGPARAMS(fPtTID(0)) |
| 229 SkDEBUGPARAMS(fSegmentID(0)) | 243 SkDEBUGPARAMS(fSegmentID(0)) |
| 230 SkDEBUGPARAMS(fSpanID(0)) | 244 SkDEBUGPARAMS(fSpanID(0)) |
| 231 SkDEBUGPARAMS(fDebugSkipAssert(debugSkipAssert)) { | 245 SkDEBUGPARAMS(fDebugSkipAssert(debugSkipAssert)) { |
| 232 if (coincidence) { | |
| 233 coincidence->debugSetGlobalState(this); | |
| 234 } | |
| 235 #if DEBUG_T_SECT_LOOP_COUNT | 246 #if DEBUG_T_SECT_LOOP_COUNT |
| 236 debugResetLoopCounts(); | 247 debugResetLoopCounts(); |
| 237 #endif | 248 #endif |
| 238 } | 249 } |
| 239 | |
| OLD | NEW |