https://gcc.gnu.org/bugzilla/show_bug.cgi?id=93871
--- Comment #5 from Thomas Henlich <thenlich at gcc dot gnu.org> --- Created attachment 47884 --> https://gcc.gnu.org/bugzilla/attachment.cgi?id=47884&action=edit Test case Output: th@THe-Surface:~$ /opt/gcc/bin/gfortran -L/opt/gcc/lib64 -Wl,-rpath -Wl,/opt/gcc/lib64 -fdec-math /mnt/c/Temp/test_cotan_complex.f90 && ./a.out Testing (+1.00000000000000002E-003,+0.0000000000000000) 1/tan= (+999.99966666664454,+0.0000000000000000) time: +0.243003011 cotan= (+999.99966666664443,-0.0000000000000000) time: +0.535230994 cos/sin= (+999.99966666664443,-0.0000000000000000) time: +0.530906022 tan(pi/2-a)= (+1000.0433799820958,+0.0000000000000000) time: +0.263039947 Testing (+1.0000000000000000,+0.0000000000000000) 1/tan= (+0.64209261593433076,+0.0000000000000000) time: +0.264049053 cotan= (+0.64209261593433076,-0.0000000000000000) time: +0.602337122 cos/sin= (+0.64209261593433076,-0.0000000000000000) time: +0.600989103 tan(pi/2-a)= (+0.64209267766718214,+0.0000000000000000) time: +0.247081995 Testing (+1.5700000000000001,+0.0000000000000000) 1/tan= (+7.96326963223192592E-004,+0.0000000000000000) time: +0.258723021 cotan= (+7.96326963223192592E-004,-0.0000000000000000) time: +0.588716030 cos/sin= (+7.96326963223192592E-004,-0.0000000000000000) time: +0.590244293 tan(pi/2-a)= (+7.96370674640914985E-004,+0.0000000000000000) time: +0.230789661 Testing (+3.0000000000000000,+0.0000000000000000) 1/tan= (-7.0152525514345339,-0.0000000000000000) time: +0.321346283 cotan= (-7.0152525514345339,-0.0000000000000000) time: +0.733272552 cos/sin= (-7.0152525514345339,-0.0000000000000000) time: +0.718843937 tan(pi/2-a)= (-7.0152503565215953,+0.0000000000000000) time: +0.266773224 t