Revision 35efca7662eafe5c0648fbb4e97a7bc1eda84322 authored by Christoph Neuhauser on 23 October 2023, 22:10:04 UTC, committed by Christoph Neuhauser on 23 October 2023, 22:10:04 UTC
1 parent ef25cab
RbcBinFileLoader.cpp
/*
* BSD 2-Clause License
*
* Copyright (c) 2022, Christoph Neuhauser
* All rights reserved.
*
* Redistribution and use in source and binary forms, with or without
* modification, are permitted provided that the following conditions are met:
*
* * Redistributions of source code must retain the above copyright notice, this
* list of conditions and the following disclaimer.
*
* * Redistributions in binary form must reproduce the above copyright notice,
* this list of conditions and the following disclaimer in the documentation
* and/or other materials provided with the distribution.
*
* THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
* AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
* IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
* DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
* FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
* DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
* SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
* CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
* OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
* OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
*/
#include <Utils/File/Logfile.hpp>
#include <Utils/File/FileLoader.hpp>
#include "Volume/VolumeData.hpp"
#include "DataSet.hpp"
#include "LoadersUtil.hpp"
#include "RbcBinFileLoader.hpp"
bool RbcBinFileLoader::setInputFiles(
VolumeData* volumeData, const std::string& filePath, const DataSetInformation& _dataSetInformation) {
dataSourceFilename = filePath;
dataSetInformation = _dataSetInformation;
xs = 1024;
ys = 32;
zs = 1024;
cellStep = 1.0f / 1023.0f;
volumeData->setGridExtent(xs, ys, zs, cellStep, cellStep, cellStep);
std::unordered_map<FieldType, std::vector<std::string>> fieldNameMap;
fieldNameMap[FieldType::SCALAR].push_back("Temperature");
fieldNameMap[FieldType::VECTOR].push_back("Velocity");
fieldNameMap[FieldType::VECTOR].push_back("Vorticity");
fieldNameMap[FieldType::SCALAR].push_back("Velocity Magnitude");
fieldNameMap[FieldType::SCALAR].push_back("Vorticity Magnitude");
fieldNameMap[FieldType::SCALAR].push_back("Helicity");
volumeData->setFieldNames(fieldNameMap);
return true;
}
bool RbcBinFileLoader::getFieldEntry(
VolumeData* volumeData, FieldType fieldType, const std::string& fieldName,
int timestepIdx, int memberIdx, HostCacheEntryType*& fieldEntry) {
uint8_t* buffer = nullptr;
size_t length = 0;
bool loaded = sgl::loadFileFromSource(dataSourceFilename, buffer, length, false);
if (!loaded) {
sgl::Logfile::get()->throwError(
"Error in RbcBinFileLoader::load: Couldn't open file \"" + dataSourceFilename + "\".");
}
auto* dataField = reinterpret_cast<float*>(buffer);
if (length != size_t(xs) * size_t(ys) * size_t(zs) * sizeof(float) * 4) {
sgl::Logfile::get()->throwError(
"Error in RbcBinFileLoader::load: Inconsistent number of bytes read from file \""
+ dataSourceFilename + "\".");
}
int vectorFieldNumEntries = xs * ys * zs * 3;
int scalarFieldNumEntries = xs * ys * zs;
auto* velocityField = new float[vectorFieldNumEntries];
auto* velocityMagnitudeField = new float[scalarFieldNumEntries];
auto* vorticityField = new float[vectorFieldNumEntries];
auto* vorticityMagnitudeField = new float[scalarFieldNumEntries];
auto* helicityField = new float[scalarFieldNumEntries];
auto* temperatureField = new float[scalarFieldNumEntries];
for (int z = 0; z < zs; z++) {
for (int y = 0; y < ys; y++) {
for (int x = 0; x < xs; x++) {
velocityField[IDXV(x, y, z, 0)] = dataField[IDXV4(x, y, z, 0)];
velocityField[IDXV(x, y, z, 1)] = dataField[IDXV4(x, y, z, 1)];
velocityField[IDXV(x, y, z, 2)] = dataField[IDXV4(x, y, z, 2)];
temperatureField[IDXS(x, y, z)] = dataField[IDXV4(x, y, z, 3)];
}
}
}
computeVectorMagnitudeField(velocityField, velocityMagnitudeField, xs, ys, zs);
computeVorticityField(velocityField, vorticityField, xs, ys, zs, cellStep, cellStep, cellStep);
computeVectorMagnitudeField(vorticityField, vorticityMagnitudeField, xs, ys, zs);
computeHelicityFieldNormalized(
velocityField, vorticityField, helicityField, xs, ys, zs,
dataSetInformation.useNormalizedVelocity,
dataSetInformation.useNormalizedVorticity);
volumeData->addField(velocityField, FieldType::VECTOR, "Velocity", timestepIdx, memberIdx);
volumeData->addField(vorticityField, FieldType::VECTOR, "Vorticity", timestepIdx, memberIdx);
volumeData->addField(helicityField, FieldType::SCALAR, "Helicity", timestepIdx, memberIdx);
volumeData->addField(velocityMagnitudeField, FieldType::SCALAR, "Velocity Magnitude", timestepIdx, memberIdx);
volumeData->addField(vorticityMagnitudeField, FieldType::SCALAR, "Vorticity Magnitude", timestepIdx, memberIdx);
volumeData->addField(temperatureField, FieldType::SCALAR, "Temperature", timestepIdx, memberIdx);
delete[] buffer;
return true;
}
![swh spinner](/static/img/swh-spinner.gif)
Computing file changes ...