Relational operations problems (real and imag parts)

Hello all!

I want to replace this line

if imag(LAMBD(m,m)) ~= 0

into

if imag(LAMBD(m,m)) == 0

because the operation "~=" (even without imag()) is testing the imaginary part of LAMBD, but it doesn't work.

I have no idea why...

Thanks,

Kate

[Matt J]

I'm skeptical that it doesn't work. And you haven't provided any illustration to prove it...

Answer 1

[Ekaterina]

It seems that in this case "~=" doesn't mean testing the imginary part of LAMBD. if f(x) ~= 0 menas that f(x) is close to zero...

[Azzi Abdelmalek]

use

 abs(y)>eps

[Azzi Abdelmalek]

Ekaterina , dont add an answer, just add a comment

Answer 2

[Azzi Abdelmalek]

Maybe it's about numerical problems, for example try this

0.3-0.2==0.4-0.3

%you will find 0

Answer 3

[Ekaterina]

It would be difficult as I have 600 lines of code.

[Matt J]

Yes, and we wouldn't be interested in seeing all 600 lines of your code. However, only a small number of those lines are actually causing trouble. Therefore, it should be easy for you to construct a separate, small, simplified example based on those lines.