Asked by Clarisha Nijman
on 2 Jan 2019

Dear all,

I have a function that returns me a PMatrix. But sometimes that matrix is full with nan's. The cause is known, but the best solution I found is to allow the function to compute such kind of "nan" matrices and to replace that "NaN matrix in the end by the identity matrix.

This is my actual code that works fine, uptill now.

Som=sum(sum(PMatrix));

Som(isnan(Som))=1;

if Som==1

PMatrix=eye(k,k);

end

I find this code clumsy and would prefer a more straight forward code saying:

if sum(sum(PMatrix))=NaN

PMatrix=eye(k,k)

end

But I can't figure out how to compose such a straight forward code. Does anyone has a suggestion for me?

Thank you in advance.

Answer by Steven Lord
on 2 Jan 2019

Accepted Answer

Which release are you using? For release R2018b or later:

if any(isnan(PMatrix), 'all')

end

For release R2018a or earlier:

if any(isnan(Pmatrix(:)))

end

Clarisha Nijman
on 3 Jan 2019

Thanks a lot!!!!,

This works for me the best! I agree really straight forward!

Steven Lord
on 3 Jan 2019

NaN == NaN % returns false

Clarisha Nijman
on 4 Jan 2019

Sign in to comment.

Answer by GT
on 2 Jan 2019

Edited by per isakson
on 3 Jan 2019

Have you tried using the function isnan?

Clarisha Nijman
on 3 Jan 2019

Yes, but I do not want to transfer the entries in ones, It would imply a matrix full of ones. I just need the identity matrix. I figured already out that sum(sum(P))==nan. So an if loop saying:

D+sum(sum(P))==nan

if D==nan

P=eye(k)

end

but it does not work. The if statement is the problem. How should it be corrected?

Sign in to comment.

Opportunities for recent engineering grads.

Apply Today
## 0 Comments

Sign in to comment.