Implementation of the reciprocal square root in MPFR

Author Paul Zimmermann

Thumbnail PDF


  • Filesize: 159 kB
  • 3 pages

Document Identifiers

Author Details

Paul Zimmermann

Cite AsGet BibTex

Paul Zimmermann. Implementation of the reciprocal square root in MPFR. In Numerical Validation in Current Hardware Architectures. Dagstuhl Seminar Proceedings, Volume 8021, pp. 1-3, Schloss Dagstuhl – Leibniz-Zentrum für Informatik (2008)


We describe the implementation of the reciprocal square root --- also called inverse square root --- as a native function in the MPFR library. The difficulty is to implement Newton's iteration for the reciprocal square root on top's of GNU MP's extsc{mpn} layer, while guaranteeing a rigorous $1/2$ ulp bound on the roundoff error.
  • Multiple precision
  • floating-point
  • inverse square root
  • correct rounding
  • MPFR library


  • Access Statistics
  • Total Accesses (updated on a weekly basis)
    PDF Downloads
Questions / Remarks / Feedback

Feedback for Dagstuhl Publishing

Thanks for your feedback!

Feedback submitted

Could not send message

Please try again later or send an E-mail