On Sun, Nov 13, 2011 at 8:45 PM, Thomas Klimpel
<Thomas.Klimpel@synopsys.com> wrote:
Bruno Abreu Calfa wrote:
> Then, I ran the code in the attachment.
> The L matrix is all zeros, but if I print the report from UMFPACK
> it works and obviously shows L and U not all zeros.
> I am suspecting that the values are not being copied
> correctly to my matrices. Does anyone have any idea why?
The "filled1_" and "filled2_" members of the compressed_matrix Lx are not set correctly. I first checked that "Lx.value_data()" contains the requested data, then I added the call "Lx.set_filled(n_row+1,lnz1);" after the call to "umf::getNumeric(Lx, Ux, P, Q, Dx, &do_recip, Rs, Numeric);". Then everything seemed fine.
Note however that I haven't checked whether "filled1" must be set to "n_row+1", "n_col+1", "std::min(n_row, n_col)+1", or ...
Regards,
Thomas