# HG changeset patch
# User William Stein <wstein@gmail.com>
# Date 1272693754 25200
# Node ID d7696684971481511b41a5d1943ecb677f27aafa
# Parent d26c205cafef57b8f3737df808ee7c9f7cb54a32
trac 8833 -- fix some numerical noise issues in riemanny.pyx
diff --git a/sage/calculus/riemann.pyx b/sage/calculus/riemann.pyx
|
a
|
b
|
|
| 488 | 488 | sage: m.riemann_map(0.25 + sqrt(-0.5)) # long time |
| 489 | 489 | (0.137514137885...+0.876696023004...j) |
| 490 | 490 | sage: m.riemann_map(1.3*I) # long time |
| 491 | | (-1.5610293963...+0.989694535737...j) |
| | 491 | (-1.561029396...+0.989694535737...j) |
| 492 | 492 | sage: I = CDF.gen() # long time |
| 493 | 493 | sage: m.riemann_map(0.4) # long time |
| 494 | | (0.733242677182...+3.5076771462...j) |
| | 494 | (0.733242677182...+3.50767714...j) |
| 495 | 495 | sage: import numpy as np # long time |
| 496 | 496 | sage: m.riemann_map(np.complex(-3, 0.0001)) # long time |
| 497 | | (1.40575713549...+1.05106...j) |
| | 497 | (1.405757135...+1.05106...j) |
| 498 | 498 | """ |
| 499 | 499 | pt1 = np.complex(pt) |
| 500 | 500 | cdef np.ndarray[double complex, ndim=1] q_vector = 1 / ( |
| … |
… |
|
| 564 | 564 | sage: m.inverse_riemann_map(0.5 + sqrt(-0.5)) # long time |
| 565 | 565 | (0.406880548363...+0.361470279816...j) |
| 566 | 566 | sage: m.inverse_riemann_map(0.95) # long time |
| 567 | | (0.486319431795...-4.900190524...j) |
| | 567 | (0.486319431795...-4.90019052...j) |
| 568 | 568 | sage: m.inverse_riemann_map(0.25 - 0.3*I) # long time |
| 569 | 569 | (0.165324498558...-0.180936785500...j) |
| 570 | 570 | sage: import numpy as np # long time |