From p.dalgaard@biostat.ku.dk  Thu Dec  7 16:33:22 2000
From: p.dalgaard@biostat.ku.dk (Peter Dalgaard BSA)
Date: 07 Dec 2000 17:33:22 +0100
Subject: [Rd] BUG: polyroot() (PR#751)
In-Reply-To: Paul Gilbert's message of "Thu, 07 Dec 2000 10:18:29 -0500"
References: <200011280525.GAA06952@pubhealth.ku.dk> <x2k89dxfto.fsf@blueberry.kubism.ku.dk> <3A2FAA45.6CBA4942@bank-banque-canada.ca>
Message-ID: <x2lmts5gz1.fsf@blueberry.kubism.ku.dk>

Paul Gilbert <pgilbert@bank-banque-canada.ca> writes:

> Peter
> 
> I noticed with the original problem that there was a difference between R 1.1.1
> and the R 1.2  snapshot from a few weeks ago. I would guess you are not running 1.2
> since you get
> 
> >> polyeval(1:5,polyroot(1:5))
> >[1] -2.220446e-16-2.664833e-16i -2.220446e-16-1.554475e-16i
> >[3] -2.220446e-16+1.436839e-16i -2.220446e-16-5.076235e-16i
> 
> and with my snapshot I get
> > polyeval(1:5,polyroot(1:5))
> [1]  3.330669e-16+5.828671e-16i  3.330669e-16+2.220446e-16i
> [3] -2.220446e-16-3.885781e-16i  3.330669e-16-5.828671e-16i
> 
> I expect you probably know more about the change than I do, but though I would warn
> you anyway.

Hm? With the current version, I'm still getting

>  polyeval(1:5,polyroot(1:5))
[1] -2.220446e-16-2.664833e-16i -2.220446e-16-1.554475e-16i
[3] -2.220446e-16+1.436839e-16i -2.220446e-16-5.076235e-16i
>

> version
         _                           
platform i586-pc-linux-gnu           
arch     i586                        
os       linux-gnu                   
system   i586, linux-gnu             
status   Under development (unstable)
major    1                           
minor    2.0                         
year     2000                        
month    12                          
day      07                          
language R                           

However, my 1.1.1 is different:

> polyeval(1:5,polyroot(1:5))
[1] 1.110223e-16-3.669889e-16i 0.000000e+00-6.722053e-16i
[3] 0.000000e+00-5.759824e-16i 6.661338e-16-4.592409e-16i
> polyeval(1:10,polyroot(1:10))
[1] -2.815017e-08-3.472720e-09i -2.815014e-08-3.472730e-09i
[3]  3.468454e+04-3.084396e+04i  5.016033e+04-6.107187e+04i
[5]  2.728165e+04+3.719825e+04i -3.455245e+04-4.233822e+03i
[7] -9.016302e+04+5.754516e+04i -7.508538e+04-7.772126e+04i
[9]  3.560675e+04+6.973712e+04i
> polyeval(1:9,polyroot(1:9))
[1] 1.127987e-13+7.833802e-13i 1.130207e-13+7.832323e-13i
[3] 1.127987e-13+7.838768e-13i 1.117995e-13+7.844319e-13i
[5] 5.270229e-13+8.071365e-13i 2.582379e-13+7.174489e-13i
[7] 1.179057e-13+6.123224e-13i 1.865175e-13+6.731772e-13i


-- 
   O__  ---- Peter Dalgaard             Blegdamsvej 3  
  c/ /'_ --- Dept. of Biostatistics     2200 Cph. N   
 (*) \(*) -- University of Copenhagen   Denmark      Ph: (+45) 35327918
~~~~~~~~~~ - (p.dalgaard@biostat.ku.dk)             FAX: (+45) 35327907
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-devel mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !)  To: r-devel-request@stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._


[Rd] BUG: polyroot() (PR#751)
Peter Dalgaard BSA