|Re: [eigen] compilation failed for not aligned tiny matrix|
[ Thread Index |
| More lists.tuxfamily.org/eigen Archives
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: Re: [eigen] compilation failed for not aligned tiny matrix
- From: Benoit Jacob <jacob.benoit.1@xxxxxxxxx>
- Date: Tue, 6 Oct 2009 09:12:17 -0400
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:mime-version:received:in-reply-to:references :date:message-id:subject:from:to:content-type; bh=XWLYxGpJjM/LLLUlHXilO6gnlgIx1vKsqHO/oO26584=; b=hYAVM4XYuwULqu53CMM+91+MW1K5hoknmzI9P340xGNR5St/0E+lge1mnCv/3NbqrE eWiZiTqYEGs9ujftdaEpdHf5FLYEarC/UYPEEmRqW6AQKEO6OjMB/E33IqCRGXUmhDXj FGN5u7ajhh4zdt+f5/QIad5O00ODvxTB4GsFM=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=mime-version:in-reply-to:references:date:message-id:subject:from:to :content-type; b=i3EzF8cTOzZggSIDHXSeBd6BF+ZLs5CuqEMx6bUChxhmICkcGvwfSSalMoprnd/rQh v+G7liqLPukBqpVoGZtJVGOIEo22LLcCeg1iwPRDytlTGyF5JLPsiullwM++xCHzcpEU xbhd/uVtnN9TzKo+QYh0zViScPAnasiDCTQoc=
2009/10/6 Gael Guennebaud <gael.guennebaud@xxxxxxxxx>:
> ok, here the bug is that computeInverse() expect a PlainMatrixType which in
> your case is equivalent to TinyRMatrix but without the DontAlign option. Yes
> that's a bug in Eigen.
> But now, I'm wondering why computeInverse is not a template function ? Think
> about submatrices, or, e.g., to perform: m.transpose().computeInverse(inv);
> - here m.transpose() is row major and so computeInverse will expect a row
> major matrix.
Yes, all these functions that take a ptr/ref to result should be
template in the result type, you're perfectly right.
> Moreover I'm wondering whether computeInverse() is still useful since we can
> make inv = a.inverse() works in place just like compute inverse.
Indeed, but that part (Inverse.h) is the trickiest of the current
solve API refactoring (btw, sorry for being slow) for reasons i'll
explain, so it's the part that i'll do last. The problems is that
these API changes go together with a reformulation of the problem of
"checking for existence of a solution" and the case of inverting tiny
matrices is special here.