Extract offset helper into target and layout specific classes

2 files changed, 64 insertions, 0 deletions

diff --git a/boltzgen/kernel/target/cl.py b/boltzgen/kernel/target/cl.py
new file mode 100644
index 0000000..a2fc819
--- /dev/null
+++ b/boltzgen/kernel/target/cl.py
@@ -0,0 +1,19 @@
+class SOA:
+    def __init__(self, descriptor, geometry):
+        self.descriptor = descriptor
+        self.geometry = geometry
+
+    def gid(self):
+        return {
+            2: 'get_global_id(1)*%d + get_global_id(0)' % self.geometry.size_x,
+            3: 'get_global_id(2)*%d + get_global_id(1)*%d + get_global_id(0)' % (self.geometry.size_x*self.geometry.size_y, self.geometry.size_x)
+        }.get(self.descriptor.d)
+
+    def pop_offset(self, i):
+        return i * self.geometry.volume
+
+    def neighbor_offset(self, c_i):
+        return {
+            2: lambda:                                                    c_i[1]*self.geometry.size_x + c_i[0],
+            3: lambda: c_i[2]*self.geometry.size_x*self.geometry.size_y + c_i[1]*self.geometry.size_x + c_i[0]
+        }.get(self.descriptor.d)()
diff --git a/boltzgen/kernel/target/cpp.py b/boltzgen/kernel/target/cpp.py
new file mode 100644
index 0000000..de79716
--- /dev/null
+++ b/boltzgen/kernel/target/cpp.py
@@ -0,0 +1,45 @@
+class AOS:
+    def __init__(self, descriptor, geometry):
+        self.descriptor = descriptor
+        self.geometry = geometry
+
+    def gid_offset(self):
+        return self.descriptor.q
+
+    def pop_offset(self, i):
+        return i
+
+    def neighbor_offset(self, c_i):
+        return self.descriptor.q * {
+            2: lambda:                                                    c_i[0]*self.geometry.size_y + c_i[1],
+            3: lambda: c_i[0]*self.geometry.size_y*self.geometry.size_z + c_i[1]*self.geometry.size_z + c_i[2]
+        }.get(self.descriptor.d)()
+
+    def padding(self):
+        return self.descriptor.q * {
+            2: lambda:                                               1*self.geometry.size_y + 1,
+            3: lambda: 1*self.geometry.size_y*self.geometry.size_z + 1*self.geometry.size_z + 1
+        }.get(self.descriptor.d)()
+
+class SOA:
+    def __init__(self, descriptor, geometry):
+        self.descriptor = descriptor
+        self.geometry = geometry
+
+    def gid_offset(self):
+        return 1
+
+    def pop_offset(self, i):
+        return i * self.geometry.volume
+
+    def neighbor_offset(self, c_i):
+        return {
+            2: lambda:                                                    c_i[0]*self.geometry.size_y + c_i[1],
+            3: lambda: c_i[0]*self.geometry.size_y*self.geometry.size_z + c_i[1]*self.geometry.size_z + c_i[2]
+        }.get(self.descriptor.d)()
+
+    def padding(self):
+        return {
+            2: lambda:                                               1*self.geometry.size_y + 1,
+            3: lambda: 1*self.geometry.size_y*self.geometry.size_z + 1*self.geometry.size_z + 1
+        }.get(self.descriptor.d)()