Changes to run, reworking meta reporting for easier estimation of

runtime parameters

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
 	;
 }
 }