Re: [ublas] Extending Numeric Bindings: UMFPACK

Subject: Re: [ublas] Extending Numeric Bindings: UMFPACK
From: Thomas Klimpel (Thomas.Klimpel_at_[hidden])
Date: 2011-11-13 20:45:36

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