reconstruction.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/math.hpp"

 
namespace shammodels::gsph {

    template<class Tscal>
    struct EffectiveFace {
        Tscal V2;    
        Tscal s_ast; 
    };

    template<class Tscal>
    inline EffectiveFace<Tscal> lin_v2_sast_ij(
        Tscal vol_a, Tscal vol_b, Tscal h_a, Tscal h_b, Tscal rab_inv) {
 
        const Tscal C  = (vol_a - vol_b) * rab_inv;
        const Tscal D  = Tscal{0.5} * (vol_a + vol_b);
        const Tscal h2 = Tscal{0.5} * (h_a * h_a + h_b * h_b);
 
        const Tscal V2    = Tscal{0.25} * h2 * C * C + D * D;
        const Tscal s_ast = Tscal{0.5} * h2 * C * D * sham::inv_sat_zero(V2);
 
        return EffectiveFace<Tscal>{V2, s_ast};
    }
 
} // namespace shammodels::gsph