A=ones(64,64,50100). It is huge matrix of double values. Unfortunately operations on such a matrix are very time consuming.
For example performing summation of all values along the third dimension
takes 235 seconds on Ultra Sparc Sun-Fire-440. That is why sometimes it maybe beneficial to create single or unsigned integer 32bit (uint32) matrix instead of doubles e.g:
A=uint32(ones(64,64,50100)). In this case performing sum over third dimension takes 211 seconds. When using single instead of double, the time is 204 sec. The difference is not so big, however when using uint8 the time is only 60s. Unfortunately using uint8 when one wants to sum up matrix with over 50000 numbers has no sense unless it is heavily sparse matrix - uint8 holds numbers not bigger than 256.
In conclusion, I think it is better to avoid such matrixes.