bad exit() from LinearAlgebra.singular_value_decomposition

Bug #195396 reported by Mark Poolman
2
Affects Status Importance Assigned to Milestone
Python Numeric
Unknown
Unknown
python-numeric (Ubuntu)
Confirmed
Medium
Unassigned

Bug Description

Certain medium size (e.g. ~50x100) matrices cause LinearAlgebra.singular_value_decomposition() to exit python completely with the message:

        ** On entry to DGESDD parameter number 12 had an illegal value

This only seems to happens when the full_matrices parameter is True.

The behaviour is the same using numarray instead of numeric.

The behaviour is the same for x86 and amd64 platforms.

This is not the same, but may concievably be related to, https://bugs.launchpad.net/scipy/+bug/134277

The bug is demonstrated in the attached file.

Any help gratefully recieved,

/best/*

Mark

Revision history for this message
Mark Poolman (mgpoolman) wrote :
Revision history for this message
Mark Poolman (mgpoolman) wrote :
Revision history for this message
Matthias Klose (doko) wrote :

please could you check if python-numpy is an alternative for you?

Changed in python-numeric:
importance: Undecided → Medium
status: New → Confirmed
Revision history for this message
Mark Poolman (mgpoolman) wrote : Re: [Bug 195396] Re: bad exit() from LinearAlgebra.singular_value_decomposition

On Wed, 2008-02-27 at 12:00 +0000, Matthias Klose wrote:
> please could you check if python-numpy is an alternative for you?
>
>

It wouldn't be an ideal solution as it's not in the dapper repositories,
and for various strategic reasons I want to avoid (or at least minimise)
the number of tarball installs connected with this project.

I'll have a quick look and see if I can replicate the problem with numpy
though.

More soon and thanks for your attention,

Mark

> ** Also affects: python-numeric via
> http://sourceforge.net/support/tracker.php?aid=923048
> Importance: Unknown
> Status: Unknown
>
> ** Changed in: python-numeric (Ubuntu)
> Importance: Undecided => Medium
> Status: New => Confirmed
>
--
Dr. Mark Poolman
Senior Research Fellow
Cell Systems Modelling Group
School of Life Science
Oxford Brookes University

Revision history for this message
Mark Poolman (mgpoolman) wrote : MORE: Re: [Bug 195396] Re: bad exit() from LinearAlgebra.singular_value_decomposition

On Wed, 2008-02-27 at 12:00 +0000, Matthias Klose wrote:
> please could you check if python-numpy is an alternative for you?
>

OK, fresh install of NumPy-1.0.4, nothing else to install, same machine
as previously, the bug is absent, BUT now affected by
https://bugs.dogfood.launchpad.net/scipy/+bug/134277 (which is why I'd
fallen back to numarray in the first place ;-) ).

Cheers,

Mark

--
Dr. Mark Poolman
Senior Research Fellow
Cell Systems Modelling Group
School of Life Science
Oxford Brookes University

Revision history for this message
Matthias Klose (doko) wrote :

does this mean you don't see the error when the atlas packages are not installed?
would it be possible to recheck on hardy (with or without atlas)? all fortran libs are now built using gfortran instead of g77.

Revision history for this message
Mark Poolman (mgpoolman) wrote : Re: [Bug 195396] Re: bad exit() from LinearAlgebra.singular_value_decomposition

On Wed, 2008-02-27 at 15:08 +0000, Matthias Klose wrote:
> does this mean you don't see the error when the atlas packages are not installed?

Assume you mean [Bug 134277] ?
The test script for that analysed two matrices, one of which was OK with
atlas but not without, the other failed with atlas but was OK without.
This is the same for the numpy version (attached).

[Bug 195396] is not affected by presence/absence of atlas.

> would it be possible to recheck on hardy (with or without atlas)? all
> fortran libs are now built using gfortran instead of g77.
>

I can have a look for diagnostic purposes (probably not 'till next week
though), but simply upgrading the hardy would not be at all a preferred
solution (I do want the LTS).

My solution to date has been to simply do the error check and fall back
to Gauss-Jordan (+/- Gramm-Schmidt) if it fails. Although this isn't the
most elegant solution, it does yield a tolerable performance, so maybe I
should just stick with that.

/best/*

Mark

--
Dr. Mark Poolman
Senior Research Fellow
Cell Systems Modelling Group
School of Life Science
Oxford Brookes University

Revision history for this message
Matthias Klose (doko) wrote :

I tend to close this report as won't fix. there is no Numeric development anymore, the next LTS is released.

To post a comment you must log in.
This report contains Public information  
Everyone can see this information.

Other bug subscribers

Remote bug watches

Bug watches keep track of this bug in other bug trackers.