@@ -144,6 +144,7 @@ submodule(stdlib_linalg) stdlib_linalg_norms
144
144
145
145
integer(ilp) :: sze,norm_request
146
146
real(${rk}$) :: rorder
147
+ intrinsic :: abs, sum, sqrt, norm2, maxval, minval, conjg
147
148
148
149
sze = size(a,kind=ilp)
149
150
@@ -171,7 +172,7 @@ submodule(stdlib_linalg) stdlib_linalg_norms
171
172
#:if rt.startswith('complex')
172
173
nrm = sqrt( real( sum( a * conjg(a) ), ${rk}$) )
173
174
#:else
174
- nrm = sqrt( sum( a ** 2 ) )
175
+ nrm = norm2( a )
175
176
#:endif
176
177
case(NORM_INF)
177
178
nrm = maxval( abs(a) )
@@ -244,6 +245,7 @@ submodule(stdlib_linalg) stdlib_linalg_norms
244
245
type(linalg_state_type) :: err_
245
246
integer(ilp) :: sze,norm_request
246
247
real(${rk}$) :: rorder
248
+ intrinsic :: abs, sum, sqrt, norm2, maxval, minval, conjg
247
249
248
250
sze = size(a,kind=ilp)
249
251
@@ -278,7 +280,7 @@ submodule(stdlib_linalg) stdlib_linalg_norms
278
280
#:if rt.startswith('complex')
279
281
nrm = sqrt( real( sum( a * conjg(a) , dim = dim ), ${rk}$) )
280
282
#:else
281
- nrm = sqrt( sum( a ** 2 , dim = dim ) )
283
+ nrm = norm2( a , dim = dim )
282
284
#:endif
283
285
case(NORM_INF)
284
286
nrm = maxval( abs(a) , dim = dim )
0 commit comments