From 5be771eb356263bf43c57056eba77df7235c78f4 Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Fri, 27 Jun 2025 21:27:20 +0200 Subject: Fixing things --- examples/cavity_2d.cpp | 26 -------------------------- 1 file changed, 26 deletions(-) (limited to 'examples/cavity_2d.cpp') diff --git a/examples/cavity_2d.cpp b/examples/cavity_2d.cpp index 6ecdb23..6332bcf 100644 --- a/examples/cavity_2d.cpp +++ b/examples/cavity_2d.cpp @@ -52,33 +52,7 @@ using MacroStruct = Struct< >; using CavityFieldD2Q9 = CellField>; -} - -/** - * Calculates the equilibrium for each direction - */ -template -std::array::type,Desc::Q> equilibrium( - typename saw::native_data_type::type rho, - const std::array::type, Desc::D>& vel -){ - using dfi = df_info; - - typename std::array::type,Desc::Q> eq; - - for(std::size_t i = 0u; i < eq.size(); ++i){ - auto vel_c = (vel[0u]*dfi::directions[i][0u] + vel[1u]*dfi::directions[i][1u]); - auto vel_c_cs2 = vel_c * dfi::inv_cs2; - eq[i] = dfi::weights[i] * rho * ( - 1.0 - + vel_c_cs2 - - dfi::inv_cs2 * 0.5 * ( vel[0u] * vel[0u] + vel[1u] * vel[1u] ) - + vel_c_cs2 * vel_c_cs2 * 0.5 - ); - } - return eq; -} /* template -- cgit v1.2.3