protected void DrawLine(Point3D point)
{
// point was in cube [0..1]x[0..1]x[0..1], we are transforming it
// to cube [0..width]x[0..height]x[0..depth]
Point3D start = point.TransformToBounds(bounds);
int maxLength = 10 * Math.Max(width, Math.Max(height, depth));
const int maxIterations = 400;
Size3D boundsSize = new Size3D(1.0 / width, 1.0 / height, 1.0 / depth);
Action<double, List<Point3D>> pointTracking = (direction, track) =>
{
var position01 = point;
double length = 0;
int i = 0;
do
{
var K1 = fieldWrapper.GetVector(position01).DecreaseLength(boundsSize);
var K2 = fieldWrapper.GetVector(position01 + (K1 / 2)).DecreaseLength(boundsSize);
var K3 = fieldWrapper.GetVector(position01 + (K2 / 2)).DecreaseLength(boundsSize);
var K4 = fieldWrapper.GetVector(position01 + K3).DecreaseLength(boundsSize);
var shift = ((K1 + 2 * K2 + 2 * K3 + K4) / 6);
if (shift.X.IsNaN() || shift.Y.IsNaN() || shift.Z.IsNaN())
break;
var next = position01 + direction * shift;
Point3D viewportPoint = position01.TransformToBounds(bounds);
track.Add(viewportPoint);
position01 = next;
length += shift.Length;
i++;
} while (length < maxLength && i < maxIterations);
};
var forwardTrack = new List<Point3D>();
forwardTrack.Add(start);
pointTracking(+1, forwardTrack);
if (forwardTrack.Count > 1)
{
var forwardLine = CreatePolyline(forwardTrack, start);
Children.Add(forwardLine);
}
var backwardTrack = new List<Point3D>();
backwardTrack.Add(start);
pointTracking(-1, backwardTrack);
if (backwardTrack.Count > 1)
{
var backwardLine = CreatePolyline(backwardTrack, start);
Children.Add(backwardLine);
}
}
