{

const int ndims=3;

typedef boost::multi_array<int,ndims> array;

boost::array<array::size_type,ndims> sma_dims = {{2,3,4}};

int data[] = {0,1,2,3,4,5,6,7,8,9,10,11,12,13,

14,15,16,17,18,19,20,21,22,23};

const int data_size = 24;

array sma(sma_dims);

sma.assign(data,data+data_size);

//

// subarray dims:

// [base,stride,bound)

// [0,1,2), [1,1,3), [0,2,4)

//

const array& csma = sma;

typedef array::index_range range;

array::index_gen indices;

array::const_array_view<ndims>::type csma2 =

csma[indices[range(0,2)][range(1,3)][range(0,4,2)]];

for (array::index i = 0; i != 2; ++i)

for (array::index j = 0; j != 2; ++j)

for (array::index k = 0; k != 2; ++k)

// FAIL! csma2 cannot be assigned to.

csma2[i][j][k] = 0;

return boost::report_errors();