Opened 12 years ago

Closed 11 years ago

#11626 closed enhancement (duplicate)

make the error function work in arbitrary precision

Reported by: zimmerma Owned by: AlexGhitza
Priority: minor Milestone: sage-duplicate/invalid/wontfix
Component: basic arithmetic Keywords:
Cc: Merged in:
Authors: Reviewers: Paul Zimmermann, Karl-Dieter Crisman
Report Upstream: N/A Work issues:
Branch: Commit:
Dependencies: Stopgaps:

GitHub link to the corresponding issue

Description 

sage: n(erf(2),54)
...
NotImplementedError: erf not implemented for precision higher than 53

However, MPFR contains a multiple precision implementation of erf. It should be straightforward to use.

Paul

PS: maybe this ticket is duplicate. It is not easy to search for erf in trac, a feature, like grep -w erf is missing, and the TracSearch help is not very useful...

Change History (5)

comment:1 Changed 12 years ago by kcrisman

Does #1173 solve this (it uses mpmath)? It's certainly possible there to do erf(2).n(54) or such, not sure about your syntax as I don't have a Sage install with this patch currently. But this should do it.

You are certainly correct about searching Trac. Sorry that neither this ticket nor #1173 can address that :)

comment:2 in reply to:  1 Changed 11 years ago by zimmerma

Replying to kcrisman:

Does #1173 solve this (it uses mpmath)?

yes it does: 

sage: n(erf(2),100)
0.99532226501895273416206925637
sage: erf(2).n(100)
0.99532226501895273416206925637

Thus as soon as #1173 is closed, we can close this ticket too.

Paul

comment:3 Changed 11 years ago by kcrisman

Milestone: sage-4.7.2 → sage-duplicate/invalid/wontfix
Reviewers: → Paul Zimmermann, Karl-Dieter Crisman
Status: new → needs_review

This is then now part of #1173. Putting this as duplicate.

To release manager: positive review as duplicate.

comment:4 Changed 11 years ago by kcrisman

Status: needs_review → positive_review

comment:5 Changed 11 years ago by jdemeyer

Resolution: → duplicate
Status: positive_review → closed
Note: See TracTickets for help on using tickets.