search for: 100i

Prof Brian Ripley wrote: > This is intentional: it aligns the numbers. E.g. > > >options(width=12) > >print(c(1+1i, 1-10i, 1+100i)) > [1] 1+ 1i > [2] 1- 10i > [3] 1+100i > > Neat, eh? > > What made you think this was a bug? Ah ok, I've misunderstood this feature probably perhaps because, at first sight, I found the display looks "strange" when some of the numbers have large imaginary pa...

...different parts of the function's domain. I am compiling a test suite for a package, and would like to verify that the different methods return approximately identical results. Toy example follows: R> f1 <- function(x){ (x-1)*(x+1)} R> f2 <- function(x){x^2-1} R> x <- pi+100i R> abs(f1(x) - f2(x)) [1] 9.298027e-12 R> stopifnot(abs(f1(x)-f2(x)) < 1e-11) Observe that f1() should be identically equal to f2(); any differences are due to rounding errors (needless to say, the real examples are more complex, with larger errors). My question is, I am unhappy about...

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't:...

Well, this is what i got... > -0.0841219200008394^(1/3) [1] -0.438163696867656 > (-0.0841219200008394)^(1/3) [1] NaN and i don't have a clue of why this happens or how to avoid it, any suggestions? thank you, Juan

...; preds = %entry %tmp4.i.i.i.i.i = tail call i64 @ContextFrame_getInt64Value(%struct.ContextFrame* %context, i32 0, i64 726) nounwind readonly ; <i64> [#uses=1] %tmp23.i.i.i.i.i = sitofp i64 %tmp4.i.i.i.i.i to double ; <double> [#uses=1] ... (gdb) x/100i 0xfffffd7ff9302528 0xfffffd7ff9302528: push %r14 0xfffffd7ff930252a: push %rbx 0xfffffd7ff930252b: sub $0x18,%rsp 0xfffffd7ff930252f: test %rdx,%rdx 0xfffffd7ff9302532: jne 0xfffffd7ff9302747 0xfffffd7ff9302538: mov %rdi,%rbx 0xfffffd7ff930253b: xor %r1...

...%entry >  %tmp4.i.i.i.i.i = tail call i64 > @ContextFrame_getInt64Value(%struct.ContextFrame* %context, i32 0, i64 > 726) nounwind readonly ; <i64> [#uses=1] >  %tmp23.i.i.i.i.i = sitofp i64 %tmp4.i.i.i.i.i to double ; <double> > [#uses=1] > ... > > > (gdb) x/100i 0xfffffd7ff9302528 > 0xfffffd7ff9302528:     push   %r14 > 0xfffffd7ff930252a:     push   %rbx > 0xfffffd7ff930252b:     sub    $0x18,%rsp > 0xfffffd7ff930252f:     test   %rdx,%rdx > 0xfffffd7ff9302532:     jne    0xfffffd7ff9302747 > 0xfffffd7ff9302538:     mov    %rdi,%rbx >...

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 */