[Insight-users] vnl_generalized_eigensystem

[Luis Ibanez]

Hi Liu,

Your observation is correct.

The VNL class

     vnl_generalized_eigensystem

has not been wrapped inside an ITK class.

Your best option at this point is probably to use the
VNL class directly.

Please note that if you have ITK matrices, you can
get their equivalent Vnl matrices by calling the method:

     matrix.GetVnlMatrix()

You can as well, assign a Vnl matrix to an itkMatrix
with an expression such as

    itk::Matrix<double, 3,3,> mm;
    mm = vnlMatrixAA;

or using the constructor, as in:

    itk::Matrix<double, 3,3,> mm( vnlMatrixAA );

of course, assuming that their dimensionalities are consistent.


For an example of the use of

    vnl_generalized_eigensystem

you may want to take a look at:

Insight/Code/Common/itkBloxCoreAtomPixel.txx

in lines 217-218 and 283-284.


     Regards,


           Luis

------------------------------
On Wed, Dec 2, 2009 at 6:18 PM, liu jianfei <flair.liu at gmail.com> wrote:
> Dear all
>
>    I want to solve a generalized eigen system AX=\lambda BX. I noticed that
> it could be computed through vnl_generalized_eigensystem call. Most
> numerical computation of ITK are supposed to use the VNL library. I wonder
> whether ITK has the similar class to solve the system. Some symmetrical
> eigen-computation classes can be found from ITK documentation, but I could
> not find the corresponding function of vnl_generalized_eigensystem. Does ITK
> provide this function or I have to install VNL library to call it directly?
> Thanks
>
>
>   Looking forward to your reply
>
> --
> Sincerely
> Jianfei Liu
>
> _____________________________________
> Powered by www.kitware.com
>
> Visit other Kitware open-source projects at
> http://www.kitware.com/opensource/opensource.html
>
> Kitware offers ITK Training Courses, for more information visit:
> http://www.kitware.com/products/protraining.html
>
> Please keep messages on-topic and check the ITK FAQ at:
> http://www.itk.org/Wiki/ITK_FAQ
>
> Follow this link to subscribe/unsubscribe:
> http://www.itk.org/mailman/listinfo/insight-users
>
>