35 C_INT a_dim1, a_offset, i__1, i__2;
36 double ret_val, d__1, d__2;
56 a_offset = 1 + a_dim1;
63 for (i__ = 1; i__ <= i__1; ++i__)
75 for (j = jlo; j <= i__2; ++j)
78 sum += (d__1 = a[i1 - j + j * a_dim1], fabs(d__1)) / w[j];
82 d__1 = an, d__2 = sum * w[i__];
double dbnorm_(C_INT *n, double *a, C_INT *nra, C_INT *ml, C_INT *mu, double *w)