Just rename from lib to modules

1 files changed, 0 insertions, 249 deletions

diff --git a/lib/core/c++/boundary.hpp b/lib/core/c++/boundary.hpp
deleted file mode 100644
index 01ae7b5..0000000
--- a/lib/core/c++/boundary.hpp
+++ /dev/null
@@ -1,249 +0,0 @@
-#pragma once
-
-#include "component.hpp"
-#include "equilibrium.hpp"
-
-namespace kel {
-namespace lbm {
-namespace cmpt {
-struct BounceBack {};
-
-template<uint64_t i>
-struct AntiBounceBack {};
-
-template<bool East>
-struct ZouHeHorizontal{};
-
-struct Equilibrium {};
-
-template<bool North>
-struct ZouHePressureY{};
-
-template<bool East>
-struct ZouHeVelocityX {};
-}
-
-/**
- * Full-Way BounceBack.
- */
-template<typename T, typename Descriptor, typename Encode>
-class component<T, Descriptor, cmpt::BounceBack, Encode> {
-private:
-public:
-	component() = default;
-
-	using Component = cmpt::BounceBack;
-
-	/**
-	 * Thoughts regarding automating order setup
-	 */
-	using access = cmpt::access_tuple<
-		cmpt::access<"dfs", 1, true, true, true>
-	>;
-
-	/**
-	 * Thoughts regarding automating order setup
-	 */
-	static constexpr saw::string_literal name = "full_way_bounce_back";
-	static constexpr saw::string_literal after = "";
-	static constexpr saw::string_literal before = "streaming";
-
-	/**
-	 * Raw setup
-	 */
-	template<typename CellFieldSchema>
-	void apply(const saw::data<CellFieldSchema, Encode>& field, const saw::data<sch::FixedArray<sch::UInt64,Descriptor::D>>& index, saw::data<sch::UInt64> time_step)const{
-		bool is_even = ((time_step.get() % 2u) == 0u);
-
-		// This is a ref
-
-		// auto& dfs_f = (is_even) ? field.template get<"dfs">() : field.template get<"dfs_old">();
-		auto& dfs_old_f = (is_even) ? field.template get<"dfs_old">() : field.template get<"dfs">();
-
-		using dfi = df_info<T,Descriptor>;
-
-		// Technically use .copy()
-		auto& dfs_old = dfs_old_f.at(index);
-		auto df_cpy = dfs_old;
-
-		for(uint64_t i = 0u; i < Descriptor::Q; ++i){
-			dfs_old.at({i}) = df_cpy.at({dfi::opposite_index.at(i)});
-		}
-	}
-};
-
-/**
- * Anti Full-Way BounceBack.
- */
-template<typename T, typename Descriptor, uint64_t D, typename Encode>
-class component<T, Descriptor, cmpt::AntiBounceBack<D>, Encode> {
-private:
-public:
-	component() = default;
-
-	using Component = cmpt::AntiBounceBack<D>;
-
-	/**
-	 * Thoughts regarding automating order setup
-	 */
-	using access = cmpt::access_tuple<
-		cmpt::access<"dfs", 1, true, true, true>
-	>;
-
-	/**
-	 * Thoughts regarding automating order setup
-	 */
-	static constexpr saw::string_literal name = "full_way_anti_bounce_back";
-	static constexpr saw::string_literal after = "";
-	static constexpr saw::string_literal before = "streaming";
-
-	/**
-	 * Raw setup
-	 */
-	template<typename CellFieldSchema>
-	void apply(const saw::data<CellFieldSchema, Encode>& field, const saw::data<sch::FixedArray<sch::UInt64,Descriptor::D>>& index, saw::data<sch::UInt64> time_step)const{
-		bool is_even = ((time_step.get() % 2u) == 0u);
-
-		// This is a ref
-
-		// auto& dfs_f = (is_even) ? field.template get<"dfs">() : field.template get<"dfs_old">();
-		auto& dfs_old_f = (is_even) ? field.template get<"dfs_old">() : field.template get<"dfs">();
-
-		using dfi = df_info<T,Descriptor>;
-
-		// Technically use .copy()
-		auto& dfs_old = dfs_old_f.at(index);
-		auto df_cpy = dfs_old;
-		
-		static_assert(Descriptor::D == 2u and Descriptor::Q == 9u, "Some parts are hard coded sadly");
-
-		dfs_old.at({0u}) = df_cpy.at({0u});
-		dfs_old.at({1u}) = df_cpy.at({1u});
-		dfs_old.at({2u}) = df_cpy.at({2u});
-		dfs_old.at({3u}) = df_cpy.at({4u});
-		dfs_old.at({4u}) = df_cpy.at({3u});
-		dfs_old.at({5u}) = df_cpy.at({7u});
-		dfs_old.at({6u}) = df_cpy.at({8u});
-		dfs_old.at({7u}) = df_cpy.at({5u});
-		dfs_old.at({8u}) = df_cpy.at({6u});
-	}
-};
-
-template<typename FP, typename Descriptor, typename Encode>
-class component<FP, Descriptor, cmpt::Equilibrium, Encode> final {
-private:
-	saw::data<sch::Scalar<FP>> density_;
-	saw::data<sch::Vector<FP,Descriptor::D>> velocity_;
-public:
-	component(
-		saw::data<sch::Scalar<FP>> density__,
-		saw::data<sch::Vector<FP,Descriptor::D>> velocity__
-	):
-		density_{density__},
-		velocity_{velocity__}
-	{}
-
-	template<typename CellFieldSchema>
-	void apply(const saw::data<CellFieldSchema, Encode>& field, const saw::data<sch::FixedArray<sch::UInt64,Descriptor::D>>& index, saw::data<sch::UInt64> time_step) const {
-		
-		bool is_even = ((time_step.get() % 2u) == 0u);
-		using dfi = df_info<FP,Descriptor>;
-		
-		auto& dfs_old_f = (is_even) ? field.template get<"dfs_old">() : field.template get<"dfs">();
-
-		auto eq = equilibrium<FP,Descriptor>(density_,velocity_);
-
-		dfs_old_f.at(index) = eq;
-	}
-};
-
-/**
- * This is massively hacky and expects a lot of conditions
- * Either this or mirrored along the horizontal line works
- *
- * 0 - 2 - 2
- * 0 - 3 - 1
- * 0 - 3 - 1
- * .........
- * 0 - 3 - 1
- * 0 - 2 - 2
- *
- */
-template<typename FP, typename Descriptor, bool Dir, typename Encode>
-class component<FP, Descriptor, cmpt::ZouHeHorizontal<Dir>, Encode> final {
-private:
-	saw::data<FP> rho_setting_;
-public:
-	component(const saw::data<FP>& rho_setting__):
-		rho_setting_{rho_setting__}
-	{}
-
-	template<typename CellFieldSchema>
-	void apply(const saw::data<CellFieldSchema, Encode>& field, saw::data<sch::FixedArray<sch::UInt64,Descriptor::D>> index, saw::data<sch::UInt64> time_step) const {
-		using dfi = df_info<FP,Descriptor>;
-
-		bool is_even = ((time_step.get() % 2) == 0);
-
-		auto& info_f = field.template get<"info">();
-		
-		// auto& dfs_f = (is_even) ? field.template get<"dfs">() : field.template get<"dfs_old">();
-		auto& dfs_old_f = (is_even) ? field.template get<"dfs_old">() : field.template get<"dfs">();
-
-		auto& dfs_old = dfs_old_f.at(index);
-
-		auto rho_vel_x = [&]() -> saw::data<FP> {
-			if constexpr (Dir){
-				auto S = dfs_old.at({0u})
-					+ dfs_old.at({3u}) + dfs_old.at({4u})
-					+ (dfs_old.at({1u}) + dfs_old.at({5u}) + dfs_old.at({7u})) * 2u;
-				return rho_setting_ - S;
-			}
-			else if constexpr (not Dir) {
-				auto S = dfs_old.at({0u})
-					+ dfs_old.at({3u}) + dfs_old.at({4u})
-					+ (dfs_old.at({2u}) + dfs_old.at({6u}) + dfs_old.at({8u})) * 2u;
-				return S - rho_setting_;
-			}
-			return {};
-		}();
-
-		static_assert(Descriptor::D == 2u and Descriptor::Q == 9u, "Some parts are hard coded sadly");
-
-		if constexpr (Dir) {
-			dfs_old.at({2u}) = dfs_old.at({1u}) + saw::data<FP>{2.0 / 3.0} * rho_vel_x;
-			dfs_old.at({6u}) = dfs_old.at({5u}) + saw::data<FP>{1.0 / 6.0} * rho_vel_x + saw::data<FP>{0.5} * (dfs_old.at({3u}) - dfs_old.at({4u}));
-			dfs_old.at({8u}) = dfs_old.at({7u}) + saw::data<FP>{1.0 / 6.0} * rho_vel_x + saw::data<FP>{0.5} * (dfs_old.at({4u}) - dfs_old.at({3u}));
-		}else if constexpr (not Dir){
-			dfs_old.at({1u}) = dfs_old.at({2u}) - saw::data<FP>{2.0 / 3.0} * rho_vel_x;
-			dfs_old.at({5u}) = dfs_old.at({6u}) - saw::data<FP>{1.0 / 6.0} * rho_vel_x + saw::data<FP>{0.5} * (dfs_old.at({4u}) - dfs_old.at({3u}));
-			dfs_old.at({7u}) = dfs_old.at({8u}) - saw::data<FP>{1.0 / 6.0} * rho_vel_x + saw::data<FP>{0.5} * (dfs_old.at({3u}) - dfs_old.at({4u}));
-		}
-	}
-};
-
-/*
-template<typename FP, typename Descriptor, bool East, typename Encode>
-class component<FP, Descriptor, cmpt::ZouHeVelocityX<East>, Encode> final {
-private:
-	saw::data<sch::Vector<FP,Descriptor::D>> velocity_;
-public:
-	component(
-		saw::data<sch::Vector<FP,Descriptor::D>> velocity__
-	):
-		velocity_{velocity__}
-	{}
-
-	template<typename CellFieldSchema>
-	void apply(const saw::data<CellFieldSchema, Encode>& field, const saw::data<sch::FixedArray<sch::UInt64,Descriptor::D>>& index, saw::data<sch::UInt64> time_step) const {
-		
-		bool is_even = ((time_step.get() % 2u) == 0u);
-		using dfi = df_info<FP,Descriptor>;
-		
-		auto& dfs_old_f = (is_even) ? field.template get<"dfs_old">() : field.template get<"dfs">();
-
-	}
-};
-*/
-
-}
-}