remove xdp restriction

Author: perazz

doc/specs/stdlib_linalg.md

@@ -104,7 +104,8 @@ end interface axpy
 Note that the 128-bit functions are only provided by `stdlib` and always point to the internal implementation. 
 Because 128-bit precision is identified as [stdlib_kinds(module):qp], initials for 128-bit procedures were 
 labelled as `q` (quadruple-precision reals) and `w` ("wide" or quadruple-precision complex numbers). 
-Extended precision ([stdlib_kinds(module):xdp]) calculations are currently not supported.
+Extended precision ([stdlib_kinds(module):xdp]) calculations are labelled as `x` (extended-precision reals).
+and `y` (extended-precision complex numbers).
 
 ### Example
 

src/stdlib_linalg.fypp

@@ -567,10 +567,8 @@ module stdlib_linalg
     !! the state flag `err` is not provided. 
     !!
     !!@note The provided functions are intended for square matrices.
-    !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
     !!       
     #:for rk,rt,ri in RC_KINDS_TYPES
-    #:if rk!="xdp" 
     module function stdlib_linalg_inverse_${ri}$(a,err) result(inva)
          !> Input matrix a[n,n]
          ${rt}$, intent(in) :: a(:,:)
@@ -579,7 +577,6 @@ module stdlib_linalg
          !> [optional] state return flag. On error if not requested, the code will stop
          type(linalg_state_type), optional, intent(out) :: err
     end function stdlib_linalg_inverse_${ri}$
-    #:endif
     #:endfor
   end interface inv
 
@@ -605,10 +602,8 @@ module stdlib_linalg
     !! work spaces are provided, no internal memory allocations take place when using this interface.             
     !!
     !!@note The provided subroutines are intended for square matrices.
-    !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
     !!      
     #:for rk,rt,ri in RC_KINDS_TYPES
-    #:if rk!="xdp" 
     module subroutine stdlib_linalg_invert_inplace_${ri}$(a,pivot,err)
          !> Input matrix a[n,n]
          ${rt}$, intent(inout) :: a(:,:)
@@ -628,7 +623,6 @@ module stdlib_linalg
          !> [optional] state return flag. On error if not requested, the code will stop
          type(linalg_state_type), optional, intent(out) :: err    
     end subroutine stdlib_linalg_invert_split_${ri}$
-    #:endif
     #:endfor
   end interface invert
 
@@ -649,17 +643,14 @@ module stdlib_linalg
     !! NaNs will be returned.
     !!
     !!@note The provided functions are intended for square matrices.
-    !!@note BLAS/LAPACK backends do not currently support extended precision (``xdp``).
     !!
     #:for rk,rt,ri in RC_KINDS_TYPES
-    #:if rk!="xdp" 
     module function stdlib_linalg_inverse_${ri}$_operator(a) result(inva)
          !> Input matrix a[n,n]
          ${rt}$, intent(in) :: a(:,:)
          !> Result matrix
          ${rt}$, allocatable :: inva(:,:)        
     end function stdlib_linalg_inverse_${ri}$_operator
-    #:endif
     #:endfor
   end interface operator(.inv.)
 

src/stdlib_linalg_inverse.fypp

@@ -33,7 +33,6 @@ submodule (stdlib_linalg) stdlib_linalg_inverse
      end subroutine handle_getri_info
 
      #:for rk,rt,ri in RC_KINDS_TYPES
-     #:if rk!="xdp" 
      ! Compute the in-place square matrix inverse of a
      module subroutine stdlib_linalg_invert_inplace_${ri}$(a,pivot,err)
          !> Input matrix a[n,n]. On return, A is destroyed and replaced by the inverse
@@ -175,7 +174,6 @@ submodule (stdlib_linalg) stdlib_linalg_inverse
 
      end function stdlib_linalg_inverse_${ri}$_operator
 
-     #:endif
      #:endfor
 
 end submodule stdlib_linalg_inverse

test/linalg/test_linalg_inverse.fypp

@@ -22,17 +22,14 @@ module test_linalg_inverse
         allocate(tests(0))
 
         #:for rk,rt,ri in RC_KINDS_TYPES
-        #:if rk!="xdp"
         tests = [tests,new_unittest("${ri}$_eye_inverse",test_${ri}$_eye_inverse), &
                        new_unittest("${ri}$_singular_inverse",test_${ri}$_singular_inverse), &
                        new_unittest("${ri}$_random_spd_inverse",test_${ri}$_random_spd_inverse)]
-        #:endif
         #:endfor
 
     end subroutine test_inverse_matrix
 
     #:for rk,rt,ri in REAL_KINDS_TYPES
-    #:if rk!="xdp"
     !> Invert real identity matrix
     subroutine test_${ri}$_eye_inverse(error)
         type(error_type), allocatable, intent(out) :: error
@@ -139,12 +136,10 @@ module test_linalg_inverse
 
     end subroutine test_${ri}$_random_spd_inverse
 
-    #:endif
     #:endfor
 
     !> Invert complex identity matrix
     #:for ck,ct,ci in CMPLX_KINDS_TYPES
-    #:if ck!="xdp"
     subroutine test_${ci}$_eye_inverse(error)
         type(error_type), allocatable, intent(out) :: error
 
@@ -295,7 +290,6 @@ module test_linalg_inverse
 
     end subroutine test_${ci}$_singular_inverse
 
-    #:endif
     #:endfor
 
 end module test_linalg_inverse