 |
 |
> REAL :: x( 100 )
> DOUBLE PRECISION :: sumx
> INTEGER i
>
> ...
>
> sum = 0
> DO i = 1, 100
> sumx = sumx + x( i )
> END DO
>
>
> One is tempted to preplace the DO loop with the F90 statement:
>
> sumx = sum( x )
>
> However, AFAIK, the summation will then be done in single precision.
What about
sumx = sum( dble(x) )
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%