public SetGlobalVariable ( string name, double value ) : Variable | ||
name | string | |
value | double | |
return | Variable |
// Use this for initialization void Awake() { xmin = -1f * minmaxVal; xmax = minmaxVal; ymin = -1f * minmaxVal; ymax = minmaxVal; zmin = -1f * minmaxVal; zmax = minmaxVal; vectors = new List <Transform>(); startPts = new Vector3[2097152]; offsets = new Vector3[2097152]; solver = new AK.ExpressionSolver(); expX = new AK.Expression(); expY = new AK.Expression(); expZ = new AK.Expression(); solver.SetGlobalVariable("x", 0); solver.SetGlobalVariable("y", 0); solver.SetGlobalVariable("z", 0); max_magnitude = 0f; textureMap = new Texture3D(128, 128, 128, TextureFormat.RGBA32, true); textureMap.wrapMode = TextureWrapMode.Clamp; numComplete = 0; drawn = true; drawQueue = 0; if (loadingMsg != null) { loadingMsgText = loadingMsg.GetComponent <TextMesh>(); } sleepInterval = 126; //CalculateVectors(); //DrawDensityPlot(); }
public Vector3 findXYZ(Vector3 uvw) { Vector3 output = new Vector3(); ExpressionSet es = calcManager.expressionSet; float scale = calcManager.paramSurface.currentScale; foreach (string key in es.expressions.Keys) { AK.ExpressionSolver solver = es.solver; if (es.ranges.ContainsKey("u")) { solver.SetGlobalVariable("u", uvw.x); } if (es.ranges.ContainsKey("v")) { solver.SetGlobalVariable("v", uvw.y); } if (es.ranges.ContainsKey("w")) { solver.SetGlobalVariable("w", uvw.z); } } if (es.IsCompiled()) { output.x = (float)es.expressions["X"].AKExpression.Evaluate(); output.y = (float)es.expressions["Y"].AKExpression.Evaluate(); output.z = (float)es.expressions["Z"].AKExpression.Evaluate(); } return(output * scale); }
// Evaluate all equations for a given x (called by a thread), and popluate arrays with the results void ThreadedEvaluate(float x_temp, int inOrderCount) { AK.ExpressionSolver subSolver = new AK.ExpressionSolver(); AK.Expression subExp = new AK.Expression(); subSolver.SetGlobalVariable("x", 0); subSolver.SetGlobalVariable("y", 0); subSolver.SetGlobalVariable("z", 0); AK.Variable subX = subSolver.GetGlobalVariable("x"); AK.Variable subY = subSolver.GetGlobalVariable("y"); AK.Variable subZ = subSolver.GetGlobalVariable("z"); subExp = subSolver.SymbolicateExpression(es.expressions["X"].expression); for (float y_temp = ymin; y_temp < ymax; y_temp += delta) { for (float z_temp = zmin; z_temp < zmax; z_temp += delta) { if ((int)((z_temp - zmin) / delta) % sleepInterval == 0) { Thread.Sleep(1); } subX.value = x_temp; subY.value = z_temp; subZ.value = y_temp; float x = (float)subExp.Evaluate(); //Mathf.Clamp(x, -Mathf.Exp(20), Mathf.Exp(20)); //float y = (float)expY.Evaluate(); //Mathf.Clamp(y, -Mathf.Exp(20), Mathf.Exp(20)); //float z = (float)expZ.Evaluate(); //Mathf.Clamp(z, -Mathf.Exp(20), Mathf.Exp(20)); Vector3 target = new Vector3(x_temp, y_temp, z_temp); Vector3 result = new Vector3(x, 0, 0); if (float.IsNaN(x) // || float.IsNaN(y) // || float.IsNaN(z) || Mathf.Abs(x) == 0) { result = new Vector3(0, 0, 0); } //Vector3 direction = result.normalized; lock (lck) { max_magnitude = (Mathf.Abs(x) > max_magnitude) ? Mathf.Abs(x) : max_magnitude; } startPts[inOrderCount] = target; offsets[inOrderCount] = result; inOrderCount++; } } // numComplete++; }
void MakeContiguousProceduralGrid() { vertices = new Vector3[(sizeX + 1) * cellsPerUnit * (sizeY + 1) * cellsPerUnit]; triangles = new int[sizeX * cellsPerUnit * sizeY * cellsPerUnit * 6]; // tracker ints int v = 0, t = 0; // v == y of first nested for-loop, but this is clearer! // some calculations float cellSize = 1 / (float)cellsPerUnit; Vector3 vertexOffset = new Vector3((float)sizeX * 0.5f, 0, (float)sizeX * 0.5f); // turns out multiplication is cheaper than division int numCellsX = sizeX * cellsPerUnit; int numCellsY = sizeY * cellsPerUnit; // <= because we need the edge vertices at the end of the grid (so gridsize + 1) for (float x = 0; x <= sizeX; x += cellSize) { for (float y = 0; y <= sizeY; y += cellSize) { solver.SetGlobalVariable("x", x); solver.SetGlobalVariable("y", y); float z = (float)zFn.Evaluate(); if (float.IsPositiveInfinity(z) || float.IsNegativeInfinity(z)) { vertices[v] = vertices[v - 1]; } else { vertices[v] = new Vector3(x, z, y) - vertexOffset; // yes, x,z,y } v++; } } v = 0; // #reset for more looping! for (int x = 0; x < numCellsX; x++) { for (int y = 0; y < numCellsY; y++) { triangles[t] = v; triangles[t + 1] = v + 1; triangles[t + 2] = v + (numCellsX + 1); triangles[t + 3] = v + (numCellsX + 1); triangles[t + 4] = v + 1; triangles[t + 5] = (v + 1) + (numCellsX + 1); v++; t += 6; } v++; // idrg why yet, but something to do with the +1 to sizeX, sizeY } }
// Use this for initialization void Start() { vectors = new List <Transform>(); startPts = new List <Vector3>(); offsets = new List <Vector3>(); solver = new AK.ExpressionSolver(); expX = new AK.Expression(); expY = new AK.Expression(); expZ = new AK.Expression(); solver.SetGlobalVariable("x", 0); solver.SetGlobalVariable("y", 0); solver.SetGlobalVariable("z", 0); varX = solver.GetGlobalVariable("x"); varY = solver.GetGlobalVariable("y"); varZ = solver.GetGlobalVariable("z"); max_magnitude = 0f; //CalculateVectors(); //DrawVectorField(); }
public static void TestGlobalConstants() { ExpressionSolver solver = new ExpressionSolver(); solver.SetGlobalVariable("test", 1); var exp1 = solver.SymbolicateExpression("test+1"); AssertSameValue(2.0, exp1.Evaluate()); solver.SetGlobalVariable("test", 2); var exp2 = solver.SymbolicateExpression("test+1"); AssertSameValue(3.0, exp2.Evaluate()); AssertSameValue(exp1.Evaluate(), exp2.Evaluate()); }
void ThreadedEvaluate(int TID, int chunkSize, int extra) { AK.ExpressionSolver solver_ = new AK.ExpressionSolver(); AK.Variable u_ = solver_.SetGlobalVariable("theta", 0); AK.Variable v_ = solver_.SetGlobalVariable("phi", 0); AK.Expression exp_ = solver_.SymbolicateExpression(expr); float umin = (float)solver_.EvaluateExpression(uminExpr); float umax = (float)solver_.EvaluateExpression(umaxExpr); float vmin = (float)solver_.EvaluateExpression(vminExpr); float vmax = (float)solver_.EvaluateExpression(vmaxExpr); System.Random rand = new System.Random(); for (int i = 0; i < chunkSize + extra; i++) { float u = (float)(chunkSize * TID + i) / (resolution - 1); u_.value = umin + (umax - umin) * u; for (int j = 0; j < resolution; j++) { float v = (float)j / (resolution - 1); v_.value = vmin + (vmax - vmin) * v; float result = (float)exp_.Evaluate(); Vector3 spherical = new Vector3((float)u_.value, (float)v_.value, result); Vector3 cartesian = SphericalToCartesian(spherical); Vector3 vel = new Vector3( (float)rand.NextDouble(), (float)rand.NextDouble(), (float)rand.NextDouble()); vel = vel.normalized * 0.1f; float sqrt2 = Mathf.Sqrt(2.0f); Color c = Color.cyan * Mathf.Sqrt(u * u + v * v) / sqrt2 + Color.yellow * (1.0f - Mathf.Sqrt(u * u + v * v) / sqrt2) + Color.red * Mathf.Sqrt(u * u + (1 - v) * (1 - v)) / sqrt2 + Color.green * (1.0f - Mathf.Sqrt(u * u + (1 - v) * (1 - v)) / sqrt2); Vector3 pos = new Vector3(cartesian.x, cartesian.z, cartesian.y); Particle p = new Particle() { position = pos, velocity = vel, color = c }; lock (insert_lck) { results.Add(p); } } } }
public void SetupSolver() { depth = parameter_names.Count; width = (int)Mathf.Pow(PARTICLE_COUNT, 1f / (float)depth); for (int i = 0; i < depth; i++) { string name = parameter_names[i]; indexedParam[i] = name; parameterMin[name] = (float)solver.EvaluateExpression(parameter_min[name]); parameterMax[name] = (float)solver.EvaluateExpression(parameter_max[name]); solver.SetGlobalVariable(name, 0); } expr1 = solver.SymbolicateExpression(expression1); expr2 = solver.SymbolicateExpression(expression2); expr3 = solver.SymbolicateExpression(expression3); for (int i = 0; i < depth; i++) { string name = parameter_names[i]; AK.Variable var = solver.GetGlobalVariable(name); //vars.Add(var); vars[name] = var; } }
public bool CompileAll() { bool isValid = true; foreach (string RO in ranges.Keys) { solver.SetGlobalVariable(RO, -666); ranges[RO].Min.compileTokens(); expValidity[RO] = ranges[RO].Min.GenerateAKSolver(solver); ranges[RO].Max.compileTokens(); expValidity[RO] &= ranges[RO].Max.GenerateAKSolver(solver); isValid &= expValidity[RO]; } foreach (string EX in expressions.Keys) { expressions[EX].compileTokens(); expValidity[EX.ToString()] = expressions[EX].GenerateAKSolver(solver); isValid &= expValidity[EX.ToString()]; } StatisticsTracking.InstantEvent("Expression Value", "Value Updated", new Dictionary <string, object>() { { "valid", isValid } }); return(isValid); }
// Use this for initialization void setArgs(params double[] values) { for (int i = 0; i < depVars.Length; i++) { solver.SetGlobalVariable(depVars[i], values[i]); } }
// Use this for initialization void Awake() { mesh = GetComponent <MeshFilter>().mesh; // aww yisss Awake > Start for variable value initialization solver = GlobalDataStore.getSolver(); zFn = GlobalDataStore.getZFn(); // incase we start from "main" scene directly. if (solver == null) { solver = new AK.ExpressionSolver(); solver.SetGlobalVariable("x", 0); solver.SetGlobalVariable("y", 0); zFn = solver.SymbolicateExpression("(x^2 + y^2) / (x*y+1)"); } }
public bool CompileAll() { coefValidity = true; foreach (string PO in eqnCoefs.Keys) { solver.SetGlobalVariable(PO, -666); eqnCoefs[PO].compileTokens(); coefValidity &= eqnCoefs[PO].GenerateAKSolver(solver); } return(coefValidity); }
public static void TestExpLocalConstants() { ExpressionSolver solver = new ExpressionSolver(); var exp1 = solver.SymbolicateExpression("test+1", new string[] { "test" }); exp1.SetVariable("test", 1.0); var exp2 = solver.SymbolicateExpression("test+1", "test"); // If you define only one variable, you don't need to use the string[] format exp2.SetVariable("test", 2.0); solver.SetGlobalVariable("test", 1000); // If there is name clash with a exp-local variable, we prefer the exp-local variable. AssertSameValue(exp1.Evaluate(), 2); AssertSameValue(exp2.Evaluate(), 3); var exp3 = solver.SymbolicateExpression("test^2"); AssertSameValue(exp3.Evaluate(), 1000 * 1000); }
public bool CompileAll() { bool isValid = true; foreach (string PO in ptCoords.Keys) { solver.SetGlobalVariable(PO, -666); ptCoords[PO].X.compileTokens(); expValidity[PO] = ptCoords[PO].X.GenerateAKSolver(solver); ptCoords[PO].Y.compileTokens(); expValidity[PO] &= ptCoords[PO].Y.GenerateAKSolver(solver); ptCoords[PO].Z.compileTokens(); expValidity[PO] &= ptCoords[PO].Z.GenerateAKSolver(solver); isValid &= expValidity[PO]; } return(isValid); }
void Start() { // config depVars = new string[2] { "x", "y" }; // object retrieval InputField = GameObject.Find("UIInputFunction"); // init solver = new AK.ExpressionSolver(); // logic foreach (string depVar in depVars) { solver.SetGlobalVariable(depVar, 0); } }
public bool CompileAll() { bool isValid = true; foreach (string RO in ranges.Keys) { solver.SetGlobalVariable(RO, -666); ranges[RO].Min.compileTokens(); expValidity[RO] = ranges[RO].Min.GenerateAKSolver(solver); ranges[RO].Max.compileTokens(); expValidity[RO] &= ranges[RO].Max.GenerateAKSolver(solver); isValid &= expValidity[RO]; } foreach (ExpOptions EX in expressions.Keys) { expressions[EX].compileTokens(); expValidity[EX.ToString()] = expressions[EX].GenerateAKSolver(solver); isValid &= expValidity[EX.ToString()]; } return(isValid); }
public Vector3 findXYZ(Vector3 uvw) { Vector3 output = new Vector3(); ExpressionSet es = calcManager.expressionSet; float scale = calcManager.paramSurface.currentScale; foreach (ExpressionSet.ExpOptions key in es.expressions.Keys) { AK.ExpressionSolver solver = es.solver; if (es.ranges.ContainsKey("t")) { solver.SetGlobalVariable("t", uvw.x); } } if (es.expressions[X].AKExpression != null && es.expressions[Y].AKExpression != null && es.expressions[Z].AKExpression != null) { output.x = (float)es.expressions[X].AKExpression.Evaluate(); output.y = (float)es.expressions[Y].AKExpression.Evaluate(); output.z = (float)es.expressions[Z].AKExpression.Evaluate(); } return(output * scale); }
public static void TestExpLocalConstants() { ExpressionSolver solver = new ExpressionSolver(); var exp1 = solver.SymbolicateExpression("test+1",new string[]{"test"}); exp1.SetVariable("test",1.0); var exp2 = solver.SymbolicateExpression("test+1","test"); // If you define only one variable, you don't need to use the string[] format exp2.SetVariable("test",2.0); solver.SetGlobalVariable("test",1000); // If there is name clash with a exp-local variable, we prefer the exp-local variable. AssertSameValue(exp1.Evaluate(),2); AssertSameValue(exp2.Evaluate(),3); var exp3 = solver.SymbolicateExpression("test^2"); AssertSameValue(exp3.Evaluate(),1000*1000); }
public static void TestStringFuncs() { ExpressionSolver solver = new ExpressionSolver(); solver.AddCustomFunction("strlen",1, delegate(object[] p) { return ((string)p[0]).Length; },true); var exp = solver.SymbolicateExpression("strlen('123')"); AssertSameValue(exp.Evaluate(),3.0); exp = solver.SymbolicateExpression("strlen('12\\'3')"); AssertSameValue(exp.Evaluate(),4.0); exp = solver.SymbolicateExpression("strlen('12\\'3 4')"); AssertSameValue(exp.Evaluate(),6.0); solver.AddCustomFunction("strlen2",2, delegate(object[] p) { return ((string)p[0]).Length*(double)p[1]; }); exp = solver.SymbolicateExpression("strlen2('12\\'3 4',2.5)"); AssertSameValue(exp.Evaluate(),6.0*2.5); string[] erroneousStrings = new string[]{"strlen(''')","strlen('''')","''"}; foreach (var errorString in erroneousStrings) { try { exp = solver.SymbolicateExpression(errorString); throw new System.Exception("ExpressionSolverTest failed"); } catch (ESSyntaxErrorException) { // Parameters were not given correctly - syntax error expected } catch (System.Exception) { throw new System.Exception("ExpressionSolverTest failed"); } } // Because strlen should be evaluated at symbolication time, the following should reduce to one real value symbol: exp = solver.SymbolicateExpression("(1+strlen('123')+1)/strlen('12345')"); Assert(exp.root.type == SymbolType.RealValue); AssertSameValue(exp.root.value,1); // But if one of the parameters is not constant, then we cant do it: exp = solver.SymbolicateExpression("strlen('123')+x", new string[]{"x"}); Assert(exp.root.type == SymbolType.SubExpression); // Test string variables. Both exp-local and global exp = solver.SymbolicateExpression("strlen(stringVariableTest)","$stringVariableTest"); exp.SetVariable("stringVariableTest","test"); AssertSameValue(exp.Evaluate(),4); try { exp.SetVariable("stringVariableTest",121); Assert(false); } catch (ESParameterTypeChangedException) { } solver.SetGlobalVariable("striva","123"); exp = solver.SymbolicateExpression("strlen( ( striva ) )/3"); AssertSameValue(exp.Evaluate(),1); try { solver.SetGlobalVariable("striva",42141.0); Assert(false); } catch (ESParameterTypeChangedException) { } }
public static void TestGlobalConstants() { ExpressionSolver solver = new ExpressionSolver(); solver.SetGlobalVariable("test",1); var exp1 = solver.SymbolicateExpression("test+1"); AssertSameValue(2.0,exp1.Evaluate()); solver.SetGlobalVariable("test",2); var exp2 = solver.SymbolicateExpression("test+1"); AssertSameValue(3.0,exp2.Evaluate()); AssertSameValue(exp1.Evaluate(),exp2.Evaluate()); }
public static void TestFuncs() { ExpressionSolver solver = new ExpressionSolver(); solver.SetGlobalVariable("zero",0); var exp1 = solver.SymbolicateExpression("sin(pi/2)-cos(zero)"); AssertSameValue(exp1.Evaluate(),0); var exp2 = solver.SymbolicateExpression("2*e^zero - exp(zero)"); AssertSameValue(exp2.Evaluate(),1); var exp3 = solver.SymbolicateExpression("log(e^6)"); AssertSameValue(exp3.Evaluate(),6); var exp4 = solver.SymbolicateExpression("sqrt(2)-2^0.5"); AssertSameValue(exp4.Evaluate(),0); var exp5 = solver.SymbolicateExpression("exp(log(6))"); AssertSameValue(exp5.Evaluate(),6); var rnd = new System.Random(); solver.AddCustomFunction("Rnd1",2, delegate(double[] p) { return p[0] + (p[1]-p[0])*(rnd.NextDouble()); },false); var exp6 = solver.SymbolicateExpression("Rnd1(0,1)"); var firstRnd = exp6.Evaluate(); int iter = 0; while (true) { var secondRnd = exp6.Evaluate(); if (firstRnd != secondRnd) { break; } iter++; if (iter==10000) { // Probability of this happening is miniscule if everything works as it should throw new System.Exception("ExpressionSolverTest failed"); } } solver.AddCustomFunction("Rnd2",2, delegate(double[] p) { return p[0] + (p[1]-p[0])*(rnd.NextDouble()); },true); var exp7 = solver.SymbolicateExpression("Rnd2(0,1)"); AssertSameValue(exp7.Evaluate(),exp7.Evaluate()); var exp8 = solver.SymbolicateExpression("cos(0)+1*2"); AssertSameValue(exp8.Evaluate(),3); solver.AddCustomFunction("dist",5, delegate(double[] p) { return System.Math.Pow( (p[2]-p[0])*(p[2]-p[0]) + (p[3]-p[1])*(p[3]-p[1]) ,p[4]); },true); var exp9 = solver.SymbolicateExpression("dist(3*x,(4*x),+6*x,-1*x,sin(x))","x"); double x = 21; exp9.SetVariable("x",x); AssertSameValue(exp9.Evaluate(),System.Math.Pow( (3*x-6*x)*(3*x-6*x)+(4*x+x)*(4*x+x),System.Math.Sin(x) )); }
public static void TestStringFuncs() { ExpressionSolver solver = new ExpressionSolver(); solver.AddCustomFunction("strlen", 1, delegate(object[] p) { return(((string)p[0]).Length); }, true); var exp = solver.SymbolicateExpression("strlen('123')"); AssertSameValue(exp.Evaluate(), 3.0); exp = solver.SymbolicateExpression("strlen('12\\'3')"); AssertSameValue(exp.Evaluate(), 4.0); exp = solver.SymbolicateExpression("strlen('12\\'3 4')"); AssertSameValue(exp.Evaluate(), 6.0); solver.AddCustomFunction("strlen2", 2, delegate(object[] p) { return(((string)p[0]).Length * (double)p[1]); }); exp = solver.SymbolicateExpression("strlen2('12\\'3 4',2.5)"); AssertSameValue(exp.Evaluate(), 6.0 * 2.5); string[] erroneousStrings = new string[] { "strlen(''')", "strlen('''')", "''" }; foreach (var errorString in erroneousStrings) { try { exp = solver.SymbolicateExpression(errorString); throw new System.Exception("ExpressionSolverTest failed"); } catch (ESSyntaxErrorException) { // Parameters were not given correctly - syntax error expected } catch (System.Exception) { throw new System.Exception("ExpressionSolverTest failed"); } } // Because strlen should be evaluated at symbolication time, the following should reduce to one real value symbol: exp = solver.SymbolicateExpression("(1+strlen('123')+1)/strlen('12345')"); Assert(exp.root.type == SymbolType.RealValue); AssertSameValue(exp.root.value, 1); // But if one of the parameters is not constant, then we cant do it: exp = solver.SymbolicateExpression("strlen('123')+x", new string[] { "x" }); Assert(exp.root.type == SymbolType.SubExpression); // Test string variables. Both exp-local and global exp = solver.SymbolicateExpression("strlen(stringVariableTest)", "$stringVariableTest"); exp.SetVariable("stringVariableTest", "test"); AssertSameValue(exp.Evaluate(), 4); try { exp.SetVariable("stringVariableTest", 121); Assert(false); } catch (ESParameterTypeChangedException) { } solver.SetGlobalVariable("striva", "123"); exp = solver.SymbolicateExpression("strlen( ( striva ) )/3"); AssertSameValue(exp.Evaluate(), 1); try { solver.SetGlobalVariable("striva", 42141.0); Assert(false); } catch (ESParameterTypeChangedException) { } }
void SamplePoints() { Vector3 start = referencePoint.lastLocalPos; //lck = new object(); //thread_num = SystemInfo.processorCount; //thread_num = 1; //thread_num = (thread_num > 2) ? 2 : thread_num; tmin = (float)solver.EvaluateExpression(t_min); tmin = (tmin > 0) ? 0 : tmin; tmax = (float)solver.EvaluateExpression(t_max); tmax = (tmax < 0) ? 0 : tmax; currLocalPos = referencePoint.lastLocalPos; currExpX = vectorField.expressionX; currExpY = vectorField.expressionY; currExpZ = vectorField.expressionZ; int lastCount = positions.Count; positions.Clear(); solver = new AK.ExpressionSolver(); solver.SetGlobalVariable("x", referencePoint.lastLocalPos.x); solver.SetGlobalVariable("y", referencePoint.lastLocalPos.y); solver.SetGlobalVariable("z", referencePoint.lastLocalPos.z); expX = solver.SymbolicateExpression(vectorField.es.expressions[ExpressionSet.ExpOptions.X].expression); expY = solver.SymbolicateExpression(vectorField.es.expressions[ExpressionSet.ExpOptions.Y].expression); expZ = solver.SymbolicateExpression(vectorField.es.expressions[ExpressionSet.ExpOptions.Z].expression); varX = solver.GetGlobalVariable("x"); varY = solver.GetGlobalVariable("y"); varZ = solver.GetGlobalVariable("z"); //solvers = new AK.ExpressionSolver[thread_num]; //expXs = new AK.Expression[thread_num]; //expYs = new AK.Expression[thread_num]; //expZs = new AK.Expression[thread_num]; //for(int i = 0; i < thread_num; i++) //{ // solvers[i] = new AK.ExpressionSolver(); // solvers[i].SetGlobalVariable("x", 0); // solvers[i].SetGlobalVariable("y", 0); // solvers[i].SetGlobalVariable("z", 0); // expXs[i] = solvers[i].SymbolicateExpression(vectorField.es.expressions[ExpressionSet.ExpOptions.X].expression); // expYs[i] = solvers[i].SymbolicateExpression(vectorField.es.expressions[ExpressionSet.ExpOptions.Y].expression); // expZs[i] = solvers[i].SymbolicateExpression(vectorField.es.expressions[ExpressionSet.ExpOptions.Z].expression); //} //Thread[] threads = new Thread[thread_num]; float positiveCount = tmax / time_step; positiveCount = (positiveCount > 50000) ? 50000 : positiveCount; float negativeCount = -tmin / time_step; negativeCount = (negativeCount > 50000) ? 50000 : negativeCount; //Vector3[] startPts = new Vector3[thread_num]; //startPts[0] = referencePoint.lastLocalPos; //for(int i = 1; i < thread_num; i++) //{ // startPts[i] = RK4(startPts[i - 1], time_step); //} //for(int i = 0; i < thread_num; i++) //{ // int index = i; // threads[i] = new Thread(() => ThreadedSampling(index, startPts[index], time_step * thread_num, positiveCount / thread_num, // negativeCount / thread_num)); // threads[i].Start(); //} ////for (int i = 0; i < 5; i++) ////{ //// yield return null; ////} //for (int i = 0; i < thread_num; i++) //{ // threads[i].Join(); //} Vector3 curr = start; for (int i = 0; i < positiveCount; i++) { curr = RK4(curr, time_step); positions.Add(i + 1, curr); } curr = start; for (int i = 0; i < negativeCount; i++) { curr = RK4(curr, -time_step); positions.Add(-i, curr); } if (positions.Count != lastCount) { InitializeParticleSystem(); } //RenderParticles(); currHighlight = 0; thread_finished = true; }
public static void TestFuncs() { ExpressionSolver solver = new ExpressionSolver(); solver.SetGlobalVariable("zero", 0); var exp1 = solver.SymbolicateExpression("sin(pi/2)-cos(zero)"); AssertSameValue(exp1.Evaluate(), 0); var exp2 = solver.SymbolicateExpression("2*e^zero - exp(zero)"); AssertSameValue(exp2.Evaluate(), 1); var exp3 = solver.SymbolicateExpression("log(e^6)"); AssertSameValue(exp3.Evaluate(), 6); var exp4 = solver.SymbolicateExpression("sqrt(2)-2^0.5"); AssertSameValue(exp4.Evaluate(), 0); var exp5 = solver.SymbolicateExpression("exp(log(6))"); AssertSameValue(exp5.Evaluate(), 6); var rnd = new System.Random(); solver.AddCustomFunction("Rnd1", 2, delegate(double[] p) { return(p[0] + (p[1] - p[0]) * (rnd.NextDouble())); }, false); var exp6 = solver.SymbolicateExpression("Rnd1(0,1)"); var firstRnd = exp6.Evaluate(); int iter = 0; while (true) { var secondRnd = exp6.Evaluate(); if (firstRnd != secondRnd) { break; } iter++; if (iter == 10000) { // Probability of this happening is miniscule if everything works as it should throw new System.Exception("ExpressionSolverTest failed"); } } solver.AddCustomFunction("Rnd2", 2, delegate(double[] p) { return(p[0] + (p[1] - p[0]) * (rnd.NextDouble())); }, true); var exp7 = solver.SymbolicateExpression("Rnd2(0,1)"); AssertSameValue(exp7.Evaluate(), exp7.Evaluate()); var exp8 = solver.SymbolicateExpression("cos(0)+1*2"); AssertSameValue(exp8.Evaluate(), 3); solver.AddCustomFunction("dist", 5, delegate(double[] p) { return(System.Math.Pow((p[2] - p[0]) * (p[2] - p[0]) + (p[3] - p[1]) * (p[3] - p[1]), p[4])); }, true); var exp9 = solver.SymbolicateExpression("dist(3*x,(4*x),+6*x,-1*x,sin(x))", "x"); double x = 21; exp9.SetVariable("x", x); AssertSameValue(exp9.Evaluate(), System.Math.Pow((3 * x - 6 * x) * (3 * x - 6 * x) + (4 * x + x) * (4 * x + x), System.Math.Sin(x))); }