3 doublereal dlapy2_(doublereal *x, doublereal *y)
5 /* System generated locals */
6 doublereal ret_val, d__1;
8 /* Builtin functions */
9 double sqrt(doublereal);
12 doublereal w, z__, xabs, yabs;
15 /* -- LAPACK auxiliary routine (version 3.1) -- */
16 /* Univ. of Tennessee, Univ. of California Berkeley and NAG Ltd.. */
19 /* .. Scalar Arguments .. */
25 /* DLAPY2 returns sqrt(x**2+y**2), taking care not to cause unnecessary */
31 /* X (input) DOUBLE PRECISION */
32 /* Y (input) DOUBLE PRECISION */
33 /* X and Y specify the values x and y. */
35 /* ===================================================================== */
37 /* .. Parameters .. */
39 /* .. Local Scalars .. */
41 /* .. Intrinsic Functions .. */
43 /* .. Executable Statements .. */
52 /* Computing 2nd power */
54 ret_val = w * sqrt(d__1 * d__1 + 1.);