scipy: BUG: scipy.linalg.blas.dnrm2 may return error result when incx < 0

Describe your issue.

scipy.linalg.blas.dnrm2 may returns error result when incx<0 : input array:[10, 10, 10, 10, 10, 10, 10, 10, 10] n:5 offx:3 incx:-1 the scipy.linalg.blas.dnrm2 function returns result 20.0

If we calculate the result follow the document in netlib:
https://netlib.org/lapack/explore-html/df/d28/group__single__blas__level1_gab5393665c8f0e7d5de9bd1dd2ff0d9d0.html

we will get result sqrt(500) -> 22.36

Reproducing Code Example

import os
import time
import numpy as np
from scipy import linalg

x = [10,10,10,10,10,10,10,10,10]
n = 5
offx = 3
incx = -1
print(f"x:{x}")
print(f"n:{n}")
print(f"offx:{offx}")
print(f"incx:{incx}") 
res = linalg.blas.dnrm2(x, n, offx, incx)
print(f"res:{res}")

Error message

x:[10, 10, 10, 10, 10, 10, 10, 10, 10]
n:5
offx:3
incx:-1
res:20.0

SciPy/NumPy/Python version information

1.7.3 1.21.5 sys.version_info(major=3, minor=7, micro=5, releaselevel=‘final’, serial=0)