
<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="">Hi,<div class=""><br class=""></div><div class="">Thanks for your test. There might be a bug in my test or the way you ran it since for me it doesn’t take a few minutes but rather less than a second. That being said since it produces nans it might be that on some OS it goes into an infinite loop...</div><div class=""><br class=""></div><div class="">Olivier</div><div class=""><br class=""><div><blockquote type="cite" class=""><div class="">On 2 Feb 2017, at 00:02, Dženan Zukić <<a href="mailto:dzenanz@gmail.com" class="">dzenanz@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Hi Olivier,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br class=""></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">I tried running you example. One thing struck me right away: ITK version is much slower (after a few minutes running in release mod I killed it) than VNL version (about a second in debug mode) for your input.</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br class=""></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">I can't comment the patch content or the numeric stability. Perhaps Brad and Hans can pitch in?</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small"><br class=""></div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Regards,</div><div class="gmail_default" style="font-family:verdana,sans-serif;font-size:small">Dženan</div></div><div class="gmail_extra"><br class=""><div class="gmail_quote">On Tue, Jan 31, 2017 at 11:24 AM, Olivier Commowick <span dir="ltr" class=""><<a href="mailto:Olivier.Commowick@inria.fr" target="_blank" class="">Olivier.Commowick@inria.fr</a>></span> wrote:<br class=""><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word" 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" target="_blank" class="">http://olivier.commowick.org/<wbr class="">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<span class="HOEnZb"><font color="#888888" class=""><br class="">Olivier<br class=""><br class=""></font></span></div><div class=""></div></div><br class=""><div style="word-wrap:break-word" class=""><div class=""></div><div class=""><br class=""></div><div class="">---<br class=""><br class="">Olivier Commowick, Ph.D.<br class="">Research Scientist<br class="">INRIA Rennes - Bretagne Atlantique, VISAGES Team<br class="">Campus de Beaulieu<br class="">35042 Rennes<br class="">FRANCE<br class=""><br class="">Phone: <a href="tel:+33%202%2099%2084%2025%2092" value="+33299842592" target="_blank" class="">+33 2 99 84 25 92</a><br class="">Email: <a href="mailto:Olivier.Commowick@inria.fr" target="_blank" class="">Olivier.Commowick@inria.fr</a><br class="">Web: <a href="http://olivier.commowick.org/" target="_blank" class="">http://olivier.commowick.org/</a><br class=""><br class=""><br class=""><br class=""><br class=""></div>

<br class=""></div><br class="">______________________________<wbr class="">_______<br class="">

Powered by <a href="http://www.kitware.com/" rel="noreferrer" target="_blank" class="">www.kitware.com</a><br class="">

<br class="">

Visit other Kitware open-source projects at<br class="">

<a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank" class="">http://www.kitware.com/<wbr class="">opensource/opensource.html</a><br class="">

<br class="">

Kitware offers ITK Training Courses, for more information visit:<br class="">

<a href="http://www.kitware.com/products/protraining.php" rel="noreferrer" target="_blank" class="">http://www.kitware.com/<wbr class="">products/protraining.php</a><br class="">

<br class="">

Please keep messages on-topic and check the ITK FAQ at:<br class="">

<a href="http://www.itk.org/Wiki/ITK_FAQ" rel="noreferrer" target="_blank" class="">http://www.itk.org/Wiki/ITK_<wbr class="">FAQ</a><br class="">

<br class="">

Follow this link to subscribe/unsubscribe:<br class="">

<a href="http://public.kitware.com/mailman/listinfo/insight-users" rel="noreferrer" target="_blank" class="">http://public.kitware.com/<wbr class="">mailman/listinfo/insight-users</a><br class="">

<br class=""></blockquote></div><br class=""></div>

</div></blockquote></div><br class="">

<br class=""></div></body></html>