LatticeParametersD3Q19.h

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.
// ==============================================================================================

# pragma once

#include <memory>

#define MFQ19 (19)                    
#define MFQ27 (27)                    
#define MFQ19_H  ((MFQ19 - 1) / 2)    

namespace MF {
    namespace GU {

        class LatticeParametersD3Q19 {
        public:
            LatticeParametersD3Q19() = default;

            ~LatticeParametersD3Q19() = default;

            static std::shared_ptr<MF::GU::LatticeParametersD3Q19> New() {
                auto LatticeParametersD3Q19_Ptr = std::make_shared<MF::GU::LatticeParametersD3Q19>();
                return LatticeParametersD3Q19_Ptr;
            };

            static constexpr double LatticeConstant_CSQ = 1.0 / 3.0;            

            static constexpr int8_t DirectionVectorComponent_Cx[MFQ27]
                 //f = 0, 1, 2,  3, 4, 5, 6, 7,  8,  9,10,11,12, 13, 14,15,16,17,18,|19,20,21, 22, 23,24,25, 26
                    = {0, 1, 0, -1, 0, 0, 0, 1, -1, -1, 1, 1, 1, -1, -1, 0, 0, 0, 0, -1, 1, 1, -1, -1, 1, 1, -1};

            static constexpr int8_t DirectionVectorComponent_Cy[MFQ27]
                 //f = 0, 1, 2, 3,  4, 5, 6, 7, 8,  9, 10,11,12,13,14,15,16, 17, 18,|19,20,21,22, 23, 24, 25, 26
                    = {0, 0, 1, 0, -1, 0, 0, 1, 1, -1, -1, 0, 0, 0, 0, 1, 1, -1, -1, 1, 1, 1, 1, -1, -1, -1, -1};

            static constexpr int8_t DirectionVectorComponent_Cz[MFQ27]
                 //f = 0, 1, 2, 3, 4, 5,  6, 7, 8, 9,10,11, 12, 13,14,15, 16, 17,18,|19, 20,21,22, 23, 24,25,26
                    = {0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 1, -1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 1, -1, -1, 1, 1};

            static constexpr uint8_t ReverseDirectionVector_Rc[MFQ27]
                    //f = 0, 1, 2, 3, 4, 5, 6, 7, 8, 9,10,11,12,13,14,15,16,17,18,|19,20,21,22,23,24,25,26
                    = {0, 3, 4, 1, 2, 6, 5, 9, 10, 7, 8, 13, 14, 11, 12, 17, 18, 15, 16, 25, 26, 23, 24, 21, 22, 19, 20};

            static constexpr uint8_t SwapDirections[MFQ19_H][2] = {{1,3},{2,4},{5,6},{7,9},{8,10},{11,13},{12,14},{15,17},{16,18}};

            static constexpr double DirectionWeight_W[MFQ19]
                    //f = 0      ,1       ,2       ,3       ,4       ,5       ,6       ,7       ,8       ,9       ,10      ,11      ,12      ,13      ,14      ,15      ,16      ,17      ,18
                    = {1.0 / 3.0, 1.0 / 18.0, 1.0 / 18.0, 1.0 / 18.0, 1.0 / 18.0, 1.0 / 18.0, 1.0 / 18.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0,
                       1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0, 1.0 / 36.0};

            static uint8_t find_f(const char nz, const char ny, const char nx) {
                unsigned char k;

                for (k = 0; k < MFQ27; k++)
                    if (DirectionVectorComponent_Cz[k] == nz && DirectionVectorComponent_Cy[k] == ny && DirectionVectorComponent_Cx[k] == nx)
                        return (k);
            return (0);
            };

        };

    } /* namespace GU */
} /* namespase MF */