| Index: experimental/Intersection/CubicParameterization_TestUtility.cpp
|
| diff --git a/experimental/Intersection/CubicParameterization_TestUtility.cpp b/experimental/Intersection/CubicParameterization_TestUtility.cpp
|
| deleted file mode 100644
|
| index 356ca88f857cf6887b4808b8dcaa2bb19cf020e0..0000000000000000000000000000000000000000
|
| --- a/experimental/Intersection/CubicParameterization_TestUtility.cpp
|
| +++ /dev/null
|
| @@ -1,53 +0,0 @@
|
| -// 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);
|
| -}
|
|
|