I have been trying to find out if the GLSL specification (as per OpenGL 4.3) specifies any precision or error bounds for fp64 operations. For fp32 (float), the errors are specified as ULPs. I am not sure if the same error bounds apply for fp64. I am hoping tighter bounds are specified for fp64 than fp32.