Skip Menu |
 

This queue is for tickets about the Math-GSL CPAN distribution.

Report information
The Basics
Id: 45044
Status: resolved
Priority: 0/
Queue: Math-GSL

People
Owner: Nobody in particular
Requestors: ibmalone [...] cantab.net
Cc:
AdminCc:

Bug Information
Severity: (no value)
Broken in: (no value)
Fixed in: (no value)



Subject: math-gsl eigenpair bug
Date: Tue, 14 Apr 2009 23:37:21 +0100
To: bug-Math-Gsl [...] rt.cpan.org
From: Ian Malone <ibmalone [...] cantab.net>
Download (untitled) / with headers
text/plain 858b
Hi, It seems the Math::GSL::Matrix eigenpair method only returns the first two eigenvectors, I've tried patching the module as below and this seems to allow larger matrices: v5.10.0 built for x86_64-linux-thread-multi on Fedora, both Math::GSL 0.18 and 0.19 --- ../Math-GSL-0.18//lib/Math/GSL/Matrix.pm 2009-03-24 04:53:10.000000000 +0000 +++ ../Math-GSL-0.18/blib/lib/Math/GSL/Matrix.pm 2009-04-14 23:30:30.000000000 +0100 @@ -1653,7 +1653,7 @@ gsl_eigen_nonsymmv($self->raw, $vector->raw, $evec->raw, $eigen); - my $eigenvectors = [ map { $evec->col($_)->as_vector } (0,1) ]; + my $eigenvectors = [ map { $evec->col($_)->as_vector } (0 .. $r-1) ]; my $x = gsl_vector_complex_real($vector->raw); # vector of real components my $y = gsl_vector_complex_imag($vector->raw); # vector of imaginary components Thanks imalone
Subject: Re: [rt.cpan.org #45044] math-gsl eigenpair bug
Date: Thu, 16 Apr 2009 01:39:53 -0700
To: bug-Math-Gsl [...] rt.cpan.org
From: Jonathan Leto <jonathan [...] leto.net>
Download (untitled) / with headers
text/plain 1.5k
Just fixed in latest git bleed branch. This will be in 0.20. Thanks for reporting this bug Ian! Duke On Tue, Apr 14, 2009 at 3:37 PM, Ian Malone via RT <bug-Math-Gsl@rt.cpan.org> wrote: Show quoted text
> Tue Apr 14 18:37:52 2009: Request 45044 was acted upon. > Transaction: Ticket created by ibmalone@cantab.net >       Queue: Math-Gsl >     Subject: math-gsl eigenpair bug >   Broken in: (no value) >    Severity: (no value) >       Owner: Nobody >  Requestors: ibmalone@cantab.net >      Status: new >  Ticket <URL: https://rt.cpan.org/Ticket/Display.html?id=45044 > > > > Hi, > > It seems the Math::GSL::Matrix eigenpair method only returns the first > two eigenvectors, I've tried patching the module as below and this seems > to allow larger matrices: > > v5.10.0 built for x86_64-linux-thread-multi on Fedora, both Math::GSL > 0.18 and 0.19 > > --- ../Math-GSL-0.18//lib/Math/GSL/Matrix.pm    2009-03-24 > 04:53:10.000000000 +0000 > +++ ../Math-GSL-0.18/blib/lib/Math/GSL/Matrix.pm        2009-04-14 > 23:30:30.000000000 +0100 > @@ -1653,7 +1653,7 @@ > >      gsl_eigen_nonsymmv($self->raw, $vector->raw, $evec->raw, $eigen); > > -    my $eigenvectors = [ map { $evec->col($_)->as_vector } (0,1) ]; > +    my $eigenvectors = [ map { $evec->col($_)->as_vector } (0 .. $r-1) ]; > >      my $x = gsl_vector_complex_real($vector->raw); # vector of real > components >      my $y = gsl_vector_complex_imag($vector->raw); # vector of > imaginary components > > Thanks > imalone > >
-- [---------------------] Jonathan Leto jaleto@gmail.com
This has been fixed in the latest release, 0.20.


This service is sponsored and maintained by Best Practical Solutions and runs on Perl.org infrastructure.

Please report any issues with rt.cpan.org to rt-cpan-admin@bestpractical.com.