I have a 3D array cart = s x 3 x p double, where s and p are variables chosen as integers by an optimiser. I'm converting this array to the array sph = s x 3 x p through cart2sph;
[sph(:,1,:), sph(:,2,:), sph(:,3,:)] = cart2sph(cart(:,1,:), cart(:,2,:), cart(:,3,:));
This works fine, unless the case happens to be the optimiser choses to set s = 1. Instead of outputting sph = 1 x 3 x p, it outputs sph = p x 3, which causes issues with my code calling the function I do this in.
It only happens with s, I've tried limiting p to 1 only it seems to run fine, can anyone help?