Re: [eigen] Reusing a sparse matrix for SparseLU

[ Thread Index | Date Index | More Archives ]

hm... if you're thinking about using coeffs() << v0, v1, v2... then your matrix is likely very small. What is its typical size and number of non-zeros? If too small and not sparse enough better use a dense matrix/solver.


On Mon, Dec 3, 2018 at 2:28 PM Matthieu Brucher <matthieu.brucher@xxxxxxxxx> wrote:

Le lun. 3 déc. 2018 à 11:43, Alberto Luaces <alberto.luaces@xxxxxxxxx> a écrit :
Matthieu Brucher writes:

> Hi Alberto,
> Thanks for the answer. Unfortunately, this still makes me create
> another array,

How so?

I need to create an array for the coifs, don't I? Seems like I could use valuePtr as well.
But the speed is very bad and profiles indicate that there may be mallocs inside the factorisation or the solve :/

> I'd rather populate the elements directly.

You can do precisely that, by using the coeffs() vector.

Oh, you mean I could do coeffs() << v0, v1, v2...? 


Mail converted by MHonArc 2.6.19+