Re: [ublas] Extending Numeric Bindings: UMFPACK

Subject: Re: [ublas] Extending Numeric Bindings: UMFPACK
From: Bruno Abreu Calfa (bacalfa_at_[hidden])
Date: 2011-11-13 21:19:46

Hi, Thomas,

Could you please attach your code? I called the "set_filled" method after
umf::getNumeric but I'm getting a runtime exception. I don't understand
what the arguments filled1 and filled2 stand for. So I'm also unsure what
to put, "n_row+1", "n_col+1", "std::min(n_row, n_col)+1", or...

Thank you very much!

Bruno

On Sun, Nov 13, 2011 at 8:45 PM, Thomas Klimpel <Thomas.Klimpel_at_[hidden]
> 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
> _______________________________________________
> ublas mailing list
> ublas_at_[hidden]
> http://lists.boost.org/mailman/listinfo.cgi/ublas
> Sent to: bacalfa_at_[hidden]
>