35 #ifndef TEMPLATE_BLAS_NRM2_HEADER
36 #define TEMPLATE_BLAS_NRM2_HEADER
49 Treal norm, scale, absxi;
61 if (*n < 1 || *incx < 1) {
70 i__1 = (*n - 1) * *incx + 1;
72 for (ix = 1; i__2 < 0 ? ix >= i__1 : ix <= i__1; ix += i__2) {
74 absxi = (d__1 = x[ix],
absMACRO(d__1));
78 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:40
#define absMACRO(x)
Definition: template_blas_common.h:45
int integer
Definition: template_blas_common.h:38
Treal template_blas_sqrt(Treal x)