1! Program to test DPROD intrinsic
2program intrinsic_dprod
3 implicit none
4 real r, s, t
5 double precision dp
6
7 ! 6d60 doesn't fit in a 4-byte real
8 r = 2e30
9 s = 4e30
10 dp = dprod (r, s)
11 if ((dp .gt. 8.001d60) .or. (dp .lt. 7.999d60)) call abort
12end program
13
14