Skip to content
Snippets Groups Projects
Commit d83ead52 authored by Hkorb's avatar Hkorb
Browse files

bugfix in probe: wrong number of threads in visit

parent a7cfa0cb
No related branches found
No related tags found
1 merge request!81fixed AL
...@@ -181,7 +181,7 @@ void ActuatorLine::visit(Parameter* para, CudaMemoryManager* cudaManager, int le ...@@ -181,7 +181,7 @@ void ActuatorLine::visit(Parameter* para, CudaMemoryManager* cudaManager, int le
cudaManager->cudaCopyBladeCoordsHtoD(this); cudaManager->cudaCopyBladeCoordsHtoD(this);
unsigned int numberOfThreads = 128; uint numberOfThreads = 128;
vf::gpu::CudaGrid bladeGrid = vf::gpu::CudaGrid(numberOfThreads, this->numberOfNodes); vf::gpu::CudaGrid bladeGrid = vf::gpu::CudaGrid(numberOfThreads, this->numberOfNodes);
interpolateVelocities<<< bladeGrid.grid, bladeGrid.threads >>>( interpolateVelocities<<< bladeGrid.grid, bladeGrid.threads >>>(
......
...@@ -109,7 +109,7 @@ void Probe::init(Parameter* para, GridProvider* gridProvider, CudaMemoryManager* ...@@ -109,7 +109,7 @@ void Probe::init(Parameter* para, GridProvider* gridProvider, CudaMemoryManager*
std::vector<real> pointCoordsY_level; std::vector<real> pointCoordsY_level;
std::vector<real> pointCoordsZ_level; std::vector<real> pointCoordsZ_level;
real dx = abs(para->getParH(level)->coordX_SP[1]-para->getParH(level)->coordX_SP[para->getParH(level)->neighborX_SP[1]]); real dx = abs(para->getParH(level)->coordX_SP[1]-para->getParH(level)->coordX_SP[para->getParH(level)->neighborX_SP[1]]);
for(uint j=0; j<para->getParH(level)->size_Mat_SP; j++ ) for(uint j=1; j<para->getParH(level)->size_Mat_SP; j++ )
{ {
for(uint point=0; point<this->nProbePoints; point++) for(uint point=0; point<this->nProbePoints; point++)
{ {
...@@ -190,7 +190,7 @@ void Probe::visit(Parameter* para, CudaMemoryManager* cudaManager, int level, un ...@@ -190,7 +190,7 @@ void Probe::visit(Parameter* para, CudaMemoryManager* cudaManager, int level, un
{ {
ProbeStruct* probeStruct = this->getProbeStruct(level); ProbeStruct* probeStruct = this->getProbeStruct(level);
vf::gpu::CudaGrid grid = vf::gpu::CudaGrid(probeStruct->nPoints, 128); vf::gpu::CudaGrid grid = vf::gpu::CudaGrid(128, probeStruct->nPoints);
interpQuantities<<<grid.grid, grid.threads>>>( probeStruct->pointIndicesD, probeStruct->nPoints, interpQuantities<<<grid.grid, grid.threads>>>( probeStruct->pointIndicesD, probeStruct->nPoints,
probeStruct->distXD, probeStruct->distYD, probeStruct->distZD, probeStruct->distXD, probeStruct->distYD, probeStruct->distZD,
...@@ -202,8 +202,6 @@ void Probe::visit(Parameter* para, CudaMemoryManager* cudaManager, int level, un ...@@ -202,8 +202,6 @@ void Probe::visit(Parameter* para, CudaMemoryManager* cudaManager, int level, un
cudaManager->cudaCopyProbeQuantityArrayDtoH(this, level); cudaManager->cudaCopyProbeQuantityArrayDtoH(this, level);
this->write(para, level, t); this->write(para, level, t);
} }
} }
void Probe::free(Parameter* para, CudaMemoryManager* cudaManager) void Probe::free(Parameter* para, CudaMemoryManager* cudaManager)
......
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