[Insight-developers] Problems with the vnl library

[Bill Hoffman]

I would like to add here, that instead of using vnl_matrix_inverse,
use vnl_svd directly which has methods that will allow you to check
for singular matrices.   I think the matrix inverse was meant to 
provide a simple interface to an inverse.   However, it uses svd anyway,
so, it won't change your code to switch.

For an online reference to vnl, you can see the
section on Numerics in the vxl book:

http://www.robots.ox.ac.uk/~vxl/book/

  
At 01:41 PM 12/15/00 -0500, Paul Hughett wrote:


> > 1. vnl_matrix_inverse does not check for singular matrices.
>
>This is indeed a problem with vnl.  They don't appear to have any
>current plans to fix it, though I have the impression they would
>be more than happy to accept someone else's fix.