Your guess is correct, you are getting this error because there is not enough local memory available. Some workarounds would be:
- Increase the number of workers to use more physical machines -- this is not going to work if you use a "local" parpool
- Split your b in slices and solve the slices each at a time. In that case, I strongly recommend to use the decomposition object and do something like:
decA = decomposition(Ad);
for i = 1:nSlices
% Get bdSlicei, for instance load it from a file
xdSlicei = decA \ bdSlicei;
% Use xdSlicei, for instance save it to a file
end
The question "how much RAM is sufficient?" is hard to answer a priori. For sparse matrices "\" is basically decomposing your matrix A into a product L*U where L and U are triangular but in general these matrices are denser than A. That's why you may very well be able to store A but not L and U. The amount of sparsity lost depends on the pattern of A and it is discovered during the decomposition. In practice, it is very problem dependent unfortunately...
