Skip to content
Snippets Groups Projects
Commit 6d7cc40b authored by Henrik Asmuth's avatar Henrik Asmuth
Browse files

add CalcMacComp27 call to PlanarAverageProbe

parent 124d4b59
No related branches found
No related tags found
1 merge request!170Kernel templetization and efficiency improvements
......@@ -15,6 +15,7 @@
#include "Parameter/Parameter.h"
#include "DataStructureInitializer/GridProvider.h"
#include "GPU/CudaMemoryManager.h"
#include "GPU/GPU_Interface.h"
#include <algorithm>
......@@ -268,6 +269,15 @@ void PlanarAverageProbe::findPoints(Parameter* para, GridProvider* gridProvider,
void PlanarAverageProbe::calculateQuantities(SPtr<ProbeStruct> probeStruct, Parameter* para, uint t_level, int level)
{
// Compute macroscopic variables in entire domain
CalcMacCompSP27(para->getParD(level)->velocityX, para->getParD(level)->velocityY, para->getParD(level)->velocityZ,
para->getParD(level)->rho, para->getParD(level)->pressure, para->getParD(level)->typeOfGridNode,
para->getParD(level)->neighborX, para->getParD(level)->neighborY,
para->getParD(level)->neighborZ, para->getParD(level)->numberOfNodes,
para->getParD(level)->numberofthreads, para->getParD(level)->distributions.f[0],
para->getParD(level)->isEvenTimestep);
getLastCudaError("In PlanarAverageProbe Kernel CalcMacSP27 execution failed");
// Definition of normal and inplane directions for moveIndices kernels
uint *neighborNormal, *neighborInplane1, *neighborInplane2;
if( this->planeNormal == 'x' )
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment