| [eigen] inverse() method for TriangularView class? |
[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen@xxxxxxxxxxxxxxxxxxx
- Subject: [eigen] inverse() method for TriangularView class?
- From: Douglas Bates <bates@xxxxxxxxxxxxx>
- Date: Mon, 31 Oct 2011 11:43:00 -0500
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=mime-version:sender:date:x-google-sender-auth:message-id:subject :from:to:content-type; bh=wMqfmWe7/VV7dQ95R62g/UzQNUu2TTjMtOB6fQqoGKc=; b=OlEGm1vF3d8qkY4jQKwKyAIpQAmadY5eDc2P8ORPhqDMB1KuuClWK6QyU0t74FjloF LNeSdZwbr0NhjfThdugZJ33R5OerPWXiPav74NHlJcd/HsBjoEI4+++xvDVmZAEbXmFJ pKQ8PP553tblV/m/3oGJTmX/f6KQQ9qCmdDYE=
The TriangularView class has solve methods but not an explicit inverse
method, which could be as simple as
solve(DenseMatrixType::Identity(cols(), cols())) (at least I think
that is how one would write the general form, but I'm not sure). If
one wanted to get fancier there may be a slight advantage in taking
account of the triangularity of the returned value and solving reduced
linear systems for each column, as is done in the Lapack *trtri
subroutines.
Has such a method been considered, or even implemented and I have
somehow managed to overlook it?