 	}
       
       dim3 grid(25, 40, 24);
-      dim3 block(2, 7, 50);
+      dim3 block(2, 1, 50);
 
       spherer.rfactor = 0.005;
       spherer.roffset = -10;

				   thrust::make_transform_iterator(thrust::make_counting_iterator(0), likelihooder),
				   thrust::make_transform_iterator(thrust::make_counting_iterator(NY * NX), likelihooder));
  target[idx] = likelihood;
}

int main()

  thrust::device_vector<float> dLambdas(NY * NX);
  thrust::device_vector<float> dIntensity(50 * 300 * 1200);
  thrust::host_vector<float> hIntensity(50 * 300 * 1200);
  thrust::device_vector<float> dIntensity2(50 * 300 * 1200);
  thrust::host_vector<float> hIntensity2(50 * 300 * 1200);

  idealsphere spherer(dPhotons, dLambdas);

