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()));
}
}