samer@4: function A=meandims(A,dims)
samer@4: % meandims - Sum over multiple dimensions
samer@4: %
samer@4: % meandims :: 
samer@4: %    [Size:[[1,E]]->A] ~'E dimensional array',
samer@4: %    Dims:[[K]->1..E]  ~'dims to sum over'
samer@4: % -> [Size1->A]        :- Size1=arrset(Size,Dims,1). 
samer@4: %
samer@4: % Result is same size as input except for summed dimensions
samer@4: % which are collapsed to a size of 1.
samer@4: for i=dims, A=mean(A,i); end