Given `a + bim`, we could just return `Base.sqrt(a + bim)` since it can handle those values without throwing domain errors. Currently, `NaNMath.sqrt(a + bim)` will just yield a missing method error.