On Jun 23, 10:06pm, Robin Vowels wrote:
> Subject: Re: Extensions of DPROD
...
> > If one wants dble(x*y), it is easy enough to just say z=x*y,
> > where z is declared as double precision.
>
> No, that's not it at all. DBLE returns a DP result, not a SP
> result masquerading as a double precision result.
Sorry, but z=x*y is IDENTICAL to z=DBLE(x*y), assuming REAL x,y;
DOUBLE PRECISION z. Or else your Fortran compiler is broken.
-P.
--
*********** How can we have ethnic music without ethnic hatred? ***********
* Peter S. Shenkin; Chemistry, Columbia U.; 3000 Broadway, Mail Code 3153 *
** NY, NY 10027; [log in to unmask]; (212)854-5143; FAX: 678-9039 ***
*MacroModel WWW page: http://www.columbia.edu/cu/chemistry/mmod/mmod.html *
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|