Frank,<br>
<br>
I have a Continuous itk build (BillsBasement) that is built with -O3
and gcc4.0. As you can see, the vnl_test_svd and vnl_test_qr tests both
fail. I used to build with -O2 and they passed.<br>
<br>
<a href="http://www.itk.org/Testing/Sites/BillsBasement/Linux-gcc40-release/20060910-1145-Continuous/Test.html">http://www.itk.org/Testing/Sites/BillsBasement/Linux-gcc40-release/20060910-1145-Continuous/Test.html</a><br>
<br>
My advice, in general, avoid -O3. It's amazing after all these years,
compiler writers still have optimization woes. I remember in 1969,
running an IBM 360/44 FORTRAN compiler. Almost all of the bug reports
were due to optimization failures. The workaround was always, turn off
the optimization.<br>
<br>
Bill<br>
<br><br><div><span class="gmail_quote">On 9/9/06, <b class="gmail_sendername">Luis Ibanez</b> <<a href="mailto:luis.ibanez@kitware.com">luis.ibanez@kitware.com</a>> wrote:</span><blockquote class="gmail_quote" style="border-left: 1px solid rgb(204, 204, 204); margin: 0pt 0pt 0pt 0.8ex; padding-left: 1ex;">
Hi Frank,<br><br>You may want to try with ITK CVS version,<br>A couple of VXL bugs were fixed recently that were related to infinte<br>loops that happen only in -O3 optimization.<br><br>See the following email in the developers list:
<br><a href="http://www.itk.org/mailman/private/insight-developers/2006-September/008389.html">http://www.itk.org/mailman/private/insight-developers/2006-September/008389.html</a><br><br>Please let us know what you find.<br>
<br><br> Thanks<br><br><br> Luis<br><br><br>--------------------<br>Frank Miller wrote:<br>> This appears to be an upstream problem. I downloaded vxl and there is no<br>> "test_svd" there either.<br>
><br>> However, on a more positive note, I recompiled itk with -O2<br>> optimization, instead of -O3, and the problem seems to be gone.<br>><br>> Optimization bugs are the worst.<br>><br>> Frank<br>>
<br>> Frank Miller wrote:<br>><br>>> Greetings,<br>>><br>>> Im having a very strange problem with the svd algorithm. My experience<br>>> is consistent with the problem mentioned in the comments at line 92 in
<br>>> the file<br>>><br>>> Utilities/vxl/core/vnl/algo/vnl_svd.txx<br>>><br>>> If this is the source of my problem, the comments state that the<br>>> "core/vnl/tests/test_svd" should have failed. However, I cant find
<br>>> this test. I turned BUILD_TESTING on and there is no "test_svd". There<br>>> is a "vnl_test_all" which when run states that<br>>><br>>> The registered test names are:<br>
>> test_bignum<br>>> test_complex<br>>> test_inverse<br>>> test_diag_matrix<br>>> test_file_matrix<br>>> test_finite<br>>> test_math<br>>> test_matrix
<br>>> test_matrix_exp<br>>> test_matrix_fixed<br>>> test_matrix_fixed_ref<br>>> test_numeric_traits<br>>> test_rational<br>>> test_real_polynomial<br>>> test_resize
<br>>> test_sample<br>>> test_sym_matrix<br>>> test_transpose<br>>> test_fastops<br>>> test_vector<br>>> test_vector_fixed_ref<br>>> test_gamma<br>
>> test_random<br>>> test_arithmetic<br>>> test_hungarian_algorithm<br>>> test_integrant<br>>> test_bessel<br>>><br>>> all of which pass on my system.<br>
>><br>>> How can I run this test on my system? Is this test not being run by<br>>> Dart? Any help would be much appreciated,<br>>><br>>><br>>> Below is some more details if you are interested.
<br>>><br>>> My system is 32bit Linux, my cpu is Intel Xeon, my compiler is<br>>> gcc-4.1.1, and Im using itk-2.8.1.<br>>><br>>> When updating any pipeline that has a particular image as an input, I
<br>>> get the following message printed to standard error:<br>>><br>>> /home/frank/opt/src/itk/2.8.1/InsightToolkit-2.8.1/Utilities/vxl/core/<br>>> vnl/algo/vnl_svd.txx: suspicious return value (3) from SVDC
<br>>> /home/frank/opt/src/itk/2.8.1/InsightToolkit-2.8.1/Utilities/vxl/core/<br>>> vnl/algo/vnl_svd.txx: M is 3x3<br>>> M = [ ...<br>>> 0.9788850881159 -0.1393774186727 -5.6885917960696<br>
>> 1.0096153244376 0.1351351488847 5.5154448566180<br>>> 0.0000000300788 -1.3927854783833 1.1043995835373 ]<br>>><br>>> The correct singular values for M are, according to octave, 8.0000,
<br>>> 1.4063, 1.4062.<br>>><br>>> Here is a trace of the image:<br>>><br>>> OrientedImage (0x83f0000)<br>>> RTTI typeinfo: itk::OrientedImage<int, 3u><br>>> Reference Count: 6
<br>>> Modified Time: 714<br>>> Debug: Off<br>>> Observers:<br>>> none<br>>> Source: (0x83eb4c8)<br>>> Source output index: 0<br>>> Release Data: On<br>
>> Data Released: False<br>>> Global Release Data: Off<br>>> PipelineMTime: 238<br>>> UpdateMTime: 715<br>>> LargestPossibleRegion:<br>>> Dimension: 3<br>>> Index: [0, 0, 0]
<br>>> Size: [256, 256, 17]<br>>> BufferedRegion:<br>>> Dimension: 3<br>>> Index: [0, 0, 0]<br>>> Size: [256, 256, 17]<br>>> RequestedRegion:<br>>> Dimension: 3
<br>>> Index: [0, 0, 0]<br>>> Size: [256, 256, 17]<br>>> Spacing: [1.40625, 1.40625, 7.99999]<br>>> Origin: [-37.7846, -241.523, 175.909]<br>>> Direction:<br>>>
0.696096 -0.0991128 -0.711074<br>>> 0.717949 0.0960961 0.689431<br>>> 2.13894e-08 -0.990425 0.13805<br>>><br>>> PixelContainer:<br>>> ImportImageContainer (0x83f01b8)<br>>> RTTI typeinfo: itk::ImportImageContainer<unsigned long, int>
<br>>> Reference Count: 1<br>>> Modified Time: 288<br>>> Debug: Off<br>>> Observers:<br>>> none<br>>> Pointer: 0x83f01e8<br>>> Container manages memory: true
<br>>> Size: 1114112<br>>> Capacity: 1114112<br>>><br>>> Thanks,<br>>><br>>> Frank<br>>> _______________________________________________<br>>> Insight-users mailing list
<br>>> <a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a><br>>> <a href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users</a><br>><br>> _______________________________________________
<br>> Insight-users mailing list<br>> <a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a><br>> <a href="http://www.itk.org/mailman/listinfo/insight-users">http://www.itk.org/mailman/listinfo/insight-users
</a><br>><br>><br><br><br>_______________________________________________<br>Insight-users mailing list<br><a href="mailto:Insight-users@itk.org">Insight-users@itk.org</a><br><a href="http://www.itk.org/mailman/listinfo/insight-users">
http://www.itk.org/mailman/listinfo/insight-users</a><br></blockquote></div><br>