Fixed error handling.

Author: loiseaujc

src/stdlib_linalg_matrix_functions.fypp

@@ -48,7 +48,7 @@ contains
                                      'invalid matrix sizes: A=',[lda,n], &
                                                           ' E=',[lde,ne])
         else
-            E(:n, :n) = A(:n, :n) ; call stdlib_linalg_${ri}$_expm_inplace(E, order, err)
+            E(:n, :n) = A(:n, :n) ; call stdlib_linalg_${ri}$_expm_inplace(E, order, err0)
         endif
 
         ! Process output and return

test/linalg/test_linalg_expm.fypp

@@ -23,7 +23,7 @@ module test_linalg_expm
 
         #:for rk,rt,ri in RC_KINDS_TYPES
         tests = [tests, new_unittest("expm_${ri}$",test_expm_${ri}$)]
-        ! tests = [tests, new_unittest("Error-handling expm_${ri}$",test_error_handling_expm_${ri}$)]
+        tests = [tests, new_unittest("Error-handling expm_${ri}$",test_error_handling_expm_${ri}$)]
         #:endfor
 
     end subroutine test_expm_computation
@@ -87,8 +87,19 @@ module test_linalg_expm
         call check(error, err%error(), "Negative Pade order")
         if (allocated(error)) return
 
+        call matrix_exp(A, order=-1, err=err)
+        ! Check result.
+        call check(error, err%error(), "Negative Pade order")
+        if (allocated(error)) return
+
+        ! Compute matrix exponential.
+        call matrix_exp(A, E(:n, :n-1), err=err)
+        ! Check result.
+        call check(error, err%error(), "Invalid matrix size")
+        if (allocated(error)) return
+
         ! Compute matrix exponential.
-        call matrix_exp(A(:n, :n-1), E, err=err)
+        call matrix_exp(A(:n, :n-1), err=err)
         ! Check result.
         call check(error, err%error(), "Invalid matrix size")
         if (allocated(error)) return