public static Params(int Nin, double Lengthin, double Timein, int Ntin) { int i; N = Nin; Length = Lengthin; Time = Timein; Nt = Ntin; tau = Timein / Nt; h = Length / N; for (i = 0; i < N / 2; i++) { Ur.Add(0); P.Add(30 * 101325); Ro.Add(30 * 1.29); E.Add(Energy(P[i], Ro[i])); } for (i = N / 2; i < N; i++) { Ur.Add(0); P.Add(101325); Ro.Add(1.29); E.Add(Energy(P[i], Ro[i])); } }
private void ParseNodes(XContainer root) { var streamElements = GetElements(root, "Stream"); StreamsCount = streamElements.Count == 0 ? 1: streamElements.Max(element => int.Parse(element.Attribute("Index").Value)) + 1; for (int stream = 0; stream < StreamsCount; stream++) { Mu.Add(new Vector <double>( ParseMu(StreamsCount == 1? GetElement(root, "Mu") : GetElements(root, "Mu") .Single(element => GetAttributeValue <int>(element.Parent, "Index") == stream)))); Lambda.Add(new Vector <double>( ParseLambda(StreamsCount == 1 ? GetElement(root, "Lambda") : GetElements(root, "Lambda") .Single(element => GetAttributeValue <int>(element.Parent, "Index") == stream), stream))); Ro.Add(Lambda.Last().DivideElementWise(Mu.Last())); } }