In my project I deal with matrices using up almost all the memory available. Logically the matrix is divided into blocks which can be computed independently from one another. Straight forward use of parfor as follows
parfor a = 1:N,
[M1] = func(a,par1, par2, par3);
N(i1(a):i2(a), i3(a):i4(a)) = M1;
does not work (the code is an example to illustrate the idea). Is there a way to establish parallel computing of these blocks?