Merge pull request #228 from JuliaMath/aa/hypotl

Fix incorrect results in `hypotl` near underflow
JuliaMath · Feb 17, 2021 · f052f42 · f052f42
2 parents aeab19f + 711654e
commit f052f42
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/e_hypotl.c b/src/e_hypotl.c
@@ -82,7 +82,7 @@ hypotl(long double x, long double y)
 	        man_t manh, manl;
 		GET_LDBL_MAN(manh,manl,b);
 		if((manh|manl)==0) return a;
-		t1=0;
+		t1=1;
 		SET_HIGH_WORD(t1,ESW(MAX_EXP-2));	/* t1=2^(MAX_EXP-2) */
 		b *= t1;
 		a *= t1;