problem with parfor in a subroutine

Dear;

I rather new using parfor computation.

The parfor execute a program in a parallel basis, when this program arrive to a routine:

    function Erz=EFieldCalIdeal(r,z,V,RThimble,RCentral,LCentral)
        for i=1:length(r)
            ...
            if r(i)<=RThimble
28              Erz(i,1)=F/r(i);
            elseif
            ...
            end
        end
    end

Gives the following error:

Error using EFieldCalIdeal (line 28)

Assignment has more non-singleton rhs dimensions than non-singleton subscripts

But in fact, the program has no problem to run using a conventional for-loop.

Any idea?

Thanks a lot,

[Walter Roberson]

Is the parfor in the calling routine, or in this routine?

It looks to me as if F is non-scalar when the problem occurs.

[Luis Isaac]

The parfor call the routine

[Walter Roberson]

Is there any randomness involved inside the parfor ? And is there any solve() involved? How is F being constructed ?