Merge branch 'determinant' of github.com:perazz/stdlib into determinant

Author: perazz

src/stdlib_linalg.fypp

@@ -229,9 +229,9 @@ module stdlib_linalg
     !!### Description
     !! 
     !! This interface provides methods for computing the determinant of a matrix.
-    !! Supported data types include real and complex.
+    !! Supported data types include `real` and `complex`.
     !! 
-    !!@note The provided functions are intended for square matrices.          
+    !!@note The provided functions are intended for square matrices only.          
     !! 
     !!### Example
     !!

src/stdlib_linalg_determinant.fypp

@@ -54,7 +54,9 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
          if (m/=n .or. .not.min(m,n)>=0) then
             err0 = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid or non-square matrix: a=[',m,',',n,']')
             det  = 0.0_${rk}$
-            goto 1
+         ! Process output and return
+         call linalg_error_handling(err0)
+         return
          end if
 
          select case (m)
@@ -159,7 +161,9 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
          if (m/=n .or. .not.min(m,n)>=0) then
             err0 = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid or non-square matrix: a=[',m,',',n,']')
             det  = 0.0_${rk}$
-            goto 1
+          ! Process output and return
+         call linalg_error_handling(err0)
+         return
          end if
 
          ! Can A be overwritten? By default, do not overwrite
@@ -186,7 +190,7 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
 
                 ! Initialize a matrix temporary
                 if (copy_a) then
-                   allocate(amat(m,n),source=a)
+                   allocate(amat, source=a)
                 else
                    amat => a
                 endif