Re: [eigen] cache-friendly matrix inverse

[ Thread Index | Date Index | More Archives ]

[ Apologies to Benoit; I sent this to him but it should have gone to the
  mailing list ]

Good point. I think you're right for very large sizes, asymptotically
they're the same.
But for a size like 10x10, the matrix-vector product will be much faster.

Can you give some proof for that?

Perhaps you're right (though I think the cross-over point is closer to 3 that 10), but it flies against what people in numerical analysis have been teaching for the last thirty years. I don't know myself, I only do matrices of larger size, and for e.g. 100 x 100 I can only repeat what Christian said:
* LU decomp instead of inverse takes the same time, AND
* LU decomp instead of inverse is more precise


Mail converted by MHonArc 2.6.19+