   thrust::device_vector<float> dLambdas(BS * NY * NX);
   thrust::device_vector<float> dExpLambdas(NY * NX);
   thrust::device_vector<float> dLogLs(NY * NX);
-  thrust::device_vector<float> dIntensity(BS * 3 * 48 * 32 * 32);
-  thrust::host_vector<float> hIntensity(BS * 3 * 48 * 32 * 32);
-  thrust::device_vector<float> dIntensity2(BS * 3 * 48 * 32 * 32);
-  thrust::host_vector<float> hIntensity2(BS * 3 * 48 * 32 * 32);
+  thrust::device_vector<float> dIntensity(BS * 1 * 48 * 32 * 32);
+  thrust::host_vector<float> hIntensity(BS * 1 * 48 * 32 * 32);
+  thrust::device_vector<float> dIntensity2(BS * 1 * 48 * 32 * 32);
+  thrust::host_vector<float> hIntensity2(BS * 1 * 48 * 32 * 32);
   thrust::device_vector<int> dpsum(BS);
   thrust::device_vector<float> dlsum(BS);
   thrust::host_vector<long long> hPhc(BS * 3);



    float factor = (/*0.01f + 0.99f **/ expf(- gr / (2 * sigmag * sigmag)));
    //factor *= sinc(qx / sigma * 0.5) * sinc(qy / sigma * 0.5);

    //float factor = 1.0f;


intens[idx] = 0;*/
intensityfactor = 0;
}
//  intensityfactor *= pow(1.1f, -1.f + by * 2.f);
  likelihood likelihooder(myspherer, intensityfactor);

  return likelihooder;
