[PYTHON MATRIX-SIG] RE: LinearAlgebra on Alpha.

[Doug Heisterkamp]

From: Rob.Hooft@embl-heidelberg.de
To: matrix-sig@python.org
Subject: [PYTHON MATRIX-SIG] LinearAlgebra on Alpha.
Mime-Version: 1.0 (generated by tm-edit 7.89)
Content-Type: text/plain; charset=US-ASCII
Sender: owner-matrix-sig@python.org
Precedence: bulk
Status: OR

Hi,

Are you using existing LAPACK libraries or the source sent with
NumPy?  

In either case, make sure that the typedef of f2c_complex and 
f2c_doublecomplex at the begining of lapack_litemodule.c matches
the form with which your LAPACK library was built.  Now that f2c.h is
distributed with NumPy, the lapack module should read it so all
the type changes are in one file.

Doug Heisterkamp
drh@oasis.unl.edu


> I just tried some LinearAlgebra routines on DEC ALPHA, and got some
> unexpected results: All zero's. What can I do to debug this?
>
> >>> import ni,Numeric
> >>> import Numeric.LinearAlgebra
> >>> LA=Numeric.LinearAlgebra
> >>> a=Numeric.array([[1.,0.,0.],[0.,1.,0.],[0.,0.,1.]])
> >>> LA.eigenvectors(a)
> ( 0.  0.  0.,  0.  0.  0.
>  0.  0.  0.
>  0.  0.  0.)
> >>> LA.singular_value_decomposition(a)
> ( 0.  0.  0.
>  0.  0.  0.
>  0.  0.  0.,  0.  0.  0.,  0.  0.  0.
>  0.  0.  0.
>  0.  0.  0.)
> >>> LA.generalized_inverse(a) 
>  0.  0.  0.
>  0.  0.  0.
>  0.  0.  0.
>
> Regards,
> -- 
> === Rob.Hooft@EMBL-Heidelberg.DE   http://www.Sander.EMBL-Heidelberg.DE/rob/ ==
> ==== In need of protein modeling?  http://www.Sander.EMBL-Heidelberg.DE/whatif/
> Validation of protein structures?  http://biotech.EMBL-Heidelberg.DE:8400/ ====
> == PGPid 0xFA19277D == Use Linux!  Free Software Rules The World! =============


=================
MATRIX-SIG  - SIG on Matrix Math for Python

send messages to: matrix-sig@python.org
administrivia to: matrix-sig-request@python.org
=================