The result
is calculated using Smith's algorithm with scaling, from
Li et al. (2002), which ensures that no unnecessary overflow or underflow occurs at intermediate stages of the computation.
Li X S, Demmel J W, Bailey D H, Henry G, Hida Y, Iskandar J, Kahan W, Kapur A, Martin M C, Tung T and Yoo D J (2002) Design, implementation and testing of extended and mixed precision BLAS ACM Trans. Math. Soft. 28(2) 152–205
None.
Background information to multithreading can be found in the
Multithreading documentation.