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))); }