Merge pull request #1178 from peastman/nancheck

Added checks for nan coordinates in CUDA and OpenCL

Merge pull request #1178 from peastman/nancheck
Added checks for nan coordinates in CUDA and OpenCL
56d3b854 · peastman · 63e95336 · dd5e8e8b · 56d3b854 · 56d3b854
Commit 56d3b854 authored Oct 09, 2015 by peastman
Show whitespace changes
Inline Side-by-side

Showing with 4 additions and 0 deletions

platforms/cuda/src/CudaContext.cpp platforms/cuda/src/CudaContext.cpp +2 -0

platforms/opencl/src/OpenCLContext.cpp platforms/opencl/src/OpenCLContext.cpp +2 -0

No files found.
--- a/platforms/cuda/src/CudaContext.cpp
+++ b/platforms/cuda/src/CudaContext.cpp
@@ -1196,6 +1196,8 @@ void CudaContext::reorderAtomsImpl() {
            molPos[i].x *= invNumAtoms;
            molPos[i].y *= invNumAtoms;
            molPos[i].z *= invNumAtoms;
+            if (molPos[i].x != molPos[i].x)
+                throw OpenMMException("Particle coordinate is nan");
        }
        if (nonbonded->getUsePeriodic()) {
            // Move each molecule position into the same box.

--- a/platforms/opencl/src/OpenCLContext.cpp
+++ b/platforms/opencl/src/OpenCLContext.cpp
@@ -1112,6 +1112,8 @@ void OpenCLContext::reorderAtomsImpl() {
            molPos[i].x *= invNumAtoms;
            molPos[i].y *= invNumAtoms;
            molPos[i].z *= invNumAtoms;
+            if (molPos[i].x != molPos[i].x)
+                throw OpenMMException("Particle coordinate is nan");
        }
        if (nonbonded->getUsePeriodic()) {
            // Move each molecule position into the same box.