BUG: Fix bug with zero-assignment with nans, see #1263

odlgroup · Jan 21, 2018 · 78fab76 · 78fab76
1 parent 4d30821
commit 78fab76
Showing 1 changed file with 4 additions and 1 deletion.
diff --git a/odl/space/npy_tensors.py b/odl/space/npy_tensors.py
@@ -1856,7 +1856,10 @@ def fallback_copy(x1, x2, n):
         _lincomb_impl(a + b, x1, 0, x1, out)
     elif out is x1 and out is x2:
         # All the vectors are aligned -> out = (a+b)*out
-        scal(a + b, out_arr, size)
+        if (a + b) != 0:
+            scal(a + b, out_arr, size)
+        else:
+            out_arr[:] = 0
     elif out is x1:
         # out is aligned with x1 -> out = a*out + b*x2
         if a != 1: