fix: linalg enable det() for xdp kinds (#1113)

Author: perazz

src/linalg/stdlib_linalg.fypp

@@ -876,8 +876,7 @@ module stdlib_linalg
     !! This interface provides methods for computing the determinant of a matrix.
     !! Supported data types include `real` and `complex`.
     !! 
-    !!@note The provided functions are intended for square matrices only.          
-    !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
+    !!@note The provided functions are intended for square matrices only.
     !! 
     !!### Example
     !!
@@ -898,10 +897,8 @@ module stdlib_linalg
     !!```
     !!     
     #:for rk,rt in RC_KINDS_TYPES
-    #:if rk!="xdp" 
     module procedure stdlib_linalg_${rt[0]}$${rk}$determinant
     module procedure stdlib_linalg_pure_${rt[0]}$${rk}$determinant
-    #:endif
     #:endfor
   end interface det
 
@@ -920,7 +917,6 @@ module stdlib_linalg
     !! Supported data types include real and complex.
     !!
     !!@note The provided functions are intended for square matrices.
-    !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
     !!
     !!### Example
     !!
@@ -935,15 +931,12 @@ module stdlib_linalg
     !!```
     !     
     #:for rk,rt in RC_KINDS_TYPES
-    #:if rk!="xdp"
     module procedure stdlib_linalg_pure_${rt[0]}$${rk}$determinant
-    #:endif
     #:endfor        
   end interface operator(.det.)
 
   interface
     #:for rk,rt in RC_KINDS_TYPES
-    #:if rk!="xdp" 
     module function stdlib_linalg_${rt[0]}$${rk}$determinant(a,overwrite_a,err) result(det)
         !> Input matrix a[m,n]
         ${rt}$, intent(inout), target :: a(:,:)
@@ -960,7 +953,6 @@ module stdlib_linalg
         !> Matrix determinant
         ${rt}$ :: det                
     end function stdlib_linalg_pure_${rt[0]}$${rk}$determinant
-    #:endif
     #:endfor
   end interface  
 
@@ -1216,7 +1208,6 @@ module stdlib_linalg
      !! Preallocated space for both eigenvalues `lambda` and the eigenvector matrices must be user-provided.      
      !! 
      !!@note The solution is based on LAPACK's general eigenproblem solvers `*GEEV`.
-     !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
      !!       
     #:for rk,rt,ri in RC_KINDS_TYPES
     #:for ep,ei in EIG_PROBLEM_LIST
@@ -1294,7 +1285,6 @@ module stdlib_linalg
      !! as an optional `type(linalg_state_type)` output flag. 
      !! 
      !!@note The solution is based on LAPACK's general eigenproblem solvers `*GEEV`.
-     !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
      !!       
     #:for rk,rt,ri in RC_KINDS_TYPES
     #:for ep,ei in EIG_PROBLEM_LIST
@@ -1349,7 +1339,6 @@ module stdlib_linalg
      !! Preallocated space for both eigenvalues `lambda` and the eigenvector matrix must be user-provided.      
      !! 
      !!@note The solution is based on LAPACK's eigenproblem solvers `*SYEV`/`*HEEV`.
-     !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
      !!      
     #:for rk,rt,ri in RC_KINDS_TYPES
     module subroutine stdlib_linalg_eigh_${ri}$(a,lambda,vectors,upper_a,overwrite_a,err)
@@ -1392,7 +1381,6 @@ module stdlib_linalg
      !! as an optional `type(linalg_state_type)` output flag. 
      !! 
      !!@note The solution is based on LAPACK's eigenproblem solvers `*SYEV`/`*HEEV`.
-     !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
      !!         
     #:for rk,rt,ri in RC_KINDS_TYPES
     module function stdlib_linalg_eigvalsh_${ri}$(a,upper_a,err) result(lambda)

src/linalg/stdlib_linalg_determinant.fypp

@@ -13,9 +13,7 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
 
      contains
 
-     ! BLAS/LAPACK backends do not currently support xdp
      #:for rk,rt in RC_KINDS_TYPES
-     #:if rk!="xdp" 
      pure module function stdlib_linalg_pure_${rt[0]}$${rk}$determinant(a) result(det)     
      !!### Summary
      !! Compute determinant of a real square matrix (pure interface).
@@ -233,7 +231,6 @@ submodule (stdlib_linalg) stdlib_linalg_determinant
 
      end function stdlib_linalg_${rt[0]}$${rk}$determinant
 
-     #:endif
      #:endfor
 
 end submodule stdlib_linalg_determinant

test/linalg/test_linalg_determinant.fypp

@@ -22,21 +22,16 @@ module test_linalg_determinant
         allocate(tests(0))
 
         #:for rk,rt in RC_KINDS_TYPES
-        #:if rk!="xdp"
         call add_test(tests,new_unittest("$eye_det_${rt[0]}$${rk}$",test_${rt[0]}$${rk}$_eye_determinant))
         call add_test(tests,new_unittest("$eye_det_multiple_${rt[0]}$${rk}$",test_${rt[0]}$${rk}$_eye_multiple))
-        #:endif
         #:endfor
         #:for ck,ct in CMPLX_KINDS_TYPES
-        #:if ck!="xdp"
         call add_test(tests,new_unittest("$complex_det_${rt[0]}$${rk}$",test_${ct[0]}$${ck}$_complex_determinant))
-        #:endif
         #: endfor
 
     end subroutine test_matrix_determinant
 
     #:for rk,rt in RC_KINDS_TYPES
-    #:if rk!="xdp"
     !> Determinant of identity matrix
     subroutine test_${rt[0]}$${rk}$_eye_determinant(error)
         type(error_type), allocatable, intent(out) :: error
@@ -96,12 +91,10 @@ module test_linalg_determinant
                          'det(0.01*eye(n))==0.01^n')
 
     end subroutine test_${rt[0]}$${rk}$_eye_multiple
-    #:endif
     #:endfor
 
     !> Determinant of complex identity matrix
     #:for ck,ct in CMPLX_KINDS_TYPES
-    #:if ck!="xdp"
     subroutine test_${ct[0]}$${ck}$_complex_determinant(error)
         type(error_type), allocatable, intent(out) :: error
 
@@ -141,7 +134,6 @@ module test_linalg_determinant
 
     end subroutine test_${ct[0]}$${ck}$_complex_determinant
 
-    #:endif
     #:endfor
 
     ! gcc-15 bugfix utility