37 #ifndef TEMPLATE_BLAS_NRM2_HEADER 38 #define TEMPLATE_BLAS_NRM2_HEADER 51 Treal norm, scale, absxi;
63 if (*n < 1 || *incx < 1) {
72 i__1 = (*n - 1) * *incx + 1;
74 for (ix = 1; i__2 < 0 ? ix >= i__1 : ix <= i__1; ix += i__2) {
76 absxi = (d__1 = x[ix],
absMACRO(d__1));
80 ssq = ssq * (d__1 * d__1) + 1.;
Treal template_blas_nrm2(const integer *n, const Treal *x, const integer *incx)
Definition: template_blas_nrm2.h:42
#define absMACRO(x)
Definition: template_blas_common.h:47
int integer
Definition: template_blas_common.h:40
Treal template_blas_sqrt(Treal x)