<html><head><meta http-equiv="Content-Type" content="text/html charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space;" class=""><div class="">Hello,<br class=""><br class="">I think I have found a bug in the behaviour of ITK’s symmetric eigensystem. In the attached link (<a href="http://olivier.commowick.org/tmp/Test.zip" class="">http://olivier.commowick.org/tmp/Test.zip</a>) you will find an example code and matrix that fails with it and not with VNL’s symmetric eigensystem (run on macOS 10.12, Xcode 8.2.1).<br class=""><br class="">It seems to me that there are subtle differences in between ITK and VNL versions, VNL being more robust (and the same as what Matlab gives, I didn’t test python). Here is what I have found so far :<br class=""><br class="">- the hypothenuse function (hypot for ITK and pythag for VNL) do not get always the exact same result up to a very far away decimal<br class="">- some loop indices were wrong I think in the ITK eigensystem. Some do vary in the VNL version, not in ITK’s equivalent<br class=""><br class="">In the patch file attached, I have made a set of corrections that make things in ITK eigensystem closer to VNL but still the ITK version fails. I have to say I am getting completely lost here so I am calling upon the community for help. If anyone knows this enough, could you try the test attached and see if you have any ideas of what is going on. That would help us a lot.<br class=""><br class="">The test attached can be compiled with cmake and run with the text file included : ./test_eigensystem data.txt <br class=""><br class="">Thanks in advance<br class="">Olivier<br class=""><br class=""></div><div class=""></div></body></html>