handle incx < 0

This commit is contained in:
Martin Kroeker 2024-03-12 20:44:01 +01:00 committed by GitHub
parent db70c7f7fb
commit 3752e73919
No known key found for this signature in database
GPG Key ID: B5690EEEBB952194
1 changed files with 2 additions and 2 deletions

View File

@ -101,7 +101,7 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x)
{ {
BLASLONG i=0; BLASLONG i=0;
if (n <= 0 || inc_x <= 0) return(0.0); if (n <= 0 || inc_x == 0) return(0.0);
if(n == 1) return (ABS(x[0])); if(n == 1) return (ABS(x[0]));
unsigned int gvl = 0; unsigned int gvl = 0;
@ -190,7 +190,7 @@ FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x)
//finish any tail using scalar ops //finish any tail using scalar ops
i*=gvl*inc_x; i*=gvl*inc_x;
n*=inc_x; n*=inc_x;
while(i < n){ while(abs(i) < abs(n)){
if ( x[i] != 0.0 ){ if ( x[i] != 0.0 ){
FLOAT absxi = ABS( x[i] ); FLOAT absxi = ABS( x[i] );
if ( scale < absxi ){ if ( scale < absxi ){