I expect that the method is the same as for other polynomial or analytic functions. If you start with exact data, the results are symbolic Root objects and can be computed with arbitrary precision.
Root