FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x) { BLASLONG i=0; BLASLONG ix=0; FLOAT maxf[2]; BLASLONG max=0; BLASLONG inc_x2; if (n < 0 || inc_x < 1 ) return(0.0); inc_x2 = 2 * inc_x; maxf[0] = ABS(x[ix]); maxf[1] = ABS(x[ix+1]); while(i < n) { if( CABS1(x,ix) > CABS1(maxf,0) ) { max = i; maxf[0] = ABS(x[ix]); maxf[1] = ABS(x[ix+1]); } ix += inc_x2; i++; } return(CABS1(maxf,0)); }
FLOAT CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x) { BLASLONG i=0; BLASLONG ix=0; FLOAT maxf; BLASLONG inc_x2; if (n <= 0 || inc_x <= 0) return(0.0); inc_x2 = 2 * inc_x; maxf = CABS1(x,0); ix += inc_x2; i++; while(i < n) { if( CABS1(x,ix) > maxf ) { maxf = CABS1(x,ix); } ix += inc_x2; i++; } return(maxf); }
BLASLONG CNAME(BLASLONG n, FLOAT *x, BLASLONG inc_x) { BLASLONG i=0; BLASLONG ix=0; FLOAT minf; BLASLONG min=0; BLASLONG inc_x2; if (n <= 0 || inc_x <= 0) return(min); inc_x2 = 2 * inc_x; minf = CABS1(x,0); ix += inc_x2; i++; while(i < n) { if( CABS1(x,ix) < minf ) { min = i; minf = CABS1(x,ix); } ix += inc_x2; i++; } return(min+1); }