                                                                              counts.e0 + counts.e1,
                                                                              counts.e0 + counts.e1 + counts.e2,
                                                                              counts.e0 + counts.e1 + counts.e2 + counts.e3);
+          }
+          if (indexed && vsum) {
+            out_vertex_level0_d[32 * 5 * blockIdx.x + 32 * y + thread] = make_uint4(vc0,
+                                                                                    vc0 + vc1,
+                                                                                    vc0 + vc1 + vc2,
+                                                                                    vc0 + vc1 + vc2 + vc3);
+
           }
           
         }

  {
    return make_float2(__fmaf_rn(float(1 << 4), z1.x, z0.x),
                       __fmaf_rn(float(1 << 4), z1.y, z0.y));
  }

  __device__ __forceinline__ uint32_t axisCountFromCase(const uint32_t c)

      hp5_result r = traverseDown(pyramid, level_offset, level_offset[15], ix);

      uint8_t index_case = index_cases[r.offset];
      uint3 p0 = decodeCellPosition(chunks, r.offset);

      uint32_t axes = piercingAxesFromCase(index_case);
      if (0 < r.remainder) axes = axes & (axes - 1);  // remove rightmost bit
      if (1 < r.remainder) axes = axes & (axes - 1);  // remove rightmost bit
      axes = axes & (-axes);                          // isolate rightmost bit

      uint3 p1 = p0;
      p1.x += (axes >> 1) & 1;
      p1.y += (axes >> 2) & 1;
      p1.z += (axes >> 4) & 1;

      emitVertexNormal(output + 6 * ix,
                       field,
