[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: "[eigen] mailing list" <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: Re: [eigen] Matrix power
- From: Dale Lukas Peterson <hazelnusse@xxxxxxxxx>
- Date: Tue, 18 Jun 2013 19:59:51 -0700
- 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 :content-type; bh=meVMcx68J5eTOMTjykfjnjsf0jvrjSk1yT9V86sb1SU=; b=jPdGHSXqqqJwxtk0pHc+r8SPQloGZ6bUzkcfPBr8ENM6MYlqJVB7mBV1V4gEW32fC3 zMrJdXEExbI4HFbmytv2o70IZR6bvSg0Tlp3OJg+YqW6HCHRLCklLTglyXacS5Da5MnD EUynE3yKMHd18VdKOZblAwnH0Q1lQA/ZMWaG4tyu5rZ77DkuEi6ZrdJDCECKjI8D2dbS dEQlF4oFfSq4bKH4vqvvuZAJ8G5s9u5uSXWZvAo/BBCtwdN3JVKnLzLHgnsZcrtx9jEm mYZsSx6TbVJIGnjhMl6lRzg6LhNEnjKAFIxMjZeeavV9mUYzYw4SWKPqo6obwqgRAOs8 vYKw==
On Mon, Jun 17, 2013 at 1:46 AM, Jitse Niesen <jitse@xxxxxxxxxxxxxxxxx> wrote:
> Re b) I need to make a pass through all the code and clean it up and then
> it's ready to go in the core. What you can do is wait till 3.3 is forked
> (which as I understand will not take too long) and then pester me until it's
> done :)
>
> Out of interest, what do you use fractional matrix powers for?
The patch you attached to the bug report fixes one of the issues I was
having. However, now I ran into another. The following minimal
working example never converges when calling A.pow(0.5):
// compile with
// g++ -Wall -O2 -g -I/usr/include/eigen3 -o testpow testpow.cc
#include <Eigen/Dense>
#include <unsupported/Eigen/MatrixFunctions>
int main()
{
Eigen::MatrixXd A(4, 4);
A << 0, 0, 0, 0,
0, 0, 0, 0,
0, 0, 0.4115557167420392, -4.285896778479871,
0, 0, -4.285896778479871, 44.63286609452607;
Eigen::MatrixXd A_sqrt = A.pow(0.5);
}
I can file a separate bug report for this or attach it to the same one
if you like.
Luke