FW: [R] Problem compiling on HP-UX 10.20
Nuzzo Art-CINT116
A.Nuzzo at motorola.com
Fri Apr 19 20:07:37 CEST 2002
This appears to have run without a problem. Here is part of the results:
> bY <-besselY(2,nu = nu <-seq(0,100,len=501))
> bY
[1] 5.103757e-01 4.203932e-01 3.010522e-01 1.662137e-01
2.736566e-02
[6] -1.070324e-01 -2.316581e-01 -3.440740e-01 -4.442857e-01
-5.342464e-01
[11] -6.174081e-01 -6.983829e-01 -7.827559e-01 -8.770708e-01
-9.890092e-01
[16] -1.127784e+00 -1.304784e+00 -1.534540e+00 -1.836098e+00
-2.234979e+00
[21] -2.765943e+00 -3.476927e+00 -4.434680e+00 -5.732883e+00
-7.503913e+00
[26] -9.935989e+00 -1.329831e+01 -1.797805e+01 -2.453517e+01
-3.378380e+01
[31] -4.691400e+01 -6.567428e+01 -9.264681e+01 -1.316640e+02
-1.884421e+02
...
[486] -3.189664e+149 -7.956669e+149 -1.985626e+150 -4.957272e+150
-1.238132e+151
[491] -3.093642e+151 -7.733056e+151 -1.933793e+152 -4.837785e+152
-1.210766e+153
[496] -3.031450e+153 -7.593065e+153 -1.902654e+154 -4.769560e+154
-1.196113e+155
[501] -3.000826e+155
> summary(bY)
Min. 1st Qu. Median Mean 3rd Qu. Max.
-3.001e+155 -1.067e+107 -1.976e+62 -9.961e+152 -2.059e+23 5.104e-01
>
Art Nuzzo
> -----Original Message-----
> From: Martin Maechler [mailto:maechler at stat.math.ethz.ch]
> Sent: Friday, April 19, 2002 12:30 PM
> To: Nuzzo Art-CINT116
> Cc: R-core at stat.math.ethz.ch
> Subject: Re: FW: [R] Problem compiling on HP-UX 10.20
>
>
> [ NB: Taken off R-help for the moment;
> We should summarize when we found the problem/solution
> ]
>
> >>>>> "Nuzzo" == Nuzzo Art-CINT116 <A.Nuzzo at motorola.com> writes:
>
> Nuzzo> Here is a copy of the last few lines in base-Ex.Rout.fail:
>
> >> x <- seq(3,500);yl <- c(-.3, .2)
> >> plot(x,x, ylim = yl, ylab="",type='n', main = "Bessel Functions
> Nuzzo> Y_nu(x)")
> >> for(nu in nus){xx <- x[x > .6*nu];
> lines(xx,besselY(xx,nu=nu), col =
> Nuzzo> nu+2)}
> >> legend(300,-.08, leg=paste("nu=",nus), col = nus+2, lwd=1)
> >>
> >> x <- seq(10,50000,by=10);yl <- c(-.1, .1)
> >> plot(x,x, ylim = yl, ylab="",type='n', main = "Bessel Functions
> Nuzzo> Y_nu(x)")
> >> for(nu in nus){xx <- x[x > .6*nu];
> lines(xx,besselY(xx,nu=nu), col =
> Nuzzo> nu+2)}
> >> summary(bY <- besselY(2,nu = nu <- seq(0,100,len=501)))
> Nuzzo> Min. 1st Qu. Median Mean 3rd
> Qu. Max.
> Nuzzo> -3.001e+155 -1.067e+107 -1.976e+62 -9.961e+152
> -2.059e+23 5.104e-01
> >> which(bY >= 0)
> Nuzzo> [1] 1 2 3 4 5
> >> summary(bY <- besselY(2,nu = nu <- seq(3,300,len=51)))
> Nuzzo> Error in vector("character", length) : negative
> length vectors are not
> Nuzzo> allowed
> Nuzzo> In addition: There were 23 warnings (use
> warnings() to see them)
> Nuzzo> Execution halted
>
> Thank you,
> now try
>
> bY <- besselY(2,nu = nu <- seq(0,100,len=501))
> bY
> summary(bY)
>
> (I guess the first (compuation) will do fine; and printing will
> be the problem; but who knows ?)
>
> On a recent 386-Linux, I get
>
> > bY <- besselY(2,nu = nu <- seq(0,100,len=501))
> > bY
> [1] 5.103757e-01 4.203932e-01 3.010522e-01
> 1.662137e-01 2.736566e-02
> [6] -1.070324e-01 -2.316581e-01 -3.440740e-01
> -4.442857e-01 -5.342464e-01
> [11] -6.174081e-01 -6.983829e-01 -7.827559e-01
> -8.770708e-01 -9.890092e-01
> [16] -1.127784e+00 -1.304784e+00 -1.534540e+00
> -1.836098e+00 -2.234979e+00
> [21] -2.765943e+00 -3.476927e+00 -4.434680e+00
> -5.732883e+00 -7.503913e+00
> [26] -9.935989e+00 -1.329831e+01 -1.797805e+01
> -2.453517e+01 -3.378380e+01
> [31] -4.691400e+01 -6.567428e+01 -9.264681e+01
> -1.316640e+02 -1.884421e+02
> .................
> .................
> [481] -3.322931e+147 -8.271870e+147 -2.060002e+148
> -5.132308e+148 -1.279200e+149
> [486] -3.189664e+149 -7.956669e+149 -1.985626e+150
> -4.957272e+150 -1.238132e+151
> [491] -3.093642e+151 -7.733056e+151 -1.933793e+152
> -4.837785e+152 -1.210766e+153
> [496] -3.031450e+153 -7.593065e+153 -1.902654e+154
> -4.769560e+154 -1.196113e+155
> [501] -3.000826e+155
> > summary(bY)
> Min. 1st Qu. Median Mean 3rd Qu.
> Max.
> -3.001e+155 -1.067e+107 -1.976e+62 -9.961e+152 -2.059e+23
> 5.104e-01
>
> ---------
>
> Martin Maechler <maechler at stat.math.ethz.ch>
http://stat.ethz.ch/~maechler/
Seminar fuer Statistik, ETH-Zentrum LEO C16 Leonhardstr. 27
ETH (Federal Inst. Technology) 8092 Zurich SWITZERLAND
phone: x-41-1-632-3408 fax: ...-1228 <><
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help 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-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list