public void Init(string constructionString)
{
dataSet = DataSet.Create(constructionString);
//Инициализируем DataSet
Variable X = dataSet.AddVariable <double>("X", "x");
Variable Y = dataSet.AddVariable <double>("Y", "y");
Variable Z = dataSet.AddVariable <double>("Z", "z");
time = dataSet.AddVariable <double>("Time", "t");
u = dataSet.AddVariable <double>("U velocity", "x", "y", "z", "t");
v = dataSet.AddVariable <double>("V velocity", "x", "y", "z", "t");
w = dataSet.AddVariable <double>("W velocity", "x", "y", "z", "t");
T = dataSet.AddVariable <double>("Temperature", "x", "y", "z", "t");
div = dataSet.AddVariable <double>("Divergence", "x", "y", "z", "t");
double[] wArr = new double[modellingParams.Nx];
for (int i = 0; i < modellingParams.Nx; i++)
{
wArr[i] = i * modellingParams.Dx;
}
X.PutData(wArr);
wArr = new double[modellingParams.Ny];
for (int i = 0; i < modellingParams.Ny; i++)
{
wArr[i] = i * modellingParams.Dy;
}
Y.PutData(wArr);
wArr = new double[modellingParams.Nz];
for (int i = 0; i < modellingParams.Nz; i++)
{
wArr[i] = i * modellingParams.Dz;
}
Z.PutData(wArr);
//Инициализируем рассчетный модуль для слоя начальными условиями
solver = new LayerSolver(prevData, modellingParams);
u.Append(prevData.U.ToArray(), "t");
v.Append(prevData.V.ToArray(), "t");
w.Append(prevData.W.ToArray(), "t");
T.Append(prevData.T.ToArray(), "t");
div.Append(prevData.Div.ToArray(), "t");
time.PutData(new double[1] {
0
});
dataSet.Commit();
}
public void Init(string constructionString)
{
dataSet = DataSet.Create(constructionString);
//Инициализируем DataSet
Variable X = dataSet.AddVariable<double>("X", "x");
Variable Y = dataSet.AddVariable<double>("Y", "y");
Variable Z = dataSet.AddVariable<double>("Z", "z");
time = dataSet.AddVariable<double>("Time", "t");
u = dataSet.AddVariable<double>("U velocity", "x", "y", "z", "t");
v = dataSet.AddVariable<double>("V velocity", "x", "y", "z", "t");
w = dataSet.AddVariable<double>("W velocity", "x", "y", "z", "t");
T = dataSet.AddVariable<double>("Temperature", "x", "y", "z", "t");
div = dataSet.AddVariable<double>("Divergence", "x", "y", "z", "t");
double[] wArr = new double[modellingParams.Nx];
for (int i = 0; i < modellingParams.Nx; i++)
{
wArr[i] = i * modellingParams.Dx;
}
X.PutData(wArr);
wArr = new double[modellingParams.Ny];
for (int i = 0; i < modellingParams.Ny; i++)
{
wArr[i] = i * modellingParams.Dy;
}
Y.PutData(wArr);
wArr = new double[modellingParams.Nz];
for (int i = 0; i < modellingParams.Nz; i++)
{
wArr[i] = i * modellingParams.Dz;
}
Z.PutData(wArr);
//Инициализируем рассчетный модуль для слоя начальными условиями
solver = new LayerSolver(prevData, modellingParams);
u.Append(prevData.U.ToArray(), "t");
v.Append(prevData.V.ToArray(), "t");
w.Append(prevData.W.ToArray(), "t");
T.Append(prevData.T.ToArray(), "t");
div.Append(prevData.Div.ToArray(), "t");
time.PutData(new double[1] { 0 });
dataSet.Commit();
}
public void PerformIteration(int i)
{
//Основной рассчет
//for (int i = 1; i < modellingParams.Nt; i++)
//{
LayerData result = solver.Solve(true);
//Кладем данные в DataSet
u.Append(result.U.ToArray(), "t");
v.Append(result.V.ToArray(), "t");
w.Append(result.W.ToArray(), "t");
T.Append(result.T.ToArray(), "t");
div.Append(result.Div.ToArray(), "t");
time.Append(new double[1] {
(double)i / modellingParams.Nt
});
dataSet.Commit();
//Переходим на следующий слой
solver = new LayerSolver(prevData, result, modellingParams);
prevData = result;
double temp = 0;
int count = result.Width * result.Height * result.Thickness;
for (int ii = 1; ii < result.Width; ii++)
{
for (int jj = 1; jj < result.Height; jj++)
{
for (int kk = 1; kk < result.Thickness; kk++)
{
temp += result.Div[ii, jj, kk];
}
}
}
temp = temp / count * modellingParams.Dx * modellingParams.Dy * modellingParams.Dz;
Console.WriteLine((double)i / modellingParams.Nt * 100 + "% Error = " + temp);
//}
//dataSet.Commit();
}
public void SolveAll(string ctorString)
{
//dataSet = ProxyDataSet.Open("msds:nc?file=../../../temp.nc");
dataSet = ProxyDataSet.Open("msds:memory");
//Инициализируем DataSet
Variable X = dataSet.AddVariable <double>("X", "x");
Variable Y = dataSet.AddVariable <double>("Y", "y");
Variable Z = dataSet.AddVariable <double>("Z", "z");
Variable time = dataSet.AddVariable <double>("Time", "t");
Variable u = dataSet.AddVariable <double>("U velocity", "x", "y", "z", "t");
Variable v = dataSet.AddVariable <double>("V velocity", "x", "y", "z", "t");
Variable w = dataSet.AddVariable <double>("W velocity", "x", "y", "z", "t");
Variable T = dataSet.AddVariable <double>("Temperature", "x", "y", "z", "t");
Variable div = dataSet.AddVariable <double>("Divergence", "x", "y", "z", "t");
dataSet.Commit();
double[] wArr = new double[modellingParams.Nx];
for (int i = 0; i < modellingParams.Nx; i++)
{
wArr[i] = i * modellingParams.Dx;
}
X.PutData(wArr);
wArr = new double[modellingParams.Ny];
for (int i = 0; i < modellingParams.Ny; i++)
{
wArr[i] = i * modellingParams.Dy;
}
Y.PutData(wArr);
wArr = new double[modellingParams.Nz];
for (int i = 0; i < modellingParams.Nz; i++)
{
wArr[i] = i * modellingParams.Dz;
}
Z.PutData(wArr);
//Инициализируем рассчетный модуль для слоя начальными условиями
LayerSolver solver = new LayerSolver(prevData, modellingParams);
u.Append(prevData.U.ToArray(), "t");
v.Append(prevData.V.ToArray(), "t");
w.Append(prevData.W.ToArray(), "t");
T.Append(prevData.T.ToArray(), "t");
div.Append(prevData.Div.ToArray(), "t");
time.PutData(new double[1] {
0
});
dataSet.Commit();
//Основной рассчет
for (int i = 1; i < modellingParams.Nt; i++)
{
LayerData result = solver.Solve(true);
//Кладем данные в DataSet
u.Append(result.U.ToArray(), "t");
v.Append(result.V.ToArray(), "t");
w.Append(result.W.ToArray(), "t");
T.Append(result.T.ToArray(), "t");
div.Append(result.Div.ToArray(), "t");
time.Append(new double[1] {
(double)i / modellingParams.Nt
});
dataSet.Commit();
//Переходим на следующий слой
solver = new LayerSolver(prevData, result, modellingParams);
prevData = result;
double temp = 0;
int count = 0;
for (int ii = 1; ii < result.Width; ii++)
{
for (int jj = 1; jj < result.Height; jj++)
{
for (int kk = 1; kk < result.Thickness; kk++)
{
temp += result.Div[ii, jj, kk];
count++;
}
}
}
temp = temp / count * modellingParams.Dx * modellingParams.Dy * modellingParams.Dz;
Console.WriteLine((double)i / modellingParams.Nt * 100 + "% Error = " + temp);
}
dataSet.Commit();
}
public static void AddVariousDataToFile(Dictionary <string, object> dataToWrite, string fileName)
{
if (!ServiceTools.CheckIfDirectoryExists(fileName))
{
return;
}
string connectionString = "msds:nc?file=";
connectionString += fileName;
NetCDFDataSet ds = null;
try
{
ds = new NetCDFDataSet(connectionString, ResourceOpenMode.OpenOrCreate);
}
catch (Exception ex)
{
throw ex;
}
foreach (KeyValuePair <string, object> keyValuePair in dataToWrite)
{
if (keyValuePair.Value.GetType() == typeof(DenseMatrix))
{
if (ds.Variables.Contains(keyValuePair.Key))
{
Variable <double> theDataVar = (Variable <double>)ds.Variables[keyValuePair.Key];
theDataVar.Append(((DenseMatrix)keyValuePair.Value).ToArray());
}
else
{
Variable <double> theDataVar = ds.AddVariable <double>(keyValuePair.Key,
((DenseMatrix)(keyValuePair.Value)).ToArray(), "y", "x");
}
}
if (keyValuePair.Value.GetType() == typeof(DenseVector))
{
if (ds.Variables.Contains(keyValuePair.Key))
{
Variable <double> theDataVar = (Variable <double>)ds.Variables[keyValuePair.Key];
theDataVar.Append(((DenseVector)keyValuePair.Value).ToArray());
}
else
{
Variable <double> theDataVar = ds.AddVariable <double>(keyValuePair.Key,
((DenseVector)(keyValuePair.Value)).ToArray(), keyValuePair.Key);
}
}
if (keyValuePair.Value.GetType() == typeof(long[]))
{
if (ds.Variables.Contains(keyValuePair.Key))
{
Variable <long> theDataVar = (Variable <long>)ds.Variables[keyValuePair.Key];
theDataVar.Append((long[])(keyValuePair.Value));
}
else
{
Variable <long> theDataVar = ds.AddVariable <long>(keyValuePair.Key, (long[])(keyValuePair.Value),
keyValuePair.Key);
}
}
if (keyValuePair.Value.GetType() == typeof(int[]))
{
if (ds.Variables.Contains(keyValuePair.Key))
{
Variable <int> theDataVar = (Variable <int>)ds.Variables[keyValuePair.Key];
theDataVar.Append((int[])(keyValuePair.Value));
}
else
{
Variable <int> theDataVar = ds.AddVariable <int>(keyValuePair.Key, (int[])(keyValuePair.Value),
keyValuePair.Key);
}
}
if (keyValuePair.Value.GetType() == typeof(string[]))
{
if (ds.Variables.Contains(keyValuePair.Key))
{
Variable <string> theDataVar = (Variable <string>)ds.Variables[keyValuePair.Key];
theDataVar.Append((string[])(keyValuePair.Value));
}
else
{
Variable <string> theDataVar = ds.AddVariable <string>(keyValuePair.Key, (string[])(keyValuePair.Value),
keyValuePair.Key);
}
}
}
try
{
ds.TryCommit();
}
catch (Exception ex)
{
throw ex;
}
ds.Dispose();
}
