Fix wrong converge check in VpSqrt (#353)

Author: tompng

ext/bigdecimal/bigdecimal.c

@@ -7261,7 +7261,7 @@ VpSqrt(Real *y, Real *x)
         VpDivd(f, r, x, y);        /* f = x/y    */
         VpAddSub(r, f, y, -1);     /* r = f - y  */
         VpMult(f, VpConstPt5, r);  /* f = 0.5*r  */
-        if (VpIsZero(f))
+        if (y_prec == y->MaxPrec && VpIsZero(f))
             goto converge;
         VpAddSub(r, f, y, 1);      /* r = y + f  */
         VpAsgn(y, r, 1);           /* y = r      */

test/bigdecimal/test_bigdecimal.rb

@@ -1225,6 +1225,8 @@ def test_sqrt_bigdecimal
     x = BigDecimal((2**200).to_s)
     assert_equal(2**100, x.sqrt(1))
 
+    assert_in_delta(BigDecimal("4.0000000000000000000125"), BigDecimal("16.0000000000000000001").sqrt(100), BigDecimal("1e-40"))
+
     BigDecimal.mode(BigDecimal::EXCEPTION_OVERFLOW, false)
     BigDecimal.mode(BigDecimal::EXCEPTION_NaN, false)
     assert_raise_with_message(FloatDomainError, "sqrt of 'NaN'(Not a Number)") { BigDecimal("NaN").sqrt(1) }