[ Thread Index |
Date Index
| More lists.tuxfamily.org/eigen Archives
]
- To: eigen <eigen@xxxxxxxxxxxxxxxxxxx>
- Subject: Re: [eigen] Matrix power
- From: Gael Guennebaud <gael.guennebaud@xxxxxxxxx>
- Date: Thu, 13 Jun 2013 22:42:08 +0200
- Cc: jdh8@xxxxxxxxxxxxxx
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=20120113; h=mime-version:in-reply-to:references:from:date:message-id:subject:to :cc:content-type:content-transfer-encoding; bh=b064q1QAuItUqb3XIz7TmWalqrn6ZbcrNJBwpNLiYS8=; b=Q0uoI1XrDOfEoFk9PdQs1lGEQ+k1zCf/xrqD3LYWQxzwdYdqZHq09MPGFqA/1QafDT Qwe9RAeN3jUpU50an5fDb4HslR6i4UggqFVnmsTG2vTUOMAleY97J8XTVefGOhPsgU8E hcxIchfzlDP3J2FWlr/zMEMzt8d1n296XDqtNKa4rR9T846sfYa4fek/IQpEVUO+qTPD kYm41bryprhyZ1g915pym7l2exzq8TmJkk5dM7GVIXR2EXQ6KyGE9+Efmwc2ezBY+RDL j90bMytZ4j2Up9a/cSBtLTJDVz5CdCsgs9U19flTbd8tN47BCrs8ZJiaAHs7vhmOMr+3 3e9A==
Hi,
I answered your bug reports/comments. Jitse and Chen-Pang He can tell
you more about the details of this matrix-functions module. On a
general level, you will find the respective unit tests in the
unsupported/test/ folder. Feel free to extend them. We usually prefer
to exploit remarkable identities and the likes rather than explicit
comparisons to precomputed reference solutions.
On this page you will find many links explaining how to run the
unit-tests, how to deal with hg and submit patches, etc.:
http://eigen.tuxfamily.org/index.php?title=Category:Developer
cheers,
gael
On Thu, Jun 13, 2013 at 8:34 PM, Dale Lukas Peterson
<hazelnusse@xxxxxxxxx> wrote:
> I recently need to compute the matrix power and tried to use
> functionality provided in the unsupported module. Some stuff worked,
> but I ran in to both segfaults as well as some incorrect results (i
> think some sort of aliasing issue, I'm not sure). I filed a bug
> report [0] for the segfaults and if I can come up with a minimal
> working example of incorrect results I will file a bug report for that
> as well.
>
> That being said, who is the best person to talk to about a) getting
> this bug/bugs fixed, and b) getting this functionality into the core
> of Eigen. Gael started a bug report requesting this functionality
> [1], but it has been a year and I haven't seen anything.
>
> I haven't looked closely, but I'm guessing that implementing a good
> set of tests would be the first step, and then figuring out which ones
> don't pass and why. Would creating tests based on the results of
> Matlab output (or perhaps some small handworked examples) be
> appropriate?
>
> I'd like to help out on this, I'm not very comfy at all with hg, I'm
> fairly comfy with C++ but probably not as comfy with the template side
> of things, but am willing to learn. My IRC nick is DLPeterson, and I
> am on #eigen frequently.
>
> Luke
>
> [0] -- http://eigen.tuxfamily.org/bz/show_bug.cgi?id=614
> [1] -- http://eigen.tuxfamily.org/bz/show_bug.cgi?id=465
> --
> “People call me a perfectionist, but I'm not. I'm a rightist. I do
> something until it's right, and then I move on to the next thing.”
> ― James Cameron
>
>