Update pycuda/gpuarray.py

Co-authored-by: Andreas Klöckner <[email protected]>

pycuda/gpuarray.py

@@ -576,7 +576,7 @@ def _new_like_me(self, dtype=None, order="C"):
     # operators ---------------------------------------------------------------
     def mul_add(self, selffac, other, otherfac, add_timer=None, stream=None, out=None):
         """Return `selffac * self + otherfac*other`."""
-        result = self._new_like_me(_get_common_dtype(self, other)) if out is None else out
+        result = out if out is not None else self._new_like_me(_get_common_dtype(self, other))
         return self._axpbyz(selffac, other, otherfac, result, add_timer, stream=stream)
 
     def __add__(self, other):