[Bug fortran/125914] -1.0**integer should be real(-1**integer)

[tkoenig at gcc dot gnu.org via Gcc-bugs]

https://gcc.gnu.org/bugzilla/show_bug.cgi?id=125914

--- Comment #1 from Thomas Koenig <tkoenig at gcc dot gnu.org> ---
We could just translate this into

 tmp(k) = tmp(k) + (1.0 - 2.0 * iand(j+1,1_real64)) / real (2*j-1)

(same trick that we use for integer).