public void Robust_FromFile_Pscico_ShouldSucceed()
{
double[] points =
File.ReadAllLines("test_data/pscico_data/insphere3d.txt")
.Select(n => n.Split()).SelectMany(x => x)
.Select(s => double.Parse(s)).ToArray();
for (int i = 0; i < NPscicoData; i += 17)
{
Assert.Equal(Math.Sign(points[i + 16]), Math.Sign(InSphere.Robust(
new double[] { points[i + 1], points[i + 2], points[i + 3] },
new double[] { points[i + 4], points[i + 5], points[i + 6] },
new double[] { points[i + 7], points[i + 8], points[i + 9] },
new double[] { points[i + 10], points[i + 11], points[i + 12] },
new double[] { points[i + 13], points[i + 14], points[i + 15] })));
}
}
public void Fast_ShouldSucceed()
{
Assert.True(InSphere.Robust(
new double[] { 0, 0, 0 },
new double[] { 1, 0, 0 },
new double[] { 0, 1, 1 },
new double[] { 0, 0, 1 },
new double[] { 0, 0, 1 }) == 0);
Assert.True(InSphere.Robust(
new double[] { 0, 0, 0 },
new double[] { 0, 1, 0 },
new double[] { 1, 0, 0 },
new double[] { 0, 0, 1 },
new double[] { 0, 0, 10 }) < 0);
Assert.True(InSphere.Robust(
new double[] { 0, 0, 0 },
new double[] { 1, 1, 0 },
new double[] { 0, 1, 1 },
new double[] { 0, 1, 0 },
new double[] { 0, 0, 0.1 }) > 0);
}
public void Fast_FromFile_ShouldSucceed()
{
double[] points =
File.ReadAllLines("test_data/simple_data/insphere-points3d.txt")
.Select(n => n.Split()).SelectMany(x => x)
.Select(s => double.Parse(s)).ToArray();
double[] results =
File.ReadAllLines("test_data/simple_data/results-insphere.txt")
.Select(n => n.Split()).SelectMany(x => x)
.Select(s => double.Parse(s)).ToArray();
int count = 0;
for (int i = 0; i < NSimpleData; i += 15)
{
var orient3DResult = InSphere.Fast(
new double[] { points[i], points[i + 1], points[i + 2] },
new double[] { points[i + 3], points[i + 4], points[i + 5] },
new double[] { points[i + 6], points[i + 7], points[i + 8] },
new double[] { points[i + 9], points[i + 10], points[i + 11] },
new double[] { points[i + 12], points[i + 13], points[i + 14] });
Assert.Equal(Math.Sign(results[count++]), Math.Sign(orient3DResult));
}
}
