From 624f5d8181a449f36950503859271e1a052e8a65 Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Tue, 27 Jan 2026 09:17:33 +0100
Subject: Double free, but constexpr things compile now

---
 lib/core/c++/macroscopic.hpp | 25 +++++++++++++++++++++++++
 1 file changed, 25 insertions(+)

(limited to 'lib/core/c++/macroscopic.hpp')

diff --git a/lib/core/c++/macroscopic.hpp b/lib/core/c++/macroscopic.hpp
index 51368e9..1532873 100644
--- a/lib/core/c++/macroscopic.hpp
+++ b/lib/core/c++/macroscopic.hpp
@@ -4,6 +4,31 @@
 
 namespace kel {
 namespace lbm {
+template<typename T, typename Desc>
+void compute_rho_u (
+		const saw::data<sch::FixedArray<T,Desc::Q>>& dfs,
+		saw::data<T>& rho,
+		saw::data<sch::FixedArray<T,Desc::D>>& vel
+	)
+{
+	using dfi = df_info<T, Desc>;
+
+	rho = 0;
+	for(size_t j = 0; j < Desc::D; ++j){
+		vel.at({{j}}).set(0);
+	}
+
+	for(size_t j = 0; j < Desc::Q; ++j){
+		rho = rho + dfs.at({j});
+		for(size_t i = 0; i < Desc::D; ++i){
+			vel.at({{i}}) = vel.at({{i}}) + dfi::directions[j][i] * dfs.at({j}).get();
+		}
+	}
+
+	for(size_t i = 0; i < Desc::D; ++i){
+		vel.at({i}) = vel.at({i}) / rho;
+	}
+}
 /**
  * Calculate the macroscopic variables rho and u in Lattice Units.
  */
-- 
cgit v1.2.3