http://gcc.gnu.org/bugzilla/show_bug.cgi?id=57871
kargl at gcc dot gnu.org changed: What |Removed |Added ---------------------------------------------------------------------------- Status|RESOLVED |REOPENED Last reconfirmed| |2013-07-10 CC| |kargl at gcc dot gnu.org Resolution|INVALID |--- Ever confirmed|0 |1 --- Comment #2 from kargl at gcc dot gnu.org --- (In reply to Dominique d'Humieres from comment #1) > > gfortran 4.8.1 -freal-4-real-16 gives real(16) kind with > > selected_real_kind(1) but it ought to give real(8) because that is > > the smallest decimal precision available with that compiler option. > > and your example does exactly that: Closing as INVALID. I think that you may have been a little too hasty in closing this PR. The issue lies in the code for selected_real_kind. implicit none integer,parameter:: p1 = selected_real_kind(1) This initialization expression is requesting the type with smallest precision that exceeds 1. With -freal-4-real-16, that type is double precision (ie kind=8). The problem lies in gfc_simplify_selected_real_kind. It has no knowledge of the -freal-*-real-* options. The loop (lines 5447-5457) for (i = 0; gfc_real_kinds[i].kind != 0; i++) { if (gfc_real_kinds[i].precision >= precision) found_precision = 1; if (gfc_real_kinds[i].range >= range) found_range = 1; if (gfc_real_kinds[i].radix >= radix) found_radix = 1; if (gfc_real_kinds[i].precision >= precision && gfc_real_kinds[i].range >= range && gfc_real_kinds[i].radix >= radix && gfc_real_kinds[i].kind < kind) kind = gfc_real_kinds[i].kind; } searchs the ordered set [24, 53, 53, 113] (on FreeBSD-i386) or [24, 53, 64, 113] (on FreeBSD-amd64). In either case, the initialization expression returns REAL(4). However, with -freal-4-real-16, one might expect the ordered set to be searched is [113, 53, 64, 113] as REAL(4) has been promoted to REAL(16).