ENH: Added np.seterr handling of errors

ovillellas · pv · commit f48c1342026f · 2013-04-10T22:47:43.000+03:00
diff --git a/numpy/core/src/umath/gufuncs_linalg.py b/numpy/core/src/umath/gufuncs_linalg.py
@@ -13,10 +13,11 @@
 appropriate kernel depending on the parameters. All wrappers forward the keyword
 parameters to the underlying generalized ufunc (the kernel).
 
-The functions are intended to be used on arrays of matrices. For those
-functions where a result may not be possible to obtain (like the inverse of
-a matrix that is not invertible) no exception is raised, but the results for
-the elements involved are set to NaN.
+The functions are intended to be used on arrays of matrices. Functions that in
+numpy.LinAlg would generate a LinAlgError (for example, inv on a non-invertible
+matrix) will just generate NaNs as result. When this happens, invalid floating
+point status will be set. Error handling can be configured for this cases
+using np.seterr.
 
 Additional functions some fused arithmetic, useful for efficient operation over
 """
diff --git a/numpy/core/src/umath/umath_linalg.c.src b/numpy/core/src/umath/umath_linalg.c.src
@@ -812,6 +812,10 @@ nan_@TYPE@_matrix(void *dst_in, const LINEARIZE_DATA_t* data)
         }
         dst += data->row_strides/sizeof(@typ@);
     }
+
+    /* if this function gets called, the module generates nans */
+    /* set the fpe status to invalid to notify that nans were generated */
+    npy_set_floatstatus_invalid();
 }
 
 static inline void

Original file line number	Diff line number	Diff line change
`@@ -812,6 +812,10 @@ nan_@TYPE@_matrix(void dst_in, const LINEARIZE_DATA_t data)`
`812`	`812`	`}`
`813`	`813`	`dst += data->row_strides/sizeof(@typ@);`
`814`	`814`	`}`
	`815`	`+`
	`816`	`+ /* if this function gets called, the module generates nans */`
	`817`	`+ /* set the fpe status to invalid to notify that nans were generated */`
	`818`	`+ npy_set_floatstatus_invalid();`
`815`	`819`	`}`
`816`	`820`
`817`	`821`	`static inline void`