From 7d094fe719923db5eb6b30edcc4fa01c6ce2cc2e Mon Sep 17 00:00:00 2001 From: "Claudius \"keldu\" Holeksa" Date: Fri, 7 Aug 2026 12:34:22 +0200 Subject: Cleaning up code --- modules/core/c++/abstract/data.hpp | 89 +++++++++++++++++++++++++ modules/core/c++/abstract/error.hpp | 2 +- modules/core/c++/particle/aabb.hpp | 1 - modules/core/c++/particle/cuboid/aabb.hpp | 40 +++++++++++ modules/core/c++/particle/cuboid/common.hpp | 10 +++ modules/core/c++/particle/cuboid/particle.hpp | 9 +++ modules/core/c++/particle/cuboid/porosity.hpp | 10 +++ modules/core/c++/particle/particle.hpp | 4 +- modules/core/c++/particle/spheroid/particle.hpp | 9 +++ 9 files changed, 170 insertions(+), 4 deletions(-) create mode 100644 modules/core/c++/particle/cuboid/aabb.hpp create mode 100644 modules/core/c++/particle/cuboid/common.hpp create mode 100644 modules/core/c++/particle/cuboid/particle.hpp create mode 100644 modules/core/c++/particle/cuboid/porosity.hpp create mode 100644 modules/core/c++/particle/spheroid/particle.hpp (limited to 'modules/core/c++') diff --git a/modules/core/c++/abstract/data.hpp b/modules/core/c++/abstract/data.hpp index f1ae5a7..327fb63 100644 --- a/modules/core/c++/abstract/data.hpp +++ b/modules/core/c++/abstract/data.hpp @@ -10,8 +10,97 @@ class data final {}; template class data,Encode> final { +public: + using Schema = sch::Primitive; + using Encode = Encode; private: + native_data_type::type value_; public: + data(): + value_{} + {} + + data(native_data_type::type value__): + value_{value__} + {} + + constexpr auto get() const { + return value_; + } + + void set(native_data_type value__){ + value_ = value__; + } + + constexpr bool operator==(const data& rhs) const { + return value_ == rhs.value_; + } + + constexpr bool operator!=(const data& rhs) const { + return value_ != rhs.value_; + } + + constexpr bool operator>(const data& rhs) const { + return value_ > rhs.value_; + } + + constexpr bool operator<(const data& rhs) const { + return value_ < rhs.value_; + } + + constexpr bool operator>=(const data& rhs) const { + return value_ >= rhs.value_; + } + + constexpr bool operator<=(const data& rhs) const { + return value_ >= rhs.value_; + } + + constexpr bool equals(const data& rhs) const { + return (*this) == rhs; + } + + data operator+(const data& rhs) const { + return {value_ + rhs.value_}; + } + + data operator-(const data& rhs) const { + return {value_ - rhs.value_}; + } + + data operator*(const data& rhs) const { + return {value_ * rhs.value_}; + } + + data operator/(const data& rhs) const { + return {value_ / rhs.value_}; + } + + data& operator+=(const data& rhs) { + value_ += rhs.value_; + return *this + } + + data& operator-=(const data& rhs) { + value_ -= rhs.value_; + return *this + } + + data& operator*=(const data& rhs) { + value_ *= rhs.value_; + return *this + } + + data& operator/=(const data& rhs) { + value_ /= rhs.value_; + return *this + } + + template + data cast_to() const { + data val{static_cast::type>(value_)}; + return val; + } }; diff --git a/modules/core/c++/abstract/error.hpp b/modules/core/c++/abstract/error.hpp index 19917de..94a2c35 100644 --- a/modules/core/c++/abstract/error.hpp +++ b/modules/core/c++/abstract/error.hpp @@ -305,7 +305,7 @@ public: /** * This tries to catch cases where error starts including itself as a type which can happen in more complicated cases. - * So this acts as a type safe guard. + * So this acts as a type safe guard in case internal handling fails to treat conversion correctly. */ template class error_or> { private: diff --git a/modules/core/c++/particle/aabb.hpp b/modules/core/c++/particle/aabb.hpp index da39ec2..e4c930b 100644 --- a/modules/core/c++/particle/aabb.hpp +++ b/modules/core/c++/particle/aabb.hpp @@ -26,7 +26,6 @@ public: static constexpr saw::error_or> calculate(const saw::data& p_grp, const saw::data>& index, const saw::data>& meta){ static_assert(PC > 0u, "Can't calculate from no particles"); if(not (index.at({{0u}}).get() < PC) ){ - std::cerr.flush(); return saw::make_error("Too large i in particle_aabb"); } diff --git a/modules/core/c++/particle/cuboid/aabb.hpp b/modules/core/c++/particle/cuboid/aabb.hpp new file mode 100644 index 0000000..502132a --- /dev/null +++ b/modules/core/c++/particle/cuboid/aabb.hpp @@ -0,0 +1,40 @@ +#pragma once + +#include "common.hpp" +#include "../aabb.hpp" + +namespace kel { +namespace lbm { + +template +class particle_aabb< + sch::ParticleGroup> +> final { +public: + using Schema = sch::ParticleGroup>; + + using AABB = sch::Struct< + sch::Member, "a">, + sch::Member, "b"> + >; + +public: + template + static constexpr saw::error_or> calculate(const saw::data& pg, const saw::data>& index, const saw::data>& meta){ + static_assert(PC > 0u, "Can't calculate from no particle"); + if(not (index.at({0u}).get() < PC)){ + return saw::make_error("Too large i in particle_aabb coll::Cuboid"); + } + + saw::data aabb; + auto& parts = pg.template get<"particles">(); + auto& pi = parts.at(index); + auto& pirb = pi.template get<"rigid_body">(); + + /// TODO + + return aabb; + } +}; +} +} diff --git a/modules/core/c++/particle/cuboid/common.hpp b/modules/core/c++/particle/cuboid/common.hpp new file mode 100644 index 0000000..a859cbe --- /dev/null +++ b/modules/core/c++/particle/cuboid/common.hpp @@ -0,0 +1,10 @@ +#pragma once + +namespace kel { +namespace lbm { +namespace coll { +template +struct Cuboid {}; +} +} +} diff --git a/modules/core/c++/particle/cuboid/particle.hpp b/modules/core/c++/particle/cuboid/particle.hpp new file mode 100644 index 0000000..e19d543 --- /dev/null +++ b/modules/core/c++/particle/cuboid/particle.hpp @@ -0,0 +1,9 @@ +#pragma once + +#include "common.hpp" +#include "porosity.hpp" + +namespace kel { +namespace lbm { +} +} diff --git a/modules/core/c++/particle/cuboid/porosity.hpp b/modules/core/c++/particle/cuboid/porosity.hpp new file mode 100644 index 0000000..bccfebf --- /dev/null +++ b/modules/core/c++/particle/cuboid/porosity.hpp @@ -0,0 +1,10 @@ +#pragma once + +#include "common.hpp" +#include "../porosity.hpp" + +namespace kel { +namespace lbm { + +} +} diff --git a/modules/core/c++/particle/particle.hpp b/modules/core/c++/particle/particle.hpp index a3669b4..21b4b34 100644 --- a/modules/core/c++/particle/particle.hpp +++ b/modules/core/c++/particle/particle.hpp @@ -82,8 +82,8 @@ saw::data>> create_sphero total_mass.at({}) = rad_d * rad_d * density.at({}) * 3.141592; }else if constexpr ( D == 3u ){ - }else if constexpr ( D== 1u ){ - + }else if constexpr ( D == 1u ){ + total_mass.at({}) = rad_d * 2.0; } std::cout<<"Total Mass: "<