Moved handle_gesv_info and signature.

Author: loiseaujc

src/lapack/stdlib_linalg_lapack_aux.fypp

@@ -43,6 +43,7 @@ module stdlib_linalg_lapack_aux
      public :: handle_potrf_info
      public :: handle_getri_info
      public :: handle_gesdd_info
+     public :: handle_gesv_info
 
      ! SELCTG is a LOGICAL FUNCTION of three DOUBLE PRECISION arguments 
      ! used to select eigenvalues to sort to the top left of the Schur form. 
@@ -1363,4 +1364,29 @@ module stdlib_linalg_lapack_aux
 
      end subroutine handle_gesdd_info
 
+     elemental subroutine handle_gesv_info(this,info,lda,n,nrhs,err)
+         character(len=*), intent(in) :: this
+         integer(ilp), intent(in) :: info,lda,n,nrhs
+         type(linalg_state_type), intent(out) :: err
+
+         ! Process output
+         select case (info)
+            case (0)
+                ! Success
+            case (-1)
+                err = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid problem size n=',n)
+            case (-2)
+                err = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid rhs size n=',nrhs)
+            case (-4)
+                err = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid matrix size a=',[lda,n])
+            case (-7)
+                err = linalg_state_type(this,LINALG_ERROR,'invalid matrix size a=',[lda,n])
+            case (1:)
+                err = linalg_state_type(this,LINALG_ERROR,'singular matrix')
+            case default
+                err = linalg_state_type(this,LINALG_INTERNAL_ERROR,'catastrophic error')
+         end select
+         
+     end subroutine handle_gesv_info
+
 end module stdlib_linalg_lapack_aux

src/stdlib_linalg_solve.fypp

@@ -8,6 +8,7 @@ submodule (stdlib_linalg) stdlib_linalg_solve
 !! Solve linear system Ax=b
      use stdlib_linalg_constants
      use stdlib_linalg_lapack, only: gesv
+     use stdlib_linalg_lapack_aux, only: handle_gesv_info
      use stdlib_linalg_state, only: linalg_state_type, linalg_error_handling, LINALG_ERROR, &
           LINALG_INTERNAL_ERROR, LINALG_VALUE_ERROR     
      implicit none
@@ -16,29 +17,6 @@ submodule (stdlib_linalg) stdlib_linalg_solve
 
      contains
 
-     elemental subroutine handle_gesv_info(info,lda,n,nrhs,err)
-         integer(ilp), intent(in) :: info,lda,n,nrhs
-         type(linalg_state_type), intent(out) :: err
-
-         ! Process output
-         select case (info)
-            case (0)
-                ! Success
-            case (-1)
-                err = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid problem size n=',n)
-            case (-2)
-                err = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid rhs size n=',nrhs)
-            case (-4)
-                err = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid matrix size a=',[lda,n])
-            case (-7)
-                err = linalg_state_type(this,LINALG_ERROR,'invalid matrix size a=',[lda,n])
-            case (1:)
-                err = linalg_state_type(this,LINALG_ERROR,'singular matrix')
-            case default
-                err = linalg_state_type(this,LINALG_INTERNAL_ERROR,'catastrophic error')
-         end select
-         
-     end subroutine handle_gesv_info
 
      #:for nd,ndsuf,nde in ALL_RHS
      #:for rk,rt,ri in RC_KINDS_TYPES
@@ -152,7 +130,7 @@ submodule (stdlib_linalg) stdlib_linalg_solve
          call gesv(n,nrhs,amat,lda,ipiv,xmat,ldb,info)
 
          ! Process output
-         call handle_gesv_info(info,lda,n,nrhs,err0)
+         call handle_gesv_info(this,info,lda,n,nrhs,err0)
 
          if (copy_a) deallocate(amat)
          if (.not.present(pivot)) deallocate(ipiv)