public static string SemiLagrangeScheme(EquationParameters parameters) { StringBuilder result = new StringBuilder(); double[] data = new double[(int)(parameters.MaxVar / parameters.VarStep)]; data = EquationMethods.InitialConditionsSetter(parameters, data.Length); for (int i = 0; i < (int)(parameters.MaxTime / parameters.TimeStep); i++) { double[] reservoir = new double[data.Length]; double[] res = new double[data.Length]; int approx; int secondapprox; for (int j = 0; j < reservoir.Length; j++) { reservoir[j] = j - parameters.Coeffitient * parameters.TimeStep / parameters.VarStep; res[j] = data[j]; } for (int j = 1; j < reservoir.Length; j++) { if (reservoir[j] > 0) { approx = (int)reservoir[j]; if (approx > reservoir[j]) { secondapprox = approx - 1; } else if (approx < reservoir[j]) { secondapprox = approx + 1; } else { secondapprox = approx; } data[j] = Math.Abs(approx - reservoir[j]) * res[approx] + Math.Abs(secondapprox - reservoir[j]) * res[secondapprox]; } else { data[j] = data[0]; } } data[0] = res[reservoir.Length - 1]; result.Append(EquationMethods.FileWriter(data, parameters, i)); } return(result.ToString()); }
public static string UpwindScheme(EquationParameters parameters) { double[] data = new double[(int)(parameters.MaxVar / parameters.VarStep)]; data = EquationMethods.InitialConditionsSetter(parameters, data.Length); StringBuilder result = new StringBuilder(); for (int i = 0; i < (int)(parameters.MaxTime / parameters.TimeStep); i++) { double[] reservoir = new double[data.Length]; for (int j = 0; j < reservoir.Length; j++) { reservoir[j] = data[j]; } for (int j = 1; j < reservoir.Length; j++) { data[j] = reservoir[j] - parameters.TimeStep * parameters.Coeffitient * (reservoir[j] - reservoir[j - 1]) / parameters.VarStep; } data[0] = reservoir[reservoir.Length - 1]; result.Append(EquationMethods.FileWriter(data, parameters, i)); } return(result.ToString()); }
public static string CentralDifferences(EquationParameters parameters) { double[] data = new double[(int)(parameters.MaxVar / parameters.VarStep)]; data = EquationMethods.InitialConditionsSetter(parameters, data.Length); StringBuilder result = new StringBuilder(); //FileStream stream = new FileStream("C:/Users/mm/Documents/Visual Studio 2017/Projects/EquationsSolver/WebGUI/DataStorage/data.json", FileMode.Create, FileAccess.Write); for (int i = 0; i < (int)(parameters.MaxTime / parameters.TimeStep); i++) { double[] reservoir = new double[data.Length]; for (int j = 0; j < reservoir.Length; j++) { reservoir[j] = data[j]; } for (int j = 1; j < reservoir.Length - 1; j++) { data[j] = reservoir[j] + parameters.TimeStep * parameters.Coeffitient * (reservoir[j + 1] - 2.0 * reservoir[j] + reservoir[j - 1]) / (parameters.VarStep * parameters.VarStep); } result.Append(EquationMethods.FileWriter(data, parameters, i)); } return(result.ToString()); //stream.Close(); }