code: def test_precision():
      f1 = Float("1.01")
      f2 = Float("1.0000000000000000000001")
      for u in [
          1,
          1e-2,
          1e-6,
          1e-13,
          1e-14,
          1e-16,
          1e-20,
          1e-100,
          1e-300,
          f1,
          f2,
      ]:
          result = construct_domain([u])
          v = float(result[1][0])
          assert abs(u - v) / u < 1e-14  # Test relative accuracy
  
      result = construct_domain([f1])
      y = result[1][0]
      assert y - 1 > 1e-50
  
      result = construct_domain([f2])
      y = result[1][0]
      assert y - 1 > 1e-50
  
file_path: sympy/sympy/polys/tests/test_constructor.py
imports: from sympy.polys.constructor import construct_domain
  from sympy.polys.domains import ZZ, QQ, ZZ_I, QQ_I, RR, CC, EX
  from sympy.polys.domains.realfield import RealField
  from sympy.polys.domains.complexfield import ComplexField
  from sympy.core import Catalan, GoldenRatio
  from sympy.core.numbers import E, Float, I, Rational, pi
  from sympy.core.singleton import S
  from sympy.functions.elementary.exponential import exp
  from sympy.functions.elementary.miscellaneous import sqrt
  from sympy.functions.elementary.trigonometric import sin
  from sympy.abc import x, y
label_desc: Use image node_id 4 for calling a global function with example usage: test_precision() without return types