protected override DataSet ServerProcessInternal(DataSet ds)
{
if (serviceUri == "(local)")
{
return(LocalProcess(ds));
}
string hash = DataSetDiskCache.ComputeHash(ds);
DataSet proxyDataSet = null;
// Creating new DataSet at the service.
// TODO: fix following:
try
{
try
{
proxyDataSet = ProxyDataSet.CreateProxySync(taskQueue, ServicePort, "msds:memory", false, 10 * 60 * 1000);
}
catch (CommunicationObjectFaultedException)
{
//Connection to server closed.
//Recreate service port and try again.
if (proxyDataSet != null && !proxyDataSet.IsDisposed)
{
proxyDataSet.Dispose();
}
this._servicePort = null;
proxyDataSet = ProxyDataSet.CreateProxySync(taskQueue, ServicePort, "msds:memory", false, 10 * 60 * 1000);
}
AutoResetEvent completed = new AutoResetEvent(false);
OnCommittedHandler onCommitted = new OnCommittedHandler(completed, OnDataSetCommitted);
proxyDataSet.Committed += onCommitted.Handler;
proxyDataSet.IsAutocommitEnabled = false;
FetchClimateRequestBuilder.CopyRequestedDataSet(ds, proxyDataSet, false);
proxyDataSet.Metadata[Namings.metadataNameHash] = hash;
proxyDataSet.Commit();
if (proxyDataSet.HasChanges)
{
proxyDataSet.Commit();
}
completed.WaitOne();
proxyDataSet.IsAutocommitEnabled = true;
return(proxyDataSet);
}
catch
{
if (proxyDataSet != null && !proxyDataSet.IsDisposed)
{
proxyDataSet.Dispose();
}
throw;
}
}
public static void Main(string[] args)
{
var port = ProxyDataSet.Create("msds:nc?file=../../../temp.nc");
NumericalParameters nPar = new NumericalParameters(0.01, 0.02, 0.02, 0.01, 40, 40, 50, 40, 150, 0.78, 1.4);
solver = new FluidCurrentSolver(nPar);
//DataSetFactory.Register(typeof(NetCDFDataSet));
solver.SolveAll("msds:nc?file=../../../temp.nc");
Console.WriteLine("Done!");
}
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();
}
