CustomGBForce was storing computed values in single precision even when double was requested

platforms/cuda/src/CudaKernels.cpp

@@ -2246,7 +2246,7 @@ void CudaCalcCustomGBForceKernel::initialize(const System& system, const CustomG
 
     int numParticles = force.getNumParticles();
     params = new CudaParameterSet(cu, force.getNumPerParticleParameters(), numParticles, "customGBParameters", true);
-    computedValues = new CudaParameterSet(cu, force.getNumComputedValues(), numParticles, "customGBComputedValues", true);
+    computedValues = new CudaParameterSet(cu, force.getNumComputedValues(), numParticles, "customGBComputedValues", true, cu.getUseDoublePrecision());
     if (force.getNumGlobalParameters() > 0)
         globals = CudaArray::create<float>(cu, force.getNumGlobalParameters(), "customGBGlobals");
     vector<vector<float> > paramVector(numParticles);

platforms/opencl/src/OpenCLContext.cpp

@@ -149,9 +149,9 @@ OpenCLContext::OpenCLContext(const System& system, int platformIndex, int device
             throw OpenMMException("The specified OpenCL device is not compatible with OpenMM");
         compilationDefines["WORK_GROUP_SIZE"] = intToString(ThreadBlockSize);
         if (platformVendor.size() >= 5 && platformVendor.substr(0, 5) == "Intel")
-			defaultOptimizationOptions = "";
-		else
-	        defaultOptimizationOptions = "-cl-fast-relaxed-math";
+            defaultOptimizationOptions = "";
+        else
+            defaultOptimizationOptions = "-cl-fast-relaxed-math";
         supports64BitGlobalAtomics = (device.getInfo<CL_DEVICE_EXTENSIONS>().find("cl_khr_int64_base_atomics") != string::npos);
         supportsDoublePrecision = (device.getInfo<CL_DEVICE_EXTENSIONS>().find("cl_khr_fp64") != string::npos);
         if ((useDoublePrecision || useMixedPrecision) && !supportsDoublePrecision)

platforms/opencl/src/OpenCLKernels.cpp

@@ -2302,7 +2302,7 @@ void OpenCLCalcCustomGBForceKernel::initialize(const System& system, const Custo
 
     int numParticles = force.getNumParticles();
     params = new OpenCLParameterSet(cl, force.getNumPerParticleParameters(), numParticles, "customGBParameters", true);
-    computedValues = new OpenCLParameterSet(cl, force.getNumComputedValues(), numParticles, "customGBComputedValues", true);
+    computedValues = new OpenCLParameterSet(cl, force.getNumComputedValues(), numParticles, "customGBComputedValues", true, cl.getUseDoublePrecision());
     if (force.getNumGlobalParameters() > 0)
         globals = OpenCLArray::create<cl_float>(cl, force.getNumGlobalParameters(), "customGBGlobals", CL_MEM_READ_ONLY);
     vector<vector<cl_float> > paramVector(numParticles);