No subject

The following code:

  vnl_matrix<double> B =
    vnl_matrix_inverse<double>( N_matrix.transpose() * N_matrix) *
N_matrix.transpose() * data_matrix;

is from

Insight/Code/Common/
itkNonUniformBSpline.txx
(line 313).
http://public.kitware.com/cgi-bin/viewcvs.cgi/Code/Common/itkNonUniformBSpline.txx?revision=1.9&root=Insight&view=markup


It is trying to solve the linear system

             X  = N * B

by doing:

        B = ( N' * N )^-1  *  N'  * X


instead of invoking one of the linear solvers
that we are carrying in ITK via VXL.


----


   Luis

--000e0ce0f792df210404942b30ba
Content-Type: text/html; charset=UTF-8
Content-Transfer-Encoding: quoted-printable