Newer
Older
//=======================================================================================
// ____ ____ __ ______ __________ __ __ __ __
// \ \ | | | | | _ \ |___ ___| | | | | / \ | |
// \ \ | | | | | |_) | | | | | | | / \ | |
// \ \ | | | | | _ / | | | | | | / /\ \ | |
// \ \ | | | | | | \ \ | | | \__/ | / ____ \ | |____
// \ \ | | |__| |__| \__\ |__| \________/ /__/ \__\ |_______|
// \ \ | | ________________________________________________________________
// \ \ | | | ______________________________________________________________|
// \ \| | | | __ __ __ __ ______ _______
// \ | | |_____ | | | | | | | | | _ \ / _____)
// \ | | _____| | | | | | | | | | | \ \ \_______
// \ | | | | |_____ | \_/ | | | | |_/ / _____ |
// \ _____| |__| |________| \_______/ |__| |______/ (_______/
//
// This file is part of VirtualFluids. VirtualFluids is free software: you can
// redistribute it and/or modify it under the terms of the GNU General Public
// License as published by the Free Software Foundation, either version 3 of
// the License, or (at your option) any later version.
//
// VirtualFluids is distributed in the hope that it will be useful, but WITHOUT
// ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or
// FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License
// for more details.
//
// You should have received a copy of the GNU General Public License along
// with VirtualFluids (see COPYING.txt). If not, see <http://www.gnu.org/licenses/>.
//
//! \file GridScalingKernelManager.h
//! \ingroup KernelManager
//! \author Martin Schoenherr
//=======================================================================================
#include "KernelManager/GridScalingKernelManager.h"
#include "GPU/CudaMemoryManager.h"
#include "GPU/GPU_Interface.h"
#include "Parameter/CudaStreamManager.h"
#include "PreCollisionInteractor/PreCollisionInteractor.h"
#include "Factories/GridScalingFactory.h"
GridScalingKernelManager::GridScalingKernelManager(SPtr<Parameter> parameter, GridScalingFactory *gridScalingFactory)
: para(parameter)
{
if(para->getMaxLevel() != 0){
if(!gridScalingFactory){
throw std::runtime_error("There is more than one level, but no scalingFactory was provided.");
}
checkScalingFunction(gridScalingFactory->getGridScalingFC(parameter->getUseTurbulentViscosity()), this->para->getParD(0)->intFC, "scalingFineToCoarse");
checkScalingFunction(gridScalingFactory->getGridScalingCF(parameter->getUseTurbulentViscosity()), this->para->getParD(0)->intCF, "scalingCoarseToFine");
this->scalingFineToCoarse = gridScalingFactory->getGridScalingFC(parameter->getUseTurbulentViscosity());
this->scalingCoarseToFine = gridScalingFactory->getGridScalingCF(parameter->getUseTurbulentViscosity());
if(this->scalingFineToCoarse == nullptr)
VF_LOG_TRACE("Function for scalingFineToCoarse is nullptr");
if(this->scalingCoarseToFine == nullptr)
VF_LOG_TRACE("Function for scalingCoarseToFine is nullptr");
void GridScalingKernelManager::runFineToCoarseKernelLB(const int level, InterpolationCellFineToCoarse *icellFC, ICellNeigh &offFC, CudaStreamIndex streamIndex) const
{
cudaStream_t stream = para->getStreamManager()->getStream(streamIndex);
this->scalingFineToCoarse(para->getParD(level).get(), para->getParD(level+1).get(), icellFC, offFC, stream);
// ScaleFC_comp_D3Q27F3(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->g6.g[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellFC->ICellFCC,
// icellFC->ICellFCF,
// icellFC->kFC,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC,
// stream);
// ScaleFC_0817_comp_27(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellFC->ICellFCC,
// icellFC->ICellFCF,
// icellFC->kFC,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC,
// stream);
// ScaleFC_RhoSq_3rdMom_comp_27(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellFC->ICellFCC,
// icellFC->ICellFCF,
// icellFC->kFC,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC,
// stream);
// ScaleFC_AA2016_comp_27(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellFC->ICellFCC,
// icellFC->ICellFCF,
// icellFC->kFC,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC,
// stream);
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
//////////////////////////////////////////////////////////////////////////
// D E P R E C A T E D
//////////////////////////////////////////////////////////////////////////
//ScaleFC27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->gridNX);
//getLastCudaError("ScaleFC27 execution failed");
//ScaleFCEff27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC);
//getLastCudaError("ScaleFC27 execution failed");
//ScaleFCLast27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC);
//getLastCudaError("ScaleFC27 execution failed");
//ScaleFCpress27(para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC);
//getLastCudaError("ScaleFC27 execution failed");
// ScaleFC_Fix_comp_27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC);
// getLastCudaError("ScaleFC27 execution failed");
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// incompressible
//ScaleFC_Fix_27(para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC);
//getLastCudaError("ScaleFC27 execution failed");
//ScaleFC_NSPress_27(para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC, para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offFC);
//getLastCudaError("ScaleFC27 execution failed");
}
void GridScalingKernelManager::runFineToCoarseKernelAD(const int level) const
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
{
//A D V E C T I O N D I F F U S I O N
if (para->getDiffMod() == 7)
{
// ScaleFCThS7(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->distributionsAD7.f[0],
// para->getParD(level+1)->distributionsAD7.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intFC.ICellFCC,
// para->getParD(level)->intFC.ICellFCF,
// para->getParD(level)->K_FC,
// para->getParD(level)->vis,
// para->getParD(level)->diffusivity,
// para->getParD(level)->numberofthreads);
ScaleFCThSMG7(
para->getParD(level)->distributions.f[0],
para->getParD(level+1)->distributions.f[0],
para->getParD(level)->distributionsAD7.f[0],
para->getParD(level+1)->distributionsAD7.f[0],
para->getParD(level)->neighborX,
para->getParD(level)->neighborY,
para->getParD(level)->neighborZ,
para->getParD(level+1)->neighborX,
para->getParD(level+1)->neighborY,
para->getParD(level+1)->neighborZ,
para->getParD(level)->numberOfNodes,
para->getParD(level+1)->numberOfNodes,
para->getParD(level)->isEvenTimestep,
para->getParD(level)->intFC.coarseCellIndices,
para->getParD(level)->intFC.fineCellIndices,
para->getParD(level)->intFC.numberOfCells,
para->getParD(level)->viscosity,
para->getParD(level)->diffusivity,
para->getParD(level)->numberofthreads,
para->getParD(level)->neighborFC);
}
else if (para->getDiffMod() == 27)
{
ScaleFCThS27(
para->getParD(level)->distributions.f[0],
para->getParD(level+1)->distributions.f[0],
Martin Schönherr
committed
para->getParD(level)->distributionsAD.f[0],
para->getParD(level+1)->distributionsAD.f[0],
para->getParD(level)->neighborX,
para->getParD(level)->neighborY,
para->getParD(level)->neighborZ,
para->getParD(level+1)->neighborX,
para->getParD(level+1)->neighborY,
para->getParD(level+1)->neighborZ,
para->getParD(level)->numberOfNodes,
para->getParD(level+1)->numberOfNodes,
para->getParD(level)->isEvenTimestep,
para->getParD(level)->intFC.coarseCellIndices,
para->getParD(level)->intFC.fineCellIndices,
para->getParD(level)->intFC.numberOfCells,
para->getParD(level)->viscosity,
para->getParD(level)->diffusivity,
para->getParD(level)->numberofthreads,
para->getParD(level)->neighborFC);
void GridScalingKernelManager::runCoarseToFineKernelLB(const int level, InterpolationCellCoarseToFine* icellCF, ICellNeigh &offCF, CudaStreamIndex streamIndex) const
cudaStream_t stream = para->getStreamManager()->getStream(streamIndex);
this->scalingCoarseToFine(para->getParD(level).get(), para->getParD(level+1).get(), icellCF, offCF, stream);
// ScaleCF_comp_D3Q27F3(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level+1)->g6.g[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellCF->ICellCFC,
// icellCF->ICellCFF,
// icellCF->kCF,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// ScaleCF_0817_comp_27(
// para->getParD(level)->distributions.f[0],
// para->getParD(level + 1)->distributions.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellCF->ICellCFC,
// icellCF->ICellCFF,
// icellCF->kCF,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// ScaleCF_RhoSq_3rdMom_comp_27(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellCF->ICellCFC,
// icellCF->ICellCFF,
// icellCF->kCF,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
// ScaleCF_AA2016_comp_27(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// icellCF->ICellCFC,
// icellCF->ICellCFF,
// icellCF->kCF,
// para->getParD(level)->omega,
// para->getParD(level+1)->omega,
// para->getParD(level)->vis,
// para->getParD(level)->nx,
// para->getParD(level)->ny,
// para->getParD(level+1)->nx,
// para->getParD(level+1)->ny,
// para->getParD(level)->numberofthreads,
//////////////////////////////////////////////////////////////////////////
// D E P R E C A T E D
//////////////////////////////////////////////////////////////////////////
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
//ScaleCF27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->gridNX);
//getLastCudaError("ScaleCF27 execution failed");
//ScaleCFEff27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offCF);
//getLastCudaError("ScaleCF27 execution failed");
//ScaleCFLast27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offCF);
//getLastCudaError("ScaleCF27 execution failed");
//ScaleCFpress27(para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offCF);
//getLastCudaError("ScaleCF27 execution failed");
// ScaleCF_Fix_comp_27( para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offCF);
// getLastCudaError("ScaleCF27 execution failed");
////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////////
// incompressible
//ScaleCF_Fix_27(para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offCF);
//getLastCudaError("ScaleCF27 execution failed");
//ScaleCF_NSPress_27(para->getParD(level)->distributions.f[0], para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->neighborX, para->getParD(level)->neighborY, para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX, para->getParD(level+1)->neighborY, para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes, para->getParD(level+1)->numberOfNodes, para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC, para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF, para->getParD(level)->omega, para->getParD(level+1)->omega,
// para->getParD(level)->vis, para->getParD(level)->nx, para->getParD(level)->ny,
// para->getParD(level+1)->nx, para->getParD(level+1)->ny, para->getParD(level)->numberofthreads,
// para->getParD(level)->offCF);
//getLastCudaError("ScaleCF27 execution failed");
}
void GridScalingKernelManager::runCoarseToFineKernelAD(const int level) const
// A D V E C T I O N D I F F U S I O N
if (para->getDiffMod() == 7)
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
// ScaleCFThS7(
// para->getParD(level)->distributions.f[0],
// para->getParD(level+1)->distributions.f[0],
// para->getParD(level)->distributionsAD7.f[0],
// para->getParD(level+1)->distributionsAD7.f[0],
// para->getParD(level)->neighborX,
// para->getParD(level)->neighborY,
// para->getParD(level)->neighborZ,
// para->getParD(level+1)->neighborX,
// para->getParD(level+1)->neighborY,
// para->getParD(level+1)->neighborZ,
// para->getParD(level)->numberOfNodes,
// para->getParD(level+1)->numberOfNodes,
// para->getParD(level)->isEvenTimestep,
// para->getParD(level)->intCF.ICellCFC,
// para->getParD(level)->intCF.ICellCFF,
// para->getParD(level)->K_CF,
// para->getParD(level)->vis,
// para->getParD(level+1)->diffusivity,
// para->getParD(level)->numberofthreads);
ScaleCFThSMG7(
para->getParD(level)->distributions.f[0],
para->getParD(level+1)->distributions.f[0],
para->getParD(level)->distributionsAD7.f[0],
para->getParD(level+1)->distributionsAD7.f[0],
para->getParD(level)->neighborX,
para->getParD(level)->neighborY,
para->getParD(level)->neighborZ,
para->getParD(level+1)->neighborX,
para->getParD(level+1)->neighborY,
para->getParD(level+1)->neighborZ,
para->getParD(level)->numberOfNodes,
para->getParD(level+1)->numberOfNodes,
para->getParD(level)->isEvenTimestep,
para->getParD(level)->intCF.coarseCellIndices,
para->getParD(level)->intCF.fineCellIndices,
para->getParD(level)->intCF.numberOfCells,
para->getParD(level)->viscosity,
para->getParD(level+1)->diffusivity,
para->getParD(level)->numberofthreads,
para->getParD(level)->neighborCF);
}
else if (para->getDiffMod() == 27)
{
ScaleCFThS27(
para->getParD(level)->distributions.f[0],
para->getParD(level+1)->distributions.f[0],
Martin Schönherr
committed
para->getParD(level)->distributionsAD.f[0],
para->getParD(level+1)->distributionsAD.f[0],
para->getParD(level)->neighborX,
para->getParD(level)->neighborY,
para->getParD(level)->neighborZ,
para->getParD(level+1)->neighborX,
para->getParD(level+1)->neighborY,
para->getParD(level+1)->neighborZ,
para->getParD(level)->numberOfNodes,
para->getParD(level+1)->numberOfNodes,
para->getParD(level)->isEvenTimestep,
para->getParD(level)->intCF.coarseCellIndices,
para->getParD(level)->intCF.fineCellIndices,
para->getParD(level)->intCF.numberOfCells,
para->getParD(level)->viscosity,
para->getParD(level+1)->diffusivity,
para->getParD(level)->numberofthreads,
para->getParD(level)->neighborCF);