1 files changed, 0 insertions, 79 deletions
diff --git a/template/kernel.mako b/template/kernel.mako
index a57c2ec..025baf5 100644
--- a/template/kernel.mako
+++ b/template/kernel.mako
@@ -98,82 +98,3 @@ __kernel void collect_moments(__global __read_only  ${float_type}* f,
     moments[${pop_offset(i)} + gid] = ${ccode(expr.rhs)};
 % endfor
 }
-
-__kernel void collect_gl_moments(__global __read_only  ${float_type}* f,
-                                 __global __read_only  int* material,
-                                 __global __write_only float4* moments)
-{
-    const unsigned int gid = ${gid()};
-
-    __global __read_only ${float_type}* preshifted_f = f + gid;
-
-% for i in range(0,descriptor.q):
-    const ${float_type} f_curr_${i} = preshifted_f[${pop_offset(i)}];
-% endfor
-
-% for i, expr in enumerate(moments_subexpr):
-    const ${float_type} ${expr[0]} = ${ccode(expr[1])};
-% endfor
-
-    float4 data;
-
-    if (material[gid] == 1) {
-% if descriptor.d == 2:
-      data.x = ${ccode(moments_assignment[0].rhs)};
-      data.y = ${ccode(moments_assignment[1].rhs)};
-      data.z = ${ccode(moments_assignment[2].rhs)};
-      data.w = sqrt(data.y*data.y + data.z*data.z);
-% elif descriptor.d == 3:
-      data.x = ${ccode(moments_assignment[0].rhs)};
-      data.y = ${ccode(moments_assignment[1].rhs)};
-      data.z = ${ccode(moments_assignment[2].rhs)};
-      data.w = ${ccode(moments_assignment[3].rhs)};
-% endif
-    } else {
-      data.x = 0.0;
-      data.y = 0.0;
-      data.z = 0.0;
-      data.w = -material[gid];
-    }
-
-    moments[gid] = data;
-}
-
-__kernel void update_particles(__global __read_only  float4* moments,
-                               __global __read_only  int*    material,
-                               __global __write_only float4* particles,
-                               __global __read_only  float4* init_particles,
-                               float aging)
-{
-  const unsigned int pid = get_global_id(0);
-
-  float4 particle = particles[pid];
-
-% if descriptor.d == 2:
-  const unsigned int gid = floor(particle.y)*${memory.size_x} + floor(particle.x);
-% elif descriptor.d == 3:
-  const unsigned int gid = floor(particle.z)*${memory.size_x*memory.size_y} + floor(particle.y)*${memory.size_x} + floor(particle.x);
-% endif
-
-  const float4 moment = moments[gid];
-
-  if (material[gid] == 1) {
-    particle.x += moment.y;
-    particle.y += moment.z;
-% if descriptor.d == 2:
-    particle.w += min(particle.x, particle.y) * aging;
-% elif descriptor.d == 3:
-    particle.z += moment.w;
-    float dy = (particle.y-${geometry.size_y/2.0});
-    float dz = (particle.z-${geometry.size_z/2.0});
-    dy *= dy;
-    dz *= dz;
-    particle.w = 10.0*sqrt(moment.y*moment.y+moment.z*moment.z+moment.w*moment.w);
-% endif
-  } else {
-    particle.xyz = init_particles[pid].xyz;
-    particle.w   = 0.0;
-  }
-
-  particles[pid] = particle;
-}