Make use of stdlib_constants to define the various zeros and ones.

Author: loiseaujc

src/stdlib_specialmatrices.fypp

@@ -9,29 +9,30 @@ module stdlib_specialmatrices
     !! of partial differential equations, as well as control and signal processing
     !! applications. ([Specifications](../page/specs/stdlib_specialmatrices.html))
     use stdlib_linalg_constants
+    use stdlib_constants
     implicit none
     private
     public :: tridiagonal
     public :: spmv
     public :: dense, transpose, hermitian
     public :: operator(*), operator(+), operator(-)
 
-    real(sp), parameter :: zero_sp = 0._sp
-    real(sp), parameter :: one_sp = 1._sp
-    real(dp), parameter :: zero_dp = 0._dp
-    real(dp), parameter :: one_dp = 1._dp
-    real(xdp), parameter :: zero_xdp = 0._xdp
-    real(xdp), parameter :: one_xdp = 1._xdp
-    real(qp), parameter :: zero_qp = 0._qp
-    real(qp), parameter :: one_qp = 1._qp
-    complex(sp), parameter :: zero_csp = (0._sp,0._sp)
-    complex(sp), parameter :: one_csp = (1._sp,1._sp)
-    complex(dp), parameter :: zero_cdp = (0._dp,0._dp)
-    complex(dp), parameter :: one_cdp = (1._dp,1._dp)
-    complex(xdp), parameter :: zero_cxdp = (0._xdp,0._xdp)
-    complex(xdp), parameter :: one_cxdp = (1._xdp,1._xdp)
-    complex(qp), parameter :: zero_cqp = (0._qp,0._qp)
-    complex(qp), parameter :: one_cqp = (1._qp,1._qp)
+    ! real(sp), parameter :: zero_sp = 0._sp
+    ! real(sp), parameter :: one_sp = 1._sp
+    ! real(dp), parameter :: zero_dp = 0._dp
+    ! real(dp), parameter :: one_dp = 1._dp
+    ! real(xdp), parameter :: zero_xdp = 0._xdp
+    ! real(xdp), parameter :: one_xdp = 1._xdp
+    ! real(qp), parameter :: zero_qp = 0._qp
+    ! real(qp), parameter :: one_qp = 1._qp
+    ! complex(sp), parameter :: zero_csp = (0._sp,0._sp)
+    ! complex(sp), parameter :: one_csp = (1._sp,1._sp)
+    ! complex(dp), parameter :: zero_cdp = (0._dp,0._dp)
+    ! complex(dp), parameter :: one_cdp = (1._dp,1._dp)
+    ! complex(xdp), parameter :: zero_cxdp = (0._xdp,0._xdp)
+    ! complex(xdp), parameter :: one_cxdp = (1._xdp,1._xdp)
+    ! complex(qp), parameter :: zero_cqp = (0._qp,0._qp)
+    ! complex(qp), parameter :: one_cqp = (1._qp,1._qp)
 
     !--------------------------------------
     !-----                           ------