Re: [Rd] print(big+small*1i) -> big + 0i

From: Prof Brian Ripley <ripley_at_stats.ox.ac.uk>
Date: Thu, 25 Mar 2010 05:58:21 +0000 (GMT)

On Wed, 24 Mar 2010, William Dunlap wrote:

> Should both parts of a complex number be printed
> to the same precision? The imaginary part of 0
> looks a bit odd when log10(real/imag) >=~ getOption("digits"),
> but I'm not sure it is awful. Some people might
> expect the same number of significant digits in the
> two parts.

'significant digits' is documented in ?signif and xref-ed in ?print. So, yes, the are 'printed to the same precision' for the definition given of 'precision'. (One can argue that precision is relative to the modulus, but this definition is a close approximation.)

>From the NEWS for 2.4.0:

     o	Printing of complex numbers could misbehave when one of the
 	parts was large (so scientific notation was used) and the
 	other was so much smaller that it had no significant digits
 	and should have been printed as zero (e.g. 1e80+3e44i).

and 2.2.0

     o	signif() on complex numbers now rounds jointly to give the
 	requested number of digits in the larger component, not
 	independently for each component.


>
>> 1e7+4i
> [1] 10000000+0i
>> 1e7+5i
> [1] 10000000+0i
>> 1e10 + 1000i
> [1] 1e+10+0e+00i
>> getOption("digits")
> [1] 7

...

-- 
Brian D. Ripley,                  ripley_at_stats.ox.ac.uk
Professor of Applied Statistics,  http://www.stats.ox.ac.uk/~ripley/
University of Oxford,             Tel:  +44 1865 272861 (self)
1 South Parks Road,                     +44 1865 272866 (PA)
Oxford OX1 3TG, UK                Fax:  +44 1865 272595

______________________________________________
R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
Received on Thu 25 Mar 2010 - 06:07:11 GMT

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Thu 25 Mar 2010 - 09:11:14 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive