LoopSmoothingLengthIter.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 "shambackends/vec.hpp"
#include "shamrock/solvergraph/IFieldRefs.hpp"
#include "shamrock/solvergraph/INode.hpp"
#include "shamrock/solvergraph/ScalarEdge.hpp"
#include <memory>
 
namespace shammodels::sph::modules {
 
    template<class Tvec>
    class LoopSmoothingLengthIter : public shamrock::solvergraph::INode {
 
        std::shared_ptr<INode> iterate_smth_length_once_ptr;
 
        using Tscal = shambase::VecComponent<Tvec>;
 
        Tscal epsilon_h;
        u32 h_iter_per_subcycles;
        bool print_info;
 
        public:
        LoopSmoothingLengthIter(
            std::shared_ptr<INode> iterate_smth_length_once_ptr,
            Tscal epsilon_h,
            u32 h_iter_per_subcycles,
            bool print_info)
            : iterate_smth_length_once_ptr(std::move(iterate_smth_length_once_ptr)),
              epsilon_h(epsilon_h), h_iter_per_subcycles(h_iter_per_subcycles),
              print_info(print_info) {}
 
        struct Edges {
            const shamrock::solvergraph::IFieldRefs<Tscal> &eps_h;
            shamrock::solvergraph::ScalarEdge<bool> &is_converged;
        };
 
        inline void set_edges(
            std::shared_ptr<shamrock::solvergraph::IFieldRefs<Tscal>> eps_h,
            std::shared_ptr<shamrock::solvergraph::ScalarEdge<bool>> is_converged) {
            __internal_set_ro_edges({eps_h});
            __internal_set_rw_edges({is_converged});
        }
 
        inline Edges get_edges() {
            return Edges{
                get_ro_edge<shamrock::solvergraph::IFieldRefs<Tscal>>(0),
                get_rw_edge<shamrock::solvergraph::ScalarEdge<bool>>(0),
            };
        }
 
        void _impl_evaluate_internal();
 
        inline virtual std::string _impl_get_label() const { return "LoopSmoothingLengthIter"; };
 
        virtual std::string _impl_get_tex() const;
    };
} // namespace shammodels::sph::modules
 
template class shammodels::sph::modules::LoopSmoothingLengthIter<f64_3>;