     float norm;
     
     // These variables currently govern the tolerance (whether it recurses on a cell)
-    float epssq = 0.05 * 0.05;
-    float itolsq = 1.0f / (0.5 * 0.5);
+    float theta = 0.5f;
     InitializationKernel<<<1, 1>>>(thrust::raw_pointer_cast(errl.data()));
     gpuErrchk(cudaDeviceSynchronize());
     

void ForceCalculationKernel(int nnodesd, 
                            int nbodiesd, 
                            volatile int * __restrict errd, 
                            float itolsqd, 
                            float epssqd, 
                            volatile int * __restrict sortd, 
                            volatile int * __restrict childd, 
                            volatile float * __restrict massd, 

        gpuErrchk(cudaDeviceSynchronize());
        
        ForceCalculationKernel<<<blocks * FACTOR5, THREADS5>>>(nnodes, N_POINTS, thrust::raw_pointer_cast(errl.data()), 
                                                                    theta,
                                                                    thrust::raw_pointer_cast(sortl.data()), 
                                                                    thrust::raw_pointer_cast(childl.data()), 
                                                                    thrust::raw_pointer_cast(massl.data()), 
