search for: 1e100

...ould be' but it is not, and we should not be making statements like that. Those in the email quoted in pgamma.c exhibit hubris. There are also at least two examples of sloppy coding that lead to numeric overflow and complete loss of accuracy. Consider > pgamma(seq(0.75, 1.25, by=0.05)*1e100, shape = 1e100, log=T) [1] -3.768207e+98 NaN NaN NaN NaN [6] -6.931472e-01 NaN NaN NaN NaN [11] 0.000000e+00 Warning message: NaNs produced in: pgamma(q, shape, scale, lower.tail, log.p) > pgamma(seq(0.75, 1.25,...

Hello everyone I am finding complex numbers to be printed oddly in R-alpha. > 1e100+1e100i [1] 1e+100+1e+100i > 1e100 + 1e44i [1] 1e+100+100000000000000008821361405306422640701865984i > dput(1e100+1e44i) 1e+100+100000000000000008821361405306422640701865984i > [same thing at home on a linux box] The first line is what I would expect, the second and third lines aren...

Hello everybody, I search a function which returns the contour of map with levels like contourLines, but I would like this function return the border of the map too, because the function contourLines cannot consider the corner of the map and it is not adapted to fill polygon after that. Thanks in advance Pierre Bruyer

...by Bill that both as.Date() and format.Date() should probably be tweaked such as to signal warnings (e.g. on integer overflow for too large numbers). ## ------------------------------------------------------------------------- xDates <- lapply(c(-Inf, Inf, NA, NaN, 1e9, 4e9, 1e100, .Machine$double.xmax), as.Date, origin = "2000-01-01") str(xDates) # --> first 4 *all* show as NA sapply(xDates, is.na) # the two +-Inf are not NA (f.D <- sapply(xDates, format))# 1..4: NA, then "negative" but all the same (?!) stopifnot(is.na(f.D)[1:4])...

Full_Name: Morten Welinder Version: 2.1.0 OS: src only Submission from: (NULL) (216.223.241.212) Now that pcauchy has been fixed, it is becoming clear that qcauchy suffers from the same problems. qcauchy(pcauchy(1e100,0,1,FALSE,TRUE),0,1,FALSE,TRUE) should yield 1e100 back, but I get 1.633178e+16. The code below does much better. Notes: 1. p need not be finite. -Inf is ok in the log_p case and R_Q_P01_check already checks things. 2. No need to disallow scale=0 and infinite location. 3. The code below u...

Here's a contribution from Ian Smith that got bounced from the list. -------- Original Message -------- Subject: Re: [Rd] pt inaccurate when x is close to 0 (PR#9945) Date: Thu, 11 Oct 2007 06:02:43 -0400 From: iandjmsmith at aol.com To: murdoch at stats.uwo.ca Duncan, I tried sending the rest of this to R-devel but it was rejected as spam, hence the personal e-mail. R calculates the pt

Hi, > 1e10+5i [1] 1e+10+0e+00i > Im(1e10+5i) [1] 5 maybe little better... --- R-3.5.1.orig/src/main/complex.c 2018-03-26 07:02:25.000000000 +0900 +++ R-3.5.1/src/main/complex.c 2018-07-10 12:50:42.523874767 +0900 @@ -381,6 +381,7 @@ r->i = fround(pow10 * x->i, digits)/pow10; } else { digits = (double)(dig); + if(digits < 1) digits=1; /* a little better */

Hi there, pmvnorm(lo=c(-Inf,-Inf), up=c(Inf,Inf), mean=c(0,0) ) should give me "1", right? But it doens't - it giver me "0". Would someone help me, please? [[alternative HTML version deleted]]

Hello all how can sample from f(x)~x^(a-1)*ind(0,min(b,-log(u)) in R? where a and b is positive constand and 0<u<1 thanks khazaei

Dear All, I am trying to solve the following maximization problem with R: find x(t) (continuous) that maximizes the integral of x(t) with t from 0 to 1, subject to the constraints dx/dt = u, |u| <= 1, x(0) = x(1) = 0. The analytical solution can be obtained easily, but I am trying to understand whether R is able to solve numerically problems like this one. I have tried to find an

...DVu0hVln2UQvNhFkZy 3DvfA1z1HOMmUlVpyYuAVgJa1FUd7VB34jmKmyrB4+0PS39Cbo9+rVZvyKNxHj/z8Kh+ l/ivk5e35L2eq1cUiomqhcus8m7vZ3edA0wEQbU4Z5p3oGBQcGb/5pe3NczjRjJ0oF1q 5OEcIXAxiREMU3DuiYkJCQE28hPMFf10/1/NtGLwlFGdLTseq8LwMM9hY9h8aKmkFAyB Xrvg==^M X-Google-DKIM-Signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=1e100.net; s=20161025; h=x-gm-message-state:from:to:subject:date:message-id:mime-version :thread-index:content-language; bh=u9iSrDj6VaHxbCoUFYZ7Zp0PHd/EPA66NRtYVeLn0+Q=; b=VvYSZ0dCKbTlD3JWsIM+c/IvKn4ODGwMBBDIlHSpKHVjLOX3bAQk5NLCdJgsVrHTH3 eX+fCHmwyIgoaa2bMZK45aURu0H/kELwUc4k/ToQ+MLikNLe7J1GEcDeG2R0Mzk6oO...

Hi Gabriel, The point is that it *visually* displays as NA, but is.na() still responds as FALSE. When I (and I am sure many people) see an NA, we then use is.na(). If we see Inf displayed, we then use is.infinite(). With as.Date() this breaks down. I'm not arguing that as.Date(Inf) should be coerced to NA. I'm arguing that as.Date(Inf) should be *visually* displayed as Inf (i.e. the

...substitute function was handling NULL in its 'rho' argument incorrectly. o The code for pgamma() introduced in 2.1.0 failed for large values of 'shape' where the previous code was perfectly acceptable, despite the claim to be uniformly better. For example, pgamma(0.9e100, 1e100) was NaN. (PR#8528) o There was no command 'ls' in browser() nor 'next' in debug(), despite the documentation (which has been corrected). Command 'where' in the browser() no longer changes to step-though mode. o factor.scope() could report incorrectly that inte...

...substitute function was handling NULL in its 'rho' argument incorrectly. o The code for pgamma() introduced in 2.1.0 failed for large values of 'shape' where the previous code was perfectly acceptable, despite the claim to be uniformly better. For example, pgamma(0.9e100, 1e100) was NaN. (PR#8528) o There was no command 'ls' in browser() nor 'next' in debug(), despite the documentation (which has been corrected). Command 'where' in the browser() no longer changes to step-though mode. o factor.scope() could report incorrectly that inte...

An evolution of: https://www.redhat.com/archives/libguestfs/2018-August/msg00155.html

If I do: script/generate controller Photo I get a controller named photo_controller. OTOH, if I do: script/generate scaffold Photo I get a controller named photos_controller. (Note that the controller name is pluralized). I realize that I''m specifying the controller name explicitly in the first case, but it seems odd to me that the scaffold command generates a pluralized