<div dir="ltr"><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></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></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></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><div class="gmail_quote">On Tue, Jan 31, 2017 at 11:24 AM, Olivier Commowick <span dir="ltr"><<a href="mailto:Olivier.Commowick@inria.fr" target="_blank">Olivier.Commowick@inria.fr</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex"><div style="word-wrap:break-word"><div>Hello,<br><br>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">http://olivier.commowick.org/<wbr>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><br>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><br>- 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>- some loop indices were wrong I think in the ITK eigensystem. Some do vary in the VNL version, not in ITK’s equivalent<br><br>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><br>The test attached can be compiled with cmake and run with the text file included : ./test_eigensystem data.txt <br><br>Thanks in advance<span class="HOEnZb"><font color="#888888"><br>Olivier<br><br></font></span></div><div></div></div><br><div style="word-wrap:break-word"><div></div><div><br></div><div>---<br><br>Olivier Commowick, Ph.D.<br>Research Scientist<br>INRIA Rennes - Bretagne Atlantique, VISAGES Team<br>Campus de Beaulieu<br>35042 Rennes<br>FRANCE<br><br>Phone: <a href="tel:+33%202%2099%2084%2025%2092" value="+33299842592" target="_blank">+33 2 99 84 25 92</a><br>Email: <a href="mailto:Olivier.Commowick@inria.fr" target="_blank">Olivier.Commowick@inria.fr</a><br>Web: <a href="http://olivier.commowick.org/" target="_blank">http://olivier.commowick.org/</a><br><br><br><br><br></div>
<br></div><br>______________________________<wbr>_______<br>
Powered by <a href="http://www.kitware.com" rel="noreferrer" target="_blank">www.kitware.com</a><br>
<br>
Visit other Kitware open-source projects at<br>
<a href="http://www.kitware.com/opensource/opensource.html" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>opensource/opensource.html</a><br>
<br>
Kitware offers ITK Training Courses, for more information visit:<br>
<a href="http://www.kitware.com/products/protraining.php" rel="noreferrer" target="_blank">http://www.kitware.com/<wbr>products/protraining.php</a><br>
<br>
Please keep messages on-topic and check the ITK FAQ at:<br>
<a href="http://www.itk.org/Wiki/ITK_FAQ" rel="noreferrer" target="_blank">http://www.itk.org/Wiki/ITK_<wbr>FAQ</a><br>
<br>
Follow this link to subscribe/unsubscribe:<br>
<a href="http://public.kitware.com/mailman/listinfo/insight-users" rel="noreferrer" target="_blank">http://public.kitware.com/<wbr>mailman/listinfo/insight-users</a><br>
<br></blockquote></div><br></div>