[Insight-users] Matrix inversion

[Tom Vercauteren]

Hey Vincent,

Not sure what the status is exactly but there has been some discussion
about it on the bug tracker:
http://www.itk.org/Bug/view.php?id=7531

Tom

On Tue, Mar 29, 2011 at 10:56, Vincent Garcia <vincent.garcia at inria.fr> wrote:
> Hi,
>
> I'm using itk::Matrix in my code and I'm not happy the way it's done.
> Let me explain, I have the following code :
>
> itk::Matrix<double, 4, 4>  A;
> ...
> itk::Matrix<double, 4, 4>  B = A>GetInverse();
>
> This doesn't work because GetInverse return a vnl_matrix.
> Why returning a vnl_matrix and not an itk::Matrix?
>
> So, I can write my code using a vnl_matrix, and then copy everything into
> and itk::Matrix.
> I don't like it, it's not clean.
> My question is : what is the best way to use itk::Matrix (e.g. invert a
> matrix) and be sure that all the memory is free afterwards?
>
> Thanks!
>
> Vincent Garcia
>
> _____________________________________
> 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
>
>