https://gcc.gnu.org/bugzilla/show_bug.cgi?id=98201
--- Comment #10 from Thomas Koenig <tkoenig at gcc dot gnu.org> --- I don't have a working mingw system myself, but I dusted off my cygwin system for this, using their cross-compiler to mingw. With $ x86_64-w64-mingw32-gfortran.exe -static -static-libgfortran csqrt.f90 this gave Argument SquareRoot SR squared (-.1000E+02 0.1000E+00) (0.1581E-01 0.3162E+01) (-.1000E+02 0.1000E+00) (-.1000E+03 0.1000E+00) (0.5000E-02 0.1000E+02) (-.1000E+03 0.1000E+00) (-.1000E+04 0.1000E+00) (0.1581E-02 0.3162E+02) (-.1000E+04 0.1000E+00) (-.1000E+05 0.1000E+00) (0.5000E-03 0.1000E+03) (-.1000E+05 0.1000E+00) (-.1000E+06 0.1000E+00) (0.1581E-03 0.3162E+03) (-.1000E+06 0.1000E+00) (-.1000E+07 0.1000E+00) (0.5000E-04 0.1000E+04) (-.1000E+07 0.1000E+00) (-.1000E+08 0.1000E+00) (0.1581E-04 0.3162E+04) (-.1000E+08 0.1000E+00) (-.1000E+09 0.1000E+00) (0.5000E-05 0.1000E+05) (-.1000E+09 0.1000E+00) (-.1000E+10 0.1000E+00) (0.1581E-05 0.3162E+05) (-.1000E+10 0.1000E+00) which is OK (I tried the program both from the cygwin command line and the cmd prompt).