but I don't like the storage allocation of 'offset'. Is there a better way to do it? I would like to avoid for-loops.
Is there a way to access the current index of mapslices (ciom) ? And do something like this: mapslices(x - > m[ x[:, ciom ], ciom], i, 1); or do i need to do a for? for ciom=1:c; v[:,ciom] =m [ i[:, ciom ], ciom] ; end other ideas?