search for: arithm

On 2 March 2015 at 09:09, Duncan Murdoch wrote: | I generally recommend that people use Rcpp, which hides a lot of the | details. It will generate your .Call calls for you, and generate the | C++ code that receives them; you just need to think about the real | problem, not the interface. It has its own learning curve, but I think | it is easier than using the low-level code that you need to

...96 a > ## i * 2L 143 164 250.6057 181 234 16863 4096 a > ## 2L * i 144 177 269.5015 193 237 16119 4096 a > ## i + i 152 183 272.6179 199 243 10434 4096 a > > plot(mb, log="y", notch=TRUE) > ## hmm, looks like even the simple arithm. differ slightly ... > ## > ## ==> zoom in: > plot(mb, log="y", notch=TRUE, ylim = c(150,300)) > > dev.copy(png, file="mbenchm-doubling.png") > dev.off() # [ <- why do I need this here for png ??? ] > ##--> see the appended *png graphic > > Th...

...216 20213 4096 a | ## i * 2L 143 164 250.6057 181 234 16863 4096 a | ## 2L * i 144 177 269.5015 193 237 16119 4096 a | ## i + i 152 183 272.6179 199 243 10434 4096 a | | plot(mb, log="y", notch=TRUE) | ## hmm, looks like even the simple arithm. differ slightly ... | ## | ## ==> zoom in: | plot(mb, log="y", notch=TRUE, ylim = c(150,300)) | | dev.copy(png, file="mbenchm-doubling.png") | dev.off() # [ <- why do I need this here for png ??? ] | ##--> see the appended *png graphic | | Those who've learnt EDA...

...707 169 216 20213 4096 a ## i * 2L 143 164 250.6057 181 234 16863 4096 a ## 2L * i 144 177 269.5015 193 237 16119 4096 a ## i + i 152 183 272.6179 199 243 10434 4096 a plot(mb, log="y", notch=TRUE) ## hmm, looks like even the simple arithm. differ slightly ... ## ## ==> zoom in: plot(mb, log="y", notch=TRUE, ylim = c(150,300)) dev.copy(png, file="mbenchm-doubling.png") dev.off() # [ <- why do I need this here for png ??? ] ##--> see the appended *png graphic Those who've learnt EDA or otherwise abou...

...64 bit int types. Problem report from Mats Fredholm <matsf at init.se> - OpenBSD CVS updates: - [atomicio.c auth-krb4.c bufaux.c channels.c compress.c fingerprint.c] [packet.h radix.c rsa.c scp.c ssh-agent.c ssh-keygen.c sshconnect.c] [sshd.c] pedantic: signed vs. unsigned, void*-arithm, etc - [ssh.1 sshd.8] Various cleanups and standardizations. - Runtime error fix for HPUX from Otmar Stahl <O.Stahl at lsw.uni-heidelberg.de> -- | "Bombay is 250ms from New York in the new world order" - Alan Cox | Damien Miller - http://www.mindrot.org/ | Email: djm at mind...

...e that the relative speed of the algorithms also depends on the size > of the input vector. i + i becomes the winner for longer vectors (e.g. i > <- 1:1e6), but a proper Rcpp version is still approximately twice as fast. The difference in speed is probably due to the fact that R does safe arithmetic. C or C++ do not: > doubleThisInt(i) [1] 2147483642 2147483644 2147483646 NA -2147483646 -2147483644 > 2L * i [1] 2147483642 2147483644 2147483646 NA NA NA Warning message: In 2L * i : NAs produced by integer overflow H. > >...

...types. Problem report from Mats Fredholm <matsf at init.se> - OpenBSD CVS updates: - [atomicio.c auth-krb4.c bufaux.c channels.c compress.c fingerprint.c] [packet.h radix.c rsa.c scp.c ssh-agent.c ssh-keygen.c sshconnect.c] [sshd.c] pedantic: signed vs. unsigned, void*-arithm, etc - [ssh.1 sshd.8] Various cleanups and standardizations. - Runtime error fix for HPUX from Otmar Stahl <O.Stahl at lsw.uni-heidelberg.de> 20000316 - Fixed configure not passing LDFLAGS to Solaris. Report from David G. Hesprich <dghespri at sprintparanet.com> - Pr...

...types. Problem report from Mats Fredholm <matsf at init.se> - OpenBSD CVS updates: - [atomicio.c auth-krb4.c bufaux.c channels.c compress.c fingerprint.c] [packet.h radix.c rsa.c scp.c ssh-agent.c ssh-keygen.c sshconnect.c] [sshd.c] pedantic: signed vs. unsigned, void*-arithm, etc - [ssh.1 sshd.8] Various cleanups and standardizations. - Runtime error fix for HPUX from Otmar Stahl <O.Stahl at lsw.uni-heidelberg.de> 20000316 - Fixed configure not passing LDFLAGS to Solaris. Report from David G. Hesprich <dghespri at sprintparanet.com> - Pr...