3d matrices in C++ using gsl_vector

Question

I have written a lot of my code using the gsl_vector type from the GSL library. I have created a 3D matrix by simply using vector< vector <gsl_vector * > >. In fact, theoretically one can use any vector class in the inner template: vector< vector < VEC > >, one doesn't necessarily need a pointer, but I would like to avoid making this change for the present.

However, I would like to avoid using the std::vector class in order to avoid the additional overhead it causes, see this answer. Can anyone suggest a convenient way to do this?

NB: I need gsl_vector also to perform BLAS operations.

SUGGESTION: Use vector<gsl_matrix*> and vector views.

Is there a specific reason you chose not to take advantage of the GLM library? — Syntactic Fructose, Jan 17 '14 at 18:10

score 0 · Answer 1 · edited May 23 '17 at 12:11

0

What about allocating a 2D array of pointers to gsl_vectors. See this answer https://stackoverflow.com/a/1768382/3146587 for an example.

edited May 23 '17 at 12:11

Community

1
1

answered Jan 17 '14 at 17:43

user3146587

4,250
1
16
25

3d matrices in C++ using gsl_vector

1 Answers1