Boost logo

Ublas :

From: Rick Harding (richard_harding_at_[hidden])
Date: 2007-09-05 22:46:36


Hi,

Using lu_factorize() on the following matrix will cause an assertion
to be thrown:

A=
[8,8](
(138.625,14.2945,.00134662,-.00731894,-.00731894,-.00731894,2.717e-018,-.015625),
(14.2945,5.24033,-0.0193141,0.00646566,0.00646566,0.00646566,0,0.125),
(0.00134662,-0.0193141,47.8208,0.0470879,0.0470879,0.0470879,0,-0.0078125),
(-0.00731894,0.00646566,0.0470879,0.477962,0,0,0,0),
(-0.00731894,0.00646566,0.0470879,0,0.477962,0,0,0),
(-0.00731894,0.00646566,0.0470879,0,0,0.477962,0,0),
(2.717e-018,0,0,0,0,0,0.477962,0),
(-0.015625,0.125,-0.0078125,0,0,0,0,0)
)

Assertion failed in file C:\boost\numeric\ublas\lu.hpp at line 148:
pm (i) == i_norm_inf

The matrix is not singular. Have I done something wrong or is there a problem
in lu_factorize()?

I'm using uBLAS found in Boost 1.33.1.

Thanks!