diff options
Diffstat (limited to 'chromium/third_party/skia/experimental/Intersection/CubicParameterization_TestUtility.cpp')
-rw-r--r-- | chromium/third_party/skia/experimental/Intersection/CubicParameterization_TestUtility.cpp | 53 |
1 files changed, 53 insertions, 0 deletions
diff --git a/chromium/third_party/skia/experimental/Intersection/CubicParameterization_TestUtility.cpp b/chromium/third_party/skia/experimental/Intersection/CubicParameterization_TestUtility.cpp new file mode 100644 index 00000000000..356ca88f857 --- /dev/null +++ b/chromium/third_party/skia/experimental/Intersection/CubicParameterization_TestUtility.cpp @@ -0,0 +1,53 @@ +// included by CubicParameterization.cpp +// accesses internal functions to validate parameterized coefficients + +#include "Parameterization_Test.h" + +static void parameter_coeffs(const Cubic& cubic, double coeffs[coeff_count]) { +#if USE_SYVESTER + double ax, bx, cx, dx; + if (try_alt) + alt_set_abcd(&cubic[0].x, ax, bx, cx, dx); + else + set_abcd(&cubic[0].x, ax, bx, cx, dx); + double ay, by, cy, dy; + if (try_alt) + alt_set_abcd(&cubic[0].y, ay, by, cy, dy); + else + set_abcd(&cubic[0].y, ay, by, cy, dy); + calc_ABCD(ax, ay, coeffs); + if (!try_alt) calc_bc(dx, bx, cx); + if (!try_alt) calc_bc(dy, by, cy); +#else + double ax = cubic[0].x; + double bx = cubic[1].x; + double cx = cubic[2].x; + double dx = cubic[3].x; + double ay = cubic[0].y; + double by = cubic[1].y; + double cy = cubic[2].y; + double dy = cubic[3].y; + calc_ABCD(ax, bx, cx, dx, ay, by, cy, dy, coeffs); +#endif + for (int index = xx_coeff; index < coeff_count; ++index) { + int procIndex = index - xx_coeff; + coeffs[index] = (*calc_proc[procIndex])(ax, bx, cx, dx, ay, by, cy, dy); + } +} + +bool point_on_parameterized_curve(const Cubic& cubic, const _Point& point) { + double coeffs[coeff_count]; + parameter_coeffs(cubic, coeffs); + double xxx = coeffs[xxx_coeff] * point.x * point.x * point.x; + double xxy = coeffs[xxy_coeff] * point.x * point.x * point.y; + double xyy = coeffs[xyy_coeff] * point.x * point.y * point.y; + double yyy = coeffs[yyy_coeff] * point.y * point.y * point.y; + double xx = coeffs[ xx_coeff] * point.x * point.x; + double xy = coeffs[ xy_coeff] * point.x * point.y; + double yy = coeffs[ yy_coeff] * point.y * point.y; + double x = coeffs[ x_coeff] * point.x; + double y = coeffs[ y_coeff] * point.y; + double c = coeffs[ c_coeff]; + double sum = xxx + xxy + xyy + yyy + xx + xy + yy + x + y + c; + return approximately_zero(sum); +} |