Refactor calculateTurbuulenceIntensity (minor changes)

src/gpu/VirtualFluids_GPU/Calculation/CalcTurbulenceIntensity.cpp

@@ -14,7 +14,8 @@ void allocTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, u
 {
     int lev = para->getFine();
 	cudaManager->cudaAllocTurbulenceIntensity(lev, size);
-    resetTurbulenceIntensity(para, cudaManager, size);
+    resetVelocityFluctuationsAndMeans(para, cudaManager, size);
+    para->getParH(lev)->turbulenceIntensity.resize(size);
 }
 
 
@@ -43,11 +44,10 @@ void calcVelocityAndFluctuations(Parameter *para, CudaMemoryManager *cudaManager
 
 
 void calcTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, uint tdiff, uint size) {
+    
     calcVelocityAndFluctuations(para, cudaManager, tdiff, size);
     int lev = para->getFine();
 
-    para->getParH(lev)->turbulenceIntensity.resize(size);
-
     real fluc_squared;
     real v_mean_squared;
     for (uint i = 0; i < size; i++) {
@@ -61,7 +61,7 @@ void calcTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, ui
 }
 
 
-void resetTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, uint size)
+void resetVelocityFluctuationsAndMeans(Parameter *para, CudaMemoryManager *cudaManager, uint size)
 {
     int lev = para->getFine();
 

src/gpu/VirtualFluids_GPU/Calculation/CalcTurbulenceIntensity.h

@@ -9,7 +9,7 @@
 extern "C" void allocTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, uint size);
 extern "C" void calcVelocityAndFluctuations(Parameter *para, CudaMemoryManager *cudaManager, uint tdiff, uint size);
 extern "C" void calcTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, uint tdiff, uint size);
-extern "C" void resetTurbulenceIntensity(Parameter *para, CudaMemoryManager *cudaManager, uint size);
+extern "C" void resetVelocityFluctuationsAndMeans(Parameter *para, CudaMemoryManager *cudaManager, uint size);
 
 void writeTurbulenceIntensityToFile(Parameter *para, uint timestep, int *vectorOfSparseIndices);
 

src/gpu/VirtualFluids_GPU/LBM/Simulation.cpp

@@ -968,7 +968,7 @@ void Simulation::run()
                 calcTurbulenceIntensity(para.get(), cudaManager.get(), t_diff, para->getParH(para->getFine())->QGeom.kQ);
 				writeTurbulenceIntensityToFile(para.get(), t, para->getParH(para->getFine())->QGeom.k);
 				t_turbulenceIntensity = t;
-                resetTurbulenceIntensity(para.get(), cudaManager.get(), para->getParH(para->getFine())->QGeom.kQ);
+                resetVelocityFluctuationsAndMeans(para.get(), cudaManager.get(), para->getParH(para->getFine())->QGeom.kQ);
             }
 			////////////////////////////////////////////////////////////////////////
 			dataWriter->writeTimestep(para, t);