In GLSL there\'s rudimentary support for double precision variables and operations which can be found here. However they also mention \"Double-precision versions of angle, t
this link seem's to be the best answer
So yes, you'll need to make your own implementation for those functions.
glibc source may be your friend.