https://gcc.gnu.org/bugzilla/show_bug.cgi?id=100561
Bug ID: 100561
Summary: Error while computing square of a real number (8 and
16 bytes)
Product: gcc
Version: og10 (devel/omp/gcc-10)
Status: UNCONFIRMED
Severity: normal
Priority: P3
Component: fortran
Assignee: unassigned at gcc dot gnu.org
Reporter: h0bb_88 at protonmail dot com
Target Milestone: ---
This simple example didn't give the same result than wolfram or any calculator
:
- result from program : 9.869604950...
- result from other calculators : 9.869604401...
The rounding is not the cause, because of error from 10° decimal, in 16-byte
floating point. This error exist too in 8-byte floating-points number.
Simple example :
program main
real(16), parameter :: pi =
3.1415926535897932384626433832795028841971693993751058209749445923
write(*,*) pi**2, pi*pi
end program
