public void GetLength()
{
BezierSegment2F b = new BezierSegment2F
{
Point1 = new Vector2F(1, 2),
ControlPoint1 = new Vector2F(4, 5),
ControlPoint2 = new Vector2F(7, 8),
Point2 = new Vector2F(10, 2),
};
float lowerBound = (b.Point2 - b.Point1).Length;
float upperBound = (b.Point2 - b.ControlPoint2).Length
+ (b.ControlPoint2 - b.ControlPoint1).Length
+ (b.ControlPoint1 - b.Point1).Length;
Assert.Less(lowerBound, b.GetLength(0, 1, 100, Numeric.EpsilonF));
Assert.Greater(upperBound, b.GetLength(0, 1, 100, Numeric.EpsilonF));
float length1 = b.GetLength(0, 1, 20, Numeric.EpsilonF);
float length2 = b.GetLengthWithDeCasteljau(20, Numeric.EpsilonF);
Assert.IsTrue(Numeric.AreEqual(length1, length2));
// Compare numerical integration method and de Casteljau method.
float approxLength = 0;
const float step = 0.0001f;
for (float u = 0; u <= 1.0f; u += step)
approxLength += (b.GetPoint(u) - b.GetPoint(u + step)).Length;
Assert.IsTrue(Numeric.AreEqual(approxLength, length1, 0.01f));
Assert.IsTrue(Numeric.AreEqual(b.GetLength(0, 1, 100, Numeric.EpsilonF), b.GetLength(0, 0.5f, 100, Numeric.EpsilonF) + b.GetLength(0.5f, 1, 100, Numeric.EpsilonF)));
Assert.IsTrue(Numeric.AreEqual(b.GetLength(0, 1, 100, Numeric.EpsilonF), b.GetLength(1, 0, 100, Numeric.EpsilonF)));
}
public void GetLength()
{
BezierSegment2F b = new BezierSegment2F
{
Point1 = new Vector2F(1, 2),
ControlPoint1 = new Vector2F(4, 5),
ControlPoint2 = new Vector2F(7, 8),
Point2 = new Vector2F(10, 2),
};
float lowerBound = (b.Point2 - b.Point1).Length;
float upperBound = (b.Point2 - b.ControlPoint2).Length
+ (b.ControlPoint2 - b.ControlPoint1).Length
+ (b.ControlPoint1 - b.Point1).Length;
Assert.Less(lowerBound, b.GetLength(0, 1, 100, Numeric.EpsilonF));
Assert.Greater(upperBound, b.GetLength(0, 1, 100, Numeric.EpsilonF));
float length1 = b.GetLength(0, 1, 20, Numeric.EpsilonF);
float length2 = b.GetLengthWithDeCasteljau(20, Numeric.EpsilonF);
Assert.IsTrue(Numeric.AreEqual(length1, length2));
// Compare numerical integration method and de Casteljau method.
float approxLength = 0;
const float step = 0.0001f;
for (float u = 0; u <= 1.0f; u += step)
{
approxLength += (b.GetPoint(u) - b.GetPoint(u + step)).Length;
}
Assert.IsTrue(Numeric.AreEqual(approxLength, length1, 0.01f));
Assert.IsTrue(Numeric.AreEqual(b.GetLength(0, 1, 100, Numeric.EpsilonF), b.GetLength(0, 0.5f, 100, Numeric.EpsilonF) + b.GetLength(0.5f, 1, 100, Numeric.EpsilonF)));
Assert.IsTrue(Numeric.AreEqual(b.GetLength(0, 1, 100, Numeric.EpsilonF), b.GetLength(1, 0, 100, Numeric.EpsilonF)));
}
