Calculation of nuclear magnetic shieldings: infinite-order Foldy–Wouthuysen transformation