 |
 |
W. J. Metzger <[log in to unmask]> wrote: ... >How about sumx = sum( real( x, kind(0d0) ) ) Or (possibly better): sumx = sum( real(x, kind(sumx) ) ) This way you can change to double, quad, or whatever by simply redeclaring SUMX. No need to rewrite code fragments. Of course, this assumes that kind(sumx) is of greater precision than kind(x). -- J. Giles %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%