@@ -13,8 +13,9 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
1313
1414 contains
1515
16+ ! BLAS/LAPACK backends do not currently support xdp
1617 #:for rk,rt in RC_KINDS_TYPES
17- #:if rk!="xdp" ! BLAS/LAPACK backends do not currently support xdp
18+ #:if rk!="xdp"
1819 module function stdlib_linalg_pure_${rt[0]}$${rk}$determinant(a) result(det)
1920 !!### Summary
2021 !! Compute determinant of a real square matrix (pure interface).
@@ -54,8 +55,8 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
5455 if (m/=n .or. .not.min(m,n)>=0) then
5556 err0 = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid or non-square matrix: a=[',m,',',n,']')
5657 det = 0.0_${rk}$
57- ! Process output and return
58- call linalg_error_handling(err0)
58+ ! Process output and return
59+ call linalg_error_handling(err0)
5960 return
6061 end if
6162
@@ -110,7 +111,7 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
110111 end select
111112
112113 ! Process output and return
113- 1 call linalg_error_handling(err0)
114+ call linalg_error_handling(err0)
114115
115116 end function stdlib_linalg_pure_${rt[0]}$${rk}$determinant
116117
@@ -161,9 +162,9 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
161162 if (m/=n .or. .not.min(m,n)>=0) then
162163 err0 = linalg_state_type(this,LINALG_VALUE_ERROR,'invalid or non-square matrix: a=[',m,',',n,']')
163164 det = 0.0_${rk}$
164- ! Process output and return
165- call linalg_error_handling(err0)
166- return
165+ ! Process output and return
166+ call linalg_error_handling(err0,err )
167+ return
167168 end if
168169
169170 ! Can A be overwritten? By default, do not overwrite
@@ -228,7 +229,7 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
228229 end select
229230
230231 ! Process output and return
231- 1 call linalg_error_handling(err0,err)
232+ call linalg_error_handling(err0,err)
232233
233234 end function stdlib_linalg_${rt[0]}$${rk}$determinant
234235
0 commit comments