Ticket #11626 (closed enhancement: duplicate)

Opened 22 months ago

Last modified 21 months ago

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: Work issues:
Report Upstream: N/A Reviewers: Paul Zimmermann, Karl-Dieter Crisman
Authors: Merged in:
Dependencies: Stopgaps:

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

comment:1 follow-up: ↓ 2 Changed 22 months 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 22 months 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 22 months ago by kcrisman

  • Status changed from new to needs_review
  • Reviewers set to Paul Zimmermann, Karl-Dieter Crisman
  • Milestone changed from sage-4.7.2 to sage-duplicate/invalid/wontfix

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

To release manager: positive review as duplicate.

comment:4 Changed 22 months ago by kcrisman

  • Status changed from needs_review to positive_review

comment:5 Changed 21 months ago by jdemeyer

  • Status changed from positive_review to closed
  • Resolution set to duplicate
Note: See TracTickets for help on using tickets.