.inv. operator: error stop on issues

Author: perazz

doc/specs/stdlib_linalg.md

@@ -924,7 +924,7 @@ This interface is equivalent to the function version of  [[stdlib_linalg(module)
 
 Returns a rank-2 square array with the same type, kind and rank as `a`, that contains the inverse of `a`.
 
-If the matrix is singular or non-square, the operator returns zeros.
+Exceptions always trigger an `error stop`.
 
 ### Example
 

src/stdlib_linalg_inverse.fypp

@@ -154,12 +154,8 @@ submodule (stdlib_linalg) stdlib_linalg_inverse
          !> Result matrix
          ${rt}$, allocatable :: inva(:,:)
 
-         type(linalg_state_type) :: err0
-
-         inva = stdlib_linalg_inverse_${ri}$(a,err0)
-
-         ! On error, return zeros
-         if (err0%error()) inva = 0.0_${rk}$
+         ! Do not provide an error handler (error stop on issues)
+         inva = stdlib_linalg_inverse_${ri}$(a)
 
      end function stdlib_linalg_inverse_${ri}$_operator