Fixed compilation error when using double precision for CustomIntegrator

f820e172 · Peter Eastman · 9e4cf579 · f820e172
Commit f820e172 authored Jan 09, 2012 by Peter Eastman
Show whitespace changes
Inline Side-by-side

Showing with 6 additions and 2 deletions

platforms/opencl/src/OpenCLKernels.cpp platforms/opencl/src/OpenCLKernels.cpp +6 -2

No files found.
--- a/platforms/opencl/src/OpenCLKernels.cpp
+++ b/platforms/opencl/src/OpenCLKernels.cpp
@@ -3677,8 +3677,12 @@ void OpenCLIntegrateCustomStepKernel::execute(ContextImpl& context, CustomIntegr
                    compute << createPerDofComputation(stepType[step] == CustomIntegrator::ComputePerDof ? variable[step] : "", expression[step], i, integrator);
                string convert = (cl.getSupportsDoublePrecision() ? "convert_float4(" : "(");
                if (variable[step] == "x") {
-                    if (storePosAsDelta[step])
+                    if (storePosAsDelta[step]) {
-                        compute << "posDelta[index] = " << convert << "position-posq[index]);\n";
+                        if (cl.getSupportsDoublePrecision())
+                            compute << "posDelta[index] = convert_float4(position-convert_double4(posq[index]));\n";
+                        else
+                            compute << "posDelta[index] = position-posq[index];\n";
+                    }
                    else
                        compute << "posq[index] = " << convert << "position);\n";
                }