Reworking structure

1 files changed, 39 insertions, 0 deletions

diff --git a/lib/core/tests/descriptor.cpp b/lib/core/tests/descriptor.cpp
new file mode 100644
index 0000000..7f743ce
--- /dev/null
+++ b/lib/core/tests/descriptor.cpp
@@ -0,0 +1,39 @@
+#include <forstio/test/suite.hpp>
+
+#include "../c++/descriptor.hpp"
+
+
+namespace {
+template<typename Descriptor>
+void check_opposite_dirs(){
+	using namespace kel;
+
+	using dfi = lbm::df_info<lbm::sch::Float32, Descriptor>;
+
+	for(uint64_t k = 0u; k < Descriptor::Q; ++k){
+		auto k_inv = dfi::opposite_index[k];
+
+		for(uint64_t i = 0u; i < Descriptor::D; ++i){
+			SAW_EXPECT(dfi::directions[k][i] == (-1*dfi::directions[k_inv][i]), "Opposites are inconsistent");
+		}
+	}
+}
+
+SAW_TEST("Opposites and Dirs D1Q3"){
+	using namespace kel;
+	check_opposite_dirs<lbm::sch::Descriptor<1,3>>();
+}
+
+SAW_TEST("Opposites and Dirs D2Q5"){
+	using namespace kel;
+	check_opposite_dirs<lbm::sch::Descriptor<2,5>>();
+}
+SAW_TEST("Opposites and Dirs D2Q9"){
+	using namespace kel;
+	check_opposite_dirs<lbm::sch::Descriptor<2,9>>();
+}
+SAW_TEST("Opposites and Dirs D3Q27"){
+	using namespace kel;
+	check_opposite_dirs<lbm::sch::Descriptor<3,27>>();
+}
+}