public double Calculate(InputDataSet_DEPRECATED ds, double dt)
{
for (int i = 0; i < points.Count; i++)
{
if (points[i] != null)
{
System.Numerics.Vector3 newPos = new System.Numerics.Vector3();
newPos.X = points[i].Pos.X + points[i].Vel.Y * (float)dt;
newPos.Y = points[i].Pos.Y + points[i].Vel.X * (float)dt;
newPos.Z = points[i].Pos.Z + points[i].Vel.Z * (float)dt;
try
{
points[i] = Interpolator_DEPRECATED.IWDInterpolatePoint(newPos, ds.GetPoints(), ds.voxelSize);
}
catch (Exception)
{
points[i] = null;
}
}
if (points[i] == null)
{
//Console.WriteLine("[WARNING] point " + i + " in null");
}
}
double A, B, C, D;
if (points[1] == null || points[2] == null)
{
B = 0;
D = 0;
}
else
{
B = (points[1].Pos.X - points[2].Pos.X) / 2 * radius;
D = (points[1].Pos.Y - points[2].Pos.Y) / 2 * radius;
}
if (points[3] == null || points[4] == null)
{
A = 0;
C = 0;
}
else
{
A = (points[4].Pos.X - points[3].Pos.X) / 2 * radius;
C = (points[4].Pos.Y - points[3].Pos.Y) / 2 * radius;
}
double[,] phi = new double[2, 2] {
{ A, B }, { C, D }
};
double[,] phiT = Matrix.Transpose(phi);
double[,] phiTphi = Matrix.Dot(phiT, phi);
EigenvalueDecomposition eig = new EigenvalueDecomposition(phiTphi, true, true);
double[] eigh = eig.RealEigenvalues;
return(Matrix.Max(eigh));
}
public Seed_DEPRECATED(Point center, float radius, InputDataSet_DEPRECATED ds) : base(center)
{
this.radius = radius;
FTLE = double.NegativeInfinity;
points.Add(center.Clone());
points.Add(ds.GetPoint(new System.Numerics.Vector3(center.Pos.X, center.Pos.Y + radius, center.Pos.Z))); //UP
points.Add(ds.GetPoint(new System.Numerics.Vector3(center.Pos.X, center.Pos.Y - radius, center.Pos.Z))); //DOWN
points.Add(ds.GetPoint(new System.Numerics.Vector3(center.Pos.X - radius, center.Pos.Y, center.Pos.Z))); //LEFT
points.Add(ds.GetPoint(new System.Numerics.Vector3(center.Pos.X + radius, center.Pos.Y, center.Pos.Z))); //RIGHT
}
public Advection_DEPRECATED(InputDataSet_DEPRECATED DS, FTLEDataSet_DEPRECATED FDS)
{
this.DS = DS;
this.FDS = FDS;
}