From b9aacd9c2fdc61721c8ca3e1b40ebc92daa3772e Mon Sep 17 00:00:00 2001
From: "Claudius \"keldu\" Holeksa" <mail@keldu.de>
Date: Wed, 5 Nov 2025 13:38:04 +0100
Subject: Reworking structure

---
 lib/core/tests/equilibrium.cpp | 40 ++++++++++++++++++++++++++++++++++++++++
 1 file changed, 40 insertions(+)
 create mode 100644 lib/core/tests/equilibrium.cpp

(limited to 'lib/core/tests/equilibrium.cpp')

diff --git a/lib/core/tests/equilibrium.cpp b/lib/core/tests/equilibrium.cpp
new file mode 100644
index 0000000..9201e55
--- /dev/null
+++ b/lib/core/tests/equilibrium.cpp
@@ -0,0 +1,40 @@
+#include <forstio/test/suite.hpp>
+
+#include "../c++/equilibrium.hpp"
+
+
+namespace {
+
+template<typename Descriptor>
+void check_equilibrium(){
+	using namespace kel;
+
+	using dfi = lbm::df_info<lbm::sch::Float64,Descriptor>;
+
+	saw::data<lbm::sch::Float64> rho{1.0};
+	saw::data<lbm::sch::FixedArray<lbm::sch::Float64,Descriptor::D>> vel;
+	for(saw::data<lbm::sch::UInt64> i{0u}; i.get() < Descriptor::D; ++i){
+		vel.at(i) = {0.0};
+	}
+	auto eq = lbm::equilibrium<lbm::sch::Float64,Descriptor>(rho,vel);
+
+	for(saw::data<lbm::sch::UInt64> i{0u}; i.get() < Descriptor::Q; ++i){
+		SAW_EXPECT(eq.at(i).get() == dfi::weights[i.get()], std::string{"No velocity and normalized rho should be exactly the weights: "} + std::to_string(eq.at(i).get()) + std::string{" "} + std::to_string(dfi::weights[i.get()]));
+	}
+}
+
+SAW_TEST("Equilibrium at rest D1Q3"){
+	using namespace kel;
+	check_equilibrium<lbm::sch::Descriptor<1,3>>();
+}
+
+SAW_TEST("Equilibrium at rest D2Q5"){
+	using namespace kel;
+	check_equilibrium<lbm::sch::Descriptor<2,5>>();
+}
+
+SAW_TEST("Equilibrium at rest D2Q9"){
+	using namespace kel;
+	check_equilibrium<lbm::sch::Descriptor<2,9>>();
+}
+}
-- 
cgit v1.2.3