[Rd] Possibly useful idea

2016-01-11 Thread ProfJCNash
I've not worked changing underlying computational infrastructure, but
developers who do might want to use ideas from FlexiBLAS. Apologies in
advance if this is well-known.

Best JN

> From: Martin Koehler koehl...@mpi-magdeburg.mpg.de
> Date: January 07, 2016
> Subject: FlexiBLAS Version 1.3.0 Release
>
> It is our pleasure to announce the new release of FlexiBLAS. We do not
> provide yet another BLAS implementation in itself, but rather present
> an easy to use method for switching between different BLAS
> implementations while debugging or benchmarking other codes.
>
> Highlights since the release of the initial version (v1.0.0) are:
> - support for BLAS-like extensions known from OpenBLAS and Intel MKL,
> - switch the BLAS backend at program runtime,
> - integration in GNU Octave,
> - and a redesigned command line tool to manage the configuration.
>
> Feel free to check the current version at
> http://www.mpi-magdeburg.mpg.de/projects/flexiblas/
> and drop us a line if you have any comments or requests about it.
>
> Further information on the realization can be found in LAWN 284
> http://www.netlib.org/lapack/lawnspdf/lawn284.pdf
> or updated version of the Manuscript (coming soon) available from the
> project web page.
>
> The library is available under GPLv3. We have tested it in production
> environments on several platforms and hope you will find it useful for
> your research, as well.
>
> Martin Köhler and Jens Saak
>

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel

Re: [Rd] coerce SEXP type to C++ matrix class and back

2016-01-11 Thread Dirk Eddelbuettel

I don't want to sound disrespectful, but why not use RcppArmadillo if
simplicity is you goal?

It is hard to beat the _fully automatic_ conversion:

R> cppFunction("arma::mat doubleUp(const arma::mat & x) { return 2*x; }", 
depends="RcppArmadillo")
R> doubleUp(matrix(1:9,3,3))
 [,1] [,2] [,3]
[1,]28   14
[2,]4   10   16
[3,]6   12   18
R> 

You can of course do all that by hand too, but why do you think both Rcpp and
RcppArmadillo have, respectively, tens of thousands of lines of code?

Dirk

-- 
http://dirk.eddelbuettel.com | @eddelbuettel | e...@debian.org

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] Possibly useful idea

2016-01-11 Thread Dirk Eddelbuettel

John,

On 11 January 2016 at 07:44, ProfJCNash wrote:
| I've not worked changing underlying computational infrastructure, but
| developers who do might want to use ideas from FlexiBLAS. Apologies in
| advance if this is well-known.

You have been able to do that at the OS level, given a suitable OS variant
like Debian, since the 1990s. My timing comparison package/vignette gcbd does
just that (adjusting before calling R).

I see nothing really new here. It's just dispatching. Sure, if your OS allows
it you can roll it up into the OS. But not all OSs will so appeal to R is
limited.

Dirk

-- 
http://dirk.eddelbuettel.com | @eddelbuettel | e...@debian.org

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


Re: [Rd] coerce SEXP type to C++ matrix class and back

2016-01-11 Thread Antonio Perez
Dirk,

I am sorry if it sounded disrespectful. On the contrary, as I said,
your packages Rcpp and RcppArmadillo work very well, and use both of
them in more "complex" applications. However, in this particular case
in which I have just one object of one type (matrix double) with its
elements stored column-major similarly to R I though that could find a
reasonable method to coerce my SEXP matrix into a Matrix double in
armadillo without using two packages for this.

I hope this is understood.



I though that there might beI ptough that might result which I insist
has similar and , a matrix with one typeit seems veryas we say
literally in Spanishother

2016-01-11 13:53 GMT+01:00 Dirk Eddelbuettel :
>
> I don't want to sound disrespectful, but why not use RcppArmadillo if
> simplicity is you goal?
>
> It is hard to beat the _fully automatic_ conversion:
>
> R> cppFunction("arma::mat doubleUp(const arma::mat & x) { return 2*x; }", 
> depends="RcppArmadillo")
> R> doubleUp(matrix(1:9,3,3))
>  [,1] [,2] [,3]
> [1,]28   14
> [2,]4   10   16
> [3,]6   12   18
> R>
>
> You can of course do all that by hand too, but why do you think both Rcpp and
> RcppArmadillo have, respectively, tens of thousands of lines of code?
>
> Dirk
>
> --
> http://dirk.eddelbuettel.com | @eddelbuettel | e...@debian.org

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] On 'R CMD INSTALL' with multiple architectures

2016-01-11 Thread Berwin A Turlach
G'day all,

I guess it is still early enough in the year to wish everybody a happy
and successful new year.

I thought I should report that the installation of the CRAN package
rstan regularly fails on my machine (a 64 bit linux box running Xubuntu
15.10).  The reason being that I have the 32-bit and the 64-bit
architecture of R installed, and my /tmp file is on a partition with
about 1Gb space.  

During the installation of rstan the compilation of the source
directory regularly fills up my /tmp directory such that at some point
the compilation for the 64-bit architecture (my main architecture)
fails as /tmp runs out of space.  At that time, the src-32 subdirectory
has a size of about 460 MB (~ half of the partition!).

The most recent version of rstan I finally managed to install by
issuing the commands:

R CMD INSTALL --no-multiarch rstan_2.9.0.tar.gz
followed by
R --arch=32 CMD INSTALL --libs-only rstan_2.9.0.tar.gz

The last command, to my surprise, actually tried to compile both
libraries (32-bit and 64-bit) again.  So the installation kept failing
until I deleted the src-32 directory while the 64 libraries where build.

By now I realise that the R-admin manual only suggests/documents that
R --arch=name CMD INSTALL --libs-only pkg1 pgk2 
only installs the library for the specified architecture if the package
has "an executable configure script or a src/Makefile file", and is
quite about its behaviour otherwise.  But I wonder whether it would be
reasonable for users to expect that 'R --arch=32 CMD INSTALL
--libs-only' installs only the library for the specified architecture
in all circumstances.

Playing around with a package that also has compiled code, but is
faster to install than rstan, I realise now that my second command
should have been:
  R --arch=32 CMD INSTALL --libs-only --no-multiarch rstan_2.9.0.tar.gz

In summary, I would like to suggest that 'R CMD INSTALL' deletes
architecture specific 'src' subdirectories as soon as they are no
longer needed and/or that 'R --arch=name CMD INSTALL --libs-only'
installs only libraries for the specified architecture (as an unwary
user might expect).

Cheers,

Berwin

== Full address 
A/Prof Berwin A Turlach   Tel.: +61 (8) 6488 3338 (secr)
School of Maths and Stats (M019)+61 (8) 6488 3383 (self)
The University of Western Australia   FAX : +61 (8) 6488 1028
35 Stirling Highway   
Crawley WA 6009 e-mail: berwin.turl...@gmail.com
Australiahttp://www.maths.uwa.edu.au/~berwin
 http://www.researcherid.com/rid/A-4995-2008

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel


[Rd] Small inaccuracy in the Writing R Extensions manual

2016-01-11 Thread Berwin A Turlach
G'day all,

In Chapter 1.4 (Writing package vignettes) the Writing R Extensions
manual states:

By default @code{R CMD build} will run @code{Sweave} on all
Sweave vignette source files in @file{vignettes}.  If
@file{Makefile} is found in the vignette source directory, then
@code{R CMD build} will try to run @command{make} after the
@code{Sweave} runs, otherwise @code{texi2pdf} is run on each
@file{.tex} file produced.

This does not seem to be quite correct as stated.  'R CMD build' seems
to run make only if there was a file in the directory vignettes that
Sweave successfully processed.   If the directory vignettes contains a
Makefile and subdirectories in which the actual vignettes are, 'R CMD
build' does not run make.

Cheers,

Berwin

__
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel