SkAnalyticEdge.cpp 26 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512
  1. /*
  2. * Copyright 2006 The Android Open Source Project
  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 "include/private/SkTo.h"
  8. #include "src/core/SkAnalyticEdge.h"
  9. #include "src/core/SkFDot6.h"
  10. #include "src/core/SkMathPriv.h"
  11. #include <utility>
  12. static const int kInverseTableSize = 1024; // SK_FDot6One * 16
  13. static inline SkFixed quick_inverse(SkFDot6 x) {
  14. SkASSERT(SkAbs32(x) < kInverseTableSize);
  15. static const int32_t table[kInverseTableSize * 2] = {
  16. -4096, -4100, -4104, -4108, -4112, -4116, -4120, -4124, -4128, -4132, -4136,
  17. -4140, -4144, -4148, -4152, -4156, -4161, -4165, -4169, -4173, -4177, -4181,
  18. -4185, -4190, -4194, -4198, -4202, -4206, -4211, -4215, -4219, -4223, -4228,
  19. -4232, -4236, -4240, -4245, -4249, -4253, -4258, -4262, -4266, -4271, -4275,
  20. -4279, -4284, -4288, -4293, -4297, -4301, -4306, -4310, -4315, -4319, -4324,
  21. -4328, -4332, -4337, -4341, -4346, -4350, -4355, -4359, -4364, -4369, -4373,
  22. -4378, -4382, -4387, -4391, -4396, -4401, -4405, -4410, -4415, -4419, -4424,
  23. -4429, -4433, -4438, -4443, -4447, -4452, -4457, -4462, -4466, -4471, -4476,
  24. -4481, -4485, -4490, -4495, -4500, -4505, -4510, -4514, -4519, -4524, -4529,
  25. -4534, -4539, -4544, -4549, -4554, -4559, -4563, -4568, -4573, -4578, -4583,
  26. -4588, -4593, -4599, -4604, -4609, -4614, -4619, -4624, -4629, -4634, -4639,
  27. -4644, -4650, -4655, -4660, -4665, -4670, -4675, -4681, -4686, -4691, -4696,
  28. -4702, -4707, -4712, -4718, -4723, -4728, -4733, -4739, -4744, -4750, -4755,
  29. -4760, -4766, -4771, -4777, -4782, -4788, -4793, -4798, -4804, -4809, -4815,
  30. -4821, -4826, -4832, -4837, -4843, -4848, -4854, -4860, -4865, -4871, -4877,
  31. -4882, -4888, -4894, -4899, -4905, -4911, -4917, -4922, -4928, -4934, -4940,
  32. -4946, -4951, -4957, -4963, -4969, -4975, -4981, -4987, -4993, -4999, -5005,
  33. -5011, -5017, -5023, -5029, -5035, -5041, -5047, -5053, -5059, -5065, -5071,
  34. -5077, -5084, -5090, -5096, -5102, -5108, -5115, -5121, -5127, -5133, -5140,
  35. -5146, -5152, -5159, -5165, -5171, -5178, -5184, -5190, -5197, -5203, -5210,
  36. -5216, -5223, -5229, -5236, -5242, -5249, -5256, -5262, -5269, -5275, -5282,
  37. -5289, -5295, -5302, -5309, -5315, -5322, -5329, -5336, -5343, -5349, -5356,
  38. -5363, -5370, -5377, -5384, -5391, -5398, -5405, -5412, -5418, -5426, -5433,
  39. -5440, -5447, -5454, -5461, -5468, -5475, -5482, -5489, -5497, -5504, -5511,
  40. -5518, -5526, -5533, -5540, -5548, -5555, -5562, -5570, -5577, -5584, -5592,
  41. -5599, -5607, -5614, -5622, -5629, -5637, -5645, -5652, -5660, -5667, -5675,
  42. -5683, -5691, -5698, -5706, -5714, -5722, -5729, -5737, -5745, -5753, -5761,
  43. -5769, -5777, -5785, -5793, -5801, -5809, -5817, -5825, -5833, -5841, -5849,
  44. -5857, -5866, -5874, -5882, -5890, -5899, -5907, -5915, -5924, -5932, -5940,
  45. -5949, -5957, -5966, -5974, -5983, -5991, -6000, -6009, -6017, -6026, -6034,
  46. -6043, -6052, -6061, -6069, -6078, -6087, -6096, -6105, -6114, -6123, -6132,
  47. -6141, -6150, -6159, -6168, -6177, -6186, -6195, -6204, -6213, -6223, -6232,
  48. -6241, -6250, -6260, -6269, -6278, -6288, -6297, -6307, -6316, -6326, -6335,
  49. -6345, -6355, -6364, -6374, -6384, -6393, -6403, -6413, -6423, -6432, -6442,
  50. -6452, -6462, -6472, -6482, -6492, -6502, -6512, -6523, -6533, -6543, -6553,
  51. -6563, -6574, -6584, -6594, -6605, -6615, -6626, -6636, -6647, -6657, -6668,
  52. -6678, -6689, -6700, -6710, -6721, -6732, -6743, -6754, -6765, -6775, -6786,
  53. -6797, -6808, -6820, -6831, -6842, -6853, -6864, -6875, -6887, -6898, -6909,
  54. -6921, -6932, -6944, -6955, -6967, -6978, -6990, -7002, -7013, -7025, -7037,
  55. -7049, -7061, -7073, -7084, -7096, -7108, -7121, -7133, -7145, -7157, -7169,
  56. -7182, -7194, -7206, -7219, -7231, -7244, -7256, -7269, -7281, -7294, -7307,
  57. -7319, -7332, -7345, -7358, -7371, -7384, -7397, -7410, -7423, -7436, -7449,
  58. -7463, -7476, -7489, -7503, -7516, -7530, -7543, -7557, -7570, -7584, -7598,
  59. -7612, -7626, -7639, -7653, -7667, -7681, -7695, -7710, -7724, -7738, -7752,
  60. -7767, -7781, -7796, -7810, -7825, -7839, -7854, -7869, -7884, -7898, -7913,
  61. -7928, -7943, -7958, -7973, -7989, -8004, -8019, -8035, -8050, -8065, -8081,
  62. -8097, -8112, -8128, -8144, -8160, -8176, -8192, -8208, -8224, -8240, -8256,
  63. -8272, -8289, -8305, -8322, -8338, -8355, -8371, -8388, -8405, -8422, -8439,
  64. -8456, -8473, -8490, -8507, -8525, -8542, -8559, -8577, -8594, -8612, -8630,
  65. -8648, -8665, -8683, -8701, -8719, -8738, -8756, -8774, -8793, -8811, -8830,
  66. -8848, -8867, -8886, -8905, -8924, -8943, -8962, -8981, -9000, -9020, -9039,
  67. -9058, -9078, -9098, -9118, -9137, -9157, -9177, -9198, -9218, -9238, -9258,
  68. -9279, -9300, -9320, -9341, -9362, -9383, -9404, -9425, -9446, -9467, -9489,
  69. -9510, -9532, -9554, -9576, -9597, -9619, -9642, -9664, -9686, -9709, -9731,
  70. -9754, -9776, -9799, -9822, -9845, -9868, -9892, -9915, -9939, -9962, -9986,
  71. -10010, -10034, -10058, -10082, -10106, -10131, -10155, -10180, -10205, -10230,
  72. -10255, -10280, -10305, -10330, -10356, -10381, -10407, -10433, -10459, -10485,
  73. -10512, -10538, -10564, -10591, -10618, -10645, -10672, -10699, -10727, -10754,
  74. -10782, -10810, -10837, -10866, -10894, -10922, -10951, -10979, -11008, -11037,
  75. -11066, -11096, -11125, -11155, -11184, -11214, -11244, -11275, -11305, -11335,
  76. -11366, -11397, -11428, -11459, -11491, -11522, -11554, -11586, -11618, -11650,
  77. -11683, -11715, -11748, -11781, -11814, -11848, -11881, -11915, -11949, -11983,
  78. -12018, -12052, -12087, -12122, -12157, -12192, -12228, -12264, -12300, -12336,
  79. -12372, -12409, -12446, -12483, -12520, -12557, -12595, -12633, -12671, -12710,
  80. -12748, -12787, -12826, -12865, -12905, -12945, -12985, -13025, -13066, -13107,
  81. -13148, -13189, -13231, -13273, -13315, -13357, -13400, -13443, -13486, -13530,
  82. -13573, -13617, -13662, -13706, -13751, -13797, -13842, -13888, -13934, -13981,
  83. -14027, -14074, -14122, -14169, -14217, -14266, -14315, -14364, -14413, -14463,
  84. -14513, -14563, -14614, -14665, -14716, -14768, -14820, -14873, -14926, -14979,
  85. -15033, -15087, -15141, -15196, -15252, -15307, -15363, -15420, -15477, -15534,
  86. -15592, -15650, -15709, -15768, -15827, -15887, -15947, -16008, -16070, -16131,
  87. -16194, -16256, -16320, -16384, -16448, -16513, -16578, -16644, -16710, -16777,
  88. -16844, -16912, -16980, -17050, -17119, -17189, -17260, -17331, -17403, -17476,
  89. -17549, -17623, -17697, -17772, -17848, -17924, -18001, -18078, -18157, -18236,
  90. -18315, -18396, -18477, -18558, -18641, -18724, -18808, -18893, -18978, -19065,
  91. -19152, -19239, -19328, -19418, -19508, -19599, -19691, -19784, -19878, -19972,
  92. -20068, -20164, -20262, -20360, -20460, -20560, -20661, -20763, -20867, -20971,
  93. -21076, -21183, -21290, -21399, -21509, -21620, -21732, -21845, -21959, -22075,
  94. -22192, -22310, -22429, -22550, -22671, -22795, -22919, -23045, -23172, -23301,
  95. -23431, -23563, -23696, -23831, -23967, -24105, -24244, -24385, -24528, -24672,
  96. -24818, -24966, -25115, -25266, -25420, -25575, -25731, -25890, -26051, -26214,
  97. -26379, -26546, -26715, -26886, -27060, -27235, -27413, -27594, -27776, -27962,
  98. -28149, -28339, -28532, -28728, -28926, -29127, -29330, -29537, -29746, -29959,
  99. -30174, -30393, -30615, -30840, -31068, -31300, -31536, -31775, -32017, -32263,
  100. -32513, -32768, -33026, -33288, -33554, -33825, -34100, -34379, -34663, -34952,
  101. -35246, -35544, -35848, -36157, -36472, -36792, -37117, -37449, -37786, -38130,
  102. -38479, -38836, -39199, -39568, -39945, -40329, -40721, -41120, -41527, -41943,
  103. -42366, -42799, -43240, -43690, -44150, -44620, -45100, -45590, -46091, -46603,
  104. -47127, -47662, -48210, -48770, -49344, -49932, -50533, -51150, -51781, -52428,
  105. -53092, -53773, -54471, -55188, -55924, -56679, -57456, -58254, -59074, -59918,
  106. -60787, -61680, -62601, -63550, -64527, -65536, -66576, -67650, -68759, -69905,
  107. -71089, -72315, -73584, -74898, -76260, -77672, -79137, -80659, -82241, -83886,
  108. -85598, -87381, -89240, -91180, -93206, -95325, -97541, -99864, -102300,
  109. -104857, -107546, -110376, -113359, -116508, -119837, -123361, -127100, -131072,
  110. -135300, -139810, -144631, -149796, -155344, -161319, -167772, -174762, -182361,
  111. -190650, -199728, -209715, -220752, -233016, -246723, -262144, -279620, -299593,
  112. -322638, -349525, -381300, -419430, -466033, -524288, -599186, -699050, -838860,
  113. -1048576, -1398101, -2097152, -4194304, 0, 4194304, 2097152, 1398101, 1048576,
  114. 838860, 699050, 599186, 524288, 466033, 419430, 381300, 349525, 322638, 299593,
  115. 279620, 262144, 246723, 233016, 220752, 209715, 199728, 190650, 182361, 174762,
  116. 167772, 161319, 155344, 149796, 144631, 139810, 135300, 131072, 127100, 123361,
  117. 119837, 116508, 113359, 110376, 107546, 104857, 102300, 99864, 97541, 95325,
  118. 93206, 91180, 89240, 87381, 85598, 83886, 82241, 80659, 79137, 77672, 76260,
  119. 74898, 73584, 72315, 71089, 69905, 68759, 67650, 66576, 65536, 64527, 63550,
  120. 62601, 61680, 60787, 59918, 59074, 58254, 57456, 56679, 55924, 55188, 54471,
  121. 53773, 53092, 52428, 51781, 51150, 50533, 49932, 49344, 48770, 48210, 47662,
  122. 47127, 46603, 46091, 45590, 45100, 44620, 44150, 43690, 43240, 42799, 42366,
  123. 41943, 41527, 41120, 40721, 40329, 39945, 39568, 39199, 38836, 38479, 38130,
  124. 37786, 37449, 37117, 36792, 36472, 36157, 35848, 35544, 35246, 34952, 34663,
  125. 34379, 34100, 33825, 33554, 33288, 33026, 32768, 32513, 32263, 32017, 31775,
  126. 31536, 31300, 31068, 30840, 30615, 30393, 30174, 29959, 29746, 29537, 29330,
  127. 29127, 28926, 28728, 28532, 28339, 28149, 27962, 27776, 27594, 27413, 27235,
  128. 27060, 26886, 26715, 26546, 26379, 26214, 26051, 25890, 25731, 25575, 25420,
  129. 25266, 25115, 24966, 24818, 24672, 24528, 24385, 24244, 24105, 23967, 23831,
  130. 23696, 23563, 23431, 23301, 23172, 23045, 22919, 22795, 22671, 22550, 22429,
  131. 22310, 22192, 22075, 21959, 21845, 21732, 21620, 21509, 21399, 21290, 21183,
  132. 21076, 20971, 20867, 20763, 20661, 20560, 20460, 20360, 20262, 20164, 20068,
  133. 19972, 19878, 19784, 19691, 19599, 19508, 19418, 19328, 19239, 19152, 19065,
  134. 18978, 18893, 18808, 18724, 18641, 18558, 18477, 18396, 18315, 18236, 18157,
  135. 18078, 18001, 17924, 17848, 17772, 17697, 17623, 17549, 17476, 17403, 17331,
  136. 17260, 17189, 17119, 17050, 16980, 16912, 16844, 16777, 16710, 16644, 16578,
  137. 16513, 16448, 16384, 16320, 16256, 16194, 16131, 16070, 16008, 15947, 15887,
  138. 15827, 15768, 15709, 15650, 15592, 15534, 15477, 15420, 15363, 15307, 15252,
  139. 15196, 15141, 15087, 15033, 14979, 14926, 14873, 14820, 14768, 14716, 14665,
  140. 14614, 14563, 14513, 14463, 14413, 14364, 14315, 14266, 14217, 14169, 14122,
  141. 14074, 14027, 13981, 13934, 13888, 13842, 13797, 13751, 13706, 13662, 13617,
  142. 13573, 13530, 13486, 13443, 13400, 13357, 13315, 13273, 13231, 13189, 13148,
  143. 13107, 13066, 13025, 12985, 12945, 12905, 12865, 12826, 12787, 12748, 12710,
  144. 12671, 12633, 12595, 12557, 12520, 12483, 12446, 12409, 12372, 12336, 12300,
  145. 12264, 12228, 12192, 12157, 12122, 12087, 12052, 12018, 11983, 11949, 11915,
  146. 11881, 11848, 11814, 11781, 11748, 11715, 11683, 11650, 11618, 11586, 11554,
  147. 11522, 11491, 11459, 11428, 11397, 11366, 11335, 11305, 11275, 11244, 11214,
  148. 11184, 11155, 11125, 11096, 11066, 11037, 11008, 10979, 10951, 10922, 10894,
  149. 10866, 10837, 10810, 10782, 10754, 10727, 10699, 10672, 10645, 10618, 10591,
  150. 10564, 10538, 10512, 10485, 10459, 10433, 10407, 10381, 10356, 10330, 10305,
  151. 10280, 10255, 10230, 10205, 10180, 10155, 10131, 10106, 10082, 10058, 10034,
  152. 10010, 9986, 9962, 9939, 9915, 9892, 9868, 9845, 9822, 9799, 9776, 9754, 9731,
  153. 9709, 9686, 9664, 9642, 9619, 9597, 9576, 9554, 9532, 9510, 9489, 9467, 9446,
  154. 9425, 9404, 9383, 9362, 9341, 9320, 9300, 9279, 9258, 9238, 9218, 9198, 9177,
  155. 9157, 9137, 9118, 9098, 9078, 9058, 9039, 9020, 9000, 8981, 8962, 8943, 8924,
  156. 8905, 8886, 8867, 8848, 8830, 8811, 8793, 8774, 8756, 8738, 8719, 8701, 8683,
  157. 8665, 8648, 8630, 8612, 8594, 8577, 8559, 8542, 8525, 8507, 8490, 8473, 8456,
  158. 8439, 8422, 8405, 8388, 8371, 8355, 8338, 8322, 8305, 8289, 8272, 8256, 8240,
  159. 8224, 8208, 8192, 8176, 8160, 8144, 8128, 8112, 8097, 8081, 8065, 8050, 8035,
  160. 8019, 8004, 7989, 7973, 7958, 7943, 7928, 7913, 7898, 7884, 7869, 7854, 7839,
  161. 7825, 7810, 7796, 7781, 7767, 7752, 7738, 7724, 7710, 7695, 7681, 7667, 7653,
  162. 7639, 7626, 7612, 7598, 7584, 7570, 7557, 7543, 7530, 7516, 7503, 7489, 7476,
  163. 7463, 7449, 7436, 7423, 7410, 7397, 7384, 7371, 7358, 7345, 7332, 7319, 7307,
  164. 7294, 7281, 7269, 7256, 7244, 7231, 7219, 7206, 7194, 7182, 7169, 7157, 7145,
  165. 7133, 7121, 7108, 7096, 7084, 7073, 7061, 7049, 7037, 7025, 7013, 7002, 6990,
  166. 6978, 6967, 6955, 6944, 6932, 6921, 6909, 6898, 6887, 6875, 6864, 6853, 6842,
  167. 6831, 6820, 6808, 6797, 6786, 6775, 6765, 6754, 6743, 6732, 6721, 6710, 6700,
  168. 6689, 6678, 6668, 6657, 6647, 6636, 6626, 6615, 6605, 6594, 6584, 6574, 6563,
  169. 6553, 6543, 6533, 6523, 6512, 6502, 6492, 6482, 6472, 6462, 6452, 6442, 6432,
  170. 6423, 6413, 6403, 6393, 6384, 6374, 6364, 6355, 6345, 6335, 6326, 6316, 6307,
  171. 6297, 6288, 6278, 6269, 6260, 6250, 6241, 6232, 6223, 6213, 6204, 6195, 6186,
  172. 6177, 6168, 6159, 6150, 6141, 6132, 6123, 6114, 6105, 6096, 6087, 6078, 6069,
  173. 6061, 6052, 6043, 6034, 6026, 6017, 6009, 6000, 5991, 5983, 5974, 5966, 5957,
  174. 5949, 5940, 5932, 5924, 5915, 5907, 5899, 5890, 5882, 5874, 5866, 5857, 5849,
  175. 5841, 5833, 5825, 5817, 5809, 5801, 5793, 5785, 5777, 5769, 5761, 5753, 5745,
  176. 5737, 5729, 5722, 5714, 5706, 5698, 5691, 5683, 5675, 5667, 5660, 5652, 5645,
  177. 5637, 5629, 5622, 5614, 5607, 5599, 5592, 5584, 5577, 5570, 5562, 5555, 5548,
  178. 5540, 5533, 5526, 5518, 5511, 5504, 5497, 5489, 5482, 5475, 5468, 5461, 5454,
  179. 5447, 5440, 5433, 5426, 5418, 5412, 5405, 5398, 5391, 5384, 5377, 5370, 5363,
  180. 5356, 5349, 5343, 5336, 5329, 5322, 5315, 5309, 5302, 5295, 5289, 5282, 5275,
  181. 5269, 5262, 5256, 5249, 5242, 5236, 5229, 5223, 5216, 5210, 5203, 5197, 5190,
  182. 5184, 5178, 5171, 5165, 5159, 5152, 5146, 5140, 5133, 5127, 5121, 5115, 5108,
  183. 5102, 5096, 5090, 5084, 5077, 5071, 5065, 5059, 5053, 5047, 5041, 5035, 5029,
  184. 5023, 5017, 5011, 5005, 4999, 4993, 4987, 4981, 4975, 4969, 4963, 4957, 4951,
  185. 4946, 4940, 4934, 4928, 4922, 4917, 4911, 4905, 4899, 4894, 4888, 4882, 4877,
  186. 4871, 4865, 4860, 4854, 4848, 4843, 4837, 4832, 4826, 4821, 4815, 4809, 4804,
  187. 4798, 4793, 4788, 4782, 4777, 4771, 4766, 4760, 4755, 4750, 4744, 4739, 4733,
  188. 4728, 4723, 4718, 4712, 4707, 4702, 4696, 4691, 4686, 4681, 4675, 4670, 4665,
  189. 4660, 4655, 4650, 4644, 4639, 4634, 4629, 4624, 4619, 4614, 4609, 4604, 4599,
  190. 4593, 4588, 4583, 4578, 4573, 4568, 4563, 4559, 4554, 4549, 4544, 4539, 4534,
  191. 4529, 4524, 4519, 4514, 4510, 4505, 4500, 4495, 4490, 4485, 4481, 4476, 4471,
  192. 4466, 4462, 4457, 4452, 4447, 4443, 4438, 4433, 4429, 4424, 4419, 4415, 4410,
  193. 4405, 4401, 4396, 4391, 4387, 4382, 4378, 4373, 4369, 4364, 4359, 4355, 4350,
  194. 4346, 4341, 4337, 4332, 4328, 4324, 4319, 4315, 4310, 4306, 4301, 4297, 4293,
  195. 4288, 4284, 4279, 4275, 4271, 4266, 4262, 4258, 4253, 4249, 4245, 4240, 4236,
  196. 4232, 4228, 4223, 4219, 4215, 4211, 4206, 4202, 4198, 4194, 4190, 4185, 4181,
  197. 4177, 4173, 4169, 4165, 4161, 4156, 4152, 4148, 4144, 4140, 4136, 4132, 4128,
  198. 4124, 4120, 4116, 4112, 4108, 4104, 4100
  199. };
  200. return table[kInverseTableSize + x];
  201. }
  202. static inline SkFixed quick_div(SkFDot6 a, SkFDot6 b) {
  203. const int kMinBits = 3; // abs(b) should be at least (1 << kMinBits) for quick division
  204. const int kMaxBits = 31; // Number of bits available in signed int
  205. // Given abs(b) <= (1 << kMinBits), the inverse of abs(b) is at most 1 << (22 - kMinBits) in
  206. // SkFixed format. Hence abs(a) should be less than kMaxAbsA
  207. const int kMaxAbsA = 1 << (kMaxBits - (22 - kMinBits));
  208. SkFDot6 abs_a = SkAbs32(a);
  209. SkFDot6 abs_b = SkAbs32(b);
  210. if (abs_b >= (1 << kMinBits) && abs_b < kInverseTableSize && abs_a < kMaxAbsA) {
  211. SkASSERT((int64_t)a * quick_inverse(b) <= SK_MaxS32
  212. && (int64_t)a * quick_inverse(b) >= SK_MinS32);
  213. SkFixed ourAnswer = (a * quick_inverse(b)) >> 6;
  214. SkASSERT(
  215. (SkFDot6Div(a,b) == 0 && ourAnswer == 0) ||
  216. SkFixedDiv(SkAbs32(SkFDot6Div(a,b) - ourAnswer), SkAbs32(SkFDot6Div(a,b))) <= 1 << 10
  217. );
  218. return ourAnswer;
  219. }
  220. return SkFDot6Div(a, b);
  221. }
  222. bool SkAnalyticEdge::setLine(const SkPoint& p0, const SkPoint& p1) {
  223. fRiteE = nullptr;
  224. // We must set X/Y using the same way (e.g., times 4, to FDot6, then to Fixed) as Quads/Cubics.
  225. // Otherwise the order of the edge might be wrong due to precision limit.
  226. const int accuracy = kDefaultAccuracy;
  227. #ifdef SK_RASTERIZE_EVEN_ROUNDING
  228. SkFixed x0 = SkFDot6ToFixed(SkScalarRoundToFDot6(p0.fX, accuracy)) >> accuracy;
  229. SkFixed y0 = SnapY(SkFDot6ToFixed(SkScalarRoundToFDot6(p0.fY, accuracy)) >> accuracy);
  230. SkFixed x1 = SkFDot6ToFixed(SkScalarRoundToFDot6(p1.fX, accuracy)) >> accuracy;
  231. SkFixed y1 = SnapY(SkFDot6ToFixed(SkScalarRoundToFDot6(p1.fY, accuracy)) >> accuracy);
  232. #else
  233. const int multiplier = (1 << kDefaultAccuracy);
  234. SkFixed x0 = SkFDot6ToFixed(SkScalarToFDot6(p0.fX * multiplier)) >> accuracy;
  235. SkFixed y0 = SnapY(SkFDot6ToFixed(SkScalarToFDot6(p0.fY * multiplier)) >> accuracy);
  236. SkFixed x1 = SkFDot6ToFixed(SkScalarToFDot6(p1.fX * multiplier)) >> accuracy;
  237. SkFixed y1 = SnapY(SkFDot6ToFixed(SkScalarToFDot6(p1.fY * multiplier)) >> accuracy);
  238. #endif
  239. int winding = 1;
  240. if (y0 > y1) {
  241. using std::swap;
  242. swap(x0, x1);
  243. swap(y0, y1);
  244. winding = -1;
  245. }
  246. // are we a zero-height line?
  247. SkFDot6 dy = SkFixedToFDot6(y1 - y0);
  248. if (dy == 0) {
  249. return false;
  250. }
  251. SkFDot6 dx = SkFixedToFDot6(x1 - x0);
  252. SkFixed slope = quick_div(dx, dy);
  253. SkFixed absSlope = SkAbs32(slope);
  254. fX = x0;
  255. fDX = slope;
  256. fUpperX = x0;
  257. fY = y0;
  258. fUpperY = y0;
  259. fLowerY = y1;
  260. fDY = dx == 0 || slope == 0 ? SK_MaxS32 : absSlope < kInverseTableSize
  261. ? quick_inverse(absSlope)
  262. : SkAbs32(quick_div(dy, dx));
  263. fCurveCount = 0;
  264. fWinding = SkToS8(winding);
  265. fCurveShift = 0;
  266. return true;
  267. }
  268. // This will become a bottleneck for small ovals rendering if we call SkFixedDiv twice here.
  269. // Therefore, we'll let the outter function compute the slope once and send in the value.
  270. // Moreover, we'll compute fDY by quickly lookup the inverse table (if possible).
  271. bool SkAnalyticEdge::updateLine(SkFixed x0, SkFixed y0, SkFixed x1, SkFixed y1, SkFixed slope) {
  272. // Since we send in the slope, we can no longer snap y inside this function.
  273. // If we don't send in the slope, or we do some more sophisticated snapping, this function
  274. // could be a performance bottleneck.
  275. SkASSERT(fWinding == 1 || fWinding == -1);
  276. SkASSERT(fCurveCount != 0);
  277. // We don't chop at y extrema for cubics so the y is not guaranteed to be increasing for them.
  278. // In that case, we have to swap x/y and negate the winding.
  279. if (y0 > y1) {
  280. using std::swap;
  281. swap(x0, x1);
  282. swap(y0, y1);
  283. fWinding = -fWinding;
  284. }
  285. SkASSERT(y0 <= y1);
  286. SkFDot6 dx = SkFixedToFDot6(x1 - x0);
  287. SkFDot6 dy = SkFixedToFDot6(y1 - y0);
  288. // are we a zero-height line?
  289. if (dy == 0) {
  290. return false;
  291. }
  292. SkASSERT(slope < SK_MaxS32);
  293. SkFDot6 absSlope = SkAbs32(SkFixedToFDot6(slope));
  294. fX = x0;
  295. fDX = slope;
  296. fUpperX = x0;
  297. fY = y0;
  298. fUpperY = y0;
  299. fLowerY = y1;
  300. fDY = (dx == 0 || slope == 0)
  301. ? SK_MaxS32
  302. : absSlope < kInverseTableSize
  303. ? quick_inverse(absSlope)
  304. : SkAbs32(quick_div(dy, dx));
  305. return true;
  306. }
  307. bool SkAnalyticEdge::update(SkFixed last_y, bool sortY) {
  308. SkASSERT(last_y >= fLowerY); // we shouldn't update edge if last_y < fLowerY
  309. if (fCurveCount < 0) {
  310. return static_cast<SkAnalyticCubicEdge*>(this)->updateCubic(sortY);
  311. } else if (fCurveCount > 0) {
  312. return static_cast<SkAnalyticQuadraticEdge*>(this)->updateQuadratic();
  313. }
  314. return false;
  315. }
  316. bool SkAnalyticQuadraticEdge::setQuadratic(const SkPoint pts[3]) {
  317. fRiteE = nullptr;
  318. if (!fQEdge.setQuadraticWithoutUpdate(pts, kDefaultAccuracy)) {
  319. return false;
  320. }
  321. fQEdge.fQx >>= kDefaultAccuracy;
  322. fQEdge.fQy >>= kDefaultAccuracy;
  323. fQEdge.fQDx >>= kDefaultAccuracy;
  324. fQEdge.fQDy >>= kDefaultAccuracy;
  325. fQEdge.fQDDx >>= kDefaultAccuracy;
  326. fQEdge.fQDDy >>= kDefaultAccuracy;
  327. fQEdge.fQLastX >>= kDefaultAccuracy;
  328. fQEdge.fQLastY >>= kDefaultAccuracy;
  329. fQEdge.fQy = SnapY(fQEdge.fQy);
  330. fQEdge.fQLastY = SnapY(fQEdge.fQLastY);
  331. fWinding = fQEdge.fWinding;
  332. fCurveCount = fQEdge.fCurveCount;
  333. fCurveShift = fQEdge.fCurveShift;
  334. fSnappedX = fQEdge.fQx;
  335. fSnappedY = fQEdge.fQy;
  336. return this->updateQuadratic();
  337. }
  338. bool SkAnalyticQuadraticEdge::updateQuadratic() {
  339. int success = 0; // initialize to fail!
  340. int count = fCurveCount;
  341. SkFixed oldx = fQEdge.fQx;
  342. SkFixed oldy = fQEdge.fQy;
  343. SkFixed dx = fQEdge.fQDx;
  344. SkFixed dy = fQEdge.fQDy;
  345. SkFixed newx, newy, newSnappedX, newSnappedY;
  346. int shift = fCurveShift;
  347. SkASSERT(count > 0);
  348. do {
  349. SkFixed slope;
  350. if (--count > 0)
  351. {
  352. newx = oldx + (dx >> shift);
  353. newy = oldy + (dy >> shift);
  354. if (SkAbs32(dy >> shift) >= SK_Fixed1 * 2) { // only snap when dy is large enough
  355. SkFDot6 diffY = SkFixedToFDot6(newy - fSnappedY);
  356. slope = diffY ? quick_div(SkFixedToFDot6(newx - fSnappedX), diffY)
  357. : SK_MaxS32;
  358. newSnappedY = SkTMin<SkFixed>(fQEdge.fQLastY, SkFixedRoundToFixed(newy));
  359. newSnappedX = newx - SkFixedMul(slope, newy - newSnappedY);
  360. } else {
  361. newSnappedY = SkTMin(fQEdge.fQLastY, SnapY(newy));
  362. newSnappedX = newx;
  363. SkFDot6 diffY = SkFixedToFDot6(newSnappedY - fSnappedY);
  364. slope = diffY ? quick_div(SkFixedToFDot6(newx - fSnappedX), diffY)
  365. : SK_MaxS32;
  366. }
  367. dx += fQEdge.fQDDx;
  368. dy += fQEdge.fQDDy;
  369. }
  370. else // last segment
  371. {
  372. newx = fQEdge.fQLastX;
  373. newy = fQEdge.fQLastY;
  374. newSnappedY = newy;
  375. newSnappedX = newx;
  376. SkFDot6 diffY = (newy - fSnappedY) >> 10;
  377. slope = diffY ? quick_div((newx - fSnappedX) >> 10, diffY) : SK_MaxS32;
  378. }
  379. if (slope < SK_MaxS32) {
  380. success = this->updateLine(fSnappedX, fSnappedY, newSnappedX, newSnappedY, slope);
  381. }
  382. oldx = newx;
  383. oldy = newy;
  384. } while (count > 0 && !success);
  385. SkASSERT(newSnappedY <= fQEdge.fQLastY);
  386. fQEdge.fQx = newx;
  387. fQEdge.fQy = newy;
  388. fQEdge.fQDx = dx;
  389. fQEdge.fQDy = dy;
  390. fSnappedX = newSnappedX;
  391. fSnappedY = newSnappedY;
  392. fCurveCount = SkToS8(count);
  393. return success;
  394. }
  395. bool SkAnalyticCubicEdge::setCubic(const SkPoint pts[4], bool sortY) {
  396. fRiteE = nullptr;
  397. if (!fCEdge.setCubicWithoutUpdate(pts, kDefaultAccuracy, sortY)) {
  398. return false;
  399. }
  400. fCEdge.fCx >>= kDefaultAccuracy;
  401. fCEdge.fCy >>= kDefaultAccuracy;
  402. fCEdge.fCDx >>= kDefaultAccuracy;
  403. fCEdge.fCDy >>= kDefaultAccuracy;
  404. fCEdge.fCDDx >>= kDefaultAccuracy;
  405. fCEdge.fCDDy >>= kDefaultAccuracy;
  406. fCEdge.fCDDDx >>= kDefaultAccuracy;
  407. fCEdge.fCDDDy >>= kDefaultAccuracy;
  408. fCEdge.fCLastX >>= kDefaultAccuracy;
  409. fCEdge.fCLastY >>= kDefaultAccuracy;
  410. fCEdge.fCy = SnapY(fCEdge.fCy);
  411. fCEdge.fCLastY = SnapY(fCEdge.fCLastY);
  412. fWinding = fCEdge.fWinding;
  413. fCurveCount = fCEdge.fCurveCount;
  414. fCurveShift = fCEdge.fCurveShift;
  415. fCubicDShift = fCEdge.fCubicDShift;
  416. fSnappedY = fCEdge.fCy;
  417. return this->updateCubic(sortY);
  418. }
  419. bool SkAnalyticCubicEdge::updateCubic(bool sortY) {
  420. int success;
  421. int count = fCurveCount;
  422. SkFixed oldx = fCEdge.fCx;
  423. SkFixed oldy = fCEdge.fCy;
  424. SkFixed newx, newy;
  425. const int ddshift = fCurveShift;
  426. const int dshift = fCubicDShift;
  427. SkASSERT(count < 0);
  428. do {
  429. if (++count < 0) {
  430. newx = oldx + (fCEdge.fCDx >> dshift);
  431. fCEdge.fCDx += fCEdge.fCDDx >> ddshift;
  432. fCEdge.fCDDx += fCEdge.fCDDDx;
  433. newy = oldy + (fCEdge.fCDy >> dshift);
  434. fCEdge.fCDy += fCEdge.fCDDy >> ddshift;
  435. fCEdge.fCDDy += fCEdge.fCDDDy;
  436. }
  437. else { // last segment
  438. newx = fCEdge.fCLastX;
  439. newy = fCEdge.fCLastY;
  440. }
  441. // we want to say SkASSERT(oldy <= newy), but our finite fixedpoint
  442. // doesn't always achieve that, so we have to explicitly pin it here.
  443. if (sortY && newy < oldy) {
  444. newy = oldy;
  445. }
  446. SkFixed newSnappedY = SnapY(newy);
  447. // we want to SkASSERT(snappedNewY <= fCEdge.fCLastY), but our finite fixedpoint
  448. // doesn't always achieve that, so we have to explicitly pin it here.
  449. if (sortY && fCEdge.fCLastY < newSnappedY) {
  450. newSnappedY = fCEdge.fCLastY;
  451. count = 0;
  452. }
  453. SkFixed slope = SkFixedToFDot6(newSnappedY - fSnappedY) == 0
  454. ? SK_MaxS32
  455. : SkFDot6Div(SkFixedToFDot6(newx - oldx),
  456. SkFixedToFDot6(newSnappedY - fSnappedY));
  457. success = this->updateLine(oldx, fSnappedY, newx, newSnappedY, slope);
  458. oldx = newx;
  459. oldy = newy;
  460. fSnappedY = newSnappedY;
  461. } while (count < 0 && !success);
  462. fCEdge.fCx = newx;
  463. fCEdge.fCy = newy;
  464. fCurveCount = SkToS8(count);
  465. return success;
  466. }