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)

History Show all quoted text
  Tue Apr 14 18:37:52 2009 ibmalone [...] cantab.net - Ticket created
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>
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
  Thu Apr 16 04:40:28 2009 jonathan [...] leto.net - Correspondence added
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>
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
  Thu Apr 16 04:40:28 2009 The RT System itself - Status changed from 'new' to 'open'
  Sat May 23 20:41:22 2009 jonathan [...] leto.net - Correspondence added
This has been fixed in the latest release, 0.20.
  Sat May 23 20:41:26 2009 jonathan [...] leto.net - Status changed from 'open' to 'resolved'