@@ -6,7 +6,7 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
66 contains
77#:for ik,it,ii in LINALG_INT_KINDS_TYPES
88
9- pure real(sp) module function stdlib${ii}$_slamch( cmach )
9+ pure real(sp) module function stdlib${ii}$_slamch( cmach )
1010 !! SLAMCH determines single precision machine parameters.
1111 ! -- lapack auxiliary routine --
1212 ! -- lapack is a software package provided by univ. of tennessee, --
@@ -17,16 +17,10 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
1717 ! =====================================================================
1818
1919 ! Local Scalars
20- real(sp) :: rnd, eps, sfmin, small, rmach
20+ real(sp) :: sfmin, small, rmach
2121 ! Intrinsic Functions
2222 ! Executable Statements
2323 ! assume rounding, not chopping. always.
24- rnd = one
25- if( one==rnd ) then
26- eps = epsilon(zero) * 0.5
27- else
28- eps = epsilon(zero)
29- end if
3024 if( stdlib_lsame( cmach, 'E' ) ) then
3125 rmach = eps
3226 else if( stdlib_lsame( cmach, 'S' ) ) then
@@ -45,7 +39,7 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
4539 else if( stdlib_lsame( cmach, 'N' ) ) then
4640 rmach = digits(zero)
4741 else if( stdlib_lsame( cmach, 'R' ) ) then
48- rmach = rnd
42+ rmach = one
4943 else if( stdlib_lsame( cmach, 'M' ) ) then
5044 rmach = minexponent(zero)
5145 else if( stdlib_lsame( cmach, 'U' ) ) then
@@ -72,16 +66,10 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
7266 ! =====================================================================
7367
7468 ! Local Scalars
75- real(dp) :: rnd, eps, sfmin, small, rmach
69+ real(dp) :: sfmin, small, rmach
7670 ! Intrinsic Functions
7771 ! Executable Statements
7872 ! assume rounding, not chopping. always.
79- rnd = one
80- if( one==rnd ) then
81- eps = epsilon(zero) * 0.5
82- else
83- eps = epsilon(zero)
84- end if
8573 if( stdlib_lsame( cmach, 'E' ) ) then
8674 rmach = eps
8775 else if( stdlib_lsame( cmach, 'S' ) ) then
@@ -100,7 +88,7 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
10088 else if( stdlib_lsame( cmach, 'N' ) ) then
10189 rmach = digits(zero)
10290 else if( stdlib_lsame( cmach, 'R' ) ) then
103- rmach = rnd
91+ rmach = one
10492 else if( stdlib_lsame( cmach, 'M' ) ) then
10593 rmach = minexponent(zero)
10694 else if( stdlib_lsame( cmach, 'U' ) ) then
@@ -129,16 +117,10 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
129117 ! =====================================================================
130118
131119 ! Local Scalars
132- real(${rk}$) :: rnd, eps, sfmin, small, rmach
120+ real(${rk}$) :: sfmin, small, rmach
133121 ! Intrinsic Functions
134122 ! Executable Statements
135123 ! assume rounding, not chopping. always.
136- rnd = one
137- if( one==rnd ) then
138- eps = epsilon(zero) * 0.5
139- else
140- eps = epsilon(zero)
141- end if
142124 if( stdlib_lsame( cmach, 'E' ) ) then
143125 rmach = eps
144126 else if( stdlib_lsame( cmach, 'S' ) ) then
@@ -157,7 +139,7 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
157139 else if( stdlib_lsame( cmach, 'N' ) ) then
158140 rmach = digits(zero)
159141 else if( stdlib_lsame( cmach, 'R' ) ) then
160- rmach = rnd
142+ rmach = one
161143 else if( stdlib_lsame( cmach, 'M' ) ) then
162144 rmach = minexponent(zero)
163145 else if( stdlib_lsame( cmach, 'U' ) ) then
@@ -182,7 +164,7 @@ submodule(stdlib_lapack_base) stdlib_lapack_auxiliary_parameters
182164
183165
184166
185- pure real(sp) module function stdlib${ii}$_slamc3( a, b )
167+ pure real(sp) module function stdlib${ii}$_slamc3( a, b )
186168 ! -- lapack auxiliary routine --
187169 ! univ. of tennessee, univ. of california berkeley and nag ltd..
188170 use stdlib_blas_constants_sp, only: negone, zero, half, one, two, three, four, eight, ten, czero, chalf, cone, cnegone
0 commit comments