|Re: [eigen] Inverse when the (dense) matrix has a known structure|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] Inverse when the (dense) matrix has a known structure
- From: Matthieu Brucher <matthieu.brucher@xxxxxxxxx>
- Date: Wed, 25 May 2016 15:44:30 +0100
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:in-reply-to:references:date:message-id:subject:from:to; bh=fA6P/qEqRn/xluuSoNZr+oMAkZXQTOI3w1NAPuBFta0=; b=RarKBOOc5Qtkzt71gmKUHjRxfoEmSM+mLUBb4MyvTnREtU1Rbj3G5PU8PbULWQk8xR dHZqhLdu41+x1bUl//coZ9PVBRRWWB3dABx9j77Ij/jDWgAkoRA0dPGphMDchylH8Xbl kYX/mhyIC1dlybOzr22Xnun6r217V1CR6l2/DFMPS5Sck+PpgdSrZuVC/qBSAODhYOCD KEeSxznTPcwk0DP8hi3FZ2U8YBu2G9vMZ3LzQJNyRS6aYSItYf4qBjDZ2VySLUZq/pCv iFUKF9vugnZtX/bTG4VS2kHnZ8FFNPaLVBvd8W7GAkPa68sxqQlETyVfhYWMUTmJpcMo QRfQ==
I tried using sympy to generate just the operations required, but it seems that even if there are far less operations that a generic inverse, the time spent int he inverse is higher than letting Eigen do it with the current pattern.
Maybe interesting if I have more holes in it though.
Thanks for the feedback!