public float sampleTrilinear(Array3D <float> noise, Vector3 pos)
{
var min = pos.ToFloored();
var f = pos - min;
/*var q000 = noise.GetTiled(new Point3(0, 0, 0) + min);
* var q100 = noise.GetTiled(new Point3(1, 0, 0) + min);
* var q010 = noise.GetTiled(new Point3(0, 1, 0) + min);
* var q001 = noise.GetTiled(new Point3(0, 0, 1) + min);
* var q110 = noise.GetTiled(new Point3(1, 1, 0) + min);
* var q011 = noise.GetTiled(new Point3(0, 1, 1) + min);
* var q101 = noise.GetTiled(new Point3(1, 0, 1) + min);
* var q111 = noise.GetTiled(new Point3(1, 1, 1) + min);*/
var x0 = TWMath.nfmod(min.X + 0, noise.Size.X);;
var y0 = TWMath.nfmod(min.Y + 0, noise.Size.Y);
var z0 = TWMath.nfmod(min.Z + 0, noise.Size.Z);
var x1 = TWMath.nfmod(x0 + 1, noise.Size.X);
var y1 = TWMath.nfmod(y0 + 1, noise.Size.Y);
var z1 = TWMath.nfmod(z0 + 1, noise.Size.Z);
var q000 = noise.GetFast(x0, y0, z0);
var q100 = noise.GetFast(x1, y0, z0);
var q010 = noise.GetFast(x0, y1, z0);
var q001 = noise.GetFast(x0, y0, z1);
var q110 = noise.GetFast(x1, y1, z0);
var q011 = noise.GetFast(x0, y1, z1);
var q101 = noise.GetFast(x1, y0, z1);
var q111 = noise.GetFast(x1, y1, z1);
var ret = TWMath.triLerp(f.dx(), q000, q100, q001, q101, q010, q110, q011, q111);
return(ret);
}
public T GetTiled(Point2 pos)
{
// Dont call the this[] directly, its about 20% slower when i tested it. Probably because of the out of bounds check
// TODO: since this is always in bounds, maybe try using 'unsafe' and pointers to avoid the C# array bounds checking?
Point2 point3 = Size; // Reading this out is a bit faster
var x = TWMath.nfmod(pos.X, point3.X);
var y = TWMath.nfmod(pos.Y, point3.Y);
return(arr[y * Size.X + x]);
}
public T GetTiledFast(int x, int y, int z)
{
// Dont call the this[] directly, its about 20% slower when i tested it. Probably because of the out of bounds check
// TODO: since this is always in bounds, maybe try using 'unsafe' and pointers to avoid the C# array bounds checking?
Point3 point3 = Size; // Reading this out is a bit faster
x = TWMath.nfmod(x, point3.X);
y = TWMath.nfmod(y, point3.Y);
z = TWMath.nfmod(z, point3.Z);
return(arr[x, y, z]);
}
public void TestDensityCalcPerformance()
{
var dens = VoxelTerrainGenerationTest.createDensityFunction5Perlin(11, 10);
var times = 10000;
var s = new Stopwatch();
s.Start();
for (int i = 0; i < times; i++)
{
float a = dens(new Vector3(1));
}
s.Stop();
var densityEval = s.Elapsed.TotalSeconds / times;
Console.WriteLine("Density evaluation: 100% - {0} ms", +densityEval * 1000);
s.Reset();
s.Start();
for (int i = 0; i < times * 5; i++)// 5 trilerps per density lookup
{
float a = TWMath.triLerp(new Vector3(0.5f), 4, 7, 8, 6, 4, 8, 9, 7);
}
s.Stop();
var trilerp = s.Elapsed.TotalSeconds / times;
Console.WriteLine("Trilerp : {0:###}% - {1}ms", trilerp / densityEval * 100, trilerp * 1000);
var ar = new Array3D <float>(new Point3(16, 16, 16));
s.Reset();
s.Start();
for (int i = 0; i < times * 5 * 8; i++)// 8 getTiled per trilerp and 5 trilerps per density lookup
{
/*var pos = new Point3(5, 5, 5);
* pos.X = TWMath.nfmod(pos.X, ar.Size.X);
* pos.Y = TWMath.nfmod(pos.Y, ar.Size.Y);
* pos.Z = TWMath.nfmod(pos.Z, ar.Size.Z);
* float a = ar.arr[pos.X, pos.Y, pos.Z];*/
float a = ar.GetTiled(new Point3(5, 5, 5));
}
s.Stop();
var getTiled = s.Elapsed.TotalSeconds / times;
Console.WriteLine("GetTiled : {0:###}% - {1}ms", getTiled / densityEval * 100, getTiled * 1000);
}
public Rect GetAABB()
{
var ll = LowerLeft;
var lr = LowerRight;
var ul = UpperLeft;
var ur = UpperRight;
var minx = TWMath.Min(ll.x, lr.x, ul.x, ur.x);
var miny = TWMath.Min(ll.y, lr.y, ul.y, ur.y);
var maxx = TWMath.Max(ll.x, lr.x, ul.x, ur.x);
var maxy = TWMath.Max(ll.y, lr.y, ul.y, ur.y);
return(new Rect(minx, miny, maxx - minx, maxy - miny));
}
private static bool SeparatesOnAxis(Vector2 axis, float aMin, float aMax, Vector2 b1, Vector2 b2, Vector2 b3, Vector2 b4)
{
var dotb1 = Vector2.Dot(axis, b1);
var dotb2 = Vector2.Dot(axis, b2);
var dotb3 = Vector2.Dot(axis, b3);
var dotb4 = Vector2.Dot(axis, b4);
var bMin = TWMath.Min(dotb1, dotb2, dotb3, dotb4);
var bMax = TWMath.Max(dotb1, dotb2, dotb3, dotb4);
if (aMin > bMax || bMin > aMax)
{
return(true);
}
return(false);
}
public void TestAccessPerformance()
{
var size = 100;
var times = 10000;
var arr = new float[size, size, size];
var multiDim = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = arr[i, i + 1, i + 2];
var val2 = arr[i + 1, i + 2, i];
}
}
});
Console.WriteLine("Multidim: " + multiDim.PrettyPrint());
var arr2 = new float[size * size * size];
var singleDim = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = arr2[(i + (i + 1) * size) * size + i + 2];
var val2 = arr2[(i + (i + 2) * size) * size + i + 1];
}
}
});
Console.WriteLine("SingleDim: " + singleDim.PrettyPrint());
Assert.Less(singleDim, multiDim);
var arr3d = new Array3D <float>(new Point3(size, size, size));
var arr3DSpeed = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = arr3d[new Point3(i, i + 1, i + 2)];
var val2 = arr3d[new Point3(-i, -i + 1, -i + 2)];
}
}
});
Console.WriteLine("arr3d[]: " + arr3DSpeed.PrettyPrint());
var arr3DFastSpeed = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = arr3d.GetFast(i, i + 1, i + 2);
var val2 = arr3d.GetFast(i, i + 2, i + 1);
}
}
});
Console.WriteLine("arr3d.GetFast:" + arr3DFastSpeed.PrettyPrint());
var arr3DTiled = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = arr3d.GetTiled(new Point3(i, i + 1, i + 2));
var val2 = arr3d.GetTiled(new Point3(i, i + 2, i + 1));
}
}
});
Console.WriteLine("arr3d.GetTiled:" + arr3DTiled.PrettyPrint());
var arr3DTiledFast = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = arr3d.GetTiledFast(i, i + 1, i + 2);
var val2 = arr3d.GetTiledFast(i, i + 2, i + 1);
}
}
});
Console.WriteLine("arr3d.GetTiledFast:" + arr3DTiledFast.PrettyPrint());
var arr3DTiledCacheSize = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var a1 = TWMath.nfmod(i, size);
var a2 = TWMath.nfmod(i + 1, size);
var a3 = TWMath.nfmod(i + 2, size);
var a4 = TWMath.nfmod(i, size);
var a5 = TWMath.nfmod(i + 2, size);
var a6 = TWMath.nfmod(i + 1, size);
var val = arr3d.GetFast(a1, a2, a3);
var val2 = arr3d.GetFast(a4, a5, a6);
}
}
});
Console.WriteLine("arr3d.GetTiledCacheSize:" + arr3DTiledCacheSize.PrettyPrint());
var internalArr = arr3d.GetInternalField <float[, , ]>("arr");
var methodGet = PerformanceHelper.Measure(() =>
{
for (int j = 0; j < times; j++)
{
for (int i = 0; i < size - 2; i++)
{
var val = Val(internalArr, i, i + 1, i + 2);
var val2 = Val(internalArr, i, i + 2, i + 1);
}
}
});
Console.WriteLine("DirectMethod: " + methodGet.PrettyPrint());
Assert.Less(arr3DFastSpeed, arr3DSpeed);
}
