Opened 6 years ago

Closed 6 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:

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 follow-up: Changed 6 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 6 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 6 years ago by kcrisman

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

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

To release manager: positive review as duplicate.

comment:4 Changed 6 years ago by kcrisman

  • Status changed from needs_review to positive_review

comment:5 Changed 6 years ago by jdemeyer

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