diff options
Diffstat (limited to 'lib/core/c++/lbm.hpp')
| -rw-r--r-- | lib/core/c++/lbm.hpp | 24 |
1 files changed, 21 insertions, 3 deletions
diff --git a/lib/core/c++/lbm.hpp b/lib/core/c++/lbm.hpp index 27a7df3..a1e088e 100644 --- a/lib/core/c++/lbm.hpp +++ b/lib/core/c++/lbm.hpp @@ -26,15 +26,33 @@ namespace kel { namespace lbm { template<typename T, typename Desc> -void print_lbm_meta(const converter<T>& conv, const saw::data<sch::SiKinematicViscosity<T>>& kin_vis_si){ +void print_lbm_meta( + const converter<T>& conv, + const saw::data<sch::SiKinematicViscosity<T>>& kin_vis_si, + const saw::data<sch::SiVelocity<T>>& char_vel, + const saw::data<sch::SiMeter<T>>& char_len +){ std::cout - <<"[LBM Meta]\n" - <<"==========\n" + <<"[Meta]\n" + <<"======\n" + <<"Re: "<<(char_vel * char_len / kin_vis_si)<<"\n" + <<"Ma: "<<(char_vel * saw::data<typename saw::unit_division<sch::Pure<T>, sch::SiVelocity<T>>::Schema>{std::sqrt(df_info<T,Desc>::inv_cs2)})<<"\n" <<"\n" + <<"[SI]\n" + <<"====\n" <<"Δx: "<<conv.delta_x()<<"\n" <<"Δt: "<<conv.delta_t()<<"\n" + <<"Δv: "<<conv.delta_v()<<"\n" + <<"Δa: "<<conv.delta_a()<<"\n" <<"KinVis: "<<kin_vis_si<<"\n" + <<"CharV: "<<char_vel<<"\n" + <<"CharL: "<<char_len<<"\n" + <<"\n" + <<"[LBM]\n" + <<"=====\n" + <<"KinVis: "<<conv.kinematic_viscosity_si_to_lbm(kin_vis_si)<<"\n" <<"τ: "<<(saw::data<typename saw::unit_division<sch::Pure<T>, sch::LbmKinematicViscosity<T>>::Schema >{df_info<T,Desc>::inv_cs2} * conv.kinematic_viscosity_si_to_lbm(kin_vis_si) + saw::data<sch::Pure<T>>{0.5})<<"\n" + <<std::endl ; } } |