Displaying 20 results from an estimated 28 matches for "bessel_y".
Did you mean:
bessel_i
2007 Jun 18
1
two bessel function bugs for nu<0
....c bessel_i_new.c
#57c57
#< bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-x))/M_PI
#---
#> bessel_k(x, -alpha, expo) * ((ize == 1)? 2. : 2.*exp(-2.0*x))/M_PI
#bug 2: besselY() for nu<0
#don't know how to check in R; a few random checks against mathematica 5.2
#fix:
#$ diff bessel_y_old.c bessel_y_new.c
#55c55
#< return(bessel_y(x, -alpha) + bessel_j(x, -alpha) * sin(-M_PI * alpha));
#---
#> return(bessel_y(x, -alpha) * cos(M_PI * alpha) - bessel_j(x,
-alpha) * sin(M_PI * alpha));
h.
--
----------------------------------
Hiroyuki Kawakatsu
Business School
Dublin City...
2010 Feb 10
1
looping problem
Hi R-users,
I have this code here:
library(numDeriv)
fprime <- function(z)
{ alp <- 2.0165;
rho <- 0.868;
# simplified expressions
a <- alp-0.5
c1 <- sqrt(pi)/(gamma(alp)*(1-rho)^alp)
c2 <- sqrt(rho)/(1-rho)
t1 <- exp(-z/(1-rho))
t2 <- (z/(2*c2))^a
bes1 <- besselI(z*c2,a)
t1bes1 <- t1*bes1
c1*t1bes1*t2
}
## Newton
2002 Jan 22
1
compile problem with bessel_i.c on IRIX64 flexor 6.5 10100655 IP35 (uname -a) (PR#1275)
...d0.c
making gamma.d from gamma.c
making lgamma.d from lgamma.c
making gamma_cody.d from gamma_cody.c
making beta.d from beta.c
making lbeta.d from lbeta.c
making polygamma.d from polygamma.c
making bessel_i.d from bessel_i.c
making bessel_j.d from bessel_j.c
making bessel_k.d from bessel_k.c
making bessel_y.d from bessel_y.c
making choose.d from choose.c
making snorm.d from snorm.c
making sexp.d from sexp.c
making dgamma.d from dgamma.c
making pgamma.d from pgamma.c
making qgamma.d from qgamma.c
making rgamma.d from rgamma.c
making dbeta.d from dbeta.c
making pbeta.d from pbeta.c
making qbeta.d from q...
1999 Mar 06
0
bessel_?.c constants
Hi All,
I've been digging around in src/nmath and have discovered that there are
two sets of "machine constants" in bessel_l.c and bessel_j.c which have
different values:
bessel_i.c: static double ensig = 1e16;
bessel_i.c: static double rtnsig = 1e-4;
bessel_i.c: static double enmten = 8.9e-308;
bessel_i.c: static double enten = 1e308;
bessel_j.c: static double
1999 Mar 06
1
bessel_?.c constants (fwd)
This doesn't seem to have made it to the list, so here it is again. Sorry
if this results in duplicates.
---------- Forwarded message ----------
Date: Fri, 5 Mar 1999 21:36:44 -0800 (PST)
From: Gregory R. Warnes <warnes@biostat.washington.edu>
To: R devel <r-devel@stat.math.ethz.ch>
Subject: bessel_?.c constants
Hi All,
I've been digging around in src/nmath and have
1999 Feb 24
2
Compiling R on Sunos
...chines in the department here. I got
to compile on Solaris 2.5.1 using gcc 2.7.2.2. On my Sunos 4.1.3 (and
I've tried it on Sunos 4.1.4 as well) machines, it crashes during the
compile. I'm using gcc 2.8.1 on there. It gets to this point:
gcc -g -O2 -I../include -I../../src/include -c bessel_y.c -o bessel_y.o
gcc -g -O2 -I../include -I../../src/include -c choose.c -o choose.o
gcc -g -O2 -I../include -I../../src/include -c snorm.c -o snorm.o
gcc -g -O2 -I../include -I../../src/include -c sunif.c -o sunif.o
sunif.c: In function `unif_rand':
sunif.c:131: `RAND_MAX' undeclared (fi...
2009 Mar 17
3
R does not compile any more on FreeBSD 8.0-CURRENT
...o
rm -rf libnmath.a
ar cr libnmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fprec.o
fround.o ftrunc.o sign.o fsign.o imax2.o imin2.o chebyshev.o log1p.o
expm1.o lgammacor.o gammalims.o stirlerr.o bd0.o gamma.o lgamma.o
gamma_cody.o beta.o lbeta.o polygamma.o bessel_i.o bessel_j.o bessel_k.o
bessel_y.o choose.o snorm.o sexp.o dgamma.o pgamma.o qgamma.o rgamma.o
dbeta.o pbeta.o qbeta.o rbeta.o dunif.o punif.o qunif.o runif.o dnorm.o
pnorm.o qnorm.o rnorm.o dlnorm.o plnorm.o qlnorm.o rlnorm.o df.o pf.o
qf.o rf.o dnf.o dt.o pt.o qt.o rt.o dnt.o dchisq.o pchisq.o qchisq.o
rchisq.o rnchisq.o dbinom...
2010 Jun 15
1
Error in nlm : non-finite value supplied by 'nlm'
Hello,
I am trying to compute MLE for non-Gaussian AR(1). The error term follows a difference poisson distribution. This distribution has one parameter (vector[2]).
So in total I want to estimate two parameters: the AR(1) paramter (vector[1]) and the distribution parameter.
My function is the negative loglikelihood derived from a mixing operator.
f=function(vector)
2002 Oct 17
1
underflow handling in besselK (PR#2179)
The besselK() function knows about overflows/underflows internally;
there is a constant xmax_BESS_K in src/nmath/bessel.h (and referred to
only in bessel_k.c), equal to 705.342, which is checked if expon.scaled is
FALSE. (The equivalent number for bessel_i.c is 709, defined as
exparg_BESS in bessel.h.) However, besselK(x) silently returns +Inf if
x>705.342. This behavior is reasonable for
1999 May 03
1
problems compiling R-0.63.3 on alpha
.../include -I../../src/include -c bessel_i.c -o
bessel_i.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c bessel_j.c -o
bessel_j.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c bessel_k.c -o
bessel_k.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c bessel_y.c -o
bessel_y.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c choose.c -o
choose.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c snorm.c -o
snorm.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c sunif.c -o
sunif.o
gcc -ieee_with_inexact -g -I../...
1999 May 03
0
compilation of R-0.63.3 on alpha (PR#183)
..../../src/include -c bessel_i.c=
-o =
bessel_i.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c bessel_j.c=
-o =
bessel_j.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c bessel_k.c=
-o =
bessel_k.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c bessel_y.c=
-o =
bessel_y.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c choose.c -=
o =
choose.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c snorm.c -o=
=
snorm.o
gcc -ieee_with_inexact -g -I../include -I../../src/include -c sunif.c -o=
=
sunif.o
gcc -ieee_with_...
1998 Dec 08
0
R compile problems with AIX
...e -c polygamma.c -o polygamma.o
gcc -g -O2 -I../include -I../../src/include -c bessel_i.c -o bessel_i.o
gcc -g -O2 -I../include -I../../src/include -c bessel_j.c -o bessel_j.o
gcc -g -O2 -I../include -I../../src/include -c bessel_k.c -o bessel_k.o
gcc -g -O2 -I../include -I../../src/include -c bessel_y.c -o bessel_y.o
gcc -g -O2 -I../include -I../../src/include -c choose.c -o choose.o
gcc -g -O2 -I../include -I../../src/include -c snorm.c -o snorm.o
gcc -g -O2 -I../include -I../../src/include -c sunif.c -o sunif.o
gcc -g -O2 -I../include -I../../src/include -c sexp.c -o sexp.o
gcc -g -O2 -I.....
1998 Dec 08
0
R compile problems with AIX
...e -c polygamma.c -o polygamma.o
gcc -g -O2 -I../include -I../../src/include -c bessel_i.c -o bessel_i.o
gcc -g -O2 -I../include -I../../src/include -c bessel_j.c -o bessel_j.o
gcc -g -O2 -I../include -I../../src/include -c bessel_k.c -o bessel_k.o
gcc -g -O2 -I../include -I../../src/include -c bessel_y.c -o bessel_y.o
gcc -g -O2 -I../include -I../../src/include -c choose.c -o choose.o
gcc -g -O2 -I../include -I../../src/include -c snorm.c -o snorm.o
gcc -g -O2 -I../include -I../../src/include -c sunif.c -o sunif.o
gcc -g -O2 -I../include -I../../src/include -c sexp.c -o sexp.o
gcc -g -O2 -I.....
2010 Jun 15
0
nlm is
Hello,
I am trying to compute MLE for non-Gaussian AR(1). The error term follows a difference poisson distribution. This distribution has one parameter (vector[2]).
So in total I want to estimate two parameters: the AR(1) paramter (vector[1]) and the distribution parameter.
My function is the negative loglikelihood derived from a mixing operator.
f=function(vector)
2017 Aug 18
4
No rule to make target all.R, needed by compiler.rdb building R 3.4.1 from source, Scientific Linux release 6.9 (Carbon)
...libnmath.a*
*ar -cr libnmath.a mlutils.o d1mach.o i1mach.o fmax2.o fmin2.o fprec.o
fround.o ftrunc.o sign.o fsign.o imax2.o imin2.o chebyshev.o log1p.o
expm1.o lgammacor.o gammalims.o stirlerr.o bd0.o gamma.o lgamma.o
gamma_cody.o beta.o lbeta.o polygamma.o cospi.o bessel_i.o bessel_j.o
bessel_k.o bessel_y.o choose.o snorm.o sexp.o dgamma.o pgamma.o qgamma.o
rgamma.o dbeta.o pbeta.o qbeta.o rbeta.o dunif.o punif.o qunif.o runif.o
dnorm.o pnorm.o qnorm.o rnorm.o dlnorm.o plnorm.o qlnorm.o rlnorm.o df.o
pf.o qf.o rf.o dnf.o dt.o pt.o qt.o rt.o dnt.o dchisq.o pchisq.o qchisq.o
rchisq.o rnchisq.o dbinom....
2001 Jul 25
2
Installation help
...d0.c
making gamma.d from gamma.c
making lgamma.d from lgamma.c
making gamma_cody.d from gamma_cody.c
making beta.d from beta.c
making lbeta.d from lbeta.c
making polygamma.d from polygamma.c
making bessel_i.d from bessel_i.c
making bessel_j.d from bessel_j.c
making bessel_k.d from bessel_k.c
making bessel_y.d from bessel_y.c
making choose.d from choose.c
making snorm.d from snorm.c
making sexp.d from sexp.c
making dgamma.d from dgamma.c
making pgamma.d from pgamma.c
making qgamma.d from qgamma.c
making rgamma.d from rgamma.c
making dbeta.d from dbeta.c
making pbeta.d from pbeta.c
making qbeta.d from q...
2018 Feb 09
1
R Compilation gets stuck on Windows 64
Please note that building R on Windows is documented in "R Installation
and Administration", including links to external software. Particularly
there is a link to texinfo which is part of Rtools. The documentation is
maintained and it is a sufficient source of information for building R
on Windows.
https://cran.r-project.org/doc/manuals/r-release/R-admin.html
2000 Oct 31
1
conflicting types for `Rf_stirlerr
...rr.c
making gamma.d from gamma.c
making lgamma.d from lgamma.c
making gamma_cody.d from gamma_cody.c
making beta.d from beta.c
making lbeta.d from lbeta.c
making polygamma.d from polygamma.c
making bessel_i.d from bessel_i.c
making bessel_j.d from bessel_j.c
making bessel_k.d from bessel_k.c
making bessel_y.d from bessel_y.c
making choose.d from choose.c
making snorm.d from snorm.c
making sexp.d from sexp.c
making dgamma.d from dgamma.c
making pgamma.d from pgamma.c
making qgamma.d from qgamma.c
making rgamma.d from rgamma.c
making dbeta.d from dbeta.c
making pbeta.d from pbeta.c
making qbeta.d from q...
2012 Aug 11
2
compile fails with x86_64-alpine-linux-uclibc-gcc
...d0.c
making gamma.d from gamma.c
making lgamma.d from lgamma.c
making gamma_cody.d from gamma_cody.c
making beta.d from beta.c
making lbeta.d from lbeta.c
making polygamma.d from polygamma.c
making bessel_i.d from bessel_i.c
making bessel_j.d from bessel_j.c
making bessel_k.d from bessel_k.c
making bessel_y.d from bessel_y.c
making choose.d from choose.c
making snorm.d from snorm.c
making dgamma.d from dgamma.c
making sexp.d from sexp.c
making pgamma.d from pgamma.c
making qgamma.d from qgamma.c
making rgamma.d from rgamma.c
making dbeta.d from dbeta.c
making pbeta.d from pbeta.c
making qbeta.d from q...
2008 Dec 22
1
Error compiling R.2.8.1 with gcc 4.4 on Mac OS 10.5.6
...d0.c
making gamma.d from gamma.c
making lgamma.d from lgamma.c
making gamma_cody.d from gamma_cody.c
making beta.d from beta.c
making lbeta.d from lbeta.c
making polygamma.d from polygamma.c
making bessel_i.d from bessel_i.c
making bessel_j.d from bessel_j.c
making bessel_k.d from bessel_k.c
making bessel_y.d from bessel_y.c
making choose.d from choose.c
making snorm.d from snorm.c
making sexp.d from sexp.c
making dgamma.d from dgamma.c
making qgamma.d from qgamma.c
making pgamma.d from pgamma.c
making rgamma.d from rgamma.c
making dbeta.d from dbeta.c
making pbeta.d from pbeta.c
making qbeta.d from q...