public static double[,] EvalMatrixD_OneThread(AsyncArg arg, IExpression expr, DependencySpace depSpace, string[] deps,
string var1, double len1, int steps1, string var2, double len2, int steps2)
{
expr = ExprUtils.GetCopy_Slow(expr);
foreach (string dep in deps)
{
Complex?val = depSpace.Get(dep);
if (val.HasValue)
{
expr = ExprSimplifier.Substitute(expr, dep, new ExprConst(val.Value.Real.ToString("f15")));
}
}
double[,] res = new double[steps1, steps2]; // t, x
for (int n = 0; n < steps1; n++)
{
if (arg.Token.IsCancellationRequested)
{
return(null);
}
double v1 = len1 * n / steps1;
IExpression substituted1 = ExprSimplifier.Substitute(expr, var1, new ExprConst(v1.ToString("f15")));
for (int j = 0; j < steps2; j++)
{
double v2 = len2 * j / steps2;
IExpression substituted2 = ExprSimplifier.Substitute(substituted1, var2, new ExprConst(v2.ToString("f15")));
substituted2 = ExprSimplifier.Simplify(substituted2);
res[n, j] = ExprDoubleSimplifier.CalcConstExpr(substituted2);
}
arg.Progress?.Report((n + 1) / (float)steps1);
}
return(res);
}
//Velocity is calculated at r'(t) = <f'(t),g'(t),h'(t)>
protected override void SetForces()
{
forces = new List <Vector3>();
float step = resolution / (tMax - tMin);
for (float i = tMin; i < tMax; i += step)
{
Vector3 newVelocity = ExprUtils.numDerivative3D(i, delegate3D);
forces.Add(newVelocity);
}
}
public static double[,] EvalMatrixD(AsyncArg arg, IExpression expr, DependencySpace depSpace, string[] deps,
string var1, double len1, int steps1, string var2, double len2, int steps2)
{
expr = ExprUtils.GetCopy_Slow(expr);
foreach (string dep in deps)
{
Complex?val = depSpace.Get(dep);
if (val.HasValue)
{
expr = ExprSimplifier.Substitute(expr, dep, new ExprConst(val.Value.Real.ToString("f15")));
}
}
double[,] res = new double[steps1, steps2]; // t, x
if (!ExprSimplifier.IsDependsOn(expr, var1))
{
double[] xArr = EvalArrayD(expr, depSpace, deps, var2, len2, steps2);
for (int j = 0; j < steps2; j++)
{
for (int i = 0; i < steps1; i++)
{
res[i, j] = xArr[j];
}
}
return(res);
}
if (!ExprSimplifier.IsDependsOn(expr, var2))
{
double[] xArr = EvalArrayD(expr, depSpace, deps, var1, len1, steps1);
for (int i = 0; i < steps1; i++)
{
for (int j = 0; j < steps2; j++)
{
res[i, j] = xArr[i];
}
}
return(res);
}
int[] progress = new int[Environment.ProcessorCount];
int interval = (int)Math.Ceiling(steps1 / (double)progress.Length);
Parallel.For(0, progress.Length, (index) => EvalMatrixD(arg, expr, res, steps1, progress, index,
var1, len1 / steps1, index * interval, Math.Min((index + 1) * interval, steps1), var2, len2 / steps2, steps2));
return(res);
}
public static double[] EvalArrayD(IExpression expr, DependencySpace depSpace, string[] deps, string varName, double length, int sections)
{
expr = ExprUtils.GetCopy_Slow(expr);
foreach (string dep in deps)
{
Complex?val = depSpace.Get(dep);
if (val.HasValue)
{
expr = ExprSimplifier.Substitute(expr, dep, new ExprConst(val.Value.Real.ToString("f15")));
}
}
double[] u0 = new double[sections];
for (int j = 0; j < u0.Length; j++)
{
double v = length * j / u0.Length;
IExpression substituted = ExprSimplifier.Substitute(expr, varName, new ExprConst(v.ToString("f15")));
substituted = ExprSimplifier.Simplify(substituted);
u0[j] = ExprDoubleSimplifier.CalcConstExpr(substituted);
}
return(u0);
}
public static Complex Eval(IExpression expr, List <DependencyNode> dependencies)
{
expr = ExprUtils.GetCopy_Slow(expr);
foreach (DependencyNode dep in dependencies)
{
Complex val = dep.Value;
IExpression constExpr;
if (val.Imaginary == 0)
{
constExpr = new ExprConst(val.Real.ToString("f15"));
}
else
{
IExpression exprImg = new ExprMult(new ExprConst(val.Real.ToString("f15")), new ExprRegFunc("i", 0));
if (val.Real == 0)
{
constExpr = exprImg;
}
else
{
constExpr = new ExprSum(new ExprConst(val.Real.ToString("f15")), exprImg);
}
}
expr = ExprSimplifier.Substitute(expr, dep.Name, constExpr);
}
expr = ExprSimplifier.Simplify(expr);
//deal with static dependencies(chi?)
//replace, considering functions
if (expr is ExprConst)
{
return((expr as ExprConst).Value.ToDouble());
}
else if (expr is ExprMult)
{
ExprMult mult = expr as ExprMult;
double img = 0;
bool found = false;
foreach (IExpression e in mult.GetArgs())
{
if (e is ExprConst)
{
img = (e as ExprConst).Value.ToDouble();
}
else if (e is ExprRegFunc && (e as ExprRegFunc).Name == "i")
{
found = true;
}
}
if (found)
{
return(new Complex(0, img));
}
}
else if (expr is ExprSum)
{
ExprSum sum = expr as ExprSum;
bool foundReal = false;
bool foundImg = false;
double real = 0;
double img = 0;
foreach (IExpression eSum in sum.GetArgs())
{
if (eSum is ExprConst)
{
real = (eSum as ExprConst).Value.ToDouble();
foundReal = true;
}
else if (eSum is ExprMult)
{
ExprMult mult = eSum as ExprMult;
foreach (IExpression e in mult.GetArgs())
{
if (e is ExprConst)
{
img = (e as ExprConst).Value.ToDouble();
}
else if (e is ExprRegFunc && (e as ExprRegFunc).Name == "i")
{
foundImg = true;
}
}
}
}
if (foundReal && foundImg)
{
return(new Complex(real, img));
}
}
return(Complex.NaN);
}
