hello  Jörn Ungermann

I use array_adapter (and it's variations) to basically accomplish that:

matrix<double, column_major, array_adapter<double> > matrix(m, n, array_adapter<double>(...));

if you are looking for immediate solution.