https://gcc.gnu.org/bugzilla/show_bug.cgi?id=70235
--- Comment #12 from Jerry DeLisle <jvdelisle at gcc dot gnu.org> --- (In reply to Dominique d'Humieres from comment #11) > > Created attachment 38075 [details] > > A patch for testing > > With the patch and using the test attached to comment 5 with y = 1.0 and > d=8, I get the following output > > -8pf18.8 y= 0.0000000. > -7pf18.8 y= 0.000000.0 ... snip ... > which is wrong for P<=-2. OK great, this is related to the conditions for the shift, thanks. I will keep at it.