blob: 3ec3920240f9bea23770b220d54770dcee82b905 (
plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
|
program main
use m1, only : pi
use m4, only : halfpi
implicit none
real :: rpi, rhalfpi
rpi = pi() / 2
rhalfpi = halfpi()
print '(a,ES15.8)', 'floating point precision loss: ', rpi - rhalfpi
end program
|