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