From 50a77ad35c1d21fb789039beff386811d05447ab Mon Sep 17 00:00:00 2001
From: Henry <henry.korb@geo.uu.se>
Date: Thu, 15 Jun 2023 18:34:41 +0200
Subject: [PATCH] remove old NumberOfStressBCNodes hack
---
.../GridReaderGenerator/GridGenerator.cpp | 1 -
.../PreCollisionInteractor/Probes/WallModelProbe.cu | 2 +-
2 files changed, 1 insertion(+), 2 deletions(-)
diff --git a/src/gpu/VirtualFluids_GPU/DataStructureInitializer/GridReaderGenerator/GridGenerator.cpp b/src/gpu/VirtualFluids_GPU/DataStructureInitializer/GridReaderGenerator/GridGenerator.cpp
index fee1bce34..e3c86317c 100644
--- a/src/gpu/VirtualFluids_GPU/DataStructureInitializer/GridReaderGenerator/GridGenerator.cpp
+++ b/src/gpu/VirtualFluids_GPU/DataStructureInitializer/GridReaderGenerator/GridGenerator.cpp
@@ -263,7 +263,6 @@ void GridGenerator::allocArrays_BoundaryValues()
cudaMemoryManager->cudaCopyWallModel(level, para->getHasWallModelMonitor());
}
para->getParD(level)->stressBC.numberOfBCnodes = para->getParH(level)->stressBC.numberOfBCnodes;
- para->getParD(level)->stressBC.kArray = para->getParH(level)->stressBC.kArray;
para->getParH(level)->numberOfStressBCnodesRead = para->getParH(level)->stressBC.numberOfBCnodes * para->getD3Qxx();
para->getParD(level)->numberOfStressBCnodesRead = para->getParH(level)->stressBC.numberOfBCnodes * para->getD3Qxx();
}
diff --git a/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/WallModelProbe.cu b/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/WallModelProbe.cu
index b2a5b9122..f52c666c9 100644
--- a/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/WallModelProbe.cu
+++ b/src/gpu/VirtualFluids_GPU/PreCollisionInteractor/Probes/WallModelProbe.cu
@@ -172,7 +172,7 @@ void temporal_average(T* quantitiesArray, T currentValue, uint currentTimestep,
void WallModelProbe::calculateQuantities(SPtr<ProbeStruct> probeStruct, Parameter* para, uint t, int level)
{
bool doTmpAveraging = (t>this->getTStartTmpAveraging());
- uint numberOfStressBCPoints = para->getParD(level)->stressBC.kArray;
+ uint numberOfStressBCPoints = para->getParD(level)->stressBC.numberOfBCnodes;
if(numberOfStressBCPoints<1) return; //Skipping levels without StressBC
uint timestep = probeStruct->timestepInTimeseries;
real inv_n = c1o1/real(probeStruct->timestepInTimeAverage+1);
--
GitLab