Re: [eigen] tri1 = tri2*tri1

[ Thread Index | Date Index | More Archives ]

>> and currently I use dense square matrices for that with
>>   temp.tri = b*a.tri
>>   a.swap(temp)
> [...]
> I'm afraid there is nothing better at the moment. What are the typical
> sizes of your matrices?

They are rather small and range from 2 to about 100. So I guess there is
not much benefit in trying to optimise this.

I'm using them for block Krylov subspace methods where at some point QR
decompositions start to become the predominant computations.


Mail converted by MHonArc 2.6.19+