[Insight-users] No test for svd? Update.

[Frank Miller]

This appears to be an upstream problem. I downloaded vxl and there is no 
"test_svd" there either.

However, on a more positive note, I recompiled itk with -O2 
optimization, instead of -O3, and the problem seems to be gone.

Optimization bugs are the worst.

Frank

Frank Miller wrote:
> Greetings,
> 
> Im having a very strange problem with the svd algorithm. My experience 
> is consistent with the problem mentioned in the comments at line 92 in 
> the file
> 
>   Utilities/vxl/core/vnl/algo/vnl_svd.txx
> 
> If this is the source of my problem, the comments state that the 
> "core/vnl/tests/test_svd" should have failed. However, I cant find this 
> test. I turned BUILD_TESTING on and there is no "test_svd". There is a 
> "vnl_test_all" which when run states that
> 
>   The registered test names are:
>      test_bignum
>      test_complex
>      test_inverse
>      test_diag_matrix
>      test_file_matrix
>      test_finite
>      test_math
>      test_matrix
>      test_matrix_exp
>      test_matrix_fixed
>      test_matrix_fixed_ref
>      test_numeric_traits
>      test_rational
>      test_real_polynomial
>      test_resize
>      test_sample
>      test_sym_matrix
>      test_transpose
>      test_fastops
>      test_vector
>      test_vector_fixed_ref
>      test_gamma
>      test_random
>      test_arithmetic
>      test_hungarian_algorithm
>      test_integrant
>      test_bessel
> 
> all of which pass on my system.
> 
> How can I run this test on my system? Is this test not being run by 
> Dart? Any help would be much appreciated,
> 
> 
> Below is some more details if you are interested.
> 
> My system is 32bit Linux, my cpu is Intel Xeon, my compiler is 
> gcc-4.1.1, and Im using itk-2.8.1.
> 
> When updating any pipeline that has a particular image as an input, I 
> get the following message printed to standard error:
> 
>   /home/frank/opt/src/itk/2.8.1/InsightToolkit-2.8.1/Utilities/vxl/core/
>   vnl/algo/vnl_svd.txx: suspicious return value (3) from SVDC
>   /home/frank/opt/src/itk/2.8.1/InsightToolkit-2.8.1/Utilities/vxl/core/
>   vnl/algo/vnl_svd.txx: M is 3x3
>   M = [ ...
>    0.9788850881159 -0.1393774186727 -5.6885917960696
>    1.0096153244376  0.1351351488847  5.5154448566180
>    0.0000000300788 -1.3927854783833  1.1043995835373  ]
> 
> The correct singular values for M are, according to octave, 8.0000, 
> 1.4063, 1.4062.
> 
> Here is a trace of the image:
> 
>   OrientedImage (0x83f0000)
>     RTTI typeinfo:   itk::OrientedImage<int, 3u>
>     Reference Count: 6
>     Modified Time: 714
>     Debug: Off
>     Observers:
>       none
>     Source: (0x83eb4c8)
>     Source output index: 0
>     Release Data: On
>     Data Released: False
>     Global Release Data: Off
>     PipelineMTime: 238
>     UpdateMTime: 715
>     LargestPossibleRegion:
>       Dimension: 3
>       Index: [0, 0, 0]
>       Size: [256, 256, 17]
>     BufferedRegion:
>       Dimension: 3
>       Index: [0, 0, 0]
>       Size: [256, 256, 17]
>     RequestedRegion:
>       Dimension: 3
>       Index: [0, 0, 0]
>       Size: [256, 256, 17]
>     Spacing: [1.40625, 1.40625, 7.99999]
>     Origin: [-37.7846, -241.523, 175.909]
>     Direction:
>   0.696096 -0.0991128 -0.711074
>   0.717949 0.0960961 0.689431
>   2.13894e-08 -0.990425 0.13805
> 
>     PixelContainer:
>       ImportImageContainer (0x83f01b8)
>         RTTI typeinfo:   itk::ImportImageContainer<unsigned long, int>
>         Reference Count: 1
>         Modified Time: 288
>         Debug: Off
>         Observers:
>           none
>         Pointer: 0x83f01e8
>         Container manages memory: true
>         Size: 1114112
>         Capacity: 1114112
> 
> Thanks,
> 
> Frank
> _______________________________________________
> Insight-users mailing list
> Insight-users at itk.org
> http://www.itk.org/mailman/listinfo/insight-users