caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 1 | /* |
| 2 | * Copyright 2012 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 | #include "PathOpsCubicIntersectionTestData.h" |
| 8 | #include "PathOpsTestCommon.h" |
| 9 | #include "SkIntersections.h" |
| 10 | #include "SkPathOpsRect.h" |
| 11 | #include "SkReduceOrder.h" |
| 12 | #include "Test.h" |
| 13 | |
| 14 | const int firstCubicIntersectionTest = 9; |
| 15 | |
| 16 | static void standardTestCases(skiatest::Reporter* reporter) { |
| 17 | for (size_t index = firstCubicIntersectionTest; index < tests_count; ++index) { |
| 18 | int iIndex = static_cast<int>(index); |
| 19 | const SkDCubic& cubic1 = tests[index][0]; |
| 20 | const SkDCubic& cubic2 = tests[index][1]; |
| 21 | SkReduceOrder reduce1, reduce2; |
| 22 | int order1 = reduce1.reduce(cubic1, SkReduceOrder::kNo_Quadratics, |
| 23 | SkReduceOrder::kFill_Style); |
| 24 | int order2 = reduce2.reduce(cubic2, SkReduceOrder::kNo_Quadratics, |
| 25 | SkReduceOrder::kFill_Style); |
| 26 | const bool showSkipped = false; |
| 27 | if (order1 < 4) { |
| 28 | if (showSkipped) { |
| 29 | SkDebugf("%s [%d] cubic1 order=%d\n", __FUNCTION__, iIndex, order1); |
| 30 | } |
| 31 | continue; |
| 32 | } |
| 33 | if (order2 < 4) { |
| 34 | if (showSkipped) { |
| 35 | SkDebugf("%s [%d] cubic2 order=%d\n", __FUNCTION__, iIndex, order2); |
| 36 | } |
| 37 | continue; |
| 38 | } |
| 39 | SkIntersections tIntersections; |
| 40 | tIntersections.intersect(cubic1, cubic2); |
| 41 | if (!tIntersections.used()) { |
| 42 | if (showSkipped) { |
| 43 | SkDebugf("%s [%d] no intersection\n", __FUNCTION__, iIndex); |
| 44 | } |
| 45 | continue; |
| 46 | } |
| 47 | if (tIntersections.isCoincident(0)) { |
| 48 | if (showSkipped) { |
| 49 | SkDebugf("%s [%d] coincident\n", __FUNCTION__, iIndex); |
| 50 | } |
| 51 | continue; |
| 52 | } |
| 53 | for (int pt = 0; pt < tIntersections.used(); ++pt) { |
| 54 | double tt1 = tIntersections[0][pt]; |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 55 | SkDPoint xy1 = cubic1.ptAtT(tt1); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 56 | double tt2 = tIntersections[1][pt]; |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 57 | SkDPoint xy2 = cubic2.ptAtT(tt2); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 58 | if (!xy1.approximatelyEqual(xy2)) { |
| 59 | SkDebugf("%s [%d,%d] x!= t1=%g (%g,%g) t2=%g (%g,%g)\n", |
| 60 | __FUNCTION__, (int)index, pt, tt1, xy1.fX, xy1.fY, tt2, xy2.fX, xy2.fY); |
| 61 | } |
| 62 | REPORTER_ASSERT(reporter, xy1.approximatelyEqual(xy2)); |
| 63 | } |
| 64 | } |
| 65 | } |
| 66 | |
| 67 | static const SkDCubic testSet[] = { |
| 68 | // FIXME: uncommenting these two will cause this to fail |
| 69 | // this results in two curves very nearly but not exactly coincident |
| 70 | #if 0 |
| 71 | {{{67.426548091427676, 37.993772624988935}, {23.483695892376684, 90.476863174921306}, |
| 72 | {35.597065061143162, 79.872482633158796}, {75.38634169631932, 18.244890038969412}}}, |
| 73 | {{{67.4265481, 37.9937726}, {23.4836959, 90.4768632}, {35.5970651, 79.8724826}, |
| 74 | {75.3863417, 18.24489}}}, |
| 75 | #endif |
| 76 | |
| 77 | {{{0, 0}, {0, 1}, {1, 1}, {1, 0}}}, |
| 78 | {{{1, 0}, {0, 0}, {0, 1}, {1, 1}}}, |
| 79 | |
| 80 | {{{0, 1}, {4, 5}, {1, 0}, {5, 3}}}, |
| 81 | {{{0, 1}, {3, 5}, {1, 0}, {5, 4}}}, |
| 82 | |
| 83 | {{{0, 1}, {1, 6}, {1, 0}, {1, 0}}}, |
| 84 | {{{0, 1}, {0, 1}, {1, 0}, {6, 1}}}, |
| 85 | |
| 86 | {{{0, 1}, {3, 4}, {1, 0}, {5, 1}}}, |
| 87 | {{{0, 1}, {1, 5}, {1, 0}, {4, 3}}}, |
| 88 | |
| 89 | {{{0, 1}, {1, 2}, {1, 0}, {6, 1}}}, |
| 90 | {{{0, 1}, {1, 6}, {1, 0}, {2, 1}}}, |
| 91 | |
| 92 | {{{0, 1}, {0, 5}, {1, 0}, {4, 0}}}, |
| 93 | {{{0, 1}, {0, 4}, {1, 0}, {5, 0}}}, |
| 94 | |
| 95 | {{{0, 1}, {3, 4}, {1, 0}, {3, 0}}}, |
| 96 | {{{0, 1}, {0, 3}, {1, 0}, {4, 3}}}, |
| 97 | |
| 98 | {{{0, 0}, {1, 2}, {3, 4}, {4, 4}}}, |
| 99 | {{{0, 0}, {1, 2}, {3, 4}, {4, 4}}}, |
| 100 | {{{4, 4}, {3, 4}, {1, 2}, {0, 0}}}, |
| 101 | |
| 102 | {{{0, 1}, {2, 3}, {1, 0}, {1, 0}}}, |
| 103 | {{{0, 1}, {0, 1}, {1, 0}, {3, 2}}}, |
| 104 | |
| 105 | {{{0, 2}, {0, 1}, {1, 0}, {1, 0}}}, |
| 106 | {{{0, 1}, {0, 1}, {2, 0}, {1, 0}}}, |
| 107 | |
| 108 | {{{0, 1}, {0, 2}, {1, 0}, {1, 0}}}, |
| 109 | {{{0, 1}, {0, 1}, {1, 0}, {2, 0}}}, |
| 110 | |
| 111 | {{{0, 1}, {1, 6}, {1, 0}, {2, 0}}}, |
| 112 | {{{0, 1}, {0, 2}, {1, 0}, {6, 1}}}, |
| 113 | |
| 114 | {{{0, 1}, {5, 6}, {1, 0}, {1, 0}}}, |
| 115 | {{{0, 1}, {0, 1}, {1, 0}, {6, 5}}}, |
| 116 | |
| 117 | {{{95.837747722788592, 45.025976907939643}, {16.564570095652982, 0.72959763963222402}, |
| 118 | {63.209855865319199, 68.047528419665767}, {57.640240647662544, 59.524565264361243}}}, |
| 119 | {{{51.593891741518817, 38.53849970667553}, {62.34752929878772, 74.924924725166022}, |
| 120 | {74.810149322641152, 34.17966562983564}, {29.368398119401373, 94.66719277886078}}}, |
| 121 | |
| 122 | {{{39.765160968417838, 33.060396198677083}, {5.1922921581157908, 66.854301452103215}, |
| 123 | {31.619281802149157, 25.269248720849514}, {81.541621071073038, 70.025341524754353}}}, |
| 124 | {{{46.078911165743556, 48.259962651999651}, {20.24450549867214, 49.403916182650214}, |
| 125 | {0.26325131778756683, 24.46489805563581}, {15.915006546264051, 83.515023059917155}}}, |
| 126 | |
| 127 | {{{65.454505973241524, 93.881892270353575}, {45.867360264932437, 92.723972719499827}, |
| 128 | {2.1464054482739447, 74.636369140183717}, {33.774068594804994, 40.770872887582925}}}, |
| 129 | {{{72.963387832494163, 95.659300729473728}, {11.809496633619768, 82.209921247423594}, |
| 130 | {13.456139067865974, 57.329313623406605}, {36.060621606214262, 70.867335643091849}}}, |
| 131 | |
| 132 | {{{32.484981432782945, 75.082940782924624}, {42.467313093350882, 48.131159948246157}, |
| 133 | {3.5963115764764657, 43.208665839959245}, {79.442476890721579, 89.709102357602262}}}, |
| 134 | {{{18.98573861410177, 93.308887208490106}, {40.405250173250792, 91.039661826118675}, |
| 135 | {8.0467721950480584, 42.100282172719147}, {40.883324221187891, 26.030185504830527}}}, |
| 136 | |
| 137 | {{{7.5374809128872498, 82.441702896003477}, {22.444346930107265, 22.138854312775123}, |
| 138 | {66.76091829629658, 50.753805856571446}, {78.193478508942519, 97.7932997968948}}}, |
| 139 | {{{97.700573130371311, 53.53260215070685}, {87.72443481149358, 84.575876772671876}, |
| 140 | {19.215031396232092, 47.032676472809484}, {11.989686410869325, 10.659507480757082}}}, |
| 141 | |
| 142 | {{{26.192053931854691, 9.8504326817814416}, {10.174241480498686, 98.476562741434464}, |
| 143 | {21.177712558385782, 33.814968789841501}, {75.329030899018534, 55.02231980442177}}}, |
| 144 | {{{56.222082700683771, 24.54395039218662}, {95.589995289030483, 81.050822735322086}, |
| 145 | {28.180450866082897, 28.837706255185282}, {60.128952916771617, 87.311672180570511}}}, |
| 146 | |
| 147 | {{{42.449716172390481, 52.379709366885805}, {27.896043159019225, 48.797373636065686}, |
| 148 | {92.770268299044233, 89.899302036454571}, {12.102066544863426, 99.43241951960718}}}, |
| 149 | {{{45.77532924980639, 45.958701495993274}, {37.458701356062065, 68.393691335056758}, |
| 150 | {37.569326692060258, 27.673713456687381}, {60.674866037757539, 62.47349659096146}}}, |
| 151 | |
| 152 | {{{67.426548091427676, 37.993772624988935}, {23.483695892376684, 90.476863174921306}, |
| 153 | {35.597065061143162, 79.872482633158796}, {75.38634169631932, 18.244890038969412}}}, |
| 154 | {{{61.336508189019057, 82.693132843213675}, {44.639380902349664, 54.074825790745592}, |
| 155 | {16.815615499771951, 20.049704667203923}, {41.866884958868326, 56.735503699973002}}}, |
| 156 | |
| 157 | {{{18.1312339, 31.6473732}, {95.5711034, 63.5350219}, {92.3283165, 62.0158945}, |
| 158 | {18.5656052, 32.1268808}}}, |
| 159 | {{{97.402018, 35.7169972}, {33.1127443, 25.8935163}, {1.13970027, 54.9424981}, |
| 160 | {56.4860195, 60.529264}}}, |
| 161 | }; |
| 162 | |
caryclark@google.com | ad65a3e | 2013-04-15 19:13:59 +0000 | [diff] [blame] | 163 | const size_t testSetCount = SK_ARRAY_COUNT(testSet); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 164 | |
| 165 | static const SkDCubic newTestSet[] = { |
caryclark@google.com | 570863f | 2013-09-16 15:55:01 +0000 | [diff] [blame^] | 166 | {{{134, 11414}, {131.990234375, 11414}, {130.32666015625, 11415.482421875}, {130.04275512695312, 11417.4130859375}}}, |
| 167 | {{{132, 11419}, {130.89543151855469, 11419}, {130, 11418.1044921875}, {130, 11417}}}, |
| 168 | |
skia.committer@gmail.com | 7f1af50 | 2013-07-24 07:01:12 +0000 | [diff] [blame] | 169 | {{{3, 4}, {1, 5}, {4, 3}, {6, 4}}}, |
| 170 | {{{3, 4}, {4, 6}, {4, 3}, {5, 1}}}, |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 171 | |
skia.committer@gmail.com | 7f1af50 | 2013-07-24 07:01:12 +0000 | [diff] [blame] | 172 | {{{130.04275512695312, 11417.413085937500 }, |
| 173 | {130.23312377929687, 11418.319335937500 }, |
| 174 | {131.03707885742187, 11419.000000000000 }, |
| 175 | {132.00000000000000, 11419.000000000000 }}}, |
| 176 | |
| 177 | {{{132.00000000000000, 11419.000000000000 }, |
| 178 | {130.89543151855469, 11419.000000000000 }, |
| 179 | {130.00000000000000, 11418.104492187500 }, |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 180 | {130.00000000000000, 11417.000000000000 }}}, |
| 181 | |
robertphillips@google.com | 8dcea7d | 2013-07-15 15:38:55 +0000 | [diff] [blame] | 182 | {{{1.0516976506771041, 2.9684399028541346 }, |
| 183 | {1.0604363140895228, 2.9633503074444141 }, |
| 184 | {1.0692548215065762, 2.9580354426587459 }, |
| 185 | {1.0781560339512140, 2.9525043684031349 }}}, |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 186 | |
robertphillips@google.com | 8dcea7d | 2013-07-15 15:38:55 +0000 | [diff] [blame] | 187 | {{{1.0523038101345104, 2.9523755204833737 }, |
| 188 | {1.0607035288264237, 2.9580853881628375 }, |
| 189 | {1.0690530472271964, 2.9633896794787749 }, |
| 190 | {1.0773566568712512, 2.9682969775000219 }}}, |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 191 | |
robertphillips@google.com | 8dcea7d | 2013-07-15 15:38:55 +0000 | [diff] [blame] | 192 | {{{1.0386522625066592, 2.9759024812329078 }, |
| 193 | {1.0559713690392631, 2.9661782500838885 }, |
| 194 | {1.0736041309019990, 2.9555348259177858 }, |
| 195 | {1.0915734362784633, 2.9440446879826569 }}}, |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 196 | |
robertphillips@google.com | 8dcea7d | 2013-07-15 15:38:55 +0000 | [diff] [blame] | 197 | {{{1.0396670794879301, 2.9435062123457261 }, |
| 198 | {1.0565690546812769, 2.9557413250983462 }, |
| 199 | {1.0732616463413533, 2.9663369676594282 }, |
| 200 | {1.0897791867435489, 2.9753618045797472 }}}, |
| 201 | |
| 202 | {{{0.8685656183311091, 3.0409266475785208 }, |
| 203 | {0.99189542936395292, 3.0212163698184424 }, |
| 204 | {1.1302108367493320, 2.9265646471747306 }, |
| 205 | {1.2952305904872474, 2.7940808546473788 }}}, |
| 206 | |
| 207 | {{{0.85437872843682727, 2.7536036928549055 }, |
| 208 | {1.0045584590592620, 2.9493041024831705 }, |
| 209 | {1.1336998329885613, 3.0248027987251747 }, |
| 210 | {1.2593809752247314, 3.0152560315809107 }}}, |
| 211 | |
| 212 | {{{0, 1}, {1, 6}, {1, 0}, {6, 2}}}, |
| 213 | {{{0, 1}, {2, 6}, {1, 0}, {6, 1}}}, |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 214 | |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 215 | {{{134,11414}, {131.990234375,11414}, {130.32666015625,11415.482421875}, {130.04275512695312,11417.4130859375}}}, |
| 216 | {{{132,11419}, {130.89543151855469,11419}, {130,11418.1044921875}, {130,11417}}}, |
| 217 | |
| 218 | {{{132,11419}, {130.89543151855469,11419}, {130,11418.1044921875}, {130,11417}}}, |
| 219 | {{{130.04275512695312,11417.4130859375}, {130.23312377929687,11418.3193359375}, {131.03707885742187,11419}, {132,11419}}}, |
| 220 | |
caryclark@google.com | a5e5592 | 2013-05-07 18:51:31 +0000 | [diff] [blame] | 221 | {{{0, 1}, {2, 3}, {5, 1}, {4, 3}}}, |
| 222 | {{{1, 5}, {3, 4}, {1, 0}, {3, 2}}}, |
| 223 | |
caryclark@google.com | b3f0921 | 2013-04-17 15:49:16 +0000 | [diff] [blame] | 224 | {{{3, 5}, {1, 6}, {5, 0}, {3, 1}}}, |
| 225 | {{{0, 5}, {1, 3}, {5, 3}, {6, 1}}}, |
| 226 | |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 227 | {{{0, 1}, {1, 5}, {1, 0}, {1, 0}}}, |
| 228 | {{{0, 1}, {0, 1}, {1, 0}, {5, 1}}}, |
| 229 | |
| 230 | {{{1, 3}, {5, 6}, {5, 3}, {5, 4}}}, |
| 231 | {{{3, 5}, {4, 5}, {3, 1}, {6, 5}}}, |
| 232 | |
| 233 | {{{0, 5}, {0, 5}, {5, 4}, {6, 4}}}, |
| 234 | {{{4, 5}, {4, 6}, {5, 0}, {5, 0}}}, |
| 235 | |
| 236 | {{{0, 4}, {1, 3}, {5, 4}, {4, 2}}}, |
| 237 | {{{4, 5}, {2, 4}, {4, 0}, {3, 1}}}, |
| 238 | |
| 239 | {{{0, 2}, {1, 5}, {3, 2}, {4, 1}}}, |
| 240 | {{{2, 3}, {1, 4}, {2, 0}, {5, 1}}}, |
| 241 | |
| 242 | {{{0, 2}, {2, 3}, {5, 1}, {3, 2}}}, |
| 243 | {{{1, 5}, {2, 3}, {2, 0}, {3, 2}}}, |
| 244 | |
| 245 | {{{2, 6}, {4, 5}, {1, 0}, {6, 1}}}, |
| 246 | {{{0, 1}, {1, 6}, {6, 2}, {5, 4}}}, |
| 247 | |
| 248 | {{{0, 1}, {1, 2}, {6, 5}, {5, 4}}}, |
| 249 | {{{5, 6}, {4, 5}, {1, 0}, {2, 1}}}, |
| 250 | |
| 251 | {{{2.5119999999999996, 1.5710000000000002}, {2.6399999999999983, 1.6599999999999997}, |
| 252 | {2.8000000000000007, 1.8000000000000003}, {3, 2}}}, |
| 253 | {{{2.4181876227114887, 1.9849772580462195}, {2.8269904869227211, 2.009330650246834}, |
| 254 | {3.2004679292461624, 1.9942047174679169}, {3.4986199496818058, 2.0035994597094731}}}, |
| 255 | |
| 256 | {{{2, 3}, {1, 4}, {1, 0}, {6, 0}}}, |
| 257 | {{{0, 1}, {0, 6}, {3, 2}, {4, 1}}}, |
| 258 | |
| 259 | {{{0, 2}, {1, 5}, {1, 0}, {6, 1}}}, |
| 260 | {{{0, 1}, {1, 6}, {2, 0}, {5, 1}}}, |
| 261 | |
| 262 | {{{0, 1}, {1, 5}, {2, 1}, {4, 0}}}, |
| 263 | {{{1, 2}, {0, 4}, {1, 0}, {5, 1}}}, |
| 264 | |
| 265 | {{{0, 1}, {3, 5}, {2, 1}, {3, 1}}}, |
| 266 | {{{1, 2}, {1, 3}, {1, 0}, {5, 3}}}, |
| 267 | |
| 268 | {{{0, 1}, {2, 5}, {6, 0}, {5, 3}}}, |
| 269 | {{{0, 6}, {3, 5}, {1, 0}, {5, 2}}}, |
| 270 | |
| 271 | {{{0, 1}, {3, 6}, {1, 0}, {5, 2}}}, |
| 272 | {{{0, 1}, {2, 5}, {1, 0}, {6, 3}}}, |
| 273 | |
| 274 | {{{1, 2}, {5, 6}, {1, 0}, {1, 0}}}, |
| 275 | {{{0, 1}, {0, 1}, {2, 1}, {6, 5}}}, |
| 276 | |
| 277 | {{{0, 6}, {1, 2}, {1, 0}, {1, 0}}}, |
| 278 | {{{0, 1}, {0, 1}, {6, 0}, {2, 1}}}, |
| 279 | |
| 280 | {{{0, 2}, {0, 1}, {3, 0}, {1, 0}}}, |
| 281 | {{{0, 3}, {0, 1}, {2, 0}, {1, 0}}}, |
| 282 | }; |
| 283 | |
caryclark@google.com | ad65a3e | 2013-04-15 19:13:59 +0000 | [diff] [blame] | 284 | const size_t newTestSetCount = SK_ARRAY_COUNT(newTestSet); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 285 | |
| 286 | static void oneOff(skiatest::Reporter* reporter, const SkDCubic& cubic1, const SkDCubic& cubic2) { |
caryclark@google.com | 8d0a524 | 2013-07-16 16:11:16 +0000 | [diff] [blame] | 287 | SkASSERT(ValidCubic(cubic1)); |
| 288 | SkASSERT(ValidCubic(cubic2)); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 289 | #if ONE_OFF_DEBUG |
| 290 | SkDebugf("computed quadratics given\n"); |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 291 | SkDebugf(" {{%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}},\n", |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 292 | cubic1[0].fX, cubic1[0].fY, cubic1[1].fX, cubic1[1].fY, |
| 293 | cubic1[2].fX, cubic1[2].fY, cubic1[3].fX, cubic1[3].fY); |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 294 | SkDebugf(" {{%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}},\n", |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 295 | cubic2[0].fX, cubic2[0].fY, cubic2[1].fX, cubic2[1].fY, |
| 296 | cubic2[2].fX, cubic2[2].fY, cubic2[3].fX, cubic2[3].fY); |
| 297 | #endif |
caryclark@google.com | d892bd8 | 2013-06-17 14:10:36 +0000 | [diff] [blame] | 298 | SkTArray<SkDQuad, true> quads1; |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 299 | CubicToQuads(cubic1, cubic1.calcPrecision(), quads1); |
| 300 | #if ONE_OFF_DEBUG |
| 301 | SkDebugf("computed quadratics set 1\n"); |
| 302 | for (int index = 0; index < quads1.count(); ++index) { |
| 303 | const SkDQuad& q = quads1[index]; |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 304 | SkDebugf(" {{%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}},\n", q[0].fX, q[0].fY, |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 305 | q[1].fX, q[1].fY, q[2].fX, q[2].fY); |
| 306 | } |
| 307 | #endif |
caryclark@google.com | d892bd8 | 2013-06-17 14:10:36 +0000 | [diff] [blame] | 308 | SkTArray<SkDQuad, true> quads2; |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 309 | CubicToQuads(cubic2, cubic2.calcPrecision(), quads2); |
| 310 | #if ONE_OFF_DEBUG |
| 311 | SkDebugf("computed quadratics set 2\n"); |
| 312 | for (int index = 0; index < quads2.count(); ++index) { |
| 313 | const SkDQuad& q = quads2[index]; |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 314 | SkDebugf(" {{%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}},\n", q[0].fX, q[0].fY, |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 315 | q[1].fX, q[1].fY, q[2].fX, q[2].fY); |
| 316 | } |
| 317 | #endif |
| 318 | SkIntersections intersections; |
| 319 | intersections.intersect(cubic1, cubic2); |
| 320 | double tt1, tt2; |
| 321 | SkDPoint xy1, xy2; |
| 322 | for (int pt3 = 0; pt3 < intersections.used(); ++pt3) { |
| 323 | tt1 = intersections[0][pt3]; |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 324 | xy1 = cubic1.ptAtT(tt1); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 325 | tt2 = intersections[1][pt3]; |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 326 | xy2 = cubic2.ptAtT(tt2); |
caryclark@google.com | 5ec2b1e | 2013-07-15 20:50:24 +0000 | [diff] [blame] | 327 | const SkDPoint& iPt = intersections.pt(pt3); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 328 | #if ONE_OFF_DEBUG |
| 329 | SkDebugf("%s t1=%1.9g (%1.9g, %1.9g) (%1.9g, %1.9g) (%1.9g, %1.9g) t2=%1.9g\n", |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 330 | __FUNCTION__, tt1, xy1.fX, xy1.fY, iPt.fX, |
| 331 | iPt.fY, xy2.fX, xy2.fY, tt2); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 332 | #endif |
caryclark@google.com | cffbcc3 | 2013-06-04 17:59:42 +0000 | [diff] [blame] | 333 | REPORTER_ASSERT(reporter, xy1.approximatelyEqual(iPt)); |
| 334 | REPORTER_ASSERT(reporter, xy2.approximatelyEqual(iPt)); |
| 335 | REPORTER_ASSERT(reporter, xy1.approximatelyEqual(xy2)); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 336 | } |
| 337 | } |
| 338 | |
| 339 | static void oneOff(skiatest::Reporter* reporter, int outer, int inner) { |
| 340 | const SkDCubic& cubic1 = testSet[outer]; |
| 341 | const SkDCubic& cubic2 = testSet[inner]; |
| 342 | oneOff(reporter, cubic1, cubic2); |
| 343 | } |
| 344 | |
| 345 | static void newOneOff(skiatest::Reporter* reporter, int outer, int inner) { |
| 346 | const SkDCubic& cubic1 = newTestSet[outer]; |
| 347 | const SkDCubic& cubic2 = newTestSet[inner]; |
| 348 | oneOff(reporter, cubic1, cubic2); |
| 349 | } |
| 350 | |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 351 | static void oneOffTests(skiatest::Reporter* reporter) { |
| 352 | for (size_t outer = 0; outer < testSetCount - 1; ++outer) { |
| 353 | for (size_t inner = outer + 1; inner < testSetCount; ++inner) { |
| 354 | oneOff(reporter, outer, inner); |
| 355 | } |
| 356 | } |
| 357 | for (size_t outer = 0; outer < newTestSetCount - 1; ++outer) { |
| 358 | for (size_t inner = outer + 1; inner < newTestSetCount; ++inner) { |
caryclark@google.com | 5ec2b1e | 2013-07-15 20:50:24 +0000 | [diff] [blame] | 359 | newOneOff(reporter, outer, inner); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 360 | } |
| 361 | } |
| 362 | } |
| 363 | |
| 364 | #define DEBUG_CRASH 0 |
| 365 | |
| 366 | static void CubicIntersection_RandTest(skiatest::Reporter* reporter) { |
| 367 | srand(0); |
| 368 | const int tests = 10000000; |
caryclark@google.com | db60de7 | 2013-04-11 12:33:23 +0000 | [diff] [blame] | 369 | #if !defined(SK_BUILD_FOR_WIN) && !defined(SK_BUILD_FOR_ANDROID) |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 370 | unsigned seed = 0; |
| 371 | #endif |
| 372 | for (int test = 0; test < tests; ++test) { |
| 373 | SkDCubic cubic1, cubic2; |
| 374 | for (int i = 0; i < 4; ++i) { |
| 375 | cubic1[i].fX = static_cast<double>(SK_RAND(seed)) / RAND_MAX * 100; |
| 376 | cubic1[i].fY = static_cast<double>(SK_RAND(seed)) / RAND_MAX * 100; |
| 377 | cubic2[i].fX = static_cast<double>(SK_RAND(seed)) / RAND_MAX * 100; |
| 378 | cubic2[i].fY = static_cast<double>(SK_RAND(seed)) / RAND_MAX * 100; |
| 379 | } |
| 380 | #if DEBUG_CRASH |
| 381 | char str[1024]; |
| 382 | snprintf(str, sizeof(str), |
| 383 | "{{{%1.9g, %1.9g}, {%1.9g, %1.9g}, {%1.9g, %1.9g}, {%1.9g, %1.9g}}},\n" |
| 384 | "{{{%1.9g, %1.9g}, {%1.9g, %1.9g}, {%1.9g, %1.9g}, {%1.9g, %1.9g}}},\n", |
| 385 | cubic1[0].fX, cubic1[0].fY, cubic1[1].fX, cubic1[1].fY, cubic1[2].fX, cubic1[2].fY, |
| 386 | cubic1[3].fX, cubic1[3].fY, |
| 387 | cubic2[0].fX, cubic2[0].fY, cubic2[1].fX, cubic2[1].fY, cubic2[2].fX, cubic2[2].fY, |
| 388 | cubic2[3].fX, cubic2[3].fY); |
| 389 | #endif |
| 390 | SkDRect rect1, rect2; |
| 391 | rect1.setBounds(cubic1); |
| 392 | rect2.setBounds(cubic2); |
| 393 | bool boundsIntersect = rect1.fLeft <= rect2.fRight && rect2.fLeft <= rect2.fRight |
| 394 | && rect1.fTop <= rect2.fBottom && rect2.fTop <= rect1.fBottom; |
| 395 | if (test == -1) { |
| 396 | SkDebugf("ready...\n"); |
| 397 | } |
| 398 | SkIntersections intersections2; |
| 399 | int newIntersects = intersections2.intersect(cubic1, cubic2); |
| 400 | if (!boundsIntersect && newIntersects) { |
| 401 | #if DEBUG_CRASH |
| 402 | SkDebugf("%s %d unexpected intersection boundsIntersect=%d " |
| 403 | " newIntersects=%d\n%s %s\n", __FUNCTION__, test, boundsIntersect, |
| 404 | newIntersects, __FUNCTION__, str); |
| 405 | #endif |
| 406 | REPORTER_ASSERT(reporter, 0); |
| 407 | } |
| 408 | for (int pt = 0; pt < intersections2.used(); ++pt) { |
| 409 | double tt1 = intersections2[0][pt]; |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 410 | SkDPoint xy1 = cubic1.ptAtT(tt1); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 411 | double tt2 = intersections2[1][pt]; |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 412 | SkDPoint xy2 = cubic2.ptAtT(tt2); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 413 | REPORTER_ASSERT(reporter, xy1.approximatelyEqual(xy2)); |
| 414 | } |
| 415 | } |
| 416 | } |
| 417 | |
| 418 | static void intersectionFinder(int index0, int index1, double t1Seed, double t2Seed, |
| 419 | double t1Step, double t2Step) { |
| 420 | const SkDCubic& cubic1 = newTestSet[index0]; |
| 421 | const SkDCubic& cubic2 = newTestSet[index1]; |
| 422 | SkDPoint t1[3], t2[3]; |
| 423 | bool toggle = true; |
| 424 | do { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 425 | t1[0] = cubic1.ptAtT(t1Seed - t1Step); |
| 426 | t1[1] = cubic1.ptAtT(t1Seed); |
| 427 | t1[2] = cubic1.ptAtT(t1Seed + t1Step); |
| 428 | t2[0] = cubic2.ptAtT(t2Seed - t2Step); |
| 429 | t2[1] = cubic2.ptAtT(t2Seed); |
| 430 | t2[2] = cubic2.ptAtT(t2Seed + t2Step); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 431 | double dist[3][3]; |
| 432 | dist[1][1] = t1[1].distance(t2[1]); |
| 433 | int best_i = 1, best_j = 1; |
| 434 | for (int i = 0; i < 3; ++i) { |
| 435 | for (int j = 0; j < 3; ++j) { |
| 436 | if (i == 1 && j == 1) { |
| 437 | continue; |
| 438 | } |
| 439 | dist[i][j] = t1[i].distance(t2[j]); |
| 440 | if (dist[best_i][best_j] > dist[i][j]) { |
| 441 | best_i = i; |
| 442 | best_j = j; |
| 443 | } |
| 444 | } |
| 445 | } |
| 446 | if (best_i == 0) { |
| 447 | t1Seed -= t1Step; |
| 448 | } else if (best_i == 2) { |
| 449 | t1Seed += t1Step; |
| 450 | } |
| 451 | if (best_j == 0) { |
| 452 | t2Seed -= t2Step; |
| 453 | } else if (best_j == 2) { |
| 454 | t2Seed += t2Step; |
| 455 | } |
| 456 | if (best_i == 1 && best_j == 1) { |
| 457 | if ((toggle ^= true)) { |
| 458 | t1Step /= 2; |
| 459 | } else { |
| 460 | t2Step /= 2; |
| 461 | } |
| 462 | } |
| 463 | } while (!t1[1].approximatelyEqual(t2[1])); |
| 464 | t1Step = t2Step = 0.1; |
| 465 | double t10 = t1Seed - t1Step * 2; |
| 466 | double t12 = t1Seed + t1Step * 2; |
| 467 | double t20 = t2Seed - t2Step * 2; |
| 468 | double t22 = t2Seed + t2Step * 2; |
| 469 | SkDPoint test; |
| 470 | while (!approximately_zero(t1Step)) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 471 | test = cubic1.ptAtT(t10); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 472 | t10 += t1[1].approximatelyEqual(test) ? -t1Step : t1Step; |
| 473 | t1Step /= 2; |
| 474 | } |
| 475 | t1Step = 0.1; |
| 476 | while (!approximately_zero(t1Step)) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 477 | test = cubic1.ptAtT(t12); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 478 | t12 -= t1[1].approximatelyEqual(test) ? -t1Step : t1Step; |
| 479 | t1Step /= 2; |
| 480 | } |
| 481 | while (!approximately_zero(t2Step)) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 482 | test = cubic2.ptAtT(t20); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 483 | t20 += t2[1].approximatelyEqual(test) ? -t2Step : t2Step; |
| 484 | t2Step /= 2; |
| 485 | } |
| 486 | t2Step = 0.1; |
| 487 | while (!approximately_zero(t2Step)) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 488 | test = cubic2.ptAtT(t22); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 489 | t22 -= t2[1].approximatelyEqual(test) ? -t2Step : t2Step; |
| 490 | t2Step /= 2; |
| 491 | } |
| 492 | #if ONE_OFF_DEBUG |
| 493 | SkDebugf("%s t1=(%1.9g<%1.9g<%1.9g) t2=(%1.9g<%1.9g<%1.9g)\n", __FUNCTION__, |
| 494 | t10, t1Seed, t12, t20, t2Seed, t22); |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 495 | SkDPoint p10 = cubic1.ptAtT(t10); |
| 496 | SkDPoint p1Seed = cubic1.ptAtT(t1Seed); |
| 497 | SkDPoint p12 = cubic1.ptAtT(t12); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 498 | SkDebugf("%s p1=(%1.9g,%1.9g)<(%1.9g,%1.9g)<(%1.9g,%1.9g)\n", __FUNCTION__, |
| 499 | p10.fX, p10.fY, p1Seed.fX, p1Seed.fY, p12.fX, p12.fY); |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 500 | SkDPoint p20 = cubic2.ptAtT(t20); |
| 501 | SkDPoint p2Seed = cubic2.ptAtT(t2Seed); |
| 502 | SkDPoint p22 = cubic2.ptAtT(t22); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 503 | SkDebugf("%s p2=(%1.9g,%1.9g)<(%1.9g,%1.9g)<(%1.9g,%1.9g)\n", __FUNCTION__, |
| 504 | p20.fX, p20.fY, p2Seed.fX, p2Seed.fY, p22.fX, p22.fY); |
| 505 | #endif |
| 506 | } |
| 507 | |
| 508 | static void CubicIntersection_IntersectionFinder() { |
| 509 | // double t1Seed = 0.87; |
| 510 | // double t2Seed = 0.87; |
| 511 | double t1Step = 0.000001; |
| 512 | double t2Step = 0.000001; |
| 513 | intersectionFinder(0, 1, 0.855895664, 0.864850875, t1Step, t2Step); |
| 514 | intersectionFinder(0, 1, 0.865207906, 0.865207887, t1Step, t2Step); |
| 515 | intersectionFinder(0, 1, 0.865213351, 0.865208087, t1Step, t2Step); |
| 516 | } |
| 517 | |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 518 | static const SkDCubic selfSet[] = { |
| 519 | {{{2, 3}, {0, 4}, {3, 2}, {5, 3}}}, |
| 520 | {{{3, 6}, {2, 3}, {4, 0}, {3, 2}}}, |
| 521 | {{{0, 2}, {2, 3}, {5, 1}, {3, 2}}}, |
| 522 | {{{0, 2}, {3, 5}, {5, 0}, {4, 2}}}, |
| 523 | {{{3.34, 8.98}, {1.95, 10.27}, {3.76, 7.65}, {4.96, 10.64}}}, |
| 524 | {{{3.13, 2.74}, {1.08, 4.62}, {3.71, 0.94}, {2.01, 3.81}}}, |
| 525 | {{{6.71, 3.14}, {7.99, 2.75}, {8.27, 1.96}, {6.35, 3.57}}}, |
| 526 | {{{12.81, 7.27}, {7.22, 6.98}, {12.49, 8.97}, {11.42, 6.18}}}, |
| 527 | }; |
| 528 | size_t selfSetCount = SK_ARRAY_COUNT(selfSet); |
| 529 | |
| 530 | static void selfOneOff(skiatest::Reporter* reporter, int index) { |
| 531 | const SkDCubic& cubic = selfSet[index]; |
| 532 | #if ONE_OFF_DEBUG |
| 533 | int idx2; |
| 534 | double max[3]; |
| 535 | int ts = cubic.findMaxCurvature(max); |
| 536 | for (idx2 = 0; idx2 < ts; ++idx2) { |
| 537 | SkDebugf("%s max[%d]=%1.9g (%1.9g, %1.9g)\n", __FUNCTION__, idx2, |
| 538 | max[idx2], cubic.ptAtT(max[idx2]).fX, cubic.ptAtT(max[idx2]).fY); |
| 539 | } |
| 540 | SkTArray<double, true> ts1; |
| 541 | SkTArray<SkDQuad, true> quads1; |
| 542 | cubic.toQuadraticTs(cubic.calcPrecision(), &ts1); |
| 543 | for (idx2 = 0; idx2 < ts1.count(); ++idx2) { |
| 544 | SkDebugf("%s t[%d]=%1.9g\n", __FUNCTION__, idx2, ts1[idx2]); |
| 545 | } |
| 546 | CubicToQuads(cubic, cubic.calcPrecision(), quads1); |
| 547 | for (idx2 = 0; idx2 < quads1.count(); ++idx2) { |
| 548 | const SkDQuad& q = quads1[idx2]; |
| 549 | SkDebugf(" {{{%1.9g,%1.9g}, {%1.9g,%1.9g}, {%1.9g,%1.9g}}},\n", |
| 550 | q[0].fX, q[0].fY, q[1].fX, q[1].fY, q[2].fX, q[2].fY); |
| 551 | } |
| 552 | SkDebugf("\n"); |
| 553 | #endif |
| 554 | SkIntersections i; |
| 555 | int result = i.intersect(cubic); |
| 556 | REPORTER_ASSERT(reporter, result == 1); |
| 557 | REPORTER_ASSERT(reporter, i.used() == 1); |
| 558 | REPORTER_ASSERT(reporter, !approximately_equal(i[0][0], i[1][0])); |
| 559 | SkDPoint pt1 = cubic.ptAtT(i[0][0]); |
| 560 | SkDPoint pt2 = cubic.ptAtT(i[1][0]); |
| 561 | REPORTER_ASSERT(reporter, pt1.approximatelyEqual(pt2)); |
| 562 | } |
| 563 | |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 564 | static void cubicIntersectionSelfTest(skiatest::Reporter* reporter) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 565 | size_t firstFail = 0; |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 566 | for (size_t index = firstFail; index < selfSetCount; ++index) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 567 | selfOneOff(reporter, index); |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 568 | } |
| 569 | } |
| 570 | |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 571 | static void PathOpsCubicIntersectionOneOffTest(skiatest::Reporter* reporter) { |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 572 | newOneOff(reporter, 0, 1); |
| 573 | } |
| 574 | |
| 575 | static void PathOpsCubicSelfOneOffTest(skiatest::Reporter* reporter) { |
| 576 | selfOneOff(reporter, 0); |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 577 | } |
| 578 | |
caryclark@google.com | ad65a3e | 2013-04-15 19:13:59 +0000 | [diff] [blame] | 579 | static void PathOpsCubicIntersectionTest(skiatest::Reporter* reporter) { |
caryclark@google.com | 9166dcb | 2013-04-08 11:50:00 +0000 | [diff] [blame] | 580 | oneOffTests(reporter); |
| 581 | cubicIntersectionSelfTest(reporter); |
| 582 | standardTestCases(reporter); |
| 583 | if (false) CubicIntersection_IntersectionFinder(); |
| 584 | if (false) CubicIntersection_RandTest(reporter); |
| 585 | } |
| 586 | |
| 587 | #include "TestClassDef.h" |
caryclark@google.com | ad65a3e | 2013-04-15 19:13:59 +0000 | [diff] [blame] | 588 | DEFINE_TESTCLASS_SHORT(PathOpsCubicIntersectionTest) |
caryclark@google.com | fa2aeee | 2013-07-15 13:29:13 +0000 | [diff] [blame] | 589 | |
| 590 | DEFINE_TESTCLASS_SHORT(PathOpsCubicIntersectionOneOffTest) |
caryclark@google.com | 4fdbb22 | 2013-07-23 15:27:41 +0000 | [diff] [blame] | 591 | |
| 592 | DEFINE_TESTCLASS_SHORT(PathOpsCubicSelfOneOffTest) |