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++ (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) bioinfo.cnio.es/~rdiaz