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