public void Test_OptimisedTTMProfiler_TwoTrianglesWitGapAtOrigin_DiagonalProfileLine(double atElevation)
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithTwoFlatUnitTrianglesWithGapAtOrigin(atElevation);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (-100, -100) to (100, 100) to bisect the two triangles with a gap
// between (0.5, 0.5) and (1.0, 1.0)
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
var profilePoints = profiler.Compute(new[] { new XYZ(-100, -100), new XYZ(100, 100) });
Assert.True(profilePoints.Count == 5, $"Profile operation returned {profilePoints.Count} intercepts instead of 5");
Assert.True(Math.Abs(profilePoints[0].X) < epsilon &&
Math.Abs(profilePoints[0].Y) < epsilon &&
Math.Abs(profilePoints[0].Z - atElevation) < epsilon, $"First profile point not at (0, 0, {atElevation}), but is at {profilePoints[0]}");
Assert.True(Math.Abs(profilePoints[1].X - 0.5) < epsilon &&
Math.Abs(profilePoints[1].Y - 0.5) < epsilon &&
Math.Abs(profilePoints[1].Z - atElevation) < epsilon, $"Second profile point not at (0.5, 0.5, {atElevation}), but is at {profilePoints[1]}");
Assert.True(Math.Abs(profilePoints[2].X - 1.0) < epsilon &&
Math.Abs(profilePoints[2].Y - 1.0) < epsilon &&
profilePoints[2].Z == TRex.Common.Consts.NullDouble, $"Gap profile point not at (1.0, 1.0, {atElevation}), but is at {profilePoints[2]}");
Assert.True(Math.Abs(profilePoints[3].X - 1.0) < epsilon &&
Math.Abs(profilePoints[3].Y - 1.0) < epsilon &&
Math.Abs(profilePoints[3].Z - atElevation) < epsilon, $"Third profile point not at (1.0, 1.0, {atElevation}), but is at {profilePoints[3]}");
Assert.True(Math.Abs(profilePoints[4].X - 1.5) < epsilon &&
Math.Abs(profilePoints[4].Y - 1.5) < epsilon &&
Math.Abs(profilePoints[4].Z - atElevation) < epsilon, $"Fourth profile point not at (1.5, 1.5, {atElevation}), but is at {profilePoints[4]}");
}
public void Test_OptimisedTTMProfiler_TwoTrianglesAtOrigin(double atElevation)
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithTwoFlatUnitTrianglesAtOrigin(atElevation);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (-100, -100) to (100, 100) to bisect the triangles
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
var profilePoints = profiler.Compute(new[] { new XYZ(-100, -100), new XYZ(100, 100) });
Assert.True(profilePoints.Count == 3, $"Profile operation returned {profilePoints.Count} intercepts instead of 3");
Assert.True(Math.Abs(profilePoints[0].X) < epsilon &&
Math.Abs(profilePoints[0].Y) < epsilon &&
Math.Abs(profilePoints[0].Z - atElevation) < epsilon, $"First profile point not at (0, 0, {atElevation}), but is at {profilePoints[0]}");
Assert.True(profilePoints[0].TriIndex == 0, "Triangle index not set");
Assert.True(Math.Abs(profilePoints[1].X - 0.5) < epsilon &&
Math.Abs(profilePoints[1].Y - 0.5) < epsilon &&
Math.Abs(profilePoints[1].Z - atElevation) < epsilon, $"Second profile point not at (0.5, 0.5, {atElevation}), but is at {profilePoints[1]}");
Assert.True(profilePoints[1].TriIndex == 0, "Triangle index not set");
Assert.True(Math.Abs(profilePoints[2].X - 1.0) < epsilon &&
Math.Abs(profilePoints[2].Y - 1.0) < epsilon &&
Math.Abs(profilePoints[2].Z - atElevation) < epsilon, $"Second profile point not at (1.0, 1.0, {atElevation}), but is at {profilePoints[2]}");
Assert.True(profilePoints[2].TriIndex == 1, "Triangle index not set");
}
public void Test_OptimisedTTMDesignBuilder_OneTriangle()
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithFlatUnitTriangleAtOrigin(0.0);
Assert.True(oneTriangleModel.Vertices.Items.Length == 3, "Invalid number of vertices for single triangle model");
Assert.True(oneTriangleModel.Triangles.Items.Length == 1, "Invalid number of triangles for single triangle model");
OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var tree, out var indices);
Assert.NotNull(tree);
Assert.NotNull(indices);
Assert.True(indices.Length == 2, $"Number of indices [{indices.Length}] incorrect, should be 2");
}
public void Test_ProfilerBuilder_OneTriangle()
{
const int expectedInterceptCount = 8;
// Create a model with a single triangle at (0, 0), (0, 1), (1, 0)
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithFlatUnitTriangleAtOrigin(0.0);
OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var tree, out var indices);
var builder = new OptimisedTTMCellProfileBuilder(tree.CellSize, true);
// Build a profile line from (-100, -100) to (100, 100) to bisect the single triangle
var result = builder.Build(new [] { new XYZ(-100, -100), new XYZ(100, 100) }, 0);
Assert.True(result, "Build() failed");
Assert.True(builder.VtHzIntercepts.Count == expectedInterceptCount, $"Intercept count [{builder.VtHzIntercepts.Count}] wrong, expected {expectedInterceptCount}");
}
public void Test_OptimisedTTMProfiler_SingleTriangleAtOrigin_OutsideToInsideWithBisection(double atElevation)
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithFlatUnitTriangleAtOrigin(atElevation);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (-100, -100) to (0.25, 0.25) to enter the single triangle at a corner and terminate within it
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
var profilePoints = profiler.Compute(new[] { new XYZ(-100, -100), new XYZ(0.25, 0.25) });
Assert.True(profilePoints.Count == 2, $"Profile operation returned {profilePoints.Count} intercepts instead of 2");
Assert.True(Math.Abs(profilePoints[0].X) < epsilon &&
Math.Abs(profilePoints[0].Y) < epsilon &&
Math.Abs(profilePoints[0].Z - atElevation) < epsilon, $"First profile point not at origin (0, 0, {atElevation}), but is at {profilePoints[0]}");
Assert.True(Math.Abs(profilePoints[1].X - 0.25) < epsilon &&
Math.Abs(profilePoints[1].Y - 0.25) < epsilon &&
Math.Abs(profilePoints[1].Z - atElevation) < epsilon, $"Second profile point not at origin (0.25, 0.25, {atElevation}), but is at {profilePoints[1]}");
}
public void Test_OptimisedTTMProfiler_TwoTrianglesAtOrigin_DiagonalColinearProfileLine(double atElevation)
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithFlatUnitTriangleAtOrigin(atElevation);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (-1, 2) to (2, -1) to be co-linear with diagonal edge of triangle
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
var profilePoints = profiler.Compute(new[] { new XYZ(-1, 2), new XYZ(2, -1) });
Assert.True(profilePoints.Count == 2, $"Profile operation returned {profilePoints.Count} intercepts instead of 2");
Assert.True(Math.Abs(profilePoints[0].X) < epsilon &&
Math.Abs(profilePoints[0].Y - 1.0) < epsilon &&
Math.Abs(profilePoints[0].Z - atElevation) < epsilon, $"First profile point not at (0, 1.0, {atElevation}), but is at {profilePoints[0]}");
Assert.True(Math.Abs(profilePoints[1].X - 1.0) < epsilon &&
Math.Abs(profilePoints[1].Y) < epsilon &&
Math.Abs(profilePoints[1].Z - atElevation) < epsilon, $"Second profile point not at (1.0, 0, {atElevation}), but is at {profilePoints[1]}");
}
public void Test_OptimisedTTMProfiler_TwoTrianglesAtOrigin_ManyTimes(int runCount)
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithTwoFlatUnitTrianglesAtOrigin(0.0);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (-100, -100) to (100, 100) to bisect the single triangle
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
DateTime startTime = DateTime.Now;
for (int i = 0; i < runCount; i++)
{
var profilePoints = profiler.Compute(new[] { new XYZ(-100, -100), new XYZ(100, 100) });
Assert.True(profilePoints.Count > 0);
}
output.WriteLine($"Times to run profile {runCount} times: {DateTime.Now - startTime}");
Assert.True(true);
}
public void Test_OptimisedTTMProfiler_1024TrianglesAtOrigin(double atElevation, double startX, double startY, double endX, double endY, int expectedPointCount)
{
var startTime = DateTime.Now;
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_With32x32FlatTrianglesAtOrigin(atElevation);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (startX, startX) to (endX, endY) to cross the patch of triangles
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
var profilePoints = profiler.Compute(new[] { new XYZ(startX, startY), new XYZ(endX, endY) });
output.WriteLine($"Total time to perform profile: {DateTime.Now - startTime}");
Assert.True(profilePoints.Count == expectedPointCount, $"Profile operation returned {profilePoints.Count} intercepts instead of {expectedPointCount}");
foreach (var pt in profilePoints)
{
Assert.True(Math.Abs(pt.Z - atElevation) < epsilon, $"Elevation {pt.Z} incorrect, should be {atElevation}");
Assert.True(Math.Abs(pt.Station - MathUtilities.Hypot(startX - pt.X, startY - pt.Y)) < epsilon, $"Station {pt.Station} incorrect, should be {MathUtilities.Hypot(startX - pt.X, startY - pt.Y)}");
}
}
public void Test_OptimisedTTMProfiler_SingleTriangleAtOrigin_MultipleProfileSegments(double atElevation)
{
var oneTriangleModel = OptimisedTTMDesignBuilder.CreateOptimisedTTM_WithFlatUnitTriangleAtOrigin(atElevation);
Assert.True(OptimisedTTMDesignBuilder.CreateOptimisedIndexForModel(oneTriangleModel, out var index, out var indices));
// Build a profile line from (-100, -100) to (0.25, 0.25) to (100, 100) to enter the single triangle at a corner
// stop within it then exit via the diagonal face
var profiler = new OptimisedTTMProfiler(oneTriangleModel, index, indices);
var profilePoints = profiler.Compute(new[] { new XYZ(-100, -100), new XYZ(0.25, 0.25), new XYZ(100, 100) });
Assert.True(profilePoints.Count == 3, $"Profile operation returned {profilePoints.Count} intercepts instead of 3");
Assert.True(Math.Abs(profilePoints[0].X) < epsilon &&
Math.Abs(profilePoints[0].Y) < epsilon &&
Math.Abs(profilePoints[0].Z - atElevation) < epsilon, $"First profile point not at origin (0, 0, {atElevation}), but is at {profilePoints[0]}");
Assert.True(Math.Abs(profilePoints[1].X - 0.25) < epsilon &&
Math.Abs(profilePoints[1].Y - 0.25) < epsilon &&
Math.Abs(profilePoints[1].Z - atElevation) < epsilon, $"Second profile point not at origin (0.25, 0.25, {atElevation}), but is at {profilePoints[1]}");
Assert.True(Math.Abs(profilePoints[2].X - 0.5) < epsilon &&
Math.Abs(profilePoints[2].Y - 0.5) < epsilon &&
Math.Abs(profilePoints[2].Z - atElevation) < epsilon, $"Second profile point not at origin (0.5, 0.5, {atElevation}), but is at {profilePoints[2]}");
}