How to fix "Error using == Arrays have incompatible sizes for this operation." in this code?

Hi, I am writing the following code:

for k = 1 : length(theFiles) % files are alreaady defined in original code
    p=1;
    t=1;
    for q=1:43 % there are 43 inputs for arr
c=corr2(B0,mask);
arr(p)=abs(c); 
if(arr(p)<1)
if(arr(p)>0.85)
    cor(t)=arr(p);  
    t=t+1;
end
end
p=p+1;
    end %end of q loop
mx=max(cor);   %conditioned arr 
ixa=find(arr==mx,1) %find the location for which mx matches arr value the first instance

But when run, the result shows following:

Error using  == 
Arrays have incompatible sizes for this operation.

Error in Untitled12 (line 264)
ixa=find(arr==mx,1);

It is noteworthy that, the error shows only for certain image inputs. For other inputs, the code works fine.

How can I fix it?

[kanika bhalla]

@Tawsif Mostafiz Can you add those images on which your code gives error so that we can try and check what is a problem?

Answer 1

[KSSV]

Error is clear, you are trying to equate arr and mx. Check do they have same dimensions? They have different dimensions that's why error.

A = rand(2) ; B = rand(2) ; 
A == B    % no error as the dimensions are same 

ans = 2×2 logical array
   0   0
   0   0

C = rand(3) ; 
A == C  % error as the dimensions are different 

Arrays have incompatible sizes for this operation.

[Tawsif Mostafiz]

I solved the problem, Thanks!

[Shreen El-Sapa]

i have the sameproblem. could u help me?

the program is

A=[  -0.00767
    0.00696
    -0.00031
    0.00005
    -0.00001
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000];
B=[ 0.12000
    0.00209
    -0.00019
    0.00001
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    .00000E+00
    .00000E+00
    .00000E+00
    .00000E+00
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000];
a = 1 ;   %RADIUS
L=.4;
c =-a/L;
b =a/L;
m =a*200; % NUMBER OF INTERVALS
[x,y]=meshgrid((c:(b-c)/m:b),(c:(b-c)/m:b)');
[I J]=find(sqrt(x.^2+y.^2)<(a-.1));
if ~isempty(I)
    x(I,J) = 0;
    y(I,J) = 0;
end
r=sqrt(x.^2+y.^2);
t=atan2(y,x);
warning off
chi=1;alpha=1;ab1=1;
k=sqrt(chi.^2+alpha.^2);   %xi=sqrt(t.^2+k.^2);
%for i=2:5
i=10;
h=@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).* (t1.*(exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))-exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i-1).*t1.^(i-1).*exp(-ab1.*t1)./factorial( i ))+((-1).^(i).*sqrt(pi.*k./2/t1.^2).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2, sqrt(t1.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))));
v=@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).* ((-t1.*exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))+sqrt(t1.^2+k.^2).*exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i).*t1.^(i-1).*exp(-ab1.*t1)./factorial( i ))+((-1).^(i-1).*sqrt(pi.*k./2/(t1.^2+k.^2)).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2, sqrt(t.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))));

h11=integral(h, 0, 10); %inf
h22=integral(v, 0, 10);%inf

Ai=A(i-1);Bi=B(i-1);
psi1=psi1+(Ai.*(r.^(-i+1)+h11)+(r.^(1./2).*besselk(i-1./2,r.*alpha)+h22).*Bi).*gegenbauerC(i,-1./2, cos(t));
%end
[DH1,hh2]=contour(x,y,psi,10,'-k');
%clabel(CH,h1,'FontSize',8);clabel(DH,h2,'FontSize',8)
hold on
m1=100;
r1=ones(1,m1+1)*a;
th=(0:2*pi/m1:2*pi);
set(polar(th,r1,'-k'),'LineWidth',1.1);
title('$\kappa=0$','Interpreter','latex','FontSize',10,'FontName','Times New Roman','FontWeight','Normal')
ylabel({'$\eta=1\quad$'},'Interpreter','latex','FontSize',10,'rot',360,'FontName','Times New Roman','FontWeight','Normal');
%title('Happel$^\prime$s model','Interpreter','latex','FontSize',10,'FontName','Times New Roman','FontWeight','Normal')
%axis square
axis off
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%

outcomes

Arrays have incompatible sizes for this operation.

Error in

test>@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).*(t1.*(exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))-exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i-1).*t1.^(i-1).*exp(-ab1.*t1)./factorial(i))+((-1).^(i).*sqrt(pi.*k./2/t1.^2).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2,sqrt(t1.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))))

(line 77)

h=@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).*

(t1.*(exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))-exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i-1).*t1.^(i-1).*exp(-ab1.*t1)./factorial( i

))+((-1).^(i).*sqrt(pi.*k./2/t1.^2).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2,

sqrt(t1.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))));

Error in integralCalc/iterateScalarValued (line 314)

fx = FUN(t);

Error in integralCalc/vadapt (line 132)

[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)

[q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 87)

Q = integralCalc(fun,a,b,opstruct);

Error in test (line 80)

h11=integral(h, 0, 10); %inf

Related documentation

>>

[Christian]

Arrays have incompatible sizes for this operation.

Error in fur_elise (line 11)

mel = t + b;

Related documentation

[Walter Roberson]

A=[  -0.00767
    0.00696
    -0.00031
    0.00005
    -0.00001
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000];
B=[ 0.12000
    0.00209
    -0.00019
    0.00001
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000
    .00000E+00
    .00000E+00
    .00000E+00
    .00000E+00
    0.00000
    0.00000
    0.00000
    0.00000
    0.00000];
a = 1 ;   %RADIUS
L=.4;
c =-a/L;
b =a/L;
m =a*200; % NUMBER OF INTERVALS
[x,y]=meshgrid((c:(b-c)/m:b),(c:(b-c)/m:b)');
[I J]=find(sqrt(x.^2+y.^2)<(a-.1));
if ~isempty(I)
    x(I,J) = 0;
    y(I,J) = 0;
end
r=sqrt(x.^2+y.^2);
t=atan2(y,x);
warning off
chi=1;alpha=1;ab1=1;
k=sqrt(chi.^2+alpha.^2);   %xi=sqrt(t.^2+k.^2);
%for i=2:5
i=10;
h=@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).* (t1.*(exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))-exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i-1).*t1.^(i-1).*exp(-ab1.*t1)./factorial( i ))+((-1).^(i).*sqrt(pi.*k./2/t1.^2).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2, sqrt(t1.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))));
v=@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).* ((-t1.*exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))+sqrt(t1.^2+k.^2).*exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i).*t1.^(i-1).*exp(-ab1.*t1)./factorial( i ))+((-1).^(i-1).*sqrt(pi.*k./2/(t1.^2+k.^2)).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2, sqrt(t.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))));

whos

  Name           Size              Bytes  Class              Attributes

  A             29x1                 232  double                       
  B             29x1                 232  double                       
  I           4053x1               32424  double                       
  J           4053x1               32424  double                       
  L              1x1                   8  double                       
  a              1x1                   8  double                       
  ab1            1x1                   8  double                       
  alpha          1x1                   8  double                       
  b              1x1                   8  double                       
  c              1x1                   8  double                       
  chi            1x1                   8  double                       
  cmdout         1x33                 66  char                         
  h              1x1                  32  function_handle              
  i              1x1                   8  double                       
  k              1x1                   8  double                       
  m              1x1                   8  double                       
  r            201x201            323208  double                       
  t            201x201            323208  double                       
  v              1x1                  32  function_handle              
  x            201x201            323208  double                       
  y            201x201            323208  double                       

h11=integral(h, 0, 10); %inf

Arrays have incompatible sizes for this operation.

Error in solution (line 77)
h=@(t1)((t1-sqrt(t1.^2+k.^2)).^(-1).* (t1.*(exp(-sqrt(t1.^2+k.^2).*(r.*cos(t)+ab1))-exp(-t1.*(r.*cos(t)+ab1))).*(((-1).^(i-1).*t1.^(i-1).*exp(-ab1.*t1)./factorial( i ))+((-1).^(i).*sqrt(pi.*k./2/t1.^2).*exp(-ab1.*sqrt(t1.^2+k.^2)).*gegenbauerC(i,-1./2, sqrt(t1.^2+k.^2)./k))).*t1.*besselj(1,t1.*r.*sin(t))));

Error in integralCalc/iterateScalarValued (line 314)
fx = FUN(t);

Error in integralCalc/vadapt (line 132)
[q,errbnd] = iterateScalarValued(u,tinterval,pathlen);

Error in integralCalc (line 75)
[q,errbnd] = vadapt(@AtoBInvTransform,interval);

Error in integral (line 87)
Q = integralCalc(fun,a,b,opstruct);

r is 201 x 201 and t is the same size, and ab1 is scalar, so r.*cos(t)+ab1 is well defined and returns 201 x 201.

t1 is the parameter passed in by integral() . integral() passes in a vector of varying length -- it is not even the same length for every invocation of the function within one integral() call. The (t1-sqrt(t1.^2+k.^2)).^(-1).* (t1.*(exp(-sqrt(t1.^2+k.^2) parts of the expression are well-formed, but of unknown size -- in particular they are probably not the same size as (r.*cos(t)+ab1)

[Walter Roberson]

If you change the calls to

h11=integral(h, 0, 10, 'arrayvalued', true); %inf
h22=integral(v, 0, 10, 'arrayvalued', true);%inf

then the integrals can proceed -- but they take a long time. Long enough that I do not know yet whether the steps after that would be compatible with the values returned by those statements.

[Walter Roberson]

psi1=psi1+(Ai.*(r.^(-i+1)+h11)+(r.^(1./2).*besselk(i-1./2,r.*alpha)+h22).*Bi).*gegenbauerC(i,-1./2, cos(t));

Variable psi1 is not defined at that point in the code. For debugging purposes we can initialize it to 0 before that line.

[DH1,hh2]=contour(x,y,psi,10,'-k');

Variable psi is not defined at that point in the code. For debugging purposes we can use psi1 instead.

psi1=psi1+(Ai.*(r.^(-i+1)+h11)+(r.^(1./2).*besselk(i-1./2,r.*alpha)+h22).*Bi).*gegenbauerC(i,-1./2, cos(t));

That comes out as all zero, because Ai and Bi are zero.