BoundaryFromThread.cpp

Go to the documentation of this file.

// ==============================================================================================
// Microflow 3D, http://www.microflow.pwr.edu.pl/
// Created by Roman Szafran on 18.05.19.
// Copyright (c) 2019 Wroclaw University of Science and Technology.
// Distributed under the Apache License, Version 2.0. You may obtain a copy of the License at
// http://www.apache.org/licenses/LICENSE-2.0 or see accompanying file license.txt.
// Redistributions of source code must retain the above copyright and license notice.
// ==============================================================================================

#include "BoundaryFromThread.h"
#include "NodeID.h"

void MF::GB::BoundaryFromThread::VDBBoundaryAdd(const openvdb::Int64Grid::Ptr& VDBGrid, const std::shared_ptr<MF::Database::ConfigData> & ConfigData_Ptr) {
    MF::GB::NodeID nodeTMP, PresentNodeID;
    nodeTMP.node_id.NodeID = {0};
    for (unsigned int iter = 0; iter < ConfigData_Ptr->getThreadsNumber(); iter++) {
        if (ConfigData_Ptr->getThreadStringParam(iter, "GridType") == "VDB") {
            if (ConfigData_Ptr->getThreadParamLength(iter, "NodeFilePath") > 0) {
                if (ConfigData_Ptr->getThreadIntParam(iter, "NodeTypeID") >= 20) { //process only threads with definition of boundary nodes
                    auto nodeID = (int16_t) ConfigData_Ptr->getThreadIntParam(iter, "NodeTypeID");
                    std::string FilePath =
                            ConfigData_Ptr->CaseFolder + ConfigData_Ptr->getProgramStringParam("GeometryFolder") +
                            ConfigData_Ptr->getThreadStringParam(iter, "NodeFilePath");
                    auto data_Ptr = std::make_unique<std::vector<double>>();
                    std::vector<std::string> VariableNameList{
                            ConfigData_Ptr->getThreadStringVectorParam(0, iter, "NodeCoordinateFileColumnNames"),
                            ConfigData_Ptr->getThreadStringVectorParam(1, iter, "NodeCoordinateFileColumnNames"),
                            ConfigData_Ptr->getThreadStringVectorParam(2, iter, "NodeCoordinateFileColumnNames")};
                    MF::RW::CSVReader::DataReadFromCSVFile(FilePath, VariableNameList, data_Ptr.get());
                    auto accessor = VDBGrid->getAccessor();
                    openvdb::Coord xyz;
                    int32_t &i = xyz[0], &j = xyz[1], &h = xyz[2];

                    for (auto it = data_Ptr->begin(); it < data_Ptr->end(); it += 3) {
                        i = int32_t(*it);
                        j = int32_t(*(it + 1));
                        h = int32_t(*(it + 2));
                        PresentNodeID.node_id.NodeID = accessor.getValue(xyz);
                        if (PresentNodeID.node_id.Node.nodeType == 3) { // only changes node type for boundary nodes
                            nodeTMP.node_id.Node.nodeType = nodeID;
                            nodeTMP.node_id.Node.uidThreadNr = (uint16_t) iter + 1; // uidThreadNr added from the thread_params.cfg starts from 1 because automatically created threads have uidThreadNr = 0;
                            accessor.setValue(xyz, nodeTMP.node_id.NodeID);
                        }
                    }
                    data_Ptr.reset();
                }
            }
        }
    }
}