FW: [R] Problem compiling on HP-UX 10.20

Nuzzo Art-CINT116
2002-04-19T18:07:37Z

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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._