diff --git a/src/gpu/VirtualFluids_GPU/Output/WriterUtilities.cpp b/src/gpu/VirtualFluids_GPU/Output/WriterUtilities.cpp index 688e55a4545251246a8190131b8dda1daae0949a..50ad3c08ad3289631ccc5f2f279cf9aaa36e1190 100644 --- a/src/gpu/VirtualFluids_GPU/Output/WriterUtilities.cpp +++ b/src/gpu/VirtualFluids_GPU/Output/WriterUtilities.cpp @@ -45,7 +45,7 @@ void WriterUtilities::getIndicesOfAllNodesInOct(std::array<uint, 8>& indices, ui void WriterUtilities::calculateRelativePositionInPart(std::array<uint, 8>& relativePositionInPart, const std::array<uint, 8>& indicesOfOct, uint startPosition) { - for (size_t i = 1; i < relativePositionInPart.size(); i++) { + for (size_t i = 0; i < relativePositionInPart.size(); i++) { relativePositionInPart[i] = indicesOfOct[i] - startPosition; } } @@ -53,11 +53,11 @@ void WriterUtilities::calculateRelativePositionInPart(std::array<uint, 8>& relat bool WriterUtilities::areAllNodesInOctValidForWriting(const std::array<uint, 8>& indicesOfOct, const LBMSimulationParameter* parH, uint endPositionOfPart) { - bool neighborsAreFluid = std::all_of(indicesOfOct.begin() + 1, indicesOfOct.end(), + bool neighborsAreFluid = std::all_of(indicesOfOct.begin(), indicesOfOct.end(), [&](uint index) { return parH->typeOfGridNode[index] == GEO_FLUID; }); bool neighborIsOutOfPart = - (std::any_of(indicesOfOct.begin() + 1, indicesOfOct.end(), [&](uint index) { return index > endPositionOfPart; })); + (std::any_of(indicesOfOct.begin(), indicesOfOct.end(), [&](uint index) { return index > endPositionOfPart; })); return neighborsAreFluid && !neighborIsOutOfPart; } \ No newline at end of file