   int rc = 0;
   for (int img = 0; img < fullsize[0]; img+=BS, rc++)
     {
-      if (rc % 64 != tid) continue;
+      if (rc % 96 != tid) continue;
       int end = min((int) fullsize[0], img + BS);
 	int imgcount = end - img;
 	count[0] = imgcount;

	      for (int dy = -0.6 * reals.sigma; dy <= 0.6 * reals.sigma; dy+= 4)
		{
		  if (dx == 0 && dy > 0) continue;
		  reals.x = dx;
		  reals.y = dy;
		  thrust::tabulate(thrust::device, d_complexSpace.data(), d_complexSpace.data() + FX * FY, reals);

      fill(&minval[0], &minval[BS], 1e30);
      fill(&maxval[0], &maxval[BS], -1e30);
      
      for (int r = 0; r < 1200; r++)
	{
	  if (r > 350) r += 4;
	  if (r > 600) r += 5;
	  if (r > 1200) r += 10;
	  if (r > 1500) r += 20;
