FindGhostLayerCandidates.hpp

Go to the documentation of this file.

// -------------------------------------------------------//
//
// SHAMROCK code for hydrodynamics
// Copyright (c) 2021-2026 Timothée David--Cléris <tim.shamrock@proton.me>
// SPDX-License-Identifier: CeCILL Free Software License Agreement v2.1
// Shamrock is licensed under the CeCILL 2.1 License, see LICENSE for more information
//
// -------------------------------------------------------//
 
#pragma once

 
#include "shammath/AABB.hpp"
#include "shammath/paving_function.hpp"
#include "shamrock/solvergraph/DDSharedScalar.hpp"
#include "shamrock/solvergraph/IFieldRefs.hpp"
#include "shamrock/solvergraph/Indexes.hpp"
#include "shamrock/solvergraph/PatchDataLayerDDShared.hpp"
#include "shamrock/solvergraph/PatchDataLayerRefs.hpp"
#include "shamrock/solvergraph/ScalarEdge.hpp"
#include "shamrock/solvergraph/ScalarsEdge.hpp"
#include "shamrock/solvergraph/SerialPatchTreeEdge.hpp"
#include "shamsolvergraph/edge/IDataEdge.hpp"
#include "shamsolvergraph/node/INode.hpp"
 
namespace shammodels::basegodunov::modules {
 
    enum class GhostType { None, Periodic, Reflective };
 
    struct GhostLayerGenMode {
        GhostType ghost_type_x;
        GhostType ghost_type_y;
        GhostType ghost_type_z;
    };
 
    template<class TgridVec>
    shammath::paving_function_general_3d<TgridVec> get_paving(
        GhostLayerGenMode mode, shammath::AABB<TgridVec> sim_box) {
 
        TgridVec box_size   = sim_box.upper - sim_box.lower;
        TgridVec box_center = (sim_box.upper + sim_box.lower) / 2;
 
        SHAM_ASSERT(sim_box.is_volume_not_null());
 
        { // check that rebuildind the AABB from size and center gives the same AABB
            shammath::AABB<TgridVec> new_box
                = {box_center - box_size / 2, box_center + box_size / 2};
            if (new_box != sim_box) {
                shambase::throw_with_loc<std::invalid_argument>(
                    "Rebuilding AABB from size and center gives a different AABB");
            }
        }
 
        return shammath::paving_function_general_3d<TgridVec>{
            box_size,
            box_center,
            mode.ghost_type_x == GhostType::Periodic,
            mode.ghost_type_y == GhostType::Periodic,
            mode.ghost_type_z == GhostType::Periodic};
    }
 
    template<class Func>
    void for_each_paving_tile(GhostLayerGenMode mode, Func &&func) {
 
        // if the ghost type is none, we do not need to repeat as there is no ghost layer
        i32 repetition_x = mode.ghost_type_x != GhostType::None;
        i32 repetition_y = mode.ghost_type_y != GhostType::None;
        i32 repetition_z = mode.ghost_type_z != GhostType::None;
 
        for (i32 xoff = -repetition_x; xoff <= repetition_x; xoff++) {
            for (i32 yoff = -repetition_y; yoff <= repetition_y; yoff++) {
                for (i32 zoff = -repetition_z; zoff <= repetition_z; zoff++) {
                    func(xoff, yoff, zoff);
                }
            }
        }
    }
 
    struct GhostLayerCandidateInfos {
        i32 xoff;
        i32 yoff;
        i32 zoff;
    };
 
#define NODE_EDGES(X_RO, X_RW)                                                                     \
    /* ------------------- inputs ------------------- */                                           \
    X_RO(shamrock::solvergraph::IDataEdge<std::vector<u64>>, ids_to_check)                         \
    X_RO(shamrock::solvergraph::ScalarEdge<shammath::AABB<TgridVec>>, sim_box)                     \
    X_RO(shamrock::solvergraph::SerialPatchTreeRefEdge<TgridVec>, patch_tree)                      \
    X_RO(shamrock::solvergraph::ScalarsEdge<shammath::AABB<TgridVec>>, patch_boxes)                \
                                                                                                   \
    /* ------------------- outputs ------------------- */                                          \
    X_RW(shamrock::solvergraph::DDSharedScalar<GhostLayerCandidateInfos>, ghost_layers_candidates)
 
    template<class TgridVec>
    class FindGhostLayerCandidates : public shamrock::solvergraph::INode {
 
        GhostLayerGenMode mode;
 
        public:
        FindGhostLayerCandidates(GhostLayerGenMode mode) : mode(mode) {}
 
        EXPAND_NODE_EDGES(NODE_EDGES)
 
        void _impl_evaluate_internal();
 
        inline virtual std::string _impl_get_label() const { return "FindGhostLayerCandidates"; };
 
        virtual std::string _impl_get_tex() const;
    };
} // namespace shammodels::basegodunov::modules
 
#undef NODE_EDGES