I did use the following test function.

def check(re, im):
    precs = [6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 31, 32, 33, 100, 127, 128, 129]
    for prec in precs:
        c = ~ComplexIntervalField(prec)(re, im)
        norm2 = re**2 + im**2

        # check real part
        inv = re / norm2
        left, right = c.real().endpoints()
        left = QQ(left.as_integer_ratio())
        right = QQ(right.as_integer_ratio())
        if left > inv or right < inv:
             print('wrong real at prec={}: left={}~{} right={}~{} value={}~{}'.format(prec, left, left.n(prec=2*prec), right, right.n(prec=2*prec), inv, inv.n(prec=2*prec)))

        # check imaginary part
        inv = - im / norm2
        left, right = c.imag().endpoints()
        left = QQ(left.as_integer_ratio())
        right = QQ(right.as_integer_ratio())
        if left > inv or right < inv:
            print('wrong imag at prec={}: left={}~{} right={}~{} value={}~{}'.format(prec, left, left.n(prec=2*prec), right, right.n(prec=2*prec), inv, inv.n(prec=2*prec)))

It seems to be the case that the tight inversion has

correct real part always
wrong imaginary part often

sage: check(1/2, 1/3)
wrong imag at prec=11: left=-473/512~-0.923828 right=-1891/2048~-0.923340 value=-12/13~-0.923077
wrong imag at prec=12: left=-1891/2048~-0.923340 right=-3781/4096~-0.923096 value=-12/13~-0.923077
wrong imag at prec=17: left=-120991/131072~-0.923088074 right=-60495/65536~-0.923080444 value=-12/13~-0.923076923
sage: check(1, 5)
wrong imag at prec=6: left=-25/128~-0.195 right=-25/128~-0.195 value=-5/26~-0.192
wrong imag at prec=7: left=-99/512~-0.193 right=-99/512~-0.193 value=-5/26~-0.192
wrong imag at prec=8: left=-197/1024~-0.1924 right=-197/1024~-0.1924 value=-5/26~-0.1923
wrong imag at prec=9: left=-197/1024~-0.19238 right=-197/1024~-0.19238 value=-5/26~-0.19231
wrong imag at prec=10: left=-197/1024~-0.19238 right=-197/1024~-0.19238 value=-5/26~-0.19231
wrong imag at prec=11: left=-197/1024~-0.192383 right=-197/1024~-0.192383 value=-5/26~-0.192308
wrong imag at prec=12: left=-3151/16384~-0.192322 right=-3151/16384~-0.192322 value=-5/26~-0.192308
wrong imag at prec=13: left=-3151/16384~-0.1923218 right=-3151/16384~-0.1923218 value=-5/26~-0.1923077
wrong imag at prec=14: left=-3151/16384~-0.19232178 right=-3151/16384~-0.19232178 value=-5/26~-0.19230769
...
sage: check(1/4, -1/2)
wrong imag at prec=6: left=13/8~1.62 right=13/8~1.62 value=8/5~1.60
wrong imag at prec=7: left=103/64~1.61 right=103/64~1.61 value=8/5~1.60
wrong imag at prec=8: left=205/128~1.602 right=205/128~1.602 value=8/5~1.600
wrong imag at prec=9: left=205/128~1.6016 right=205/128~1.6016 value=8/5~1.6000
wrong imag at prec=10: left=205/128~1.6016 right=205/128~1.6016 value=8/5~1.6000
...

Uh oh!

Reintroduce tight enclosure of interval inversion #37963

Description

Metadata

Metadata

Assignees

Labels

Type

Projects

Milestone

Relationships

Development

Issue actions