PropagationD3Q19.cpp

Go to the documentation of this file.

// ==============================================================================================
// Microflow 3D, http://www.microflow.pwr.edu.pl/
// Created by Roman Szafran on 01.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 "PropagationD3Q19.h"

void MF::GU::PropagationD3Q19::PropagationTest() {
    unsigned char k;
    openvdb::Coord xyz, XYZ;
    MF::GU::DirectionsD3Q19 value;
    openvdb::Int32Grid::Accessor accessor = m_VDBPropagationGrid_Ptr->getAccessor();

    // Iterate over all active values.
    for (openvdb::Int32Grid::ValueOnIter iter = m_VDBPropagationGrid_Ptr->beginValueOn(); iter.test(); ++iter) {
        xyz = iter.getCoord();
        value.node_t.FF = 0;
        for (k = 0; k < MFQ27; k++) {
            XYZ.z() = xyz.z() + m_LatticeParameters_Ptr->DirectionVectorComponent_Cz[k];
            XYZ.y() = xyz.y() + m_LatticeParameters_Ptr->DirectionVectorComponent_Cy[k];
            XYZ.x() = xyz.x() + m_LatticeParameters_Ptr->DirectionVectorComponent_Cx[k];
            if (accessor.isValueOn(XYZ)) {
                switch (k) {
                    case 0: {
                        value.node_t.ff.f0 = 1;
                        break;
                    }
                    case 1: {
                        value.node_t.ff.f1 = 1;
                        break;
                    }
                    case 2: {
                        value.node_t.ff.f2 = 1;
                        break;
                    }
                    case 3: {
                        value.node_t.ff.f3 = 1;
                        break;
                    }
                    case 4: {
                        value.node_t.ff.f4 = 1;
                        break;
                    }
                    case 5: {
                        value.node_t.ff.f5 = 1;
                        break;
                    }
                    case 6: {
                        value.node_t.ff.f6 = 1;
                        break;
                    }
                    case 7: {
                        value.node_t.ff.f7 = 1;
                        break;
                    }
                    case 8: {
                        value.node_t.ff.f8 = 1;
                        break;
                    }
                    case 9: {
                        value.node_t.ff.f9 = 1;
                        break;
                    }
                    case 10: {
                        value.node_t.ff.f10 = 1;
                        break;
                    }
                    case 11: {
                        value.node_t.ff.f11 = 1;
                        break;
                    }
                    case 12: {
                        value.node_t.ff.f12 = 1;
                        break;
                    }
                    case 13: {
                        value.node_t.ff.f13 = 1;
                        break;
                    }
                    case 14: {
                        value.node_t.ff.f14 = 1;
                        break;
                    }
                    case 15: {
                        value.node_t.ff.f15 = 1;
                        break;
                    }
                    case 16: {
                        value.node_t.ff.f16 = 1;
                        break;
                    }
                    case 17: {
                        value.node_t.ff.f17 = 1;
                        break;
                    }
                    case 18: {
                        value.node_t.ff.f18 = 1;
                        break;
                    }
                    case 19: {
                        value.node_t.ff.f19 = 1;
                        break;
                    }
                    case 20: {
                        value.node_t.ff.f20 = 1;
                        break;
                    }
                    case 21: {
                        value.node_t.ff.f21 = 1;
                        break;
                    }
                    case 22: {
                        value.node_t.ff.f22 = 1;
                        break;
                    }
                    case 23: {
                        value.node_t.ff.f23 = 1;
                        break;
                    }
                    case 24: {
                        value.node_t.ff.f24 = 1;
                        break;
                    }
                    case 25: {
                        value.node_t.ff.f25 = 1;
                        break;
                    }
                    case 26: {
                        value.node_t.ff.f26 = 1;
                        break;
                    }
                    default: {
                        std::cout << "Error in propagationD3Q19 test" << std::endl;
                        exit(EXIT_FAILURE);
                    }
                }
            }
        }
        accessor.setValue(xyz, value.node_t.FF);
    }
}