|
Server : Apache/2.2.17 (Unix) mod_ssl/2.2.17 OpenSSL/0.9.8e-fips-rhel5 DAV/2 PHP/5.2.17 System : Linux localhost 2.6.18-419.el5 #1 SMP Fri Feb 24 22:47:42 UTC 2017 x86_64 User : nobody ( 99) PHP Version : 5.2.17 Disable Function : NONE Directory : /proc/21572/root/usr/lib64/python2.4/site-packages/Numeric/ |
Upload File : |
mò
C1?c @ sC d k Z d k Z d k Z d k Z d k Z d k Z d e f d „ ƒ YZ h d d <d d <d d <d d
<Z d d d d
g Z h d d <d d <d d <d d <d d <d d <Z
h d d <d d <d d <d d <d d <d d <Z d d g d d g g Z d „ Z
d „ Z e i Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d d „ Z d „ Z d d „ Z d d „ Z d d „ Z d „ Z d d! „ Z e d" j or d# k Td$ „ Z e d% d& g d' d( g g ƒ Z! e d& d% g ƒ Z" e e! e" ƒ e! d) Z! e" d) Z" e e! e" ƒ n d S(* Nt LinAlgErrorc B s t Z RS( N( t __name__t
__module__( ( ( t; /usr/lib64/python2.4/site-packages/Numeric/LinearAlgebra.pyR s t fi t di t Fi t Di t st ct zt it lc G s[ d } d } x@ | D]8 } | i ƒ } t | t | ƒ } t | t | ƒ } q Wt | | S( Ni i (
t kindt precisiont arrayst at typecodet tt maxt _array_kindt _array_precisiont _array_type( R R R
R R ( ( R t _commonType s c G sœ f } xl | D]d } | i ƒ | j o# | t i t i | ƒ ƒ f } q
| t i t i | ƒ i | ƒ ƒ f } q
Wt | ƒ d j o | d Sn | Sd S( Ni i (
t cast_arraysR R R t typet copyt Numerict transposet astypet len( R R R R ( ( R t _castCopyAndTranspose# s #, c G s„ f } xT | D]L } | i ƒ | j o | t | ƒ f } q
| t | i | ƒ ƒ f } q
Wt | ƒ d j o | d Sn | Sd S( Ni i ( R R R R R t _fastCTR R ( R R R R ( ( R t _fastCopyAndTranspose5 s c G s8 x1 | D]) } t | i ƒ d j o
t d ‚ q q Wd S( Ni s Array must be two-dimensional( R R R t shapeR ( R R ( ( R t _assertRank2A s c G sA x: | D]2 } t | i ƒ t | i ƒ j o
t d ‚ q q Wd S( Ns Array must be square( R R R R" t minR ( R R ( ( R t _assertSquarenessF s c C sS t | i ƒ d j } | o | d d … t i f } n t | | ƒ t | ƒ | i d } | i d } | | i d j o
t d ‚ n t | | ƒ } t | d j o
t i } n
t i } t | | | ƒ \ } } t i | d ƒ } | | | | | | | | d ƒ } | d d j o
t d ‚ n | o t i | ƒ Sn t i | ƒ Sd S( Ni i s Incompatible dimensionsR t infos Singular matrix( R t bR" t one_eqR t NewAxisR# R R% t n_eqt n_rhsR R R R t lapack_litet zgesvt lapack_routinet dgesvR! t zerost pivotst resultst ravelt
multiarrayR ( R R' R1 R( R2 R+ R. R* R ( ( R t solve_linear_equationsN s*
!
c C s t | t i | i d ƒ ƒ S( Ni ( R5 R R t identityR" ( R ( ( R t inversek s c C sË t | ƒ t | ƒ t | ƒ } t | | ƒ } | i d } | i d } t | d j o
t
i } n
t
i
} | d | | | d ƒ } | d d j o
t d ‚ n t i t i t i | d d ƒƒ ƒ S( Ni i t LR& sK Matrix is not positive definite - Cholesky decomposition cannot be computedt k( R# R R% R R R R" t mt nR R, t zpotrfR. t dpotrfR2 R R R R t MLabt triu( R R: R2 R; R. R ( ( R t cholesky_decompositionp s
c
C s† t | ƒ t | ƒ t | ƒ } t d t | } t | | ƒ } | i d }
t i d | ƒ } t | d j oã t i } t i |
f | ƒ } t i |
f | ƒ } d } t i | f | ƒ } | d d |
| |
| | d | d | d | d ƒ } t t | d ƒ ƒ } t i | f | ƒ } | d d |
| |
| | d | d | | | d ƒ } nt i } t i |
f | ƒ } t i |
f | ƒ } d } t i | f | ƒ } | d d |
| |
| | | d | d | d d ƒ } t | d ƒ } t i | f | ƒ } | d d |
| |
| | | d | d | | d ƒ } t i i t i | d ƒ ƒ o
| } n | d | } | d d j o
t d ‚ n | S(
Ni i t Niÿÿÿÿf0.0x0.01.0R& s Eigenvalues did not converge( i ( R# R R% R R R R t real_tR! R" R; R R0 t dummyR R, t zgeevR. t wt rworkt lworkt workR2 t intt abst dgeevt wrt wit logical_andt reducet equalR (
R RB R2 R. RG RF RM RL RC RH R; R RE ( ( R t eigenvaluesƒ sH
"
R8 c C s^ t | ƒ t | ƒ t | ƒ } t d t | } t | | ƒ } | i d } d | d } t i
| f d ƒ } t | d j ot i } t i
| f | ƒ }
d } t i
| f | ƒ }
d } t i
| f | ƒ } | d | | | | |
|
d | d | | d ƒ
} t t |
d ƒ ƒ } t i
| f | ƒ }
t | d ƒ } t i
| f | ƒ } | d | | | | |
|
| | | | | d ƒ
} n³ t i } t i
| f | ƒ }
d } t i
| f | ƒ }
| d | | | | |
|
d | | d ƒ } t |
d ƒ } t i
| f | ƒ }
| d | | | | |
|
| | | d ƒ } | d d j o
t d ‚ n |
S(
Ni i i R i RA iÿÿÿÿR& s Eigenvalues did not converge( R# R R% R R R R RB R R" R; t liworkR R0 t iworkR R, t zheevdR. RE RG RH t lrworkRF t UPLOR2 RI RJ t dsyevdR ( R RV RU RB R2 R. RG RF RS RR RH R; R RE ( ( R t Heigenvalues¬ s@
04 **
c C s“ t | ƒ t | ƒ t | ƒ } t d t | } t | | ƒ } | i d } t i d | ƒ }
t | d j oÿ t i } t i | f | ƒ } t i | | f | ƒ } d } t i | f | ƒ } t i d | f | ƒ } | d d | | | | |
d | | | d | d ƒ } t t | d ƒ ƒ } t i | f | ƒ } | d d | | | | |
d | | | | | d ƒ } nüt i } t i | f | ƒ } t i | f | ƒ } t i | | f | ƒ }
d } t i | f | ƒ } | d d | | | | | |
d |
| | d d ƒ } t | d ƒ } t i | f | ƒ } | d d | | | | | |
d |
| | | d ƒ } t i i t i | d ƒ ƒ o | } |
} nÜ | d | } t i! |
t i" ƒ } t i# t i t i | d ƒ d ƒ ƒ } xŽ t% t& | ƒ d ƒ D]v } |
| d | d |
| d | d | | d | <|
| d | d |
| d | d | | d | d <qñW| d d j o
t( d
‚ n | | f S(
s½ eigenvectors(a) returns u,v where u is the eigenvalues and
v is a matrix of eigenvectors with vector v[i] corresponds to
eigenvalue u[i]. Satisfies the equation dot(a, v[i]) = u[i]*v[i]
i i i RA t Viÿÿÿÿf0.0x0.01.0R& s Eigenvalues did not convergeN( i () R# R R% R R R R RB R! R" R; R R0 RC R R, RD R. RE t vRG RH RF R2 RI RJ RK RL RM t vrRN RO RP t arrayt Complext nonzerot indt rangeR R R ( R RB R2 R. R_ RG RF RZ RM RL RC R RH R[ R; R RE ( ( R t eigenvectorsÑ s^
"
6>
c C sd t | ƒ t | ƒ t | ƒ } t d t | } t | | ƒ } | i d } d | d } t i
| f d ƒ } t | d j ot i } t i
| f | ƒ }
d } t i
| f | ƒ }
d } t i
| f | ƒ } | d | | | | |
|
d | d | | d ƒ
} t t |
d ƒ ƒ } t i
| f | ƒ }
t | d ƒ } t i
| f | ƒ } | d | | | | |
|
| | | | | d ƒ
} n³ t i } t i
| f | ƒ }
d } t i
| f | ƒ }
| d | | | | |
|
d | | d ƒ } t |
d ƒ } t i
| f | ƒ }
| d | | | | |
|
| | | d ƒ } | d d j o
t d ‚ n |
| f S(
Ni i i R i RY iÿÿÿÿR& s Eigenvalues did not converge( R# R R% R R R R RB R R" R; RR R R0 RS R R, RT R. RE RG RH RU RF RV R2 RI RJ RW R ( R RV RU RB R2 R. RG RF RS RR RH R; R RE ( ( R t
Heigenvectors
s@
04 **
c C sþ t | ƒ | i d } | i d }
t | ƒ } t d t | } t
| | ƒ } | o |
} | }
d } n% t | |
ƒ } t | |
ƒ }
d } t i t | |
ƒ f | ƒ } t i | |
f | ƒ } t i | |
f | ƒ } t i d t |
| ƒ f d ƒ } t | d j oÿ t i } t i d t |
| ƒ t |
| ƒ d t |
| ƒ f | ƒ } d } t i | f | ƒ } | | |
| | |
| | |
| |
| d | | d ƒ } t t | d ƒ ƒ } t i | f | ƒ } | | |
| | |
| | |
| |
| | | | d ƒ } n° t i } d } t i | f | ƒ } | | |
| | |
| | |
| |
| d | d ƒ } t | d ƒ } t i | f | ƒ } | | |
| | |
| | |
| |
| | | d ƒ } | d d j o
t! d
‚ n t" i# | ƒ | t" i# | ƒ f S( Ni i t At Si R i iÿÿÿÿR& s SVD did not converge($ R# R R" R; R: R R R R RB R! t
full_matricest nut nvtt optionR$ R R0 R t ut vtRS R R, t zgesddR. RF RG RH R2 RI RJ t dgesddR R4 R ( R Re RB R2 R. Rj RG RF RS Rf Rg Rh RH R: R; R Ri R ( ( R t singular_value_decomposition0 sP
" @!! !!
f1e-10c C s t i | d d ƒ} | i ƒ t i d j o t i | ƒ } n t | d ƒ \ } } } | i
d } | i
d } | t i i | ƒ } xN t t | | ƒ ƒ D]7 } | | | j o d | | | | <q d | | <q Wt i t i | ƒ | d d … t i f t i | ƒ ƒ S( NR i R] i f1.0f0.0( R R\ R R t typecodest conjugateRm Ri R Rj R" R: R; t rcondt maximumRO t cutoffR` R$ R t dotR R) ( R Rp Rr R R: R; R Rj Ri ( ( R t generalized_inverse_ s
c C sê t | ƒ t | ƒ t | ƒ } t | | ƒ } | i d } t | d j o
t i
} n
t i } t
i | f d ƒ } | | | | | | d ƒ } t
i i t
i | t
i d | d ƒ ƒ ƒ d } d d | t
i i t
i | ƒ ƒ S( Ni i R i f1.0f2.0( R# R R% R R R! R" R; R R, t zgetrfR. t dgetrfR R0 R1 R2 t addRO t not_equalt
arrayranget signt multiplyt diagonal( R R1 R2 R; R. R Rz ( ( R t determinantq s
c C s2 t | i ƒ d j } | o | d d … t i f } n t | | ƒ | i d } | i d } | i d }
t | | ƒ } | | i d j o
t
d ‚ n t | | ƒ } t d t | } t i | |
f | ƒ } t i | ƒ | d | i d … d |
… f <t | | | ƒ \ } } t i t | | ƒ f | ƒ } t d t t i t t | | ƒ ƒ d ƒ ƒ d ƒ } t i d t | | ƒ | d t | | ƒ f d ƒ }
t | d j ott i! } d } t i | f | ƒ } t i | f | ƒ } | | | |
| | | | | | d | d | |
d ƒ } t t( | d ƒ ƒ } t i | f | ƒ } t i | | f | ƒ } t i | |
f | ƒ } t i+ | | |
| | | | | | d | d |
d ƒ } t | d ƒ } t i | f | ƒ } t i | f | ƒ } | | | |
| | | | | | d | | | |
d ƒ } n° t i+ } d } t i | f | ƒ } | | | |
| | | | | | d | d |
d ƒ } t | d ƒ } t i | f | ƒ } | | | |
| | | | | | d | | |
d ƒ } | d
d j o
t
d ‚ n t i- g | ƒ } | on t i t i/ | ƒ | ƒ } | d | j o= | | j o0 t i- t i1 t i/ | ƒ | d
ƒ g ƒ } qnŒ t i t i2 | ƒ d | … d d … f ƒ } | d | j oL | | j o? t i t i1 t i2 | ƒ | d … d d … f d
ƒ ƒ } n | | | d t i | t | | ƒ ƒ f S( s solveLinearLeastSquares(a,b) returns x,resids,rank,s
where x minimizes 2-norm(|b - Ax|)
resids is the sum square residuals
rank is the rank of A
s is an rank of the singual values of A in desending order
If b is a matrix then x is also a matrix with corresponding columns.
If the rank of A is less than the number of columns of A or greater than
the numer of rows, then residuals will be returned as an empty array
otherwise resids = sum((b-dot(A,x)**2).
Singular values less than s[0]*rcond are treated as zero.
i Ni s Incompatible dimensionsf2.0i i R iÿÿÿÿR& s, SVD did not converge in Linear Least Squarest ranki (3 R R' R" R( R R) R# R R: R; R+ R t ldbR R R R R RB R0 t bstarR R R$ R RI t matht logt floatt nlvlRS R R, t zgelsdR. RG RF RH Rp R2 RJ t a_realt
bstar_realt dgelsdRU R\ t residsR3 t xt sumR ( R R' Rp R RU R( R2 R. R‰ R‡ R+ RG RF RS R€ RB R„ R† RH R: R; R R RŠ ( ( R t linear_least_squaresƒ sr
, 57 !
4. ?t __main__( t *c C s‹ d GHt | | ƒ } | t | | ƒ } | GHt | ƒ } t | | ƒ t | i d ƒ } | GHt
| ƒ } t | ƒ \ } } | | } | GHt | ƒ } t | | ƒ | | } | GHt | ƒ \ } }
} | t i | |
| ƒ } | GHt | ƒ } t | | ƒ t | i d ƒ } | GHt | ƒ }
|
t i | ƒ } | GHt | | ƒ \ } } } } | t | | ƒ } | GH| | i d GH| |
GHd S( Ns, All numbers printed should be (almost) zero:i ( R5 R R' RŠ t matrixmultiplyt checkR7 t a_invR6 R" RQ t evRa t evaluest evectorsR Rm Ri R Rj R Rt t a_ginvR} t detR{ RO RŒ t residualsR~ t sv( R R' R~ Rj R‘ R’ R R• RŠ R” R– R˜ R“ R Ri R— ( ( R t testÏ s8
f1.0f2.0f3.0f4.0x0.00.0(# R R R, R R> R4 t ExceptionR t _lapack_typet _lapack_letterR R R R R R! R R# R% R5 R7 R@ RQ RX Ra Rb Rm Rt R} RŒ R R™ R\ R R' ( R R} R R R5 R7 R Rm Rœ RŒ R™ R R R› R> R# R4 Rb Rt R' RX R R Ra R, R! R% R@ R R R RQ ( ( R t ? sJ *<< ) % 9 & / I
-