[R] arrays when calling C++ from R

[Ramon Diaz]

Dear All,

I am writing some code in C++ that I call from R (using ".C"). I would like to 
have support for 2D and 3D arrays, including reading, assignment and 
operations over whole columns, rows, and slices, etc (with linear algebra a 
plus).  In the past I have used Blitz++ (http://www.oonumerics.org/blitz/).  
However, I would eventually want to make that code a contributed package, and 
thus Blitz++ does not seem a good option (and, for the same reasons, probably 
neither is MTL). I am considering instead TNT (its "ancestor", Lapack++, is 
used in the Matrix package), or GSL. Does anybody have any suggestions as to 
the best way to go?

Thanks,


-- 
Ramón Díaz-Uriarte
Bioinformatics Unit
Centro Nacional de Investigaciones Oncológicas (CNIO)
(Spanish National Cancer Center)
Melchor Fernández Almagro, 3
28029 Madrid (Spain)
http://bioinfo.cnio.es/~rdiaz