4 #include "gromacs/utility/real.h"
5 #include "../gmx_blas.h"
8 F77_FUNC(snrm2,SNRM2)(int * n__,
13 float ssq,scale,absxi,t;
31 max_ix = 1+(n-1)*(incx);
32 for(ix=1;ix<=max_ix;ix+=incx) {
34 if(fabs(t)>GMX_FLOAT_MIN) {
35 absxi = (t>=0) ? t : (-t);
47 return scale*sqrt(ssq);