Improving the Accuracy of the Fast Inverse Square Root by Modifying Newton–Raphson Corrections

Walczyk, Cezary J.; Moroz, Leonid; Cieśliński, Jan L.

doi:10.3390/e23010086

Cited by 11 publications

(14 citation statements)

References 25 publications

Supporting

Mentioning

Contrasting

Order By: Relevance

“…Here, we briefly outline the main concepts of the basic FISR method used in the paper, which are defined in [28,29]. Suppose that we have a positive normalized FP number…”

Section: Brief Theory Of the Fisr Methodsmentioning

confidence: 99%

“…As proved in [28,29], in order to find the behavior of the relative error when calculating 0 y over the whole range of normalized FP numbers, it is sufficient to describe it in the range…”

Section: Brief Theory Of the Fisr Methodsmentioning

confidence: 99%

“…The accuracy of the RcpSqrt2 algorithm after the first iteration is 10.15 bits. Another similar algorithm from [29] has even better accuracy (20.97 bits) but requires one extra multiplication compared to Algorithms 1 and 2. return y; 10: } Our other recent modifications of the FISR algorithm are discussed in [26] and [32]. In both papers, besides the NR method, a modified second-order Householder method is used to improve the accuracy of the algorithms-in the former case, in the last iteration of the algorithm, and in the latter case, in the first iteration.…”

Section: State-of-the-art Fisr Algorithmsmentioning

confidence: 99%

“…Among the better-known methods of calculating the square root and reciprocal square root [1,2,15,20], the FISR method [20,[25][26][27][28][29] has recently gained increasing popularity in SW [8,20,27,[29][30][31][32] and HW [3,[33][34][35][36][37] applications. The algorithm was proposed for the first time in [24] but gained wider popularity through its use in the computer game Quake III Arena [27].…”

Section: Introductionmentioning

confidence: 99%

See 3 more Smart Citations

Modified Fast Inverse Square Root and Square Root Approximation Algorithms: The Method of Switching Magic Constants

2021

Self Cite

View full text Add to dashboard Cite

Many low-cost platforms that support floating-point arithmetic, such as microcontrollers and field-programmable gate arrays, do not include fast hardware or software methods for calculating the square root and/or reciprocal square root. Typically, such functions are implemented using direct lookup tables or polynomial approximations, with a subsequent application of the Newton–Raphson method. Other, more complex solutions include high-radix digit-recurrence and bipartite or multipartite table-based methods. In contrast, this article proposes a simple modification of the fast inverse square root method that has high accuracy and relatively low latency. Algorithms are given in C/C++ for single- and double-precision numbers in the IEEE 754 format for both square root and reciprocal square root functions. These are based on the switching of magic constants in the initial approximation, depending on the input interval of the normalized floating-point numbers, in order to minimize the maximum relative error on each subinterval after the first iteration—giving 13 correct bits of the result. Our experimental results show that the proposed algorithms provide a fairly good trade-off between accuracy and latency after two iterations for numbers of type float, and after three iterations for numbers of type double when using fused multiply–add instructions—giving almost complete accuracy.

show abstract

“…Here, we briefly outline the main concepts of the basic FISR method used in the paper, which are defined in [28,29]. Suppose that we have a positive normalized FP number…”

Section: Brief Theory Of the Fisr Methodsmentioning

confidence: 99%

“…As proved in [28,29], in order to find the behavior of the relative error when calculating 0 y over the whole range of normalized FP numbers, it is sufficient to describe it in the range…”

Section: Brief Theory Of the Fisr Methodsmentioning

confidence: 99%

Section: State-of-the-art Fisr Algorithmsmentioning

confidence: 99%

Section: Introductionmentioning

confidence: 99%

See 2 more Smart Citations

Modified Fast Inverse Square Root and Square Root Approximation Algorithms: The Method of Switching Magic Constants

2021

Self Cite

View full text Add to dashboard Cite

show abstract

“…To further improve the accuracy of the Algorithm 1, we can apply the additive correction of the parameter k 2 for each iteration, as proposed in [10,20,21]. In other words, the first iteration should be performed according to the modified Newton-Raphson formulas…”

Section: The Theory Of the Fast Inverse Cube Root Methods And Proposed Algorithmsmentioning

confidence: 99%

Fast Calculation of Cube and Inverse Cube Roots Using a Magic Constant and Its Implementation on Microcontrollers

et al. 2021

Self Cite

View full text Add to dashboard Cite

We develop a bit manipulation technique for single precision floating point numbers which leads to new algorithms for fast computation of the cube root and inverse cube root. It uses the modified iterative Newton–Raphson method (the first order of convergence) and Householder method (the second order of convergence) to increase the accuracy of the results. The proposed algorithms demonstrate high efficiency and reduce error several times in the first iteration in comparison with known algorithms. After two iterations 22.84 correct bits were obtained for single precision. Experimental tests showed that our novel algorithm is faster and more accurate than library functions for microcontrollers.

show abstract

A class of high-order improved fast weighted essentially non-oscillatory schemes for achieving optimal order at any critical points

et al. 2022

View full text Add to dashboard Cite

The weighted essentially non-oscillatory (WENO) scheme is one of the most popular shock-capturing schemes, and constructing a more efficient and higher-order WENO scheme has always been an intention of optimization design. In the general WENO reconstruction framework, the smoothness indicator plays an important role in identifying whether the sub-stencils are in discontinuous or smooth regions. However, the classical smoothness indicator is the most expensive one in the whole reconstruction algorithm, and its computational complexity increases sharply with the improvement of the accuracy order. Therefore, a class of efficient and superior WENO schemes called improved fast WENO (IFWENO) are proposed based on the fast WENO (FWENO). To improve efficiency, the smoothness indicator of the IFWENO scheme is simplified from the traditional version, and the nonlinear weight calculation method is modified. The parameter ε is carefully designed to obtain the superior property that the accuracy of the spatial derivatives will not degrade at any order critical point in smooth regions. The reason for the instability occurring in the high-order FWENO is revealed, and the parameter p is likewise specifically selected to improve robustness at discontinuities. The excellent multi-scale resolution of the proposed IFWENO scheme is proven by the theoretical analyses and the numerical experiments. Through several typical examples, the consistently high accuracy and efficiency of the designed scheme in both smooth and discontinuous regions are verified.

show abstract

Improving the Accuracy of the Fast Inverse Square Root by Modifying Newton–Raphson Corrections

Cited by 11 publications

References 25 publications

Modified Fast Inverse Square Root and Square Root Approximation Algorithms: The Method of Switching Magic Constants

Modified Fast Inverse Square Root and Square Root Approximation Algorithms: The Method of Switching Magic Constants

Fast Calculation of Cube and Inverse Cube Roots Using a Magic Constant and Its Implementation on Microcontrollers

A class of high-order improved fast weighted essentially non-oscillatory schemes for achieving optimal order at any critical points

Contact Info

Product

Resources

About