1 files changed, 43 insertions, 0 deletions
diff --git a/c++/particle/geometry/circle.hpp b/c++/particle/geometry/circle.hpp
new file mode 100644
index 0000000..65b8966
--- /dev/null
+++ b/c++/particle/geometry/circle.hpp
@@ -0,0 +1,43 @@
+#pragma once
+
+namespace kel {
+namespace lbm {
+
+template<typename T>
+class particle_circle_geometry {
+private:
+	saw::data<T> radius_;
+public:
+	particle_circle_geometry(saw::data<T> radius__):
+		radius_{radius__}
+	{}
+
+	template<typename MT>
+	saw::data<sch::ParticleMask<MT,2>> generate_mask(uint64_t resolution, uint64_t boundary_nodes = 0) const {
+		saw::data<sch::ParticleMask<MT,2>> mask;
+
+		auto& grid = mask.template get<"grid">();
+
+		uint64_t size = resolution + 2*boundary_nodes;
+		grid = {{size,size}};
+
+		saw::data<T> radius_squared = radius_ * radius_;
+
+
+
+		for(uint64_t i = 0; i < size; ++i){
+			for(uint64_t j = 0; j < size; ++j){
+				if(i < boundary_nodes || j < boundary_nodes || i >= resolution+boundary_nodes || j >= resolution + boundary_nodes ){
+					grid.at({i,j}).set(0);
+				}else{
+
+				}
+			}
+		}
+
+		return mask;
+	}
+};
+
+}
+}