I wrote some testcases for catanf(), and found that catanf(1+0*I) return 99999996802856924650656260769173209088+99999996802856924650656260769173209088*I, which should return pi/4 in expected. catanl, catanhf, and catanhl have the same problem. Is this a bug? Any patches existed?