public inbox for gcc-bugs@sourceware.org
help / color / mirror / Atom feed
From: "fxcoudert at gcc dot gnu dot org" <gcc-bugzilla@gcc.gnu.org>
To: gcc-bugs@gcc.gnu.org
Subject: [Bug libfortran/24685] real(16) formatted input is broken for huge values
Date: Wed, 24 May 2006 18:58:00 -0000 [thread overview]
Message-ID: <20060524185836.17305.qmail@sourceware.org> (raw)
In-Reply-To: <bug-24685-8500@http.gcc.gnu.org/bugzilla/>
------- Comment #19 from fxcoudert at gcc dot gnu dot org 2006-05-24 18:58 -------
(In reply to comment #18)
> it is still broken on powerpc{,64}
Hi Jakub, I'm not sure exactly what is still broken. On
powerpc-apple-darwin7.9.0, with mainline gfortran 20060512:
$ cat foo.f90
! { dg-do run }
! { dg-require-effective-target fortran_large_real }
program huge_real10_formatted
! This should be kind=10 on systems that support it
integer, parameter :: k = selected_real_kind (precision (0.0_8) + 1)
real(kind=k) :: a,b(2), c
character(len=180) :: tmp
! Test real(k) scalar and array formatted IO with big value
b(:) = huge (1.0_k)/2
print *, 'real(10) big value: ', b(1)
write (tmp, *) b
read (tmp, *) a, c
print *, 'same value read again: ', a
print *, 'difference: ', a-b(1)
! Test with really big value
b(:) = huge (1.0_k)
print *, 'huge value: ', b(1)
write (tmp, *) b
read (tmp, *) a, c ! This is line 19
print *, "We don't get this far!"
if (a /= b(1)) call abort ()
if (c /= b(2)) call abort ()
end program huge_real10_formatted
$ gfortran foo.f90 && ./a.out
real(10) big value: 4.4942328371557897693232629769725618E+0307
same value read again: 4.4942328371557897693232629769725618E+0307
difference: 0.000000000000000000000000000000000
huge value: 8.9884656743115795386465259539451237E+0307
We don't get this far!
$ gfortran foo.f90 -mlong-double-128 && ./a.out
real(10) big value: 4.4942328371557897693232629769725618E+0307
same value read again: 4.4942328371557897693232629769725618E+0307
difference: 0.000000000000000000000000000000000
huge value: 8.9884656743115795386465259539451237E+0307
We don't get this far!
$ cat real16.f90
character(len=100) :: a
real(kind=16) :: x, y, z
x = huge(x)
write(a,*) x
read(a,*) y
if (x /= y) print *, x, y, x-y
end
$ gfortran real16.f90 && ./a.out
$ gfortran real16.f90 -mlong-double-128 && ./a.out
All this looks like it's working fine...
--
fxcoudert at gcc dot gnu dot org changed:
What |Removed |Added
----------------------------------------------------------------------------
CC| |jakub at gcc dot gnu dot org
AssignedTo|fxcoudert at gcc dot gnu dot|unassigned at gcc dot gnu
|org |dot org
Status|ASSIGNED |NEW
http://gcc.gnu.org/bugzilla/show_bug.cgi?id=24685
next prev parent reply other threads:[~2006-05-24 18:58 UTC|newest]
Thread overview: 40+ messages / expand[flat|nested] mbox.gz Atom feed top
2005-11-05 19:34 [Bug libfortran/24685] New: real(10) " jblomqvi at cc dot hut dot fi
2005-11-05 23:24 ` [Bug libfortran/24685] " fxcoudert at gcc dot gnu dot org
2006-02-05 22:20 ` [Bug libfortran/24685] [4.1 only] " fxcoudert at gcc dot gnu dot org
2006-02-05 23:04 ` fxcoudert at gcc dot gnu dot org
2006-02-14 15:48 ` fxcoudert at gcc dot gnu dot org
2006-02-14 15:50 ` fxcoudert at gcc dot gnu dot org
2006-02-16 10:55 ` ebotcazou at gcc dot gnu dot org
2006-02-16 11:12 ` [Bug libfortran/24685] real(16) " fxcoudert at gcc dot gnu dot org
2006-02-16 17:44 ` pinskia at gcc dot gnu dot org
2006-02-23 17:44 ` sje at cup dot hp dot com
2006-03-15 7:29 ` ebotcazou at gcc dot gnu dot org
2006-03-15 16:22 ` jb at gcc dot gnu dot org
2006-03-15 16:44 ` sje at cup dot hp dot com
2006-03-16 8:30 ` jb at gcc dot gnu dot org
2006-04-01 21:34 ` ebotcazou at gcc dot gnu dot org
2006-04-01 21:35 ` ebotcazou at gcc dot gnu dot org
2006-04-10 12:03 ` jakub at gcc dot gnu dot org
2006-04-21 17:04 ` jakub at gcc dot gnu dot org
2006-04-21 17:22 ` jakub at gcc dot gnu dot org
2006-05-24 18:58 ` fxcoudert at gcc dot gnu dot org [this message]
2006-05-24 19:01 ` Andrew Pinski
2006-05-24 19:01 ` pinskia at physics dot uc dot edu
2006-05-24 20:35 ` fxcoudert at gcc dot gnu dot org
2007-05-25 19:50 ` jvdelisle at gcc dot gnu dot org
2007-05-26 17:51 ` dominiq at lps dot ens dot fr
2007-11-03 16:48 ` jvdelisle at gcc dot gnu dot org
2007-11-12 5:44 ` jvdelisle at gcc dot gnu dot org
2008-01-17 11:24 ` [Bug libfortran/24685] real(16) formatted input is broken for huge values (gfortran.dg/default_format_2.f90) burnus at gcc dot gnu dot org
2008-01-17 15:33 ` ebotcazou at gcc dot gnu dot org
2008-01-24 15:09 ` dominiq at lps dot ens dot fr
2008-01-25 15:47 ` dominiq at lps dot ens dot fr
2008-02-23 20:11 ` fxcoudert at gcc dot gnu dot org
2008-02-23 20:11 ` fxcoudert at gcc dot gnu dot org
2008-02-28 7:10 ` uros at gcc dot gnu dot org
2008-12-12 22:20 ` janis at gcc dot gnu dot org
2008-12-12 22:23 ` janis at gcc dot gnu dot org
2008-12-19 1:39 ` howarth at nitro dot med dot uc dot edu
2008-12-19 17:51 ` janis at gcc dot gnu dot org
2008-12-19 18:15 ` janis at gcc dot gnu dot org
2008-12-19 18:23 ` janis at gcc dot gnu dot org
Reply instructions:
You may reply publicly to this message via plain-text email
using any one of the following methods:
* Save the following mbox file, import it into your mail client,
and reply-to-all from there: mbox
Avoid top-posting and favor interleaved quoting:
https://en.wikipedia.org/wiki/Posting_style#Interleaved_style
* Reply using the --to, --cc, and --in-reply-to
switches of git-send-email(1):
git send-email \
--in-reply-to=20060524185836.17305.qmail@sourceware.org \
--to=gcc-bugzilla@gcc.gnu.org \
--cc=gcc-bugs@gcc.gnu.org \
/path/to/YOUR_REPLY
https://kernel.org/pub/software/scm/git/docs/git-send-email.html
* If your mail client supports setting the In-Reply-To header
via mailto: links, try the mailto: link
Be sure your reply has a Subject: header at the top and a blank line
before the message body.
This is a public inbox, see mirroring instructions
for how to clone and mirror all data and code used for this inbox;
as well as URLs for read-only IMAP folder(s) and NNTP newsgroup(s).