|
Server : Apache/2.2.2 (Fedora) System : Linux App1.pathumtani.go.th 2.6.20-1.2320.fc5smp #1 SMP Tue Jun 12 19:40:16 EDT 2007 i686 User : apache ( 48) PHP Version : 5.2.9 Disable Function : NONE Directory : /proc/self/root/usr/lib/python2.4/site-packages/Numeric/MA/ |
Upload File : |
mò
C1?c @ sã d Z d k Z d k Z d k Z d k Z d k Td k l Z l Z l Z e Z
d Z d e f d „ ƒ YZ
d f d „ ƒ YZ e d ƒ Z e i d
g ƒ i e ƒ Z e i d
d g ƒ i e ƒ Z d Z e i d
g ƒ i e ƒ Z d Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d
d
d „ Z! d „ Z" e" Z# d „ Z$ e% d „ Z& d „ Z' d „ Z( d f d „ ƒ YZ) d f d „ ƒ YZ* d f d „ ƒ YZ+ d! f d" „ ƒ YZ, d# f d$ „ ƒ YZ- d% f d& „ ƒ YZ. d' f d( „ ƒ YZ/ d) f d* „ ƒ YZ0 e- e i1 d+ e, d+ ƒ ƒ Z1 e- e i2 d, e+ d+ ƒ ƒ Z2 e- e i3 d, e+ d+ ƒ ƒ Z3 e- e i4 ƒ Z4 e- e i5 ƒ Z5 e- e i6 ƒ Z6 e- e i7 ƒ Z7 e- e i8 d+ e* d ƒ ƒ Z8 e- e i9 d+ e) d- d, ƒ ƒ Z9 e- e i: d+ e) d- d, ƒ ƒ Z: e- e i; ƒ Z; e- e i< ƒ Z< e- e i= ƒ Z= e- e i> ƒ Z> e- e i? ƒ Z? e- e i@ ƒ Z@ e- e iA ƒ ZA e- e iB ƒ ZB e- e iC ƒ ZC e- e iD ƒ ZD e- e iE ƒ ZE e- e iF ƒ ZF e- e iG d. ƒ ZG e- e iH ƒ ZH e0 e iI ƒ ZI e0 e iJ ƒ ZJ e% eJ _K e0 e iL d. d. ƒ ZL e/ e iM e. ƒ d
d. ƒ ZM e/ e iN e. ƒ d
d. ƒ ZN e/ e iO e. ƒ d
d. ƒ ZO e/ e iP e. ƒ d
d. ƒ ZP e/ e iQ e. ƒ d
d. ƒ ZQ e0 e iR ƒ ZR e0 e iS d+ d, ƒ ZS e% eS _K e0 e iT ƒ ZT e% eT _K e0 e iU ƒ ZU e% eU _K e0 e iV ƒ ZV e% eV _K e0 e iW ƒ ZW e% eW _K e0 e iX ƒ ZX e% eX _K e0 e iY ƒ ZY e% eY _K e0 e iZ ƒ ZZ e0 e i[ ƒ Z[ e0 e i\ ƒ Z\ e0 e i] ƒ Z] e0 e i^ ƒ Z^ e0 e i_ ƒ Z_ e i` Z` e ia Za e ib Zb d/ ec f d0 „ ƒ YZd ed Z d1 ed f d2 „ ƒ YZe ee ƒ Zf d3 „ Zg eg Zh eg Zi d. d4 d5 d6 „ Zj d. d7 „ Zk d4 d5 d. d
d8 „ Zl d. d
d9 „ Zm e% d. e% d: „ Zn en Zo d; „ Zp d< „ Zq d= „ Zr d> „ Zs d
d? „ Zt d@ „ Zu e% dA „ Zv ew d
dB „ Zx ew d
dC „ Zy e% dD „ Zz e% dE „ Z{ e% e% dF „ Z| eI iK Z} eL iK Z~ d
e% d
dG „ Z dH „ Z€ dI „ Z d. dJ „ Z‚ d. dK „ Zƒ d. dL „ Z„ d. dM „ Z… d. dN „ Z† d. dO „ Z‡ d. dP „ Zˆ d. dQ „ Z‰ d. dR „ ZŠ dS „ Z‹ dT „ ZŒ d
dU „ Z d
dV „ ZŽ e% dW „ Z dX „ Z dY „ Z‘ dZ „ Z’ d[ „ Z“ d\ „ Z” d] d^ „ Z• d_ f d` „ ƒ YZ– e– ƒ Z— da f db „ ƒ YZ˜ e˜ ƒ Z™ d] e% dc „ Zš d
d
d. dd „ Z› d] e% de „ Zœ d] e% df „ Z d] e% dg „ Zž dh „ ZŸ e% di „ Z dj a¡ e% e% e% dk „ Z¢ e i£ Z¤ e i¥ Z¦ d
dl „ Z§ dm „ Z¨ e§ t¡ ƒ d S(n sþ MA: a facility for dealing with missing observations
MA is generally used as a Numeric.array look-alike.
There are some differences in semantics, see manual.
In particular note that slices are copies, not references.
by Paul F. Dubois
L-264
Lawrence Livermore National Laboratory
dubois@users.sourceforge.net
Copyright 1999, 2000, 2001 Regents of the University of California.
Released for unlimited redistribution; see file Legal.htm
Documentation is in the Numeric manual; see numpy.sourceforge.net
N( t *( t es pis NewAxisf1e-35t MAErrorc B s# t Z d d „ Z d „ Z e Z RS( Nc C s
| | _ d S( s Create an exceptionN( t argst self( R R ( ( t1 /usr/lib/python2.4/site-packages/Numeric/MA/MA.pyt __init__ s c C s
t | i ƒ S( s# Calculate the string representationN( t strR R ( R ( ( R t __str__ s ( t __name__t
__module__t NoneR R t __repr__( ( ( R R s t _MaskedPrintOptionc B sG t Z d Z d „ Z d „ Z d „ Z d „ Z d d „ Z d „ Z RS( s< One instance of this class, masked_print_option, is created.c C s | i | ƒ d | _ d S( s& Create the masked print option object.i N( R t set_displayt displayt _enabled( R R ( ( R R s
c C s | i S( s# Show what prints for masked values.N( R t _display( R ( ( R R $ s c C s
| | _ d S( s2 set_display(s) sets what prints for masked values.N( t sR R ( R R ( ( R R ( s c C s | i S( s( Is the use of the display value enabled?N( R R ( R ( ( R t enabled, s i c C s
| | _ d S( s Set the enabling flag to flag.N( t flagR R ( R R ( ( R t enable0 s c C s
t | i ƒ S( N( R R R ( R ( ( R R 4 s ( R R
t __doc__R R R R R R ( ( ( R R
s s --f1e+20x0.00.0t ?i c C sD t | t i ƒ o t Sn&t | t i ƒ p t | t i ƒ o t Snø t | t i ƒ o t SnÝ t | t i
ƒ o t Sn t | t ƒ p t | t
i ƒ oš | i ƒ } | t d j o t Sn | t d j o t Sn | t d j o t Sn | t d j o t Sn | t d j o t
i t ƒ Sn t Sn t Sd S( s7 Function to calculate default fill value for an object.t Floatt Integert Complext Charactert UnsignedIntegerN( t
isinstancet objt typest FloatTypet default_real_fill_valuet IntTypet LongTypet default_integer_fill_valuet
StringTypet default_character_fill_valuet ComplexTypet default_complex_fill_valuet MaskedArrayt Numerict arraytypet typecodet xt typecodest absolutet default_object_fill_value( R R- ( ( R t default_fill_valueA s. &# c C sÝ t | t i ƒ o t Sn¿ t | t i ƒ p t | t i ƒ o t i SnŽ t | t ƒ p t | t
i ƒ oa | i ƒ } | t d j o t Sn | t d j o t i Sn | t d j o t i SqÙ n
t d ‚ d S( sD Function to calculate default fill value suitable for taking minima.R R R s( Unsuitable type for calculating minimum.N( R R R R R! R" R# t syst maxintR) R* R+ R, R- R. t TypeError( R R- ( ( R t minimum_fill_value[ s & # c C sÛ t | t i ƒ o t Sn¼ t | t i ƒ p t | t i ƒ o t Sn t | t ƒ p t | t i
ƒ o` | i ƒ } | t
d j o t Sn | t
d j o t i Sn | t
d j o d Sq× n
t d ‚ d S( sD Function to calculate default fill value suitable for taking maxima.R R R i s( Unsuitable type for calculating maximum.N( R R R R R! R" R# R$ R) R* R+ R, R- R. R2 R3 R4 ( R R- ( ( R t maximum_fill_valuel s & # c C s" t | ƒ o | i | ƒ n d S( s, Set fill value of a if it is a masked array.N( t
isMaskedArrayt at set_fill_valuet
fill_value( R8 R: ( ( R R9 } s
c C s& t | t ƒ o | i ƒ Sn d Sd S( s… Mask of values in a; could be None.
Returns None if a is not a masked array.
To get an array for sure use getmaskarray.N( R R8 R) t raw_maskR ( R8 ( ( R t getmask‚ s c C s5 t | ƒ } | d j o t t | ƒ ƒ Sn | Sd S( sË Mask of values in a; an array of zeros if mask is None
or not a masked array. Caution: has savespace attribute,
and is a byte-sized integer.
Do not try to add up entries, for example.
N( R< R8 t mR t make_mask_nonet shape( R8 R= ( ( R t getmaskarray‹ s
c C sC | d j p& t | t i ƒ o | i ƒ t j o d Sn d Sd S( s; Is m a legal mask? Does not check contents, only type.
i i N( R= R R R* t ArrayTypeR, t MaskType( R= ( ( R t is_mask— s 3c C sì | d j o d Sn§ t | t i ƒ or | i ƒ t j o? | o t i | d d ƒ} q– t i | d d d d ƒ} q» | i
t ƒ } | i d ƒ n" t i t | d ƒ t d d ƒ} | o" t i t i | ƒ ƒ o d Sn | Sd S( s make_mask(m, copy=0, flag=0)
return m as a mask, creating a copy if necessary or requested.
Can accept any sequence of integers or None. Does not check
that contents must be 0s and 1s.
if flag, return None if m contains no true elements.
t savespacei t copyi N( R= R R R* RA R, RB RE t arrayt resultt astypeRD t filledR t sometruet ravel( R= RE R RG ( ( R t make_mask s
!!c C s, t i | t ƒ } | i d ƒ | | _ | S( s& Return a mask of all zeros of shape s.i N( R* t zerosR RB RG RD R? ( R RG ( ( R R> º s
c C sn | d j o t | ƒ Sn | d j o t | ƒ Sn | | j o t | ƒ o | Sn t t i | | ƒ ƒ S( s Logical or of the mask candidates m1 and m2, treating None as false.
Result may equal m1 or m2 if the other is None.
N( t m1R RL t m2RC R* t
logical_or( RN RO ( ( R t mask_orà s
c C s t | t ƒ o | i | ƒ Sn] t | t i ƒ o | i ƒ o | Sn5 t | t i ƒ o t i
| d ƒ Sn t i
| ƒ Sd S( s“ a as a contiguous Numeric array with any masked areas replaced by value
if value is None or the special element "masked", fill_value(a)
is used instead.
If a is already a contiguous Numeric array, a itself is returned.
filled(a) can be used to be sure that the result is Numeric when
passing an object a to other software ignorant of MA, in particular to
Numeric itself.
t ON( R R8 R) RI t valueR* RA t iscontiguousR t DictTypeRF ( R8 RS ( ( R RI Ì s
c C s- t | ƒ o | i ƒ } n
t | ƒ } | S( sb
The fill value of a, if it has one; otherwise, the default fill value
for that type.
N( R7 R8 R: RG R1 ( R8 RG ( ( R R: à s
c C s1 t | ƒ } t | ƒ } | | j o | Sn d S( s: The common fill_value of a and b, if there is one, or NoneN( R: R8 t t1t bt t2R ( R8 RW RX RV ( ( R t common_fill_valueë s
t domain_check_intervalc B s t Z d Z d „ Z d „ Z RS( s9 domain_check_interval(a,b)(x) = true where x < a or y > bc C s | | _ | | _ d S( s9 domain_check_interval(a,b)(x) = true where x < a or y > bN( t y1R t y2( R R[ R\ ( ( R R õ s c C s. t i t i | | i ƒ t i | | i ƒ ƒ S( s Execute the call behavior.N( R* RP t greaterR- R R\ t lessR[ ( R R- ( ( R t __call__ú s ( R R
R R R_ ( ( ( R RZ ó s t
domain_tanc B s t Z d Z d „ Z d „ Z RS( s/ domain_tan(eps) = true where abs(cos(x)) < eps)c C s
| | _ d S( s/ domain_tan(eps) = true where abs(cos(x)) < eps)N( t epsR ( R Ra ( ( R R s c C s% t i t i t i | ƒ ƒ | i ƒ S( s Execute the call behavior.N( R* R^ R/ t cosR- R Ra ( R R- ( ( R R_ s ( R R
R R R_ ( ( ( R R` s t domain_greaterc B s t Z d Z d „ Z d „ Z RS( s( domain_greater(v)(x) = true where x <= vc C s
| | _ d S( s( domain_greater(v)(x) = true where x <= vN( t critical_valueR ( R Rd ( ( R R s c C s t i | | i ƒ S( s Execute the call behavior.N( R* t
less_equalR- R Rd ( R R- ( ( R R_ s ( R R
R R R_ ( ( ( R Rc
s t domain_greater_equalc B s t Z d Z d „ Z d „ Z RS( s- domain_greater_equal(v)(x) = true where x < vc C s
| | _ d S( s- domain_greater_equal(v)(x) = true where x < vN( Rd R ( R Rd ( ( R R s c C s t i | | i ƒ S( s Execute the call behavior.N( R* R^ R- R Rd ( R R- ( ( R R_ s ( R R
R R R_ ( ( ( R Rf s t masked_unary_operationc B s) t Z d d d „ Z d „ Z d „ Z RS( Ni c C s: | | _ | | _ | | _ t | d t | ƒ ƒ | _ d S( sÜ masked_unary_operation(aufunc, fill=0, domain=None)
aufunc(fill) must be defined
self(x) returns aufunc(x)
with masked values where domain(x) is true or getmask(x) is true.
R N( t aufuncR t ft fillt domaint getattrR R ( R Rh Rj Rk ( ( R R s
c O s
t | ƒ } t | | i ƒ } | i d j o t | | i | ƒ ƒ } n | d j oA | i
| | | Ž } t | ƒ t i j o t | ƒ Sq| Sno t | | ƒ } | i
t | | i ƒ | | Ž } t | ƒ t i j o t | | ƒ Sn | d o t Sn | Sd S( s Execute the call behavior..N( R< R8 R= RI R Rj t d1Rk R RQ Ri R t kwargsRG t typeR* RA t masked_arrayt dxt masked( R R8 R Rn R= RG Rq Rm ( ( R R_ * s"
! c C s d t | i ƒ S( Ns Masked version of ( R R Ri ( R ( ( R R A s ( R R
R R R_ R ( ( ( R Rg s t domain_safe_dividec B s t Z e d „ Z d „ Z RS( Nc C s
| | _ d S( N( t toleranceR ( R Rt ( ( R R F s c C s# t i | ƒ | i t i | ƒ j S( N( R* R/ R8 R Rt RW ( R R8 RW ( ( R R_ H s ( R R
t divide_toleranceR R_ ( ( ( R Rs E s t domained_binary_operationc B s/ t Z d Z d d d „ Z d „ Z d „ Z RS( sœ Binary operations that have a domain, like divide. These are complicated so they
are a separate class. They have no reduce, outer or accumulate.
i c C sC | | _ | | _ | | _ | | _ t | d t | ƒ ƒ | _ d S( sj abfunc(fillx, filly) must be defined.
abfunc(x, filly) = x for all x to enable reduce.
R N( t abfuncR Ri Rk t fillxt fillyRl R R ( R Rw Rk Rx Ry ( ( R R O s c
C sf t | ƒ } t | ƒ } t | | i ƒ } t | | i ƒ } | i | | ƒ } t
i | ƒ o( t | | i | ƒ } t | | ƒ } n t | | ƒ } | d j o> | i | | ƒ } t | ƒ t
i j o t | ƒ Sqà | Sn | i | | ƒ } t | ƒ t
i j oC | i | i j o t t | ƒ t | ƒ ƒ } n t | | ƒ Sn | d o t Sn | Sd S( s Execute the call behavior..N( R< R8 t maRW t mbRI R Rx Rm Ry t d2Rk t tR* RJ t whereRQ R= R Ri RG Ro RA Rp R? R@ Rr (
R R8 RW Rz RG R{ R= R} R| Rm ( ( R R_ Y s.
c C s d t | i ƒ S( Ns Masked version of ( R R Ri ( R ( ( R R t s ( R R
R R R_ R ( ( ( R Rv K s
t masked_binary_operationc B sJ t Z d d d „ Z d „ Z d d „ Z d „ Z d d „ Z d „ Z RS( Ni c C s: | | _ | | _ | | _ t | d t | ƒ ƒ | _ d S( sj abfunc(fillx, filly) must be defined.
abfunc(x, filly) = x for all x to enable reduce.
R N( Rw R Ri Rx Ry Rl R R ( R Rw Rx Ry ( ( R R x s
c O s@ t t | ƒ t | ƒ ƒ } | d j oh t | | i ƒ } t | | i
ƒ } | i | | | | Ž } t | ƒ t i j o t | ƒ Sq | Sn t | | i ƒ } t | | i
ƒ } | i | | | | Ž } t | ƒ t i j oC | i | i j o t t | ƒ t | ƒ ƒ } n t | | ƒ Sn | d o t Sn | Sd S( s Execute the call behavior..N( RQ R< R8 RW R= R RI R Rx Rm Ry R| Ri R Rn RG Ro R* RA Rp R? R@ Rr ( R R8 RW R Rn R= RG R| Rm ( ( R R_ s&
c C s t | ƒ } t | | i ƒ } | i f j o9 d | _ | d j o t | d d ƒ} d | _ qg n | d j o t
| i i | | ƒ ƒ Sn† t
| | ƒ } | i i t | | i ƒ | ƒ } t i i | | ƒ } t | t i ƒ o t
| | t | ƒ ƒ Sn | o t Sn | Sd S( s6 Reduce target along the given axis with this function.i RE N( i ( i ( R< t targetR= RI R Ry R} R? R RL Rp Ri t reducet axisR* t logical_andR RA R: Rr ( R R€ R‚ R= R} ( ( R R ˜ s$
!c C s t | ƒ } t | ƒ } | d j o | d j o
d } n+ t | ƒ } t | ƒ } t i | | ƒ } | i i t | | i
ƒ t | | i ƒ ƒ } t | | ƒ S( s< Return the function applied to the outer product of a and b.N( R< R8 Rz RW R{ R R= R@ RP t outerR Ri RI Rx Ry t dRp ( R R8 RW Rz R… R{ R= ( ( R R„ ® s
-c C s+ t | | i ƒ } t | i i | | ƒ ƒ S( s= Accumulate target along axis after filling with y fill value.N( RI R€ R Ry R} Rp Ri t
accumulateR‚ ( R R€ R‚ R} ( ( R R† » s c C s d t | i ƒ S( Ns Masked version of ( R R Ri ( R ( ( R R ¿ s ( R R
R R_ R R„ R† R ( ( ( R R w s
f0.0f1.0f-1.0i R) c B sh t Z d Z d Z e d e e e d „ Z e d „ Z d „ Z d „ Z d „ Z d „ Z
d „ Z d
„ Z d „ Z
d „ Z d
„ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z d „ Z e Z e Z e Z d „ Z d „ Z d „ Z d „ Z! e! Z" d „ Z# d „ Z$ d „ Z% d „ Z& d! „ Z' d" „ Z( d# „ Z) e) Z* d$ „ Z+ d% „ Z, d& „ Z- d' „ Z. d( „ Z/ d) „ Z0 e d* „ Z1 d+ „ Z2 d, „ Z3 d- „ Z4 d. „ Z5 d/ „ Z6 d0 „ Z7 d1 „ Z8 d2 „ Z9 d3 „ Z: d4 „ Z; d5 „ Z< d6 „ Z= d7 „ Z> d8 „ Z? e d9 „ Z@ d: „ ZA d; „ ZB e d< „ ZC d= „ ZD d> „ ZE d? „ ZF d@ „ ZG dA „ ZH dB „ ZI dC „ ZJ dD „ ZK dE „ ZL dF „ ZM dG „ ZN e dH „ ZO e dI „ ZP dJ „ ZM dK „ ZQ e dL „ ZR e dM „ ZS dN „ ZT dO „ ZU eV e e dP dQ ƒZW eV e e
dP dR ƒZX eV e e dP dS ƒZY eV e
e dP dT ƒZZ eZ Z[ RS(U s€ Arrays with possibly masked values.
Masked values of 1 exclude element from the computation.
Construction:
x = array(data, typecode=None, copy=1, savespace=0,
mask = None, fill_value=None)
If copy=0, every effort is made not to copy the data:
If data is a MaskedArray, and argument mask=None,
then the candidate data is data.raw_data() and the
mask used is data.mask(). If data is a Numeric array,
it is used as the candidate raw data.
If savespace != data.spacesaver() or typecode is not None and
is != data.typecode() then a data copy is required.
Otherwise, the candidate is used.
If a data copy is required, raw data stored is the result of:
Numeric.array(data, typecode=typecode, copy=copy, savespace=savespace)
If mask is None there are no masked values. Otherwise mask must
be convertible to an array of integers of typecode MaskType,
with values 1 or 0, and of the same shape as x.
fill_value is used to fill in masked values when necessary,
such as when printing and in method/function filled().
The fill_value is not used for computation within this module.
If savespace is 1, the data is given the spacesaver property, and
the mask is replaced by None if all its elements are true.
s MA.MaskedArrayi c C s€ | }
| }
| } t | t ƒ o¾ | i ƒ } | i ƒ } |
d j o
| }
n |
| j o
d } n | i
ƒ } |
d j o
| }
n |
| j o
d } n d }
| d j o | i ƒ } q¦| d j o
d } q¦nÇ t | t i ƒ o~ | } | i ƒ } |
d j o
| }
n |
| j o
d } n | i
ƒ } |
d j o
| }
q¦|
| j o
d } q¦n6 d } |
d j o
d }
n t i | |
d |
ƒ} | oU |
| j o" t i | d d d |
ƒ| _ q | i |
ƒ | _ | i i |
ƒ n
| | _ | d j o d | _ d | _ nBt | d |
ƒ| _ | i d j o
d | _ n| i | j | _ t | i ƒ } t | i ƒ } | | j o† | d j o+ t i | i | i i ƒ | _ d | _ qo| d j o4 t i | i | i i ƒ | _ | i i | i _ qot d ‚ nM | d j o? t | i ƒ t | i ƒ j o | i ƒ | i i | i _ n | i | ƒ d S( sõ array(data, typecode=None,copy=1, savespace=None,
mask=None, fill_value=None)
If data already a Numeric array, its typecode and spacesaver()
become the default values for typecode and savespace.
i i RD RE R s Mask and data not compatible.N(" R, t tcRD t ssRE t need_data_copiedR t dataR) t raw_datat ct ctcR t
spacesavert csst maskR* RA RF R t _dataRH t _maskt _shared_maskRL t sizet nmt ndt resizeR? R t unshare_maskR9 R: ( R RŠ R, RE RD R R: R‰ R• R– R‡ R RŒ Rˆ R ( ( R R ! s€
"
,
c C su | i d j o? t i t i | i ƒ ƒ o
t d ‚ qO d | _ d | _ n | o | i i
| ƒ Sn | i Sd S( s: Special hook for Numeric. Converts to Numeric if possible.sŠ Cannot automatically convert masked array to Numeric because data
is masked in one or more locations.
i N( R R’ R R* RJ RK R R“ R} R‘ RH ( R R} ( ( R t __array__t s
c C s
| i i S( s Return the current shape.N( R R‘ R? ( R ( ( R t
_get_shape… s c C s° | i i ƒ p4 t i | i | i i ƒ d | i i ƒ ƒ | _ n | | i _ | i d j oL | i ƒ | i i ƒ p" t i | i t d d ƒ | _ n | | i _ n d S( s Set the array's shape.i N(
R R‘ RT R* RF R, RŽ t newshapeR? R’ R R˜ RB ( R R› ( ( R t
_set_shape‰ s
"c C sb | i d j o& t | i i d d d | i ƒ ƒSn) t | i i d | i i d | i ƒ ƒSd S( s# Calculate the flat value.
R R: N( R R’ R Rp R‘ t flatR: ( R ( ( R t _get_flat• s c C s | i } | | (d S( s x.flat = valueN( R R t yRS ( R RS RŸ ( ( R t _set_flat s c C se | i d j o& t | i i d d d | i ƒ ƒSn, t | i i d | i ƒ i d | i ƒ ƒSd S( s% Get the real part of a complex array.R R: N( R R’ R Rp R‘ t realR: R R ( R ( ( R t _get_real¥ s c C s | i } | | d <d S( s x.real = value.N( R R¡ RŸ RS ( R RS RŸ ( ( R t _set_real® s c C se | i d j o& t | i i d d d | i ƒ ƒSn, t | i i d | i ƒ i d | i ƒ ƒSd S( s* Get the imaginary part of a complex array.R R: N( R R’ R Rp R‘ t imaginaryR: R R ( R ( ( R t _get_imaginary³ s c C s | i } | | d <d S( s x.imaginary = value.N( R R¤ RŸ RS ( R RS RŸ ( ( R t _set_imaginary¼ s c C s6 t i ƒ o
t } n
| i ƒ } t t | | ƒ ƒ S( s~ Calculate the str representation, using masked for fill if
it is enabled. Otherwise fill with fill value.
N( t masked_print_optionR Ri R R: R RI ( R Ri ( ( R R Á s
c C s§ d } d } d } d } t | i ƒ } | i d j oO | d j o! | h d t
| i ƒ ƒ <Sn | h d t
| i ƒ ƒ <Sn | d j oK | h d t
| i ƒ ƒ <d t
| i ƒ ƒ <d t
| i
ƒ ƒ <Sn | h d t
| i ƒ ƒ <d t
| i ƒ ƒ <d t
| i
ƒ ƒ <Sd } | i d j o! | h d t
| i ƒ ƒ <SnH | h d t
| i ƒ ƒ <d t
| i ƒ ƒ <d t
| i
ƒ ƒ <Sd S(
s Calculate the repr representation, using masked for fill if
it is enabled. Otherwise fill with fill value.
sM array(data =
%(data)s,
mask =
%(mask)s,
fill_value=%(fill)s)
sI array(data = %(data)s,
mask = %(mask)s,
fill_value=%(fill)s)
s array(
%(data)s)s array(%(data)s)i RŠ R Rj N( t with_maskt
with_mask1t without_maskt
without_mask1t lenR R? t nR’ R R RI R R: ( R R R© R« Rª R¨ ( ( R R Ë s"
!!
KG!c C s; | i ƒ | i d j o
t d ‚ n t | i ƒ d ƒ S( s Convert self to float.s0 Cannot convert masked element to a Python float..N( R t unmaskR’ R R t floatR‹ ( R ( ( R t __float__ú s
c C s; | i ƒ | i d j o
t d ‚ n t | i ƒ d ƒ S( s Convert self to int.s. Cannot convert masked element to a Python int..N( R R® R’ R R t intR‹ ( R ( ( R t __int__ s
c C sÞ | i } | i | } | i i ƒ } | i i ƒ } t
| ƒ t i j ok | d j o" t | d | d d d | ƒ} n5 t | d | d d d | d | | d | i ƒ ƒ} | Sn& | d j p | | o | Sn t Sd S( s Get copy of item described by i.R, RE i RD R R: N( R R’ R= R‘ t it doutRŽ Rˆ R, R‡ Ro R* RA R RF RG R: Rr ( R R³ R´ R= Rˆ RG R‡ ( ( R t __getitem__ s
" c C s | i } | i | | !} | i i ƒ } | i i ƒ } | d j o t | d | d d d | ƒSn6 t | d | d d d | d | | | !d | i
ƒ ƒSd S( s# Get copy of slice described by i, jR, RE i RD R R: N( R R’ R= R‘ R³ t jR´ RŽ Rˆ R, R‡ R RF R: ( R R³ R¶ R´ R= Rˆ R‡ ( ( R t __getslice__ s
c C s8 | t j o
t d ‚ n | t j oQ | i d j o" t | i i ƒ | _ d | _ n | i
ƒ d | i | <d Sn t | ƒ } t | ƒ i | i i ƒ ƒ } | | i | <| d j o/ | i d j o | i
ƒ d | i | <q4nJ | i d j o" t | i i ƒ | _ d | _ n | i
ƒ | | i | <d S( sF Set item described by index. If value is masked, mask those locations.s Cannot alter the masked element.i i N( R Rr R RS R’ R R> R‘ R? R“ R˜ t indexR< R= RI RH R, ( R R¸ RS R= ( ( R t __setitem__+ s,
c C s0 | t j o
t d ‚ n | t j oJ | i d j o" t | i i ƒ | _ d | _ n d | i | | +d Sn t | ƒ } t | ƒ i | i i ƒ ƒ } | | i | | +| d j o2 | i d j o | i ƒ d | i | | +q,nC | i d j o" t | i i ƒ | _ d | _ n | | i | | +d S( s8 Set slice i:j; if value is masked, mask those locations.s Cannot alter the masked element.i i N( R Rr R RS R’ R R> R‘ R? R“ R³ R¶ R< R= RI RH R, R˜ ( R R³ R¶ RS R= ( ( R t __setslice__F s(
c C s
t | i ƒ S( se Return length of first dimension. This is weird but Python's
slicing behavior depends on it.N( R¬ R R‘ ( R ( ( R t __len__] s c C s
t | | ƒ S( s Return bitwise_andN( t bitwise_andR t other( R R½ ( ( R t __and__b s c C s
t | | ƒ S( s Return bitwise_orN( t
bitwise_orR R½ ( R R½ ( ( R t __or__f s c C s
t | | ƒ S( s Return bitwise_xorN( t bitwise_xorR R½ ( R R½ ( ( R t __xor__j s c C s
t | ƒ S( s Return absolute(self)N( R/ R ( R ( ( R t __abs__r s c C s
t | ƒ S( s Return negative(self)N( t negativeR ( R ( ( R t __neg__v s c C s
t | ƒ S( s Return array(self)N( RF R ( R ( ( R t __pos__z s c C s
t | | ƒ S( s Return add(self, other)N( t addR R½ ( R R½ ( ( R t __add__~ s c C s
t | | ƒ S( s Return remainder(self, other)N( t remainderR R½ ( R R½ ( ( R t __mod__„ s c C s
t | | ƒ S( s Return remainder(other, self)N( RÉ R½ R ( R R½ ( ( R t __rmod__ˆ s c C s
t | | ƒ S( N( t
left_shiftR R ( R R ( ( R t
__lshift__Πs c C s
t | | ƒ S( N( t right_shiftR R ( R R ( ( R t
__rshift__ s c C s
t | | ƒ S( s Return subtract(self, other)N( t subtractR R½ ( R R½ ( ( R t __sub__’ s c C s
t | | ƒ S( s Return subtract(other, self)N( RÐ R½ R ( R R½ ( ( R t __rsub__– s c C s
t | | ƒ S( s Return multiply(self, other)N( t multiplyR R½ ( R R½ ( ( R t __mul__š s c C s
t | | ƒ S( s Return divide(self, other)N( t divideR R½ ( R R½ ( ( R t __div__ s c C s
t | | ƒ S( s Return divide(other, self)N( RÕ R½ R ( R R½ ( ( R t __rdiv__¤ s c C s
t | | ƒ S( s Return divide(self, other)N( t true_divideR R½ ( R R½ ( ( R t __truediv__¨ s c C s
t | | ƒ S( s Return divide(other, self)N( RØ R½ R ( R R½ ( ( R t __rtruediv__¬ s c C s
t | | ƒ S( s Return divide(self, other)N( t floor_divideR R½ ( R R½ ( ( R t __floordiv__° s c C s
t | | ƒ S( s Return divide(other, self)N( RÛ R½ R ( R R½ ( ( R t
__rfloordiv__´ s c C s t | | | ƒ S( s Return power(self, other, third)N( t powerR R½ t third( R R½ Rß ( ( R t __pow__¸ s c C s
t | ƒ S( s Return sqrt(self)N( t sqrtR ( R ( ( R t __sqrt__¼ s c C s | i i ƒ } t | d ƒ } | i ƒ } | | j o n<| t d j o1 | t d j o | i | ƒ } qvt
d ‚ nú | t d j oU | t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n” | t d j oy | t d j o | i | ƒ } qv| t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n
t
d ‚ | i d j o7 | i | 7_ t
| ƒ } | | _ | d j | _ nI t | t | d t
| ƒ ƒƒ } | i ƒ | _ | i ƒ | _ d | _ | S( s Add other to self in place.i R s& Incorrect type for in-place operation.R R R i N( R R‘ R, R} RI R½ Ri RV R. RH R4 R’ R R< R= R“ RÇ Rp RG R‹ R; ( R R½ Ri R= RV R} RG ( ( R t __iadd__À sF
! c C s | i i ƒ } t | d ƒ } | i ƒ } | | j o n<| t d j o1 | t d j o | i | ƒ } qvt
d ‚ nú | t d j oU | t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n” | t d j oy | t d j o | i | ƒ } qv| t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n
t
d ‚ | i d j o7 | i | 9_ t
| ƒ } | | _ | d j | _ nI t | t | d t
| ƒ ƒƒ } | i ƒ | _ | i ƒ | _ d | _ | S( s Add other to self in place.i R s& Incorrect type for in-place operation.R R R i N( R R‘ R, R} RI R½ Ri RV R. RH R4 R’ R R< R= R“ RÓ Rp RG R‹ R; ( R R½ Ri R= RV R} RG ( ( R t __imul__ë sF
! c C s | i i ƒ } t | d ƒ } | i ƒ } | | j o n<| t d j o1 | t d j o | i | ƒ } qvt
d ‚ nú | t d j oU | t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n” | t d j oy | t d j o | i | ƒ } qv| t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n
t
d ‚ | i d j o7 | i | 8_ t
| ƒ } | | _ | d j | _ nI t | t | d t
| ƒ ƒƒ } | i ƒ | _ | i ƒ | _ d | _ | S( s" Subtract other from self in place.i R s& Incorrect type for in-place operation.R R R i N( R R‘ R, R} RI R½ Ri RV R. RH R4 R’ R R< R= R“ RÐ Rp RG R‹ R; ( R R½ Ri R= RV R} RG ( ( R t __isub__ sF
! c C sâ | i i ƒ } t | d ƒ } | i ƒ } | | j o n<| t d j o1 | t d j o | i | ƒ } qvt
d ‚ nú | t d j oU | t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n” | t d j oy | t d j o | i | ƒ } qv| t d j o | i | ƒ } qv| t d j o | i | ƒ } qvt
d ‚ n
t
d ‚ t | ƒ } t
| t | d | ƒƒ } | i ƒ | _ | i ƒ } | | i j o | | _ d | _ n | S( s Divide self by other in place.i R s& Incorrect type for in-place operation.R R R i N( R R‘ R, R} RI R½ Ri RV R. RH R4 R< t moRÕ Rp RG R‹ R; t dmR’ R“ ( R R½ Rç Ri RV R} Ræ RG ( ( R t __idiv__C sB
c C s
t | | ƒ S( N( t equalR R½ ( R R½ ( ( R t __eq__j s c C s
t | | ƒ S( N( t not_equalR R½ ( R R½ ( ( R t __ne__m s c C s
t | | ƒ S( N( R^ R R½ ( R R½ ( ( R t __lt__p s c C s
t | | ƒ S( N( Re R R½ ( R R½ ( ( R t __le__s s c C s
t | | ƒ S( N( R] R R½ ( R R½ ( ( R t __gt__v s c C s
t | | ƒ S( N( t
greater_equalR R½ ( R R½ ( ( R t __ge__y s c C s; | i i | ƒ } | i | i i ƒ ƒ t | d | i ƒS( s# return self as array of given type.R N( R R‘ RH R‡ R… RD RŽ RF R’ ( R R‡ R… ( ( R RH | s c C s
| i i ƒ S( s‰ Returns the raw data field, byte_swapped. Included for consistency
with Numeric but doesn't make sense in this context.
N( R R‘ t byte_swapped( R ( ( R Rò ‚ s c C sl t i | i ƒ } | i d j o t | ƒ Sn9 d t i | i ƒ } t i | | ƒ } t | d d ƒSd S( s' A 1-D array of all the non-masked data.i RE i N( R* RK R R‘ R… R’ R RF R= t compressRŒ ( R R… R= RŒ ( ( R t
compressedˆ s c
C sI | i } | i i } t | ƒ } | d j o~ | d j o d Sn | d j o | d Sn | d j o t
d „ | ƒ Sq¬ | | } t | ƒ } | | =t | ƒ | Sn | d j o] t i | ƒ i t ƒ } t | ƒ } | d j o | d } n t i i
| ƒ } | | Sn0 t | | ƒ } t | i t ƒ | ƒ } | | Sd S( s? Count of the non-masked elements in a, or along a certain axis.i i c C s | | S( N( R- RŸ ( R- RŸ ( ( R t <lambda> s N( R R’ R= R‘ R? R R¬ t lsR R‚ R R t listR} t onesR* RK RH t Intt wR” t n1t n2RÇ t sum(
R R‚ R= R R Rö Rú R} Rû Rü ( ( R t count’ s2
c C s
t | | ƒ S( s% s.dot(other) = innerproduct(s, other)N( t innerproductR R½ ( R R½ ( ( R t dot° s c C s | i S( s Get the current fill value.N( R t _fill_value( R ( ( R R: ´ s c C s | i } | i } | d j oD | i ƒ o | Sqc t i | d | i ƒ d d d | i
ƒ ƒSn | } | d j o
| i
} n | t j o t i | ƒ } | i | _ ng yD t i | d | i ƒ d d d | i
ƒ ƒ} t i | | | ƒ Wn t i | | | f ƒ } n X| S( sð A Numeric array with masked values filled. If fill_value is None,
use self.fill_value().
If mask is None, copy data only if not contiguous.
Result is always a contiguous, Numeric array.
R, RE i RD N( R R‘ R… R’ R= R RT R* RF R, RŽ R: RS R Rr RG R? t putmaskt choose( R R: RG R… R= RS ( ( R RI ¸ s*
c C s t | i ƒ t | i ƒ f S( s) Return the ids of the data and mask areasN( t idR R‘ R’ ( R ( ( R t idsÖ s c C s
| i i ƒ S( s Is the data contiguous?N( R R‘ RT ( R ( ( R RT Ú s c C s
| i i ƒ S( s Item size of each data item.N( R R‘ t itemsize( R ( ( R R Þ s c C sG | i } | d j o | Sn& | i ƒ o | Sn t i | i ƒ Sd S( s1 Return the data mask, or None. Result contiguous.N( R R’ R= R RT R* RF ( R R= ( ( R R â s
c C s
t | | ƒ S( s' s.outer(other) = outerproduct(s, other)N( t outerproductR R½ ( R R½ ( ( R R„ ì s c C sŽ t i | i ƒ ƒ } | i d j o
| } n t i d | i | ƒ } t | ƒ t | ƒ j o
t d ‚ n t i | i
| t | ƒ ƒ d S( s\ Set the non-masked entries of self to filled(values).
No change to mask
i s( x.put(values) incorrect count of values.N( R* t arangeR R” t iotaR’ R t indRó R¬ t valuesR t putR‘ RI ( R R R
R ( ( R R ð s
c C s– | i d j o‚ t i | i ƒ ƒ } t i | i | ƒ } t | ƒ t | ƒ j o
t d ‚ n t i | i
| t | ƒ ƒ d | _ d | _ n d S( s\ Set the masked entries of self to filled(values).
Mask changed to None.
s( x.put(values) incorrect count of values.i N( R R’ R R* R R” R Ró R
R¬ R R R R‘ RI R“ ( R R R
R ( ( R R ý s
c C s | i S( s_ The raw data; portions may be meaningless.
May be noncontiguous. Expert use only.N( R R‘ ( R ( ( R R‹
s c C s | i S( sh The raw mask; portions may be meaningless.
May be noncontiguous. Expert use only.
N( R R’ ( R ( ( R R; s c C s
| i i ƒ S( s Get the spacesaver attribute.N( R R‘ RŽ ( R ( ( R RŽ s c C s | i i | ƒ d S( s% Set the spacesaver attribute to valueN( R R‘ RD RS ( R RS ( ( R RD s c C s0 | d j o t | i ƒ ƒ } n | | _ d S( s3 Set the fill value to v. Omit v to restore default.N( t vR R1 R R‹ R ( R R
( ( R R9 s
c C sT | i i } | d j o/ t | ƒ d j o d SqP t d „ | ƒ Sn | | Sd S( s5 Number of elements in array, or in a particular axis.i i c C s | | S( N( R- RŸ ( R- RŸ ( ( R Rõ * s N( R R‘ R? R R‚ R R¬ R ( R R‚ R ( ( R R” # s
c C s
| i i ƒ S( s. spacesaver() queries the spacesaver attribute.N( R R‘ RŽ ( R ( ( R RŽ . s c C s
| i i ƒ S( N( R R‘ R, ( R ( ( R R, 2 s c C s | i | ƒ i ƒ S( s Convert to listN( R RI R: t tolist( R R: ( ( R R 5 s c C s | i | ƒ i ƒ S( s Convert to stringN( R RI R: t tostring( R R: ( ( R R 9 s c C sT | i d j o d Sn t | i d d ƒ} | d j o d | _ d | _ n d S( s% Replace the mask by None if possible.NR i i ( R R’ R RL R= R“ ( R R= ( ( R R® = s
c C s9 | i o+ t | i d d d d ƒ| _ d | _ n d S( s' If currently sharing mask, make a copy.RE i R i N( R R“ RL R’ ( R ( ( R R˜ E s
t docs# tuple giving the shape of the arrays Access array in flat form.s! Access the real part of the arrays& Access the imaginary part of the array(\ R R
R t handler_cache_keyR R R™ Rš Rœ Rž R R¢ R£ R¥ R¦ R R R° R² Rµ R· R¹ Rº R» R¾ RÀ R t __rand__t __ror__t __rxor__Rà RÅ RÆ RÈ t __radd__RÊ RË RÍ RÏ RÑ RÒ RÔ t __rmul__RÖ R× RÙ RÚ RÜ RÝ Rà Râ Rã Rä Rå Rè Rê Rì Rí Rî Rï Rñ RH Rò Rô Rþ R R: RI R RT R R R„ R R R‹ R; RŽ RD R9 R” R, R R R® R˜ t propertyR? R R¡ R¤ t imag( ( ( R R) sº S
/ + + - '
t MaskedScalarc B s t Z d „ Z e e i ƒ Z RS( Nc C s8 t i | d g d d g ƒf | i _ f | i _ d S( Ni R i ( R) R R R‘ R? R’ ( R ( ( R R \ s ( R R
R R R) Rš R? ( ( ( R R [ s c C s
t | t ƒ S( s9 Is x a masked array, that is, an instance of MaskedArray?N( R R- R) ( R- ( ( R R7 e s f1.0000000000000001e-05f1e-08c C sÍ t t | ƒ t | ƒ ƒ } t | ƒ }
t | ƒ } t t |
d d d | ƒ| ƒ i
t ƒ } t t | d d d | ƒd ƒ i
t ƒ } t i t i | | ƒ | | t i | ƒ ƒ } t i t i | ƒ ƒ S( s¥ Returns true if all components of a and b are equal
subject to given tolerances.
If fill_value is 1, masked values considered equal.
If fill_value is 0, masked values considered unequal.
The relative error rtol should be positive and << 1.0
The absolute error atol comes into play for those elements
of b that are very small or zero; it says how small a must be also.
RE i R i N( RQ R< R8 RW R= RI Rm R| RF R: RH R R- RŸ R* Re R/ t atolt rtolR… t alltrueRK ( R8 RW R: R R R… R= RŸ R- R| Rm ( ( R t allclosel s **0c C sà t t | ƒ t | ƒ ƒ } | d j oD t | ƒ } t | ƒ } t i
| | ƒ } t i t i
| ƒ ƒ Snq | oe t | ƒ } t | ƒ } t i
| | ƒ } t | d | d d ƒ} t i t i
t | d ƒ ƒ ƒ Sn d Sd S( s„
True if all entries of a and b are equal, using
fill_value as a truth value where either or both are masked.
R RE i i N( RQ R< R8 RW R= R RI R- RŸ R* Ré R… R RK R: RF Rç ( R8 RW R: R… R= Rç RŸ R- ( ( R t allequal} s
#c
C s¨ t i } t | | ƒ } | i ƒ t d j ob t i | | | ƒ | | | | ƒ ƒ } t
| d d ƒ} t | d | d | d | d | ƒSn t | | | | ƒ Sd S( s´
masked_values(data, value, rtol=1.e-5, atol=1.e-8)
Create a masked array; mask is None if possible.
If copy==0, and otherwise possible, result
may share data values with original array.
Let d = filled(data, value). Returns d
masked where abs(data-value)<= atol + rtol * abs(value)
if d is of a floating point type. Otherwise returns
masked_object(d, value, copy, savespace)
R R i R RD RE R: N( R* R/ t absRI RŠ RS R… R, R. Re R R R= RL RF RD RE t
masked_object( RŠ RS R R RE RD R R= R… ( ( R t
masked_values‘ s *c
C sO t | | ƒ } t t i | | ƒ d d ƒ} t | d | d | d | d | ƒS( s5 Create array masked where exactly data equal to valueR i R RE RD R: N( RI RŠ RS R… RL R* Ré Rç RF RE RD ( RŠ RS RE RD Rç R… ( ( R R § s
c C s t t i | | | | ƒ ƒ S( sq Just like range() except it returns a array whose type can be specfied
by the keyword argument typecode.
N( RF R* t
arrayranget startt stopt stepR, ( R# R$ R% R, ( ( R R" ® s c C s t t i | | ƒ ƒ S( sA Construct a masked array from a string. Result will have no mask.N( Rp R* t
fromstringR R} ( R R} ( ( R R& ¶ s c C sk t | ƒ } | d j o& t i t | ƒ | ƒ } t | ƒ Sn) t i t | d ƒ | ƒ } t | | ƒ Sd S( s Left shift n bitsi N(
R< R8 R= R R* RÌ RI R R… Rp ( R8 R R… R= ( ( R RÌ º s
c C sk t | ƒ } | d j o& t i t | ƒ | ƒ } t | ƒ Sn) t i t | d ƒ | ƒ } t | | ƒ Sd S( s Right shift n bitsi N(
R< R8 R= R R* RÎ RI R R… Rp ( R8 R R… R= ( ( R RÎ Ä s
c C sj t | ƒ } | d j o t i | | ƒ } n t t i t | ƒ | ƒ d | ƒ} | i
t | ƒ ƒ | S( sw resize(a, new_shape) returns a new array with the specified shape.
The original array's total size can be any size.R N( R< R8 R= R R* R— t new_shapeRF RI RG R9 R: ( R8 R' R= RG ( ( R R— Î s
$c C s t | ƒ } t | t i ƒ o! t | g t | ƒ | ƒ } n t
| ƒ } | d j o t
i | | | ƒ } n t
i | | | ƒ } t | | ƒ } | i t | ƒ ƒ | S( s repeat elements of a repeats times along axis
repeats is a sequence of length a.shape[axis]
telling how many times to repeat each element.
N( RI R8 t afR t repeatsR R" t tupleR? R‚ R< R= R R* t repeatR… Rp RG R9 R: ( R8 R) R‚ R… RG R( R= ( ( R R+ Ø s !
c C s t t i | ƒ ƒ S( s< identity(n) returns the identity matrix of shape n x n.
N( RF R* t identityR ( R ( ( R R, é s c C s t t i | | ƒ ƒ S( s„ indices(dimensions,typecode=None) returns an array representing a grid
of indices with row-only, and column-only variation.
N( RF R* t indicest
dimensionsR, ( R. R, ( ( R R- î s c C s t t i | | | ƒ ƒ S( s_ zeros(n, typecode=Int, savespace=0) =
an array of all zeros of the given length or shape.N( RF R* RM R? R, RD ( R? R, RD ( ( R RM ô s c C s t t i | | | ƒ ƒ S( s] ones(n, typecode=Int, savespace=0) =
an array of all ones of the given length or shape.N( RF R* Rø R? R, RD ( R? R, RD ( ( R Rø ù s c C s t | ƒ } | i | ƒ S( s? Count of the non-masked elements in a, or along a certain axis.N( Rp R8 Rþ R‚ ( R8 R‚ ( ( R Rþ ÿ s c C s | d j o
t d ‚ n t | ƒ } t | ƒ } t | | ƒ } t
| d ƒ } t
| d ƒ } | i
ƒ t d j o t t i | | ƒ | ƒ Sn t t i | d ƒ d d ƒ} t | | ƒ } | d j o t t i | | ƒ ƒ Sn/ t i | d | ƒ } t t i | | ƒ | ƒ Sd S( s a**bs 3-argument power not supported.i R i R N( Rß R R R< R8 Rz RW R{ RQ R= RI t fat fbR, R. Rp R* RÞ RL Re t mdR~ ( R8 RW Rß R1 Rz R{ R= R/ R0 ( ( R RÞ s
c C sV | d j o- t | ƒ o | i ƒ } q: t | ƒ } n t | d | d d d | ƒS( s~ masked_array(a, mask=None) =
array(a, mask=mask, copy=0, fill_value=fill_value)
Use fill_value(a) if None.
R RE i R: N( R: R R7 R8 R1 RF R ( R8 R R: ( ( R Rp s
c
C s t | ƒ } | i ƒ } | i } | f j o
d
} n | d j o2| d j oˆ | d j o1 t i | i
ƒ i ƒ }
t d „ | d ƒ } qst | d ƒ i } t i i | i
ƒ i | ƒ }
t i i | ƒ } ~ qV| d j o@ t i | i ƒ }
t i | d ƒ i } t i i | ƒ } ~ qVt t | d ƒ t i d | ƒi } t i | i | ƒ }
t i | ƒ } ~ nà| d j oc| d j o- | | d } t i i | i
ƒ | ƒ }
qVt | d ƒ } | i } | f j o
d } n | | j oG t i | t d d ƒ} t i | | | ƒ }
t i | | ƒ } ~ qV| | | f j o… | | } t g t | ƒ } t d d d ƒ | | <t d t t | ƒ ƒ d ƒ } t i | | | ƒ }
t i | | ƒ } ~ ~ qVt d
‚ np| d j o@ t i | | ƒ }
t i | d ƒ } t i i | | ƒ } ~ n#t | d ƒ } | i } | f j o
d } n | | j oG t | t d | d d ƒ} t i | | | ƒ }
t i | | ƒ } n | | | f j o | | } t g t | ƒ } t d d d ƒ | | <t d t t | ƒ ƒ d ƒ } t i | | | ƒ }
t i | | ƒ } n
t d
‚ ~ t |
| ƒ } ~
t | t! ƒ op | i" ƒ | o[ t | t! ƒ p t | ƒ } n | i | i j p t$ | i t ƒ | } n | i" ƒ qèn | o | | f Sn | Sd S( sJ average(a, axis=0, weights=None)
Computes average along indicated axis.
If axis is None, average over the entire array
Inputs can be integer or floating types; result is of type Float.
If weights are given, result is sum(a*weights)/(sum(weights)*1.0)
weights must have a's shape or be the 1-d with length the size
of a in the given axis.
If returned, return a tuple: the result and the sum of the weights
or count of values. Results will have the same shape.
masked values in the weights will be set to 0.0
i c C s | | S( N( R- RŸ ( R- RŸ ( ( R Rõ B s f1.0f0.0R RE i s w[s ] * ones(ash, Float)s average: weights wrong shape.s( ] * masked_array(ones(ash, Float), mask)N( i ( f1.0f0.0( i ( f1.0f0.0( i (% Rp R8 R R? t ashR‚ R t weightsRÇ R R‹ R R R… RI Rú R* R RF R t wsht nit NewAxisR¬ t rt slicet evalt reprR* t
ValueErrorRÕ RG R R) R® t returnedRø (
R8 R‚ R3 R< RG R5 R2 R… R4 R R R7 Rú ( ( R t average* s¢
$
c C s t t | d ƒ d ƒ } | t j o d }
d } n0 t | ƒ }
t | ƒ } | d j o
d } n | t j o d } d } n0 t | ƒ } t | ƒ } | d j o
d } n t
i | | |
f ƒ } t
i | | | f ƒ } t | ƒ } t t | | ƒ d d d d ƒ} t | | ƒ S( s where(condition, x, y) is x where condition is nonzero, y otherwise.
condition must be convertible to an integer array.
Answer is always the shape of condition.
The type depends on x and y. It is integer if both x and y are
the value masked.
i i RE R N( RI Rë t conditiont fcR- Rr t xvt xmR< R RŸ t yvt ymR* R R… R1 R= RL RQ Rp ( R> R- RŸ R1 RC R… RA R= R? RB R@ ( ( R R~ – s*
!c C sÇ d „ }
d „ } t | d ƒ } g } | D] } | | | ƒ q, ~ } g } | D] } | |
| ƒ qS ~ } t
i | | ƒ } t
i | | ƒ } t t | t | ƒ ƒ d d d d ƒ} t | | ƒ S( s= Returns array shaped like indices with elements chosen from tc C s | t j o d Sn t | ƒ S( Ni ( R- Rr RI ( R- ( ( R t fmask´ s
c C s: | t j o d Sn t | ƒ } | d j o d Sn | S( Ni i ( R- Rr R< R= R ( R- R= ( ( R t nmask· s
i RE R i N( RD RE RI R- RŒ t _[1]R} R- t masksR8 R* R R… R= RL RQ R< Rp ( R- R} R8 RŒ RE R… R= RG RF R- RD ( ( R R ² s '''c C s@ t | d ƒ } t t | ƒ | ƒ } t t | ƒ d | d | ƒS( sj Return x as an array masked where condition is true.
Also masked where x or condition masked.
i RE R N( RI R> t cmRQ R< R- R= RF RE ( R> R- RE RH R= ( ( R t masked_whereÄ s c C s t t | | ƒ | | ƒ S( s3 masked_greater(x, value) = x masked where x > valueN( RI R] R- RS RE ( R- RS RE ( ( R t masked_greaterÌ s c C s t t | | ƒ | | ƒ S( s: masked_greater_equal(x, value) = x masked where x >= valueN( RI Rð R- RS RE ( R- RS RE ( ( R t masked_greater_equalÐ s c C s t t | | ƒ | | ƒ S( s0 masked_less(x, value) = x masked where x < valueN( RI R^ R- RS RE ( R- RS RE ( ( R t masked_lessÔ s c C s t t | | ƒ | | ƒ S( s7 masked_less_equal(x, value) = x masked where x <= valueN( RI Re R- RS RE ( R- RS RE ( ( R t masked_less_equalØ s c C sL t | d ƒ } t i | | ƒ } t | t | ƒ ƒ } t
| d | d | ƒS( s6 masked_not_equal(x, value) = x masked where x != valuei R RE N( RI R- R… R* Rë RS RŒ RQ R< R= RF RE ( R- RS RE RŒ R… R= ( ( R t masked_not_equalÜ s
c C sL t | d ƒ } t i | | ƒ } t | t | ƒ ƒ } t
| d | d | ƒS( s{ masked_equal(x, value) = x masked where x == value
For floating point consider masked_values(x, value) instead.
i R RE N( RI R- R… R* Ré RS RŒ RQ R< R= RF RE ( R- RS RE RŒ R… R= ( ( R t masked_equalã s
c C s‡ | | j o | } | } | } n t | d ƒ } t i t i | | ƒ t i | | ƒ ƒ } t | t | ƒ ƒ } t | d | d | ƒS( sj x with mask of all values of x that are inside [v1,v2]
v1 and v2 can be given in either order.
i R RE N( t v2t v1R} RI R- R… R* Rƒ Re Rð RŒ RQ R< R= RF RE ( R- RQ RP RE RŒ R… R= R} ( ( R t
masked_insideì s
*c C s‡ | | j o | } | } | } n t | d ƒ } t i t i | | ƒ t i | | ƒ ƒ } t | t | ƒ ƒ } t | d | d | ƒS( sk x with mask of all values of x that are outside [v1,v2]
v1 and v2 can be given in either order.
i R RE N( RP RQ R} RI R- R… R* RP R^ R] RŒ RQ R< R= RF RE ( R- RQ RP RE RŒ R… R= R} ( ( R t masked_outsideù s
*c C s_ t | ƒ } t i t | ƒ | ƒ } | d j o t | ƒ Sn t | d t i | | ƒ ƒSd S( s Copy of a with a new shape.R N(
R< R8 R= R* t reshapeRI R› R… R Rp ( R8 R› R… R= ( ( R RT s
c C sY t | ƒ } t i t | ƒ ƒ } | d j o t | ƒ Sn t | d t i | ƒ ƒSd S( s- a as one-dimensional, may share data and maskR N( R< R8 R= R* RK RI R… R Rp ( R8 R… R= ( ( R RK s
c C s» g } x! | D] } | i t | ƒ ƒ q
Wt i | | ƒ } x0 | D] } t | ƒ d j o PqC qC Wt
| ƒ Sg } x! | D] } | i t | ƒ ƒ q| Wt i | | ƒ } t
| d | ƒS( s+ Concatenate the arrays along the given axisR N(
R… t arraysR- t appendRI R* t concatenateR‚ R< R Rp Rç R@ ( RU R‚ Rç R… R- ( ( R RW s"
c C sz t | ƒ } t | ƒ i ƒ } | d j o t t i | | | ƒ ƒ Sn/ t t i | | | ƒ d t i | | | ƒ ƒSd S( s; take(a, indices, axis=0) returns selection of items from a.R N( R< R8 R= Rp R‹ R… R R* t takeR- R‚ ( R8 R- R‚ R… R= ( ( R RX ( s
c C sk t | ƒ } t | ƒ } | d j o t t i | | ƒ ƒ Sn) t t i | | ƒ d t i | | ƒ ƒSd S( s9 transpose(a, axes=None) reorder dimensions per tuple axesR N(
R< R8 R= RI R… R Rp R* t transposet axes( R8 RZ R… R= ( ( R RY 2 s
c C s{ | i ƒ } t | ƒ } t | ƒ } t i | | | ƒ t
| ƒ } | d j o' | i
ƒ t i | i ƒ | d ƒ n d S( s{ put(a, indices, values) sets storage-indexed locations to corresponding values. values and indices are filled if necessary.i N( R8 R‹ R… RI R- R
R R
R* R R< R= R R˜ R; ( R8 R- R R… R= R
R
( ( R R < s
c C sv | d j o d Sn t i | i ƒ | | ƒ t | ƒ } | d j o d Sn | i ƒ t i | i
ƒ | d ƒ d S( s0 put (a, mask, values) sets a where mask is true.Ni ( R R R* R R8 R‹ R R< R= R˜ R; ( R8 R R R= ( ( R R G s
c C sz t | d ƒ } t | d ƒ } t | i ƒ d j o
d | _ n t | i ƒ d j o
d | _ n t t i | | ƒ ƒ S( sþ innerproduct(a,b) returns the dot product of two arrays, which has
shape a.shape[:-1] + b.shape[:-1] with elements computed by summing the
product of the elements from the last dimensions of a and b.
Masked elements are replace by zeros.
i i N( i ( i (
RI R8 R/ RW R0 R¬ R? Rp R* Rÿ ( R8 RW R/ R0 ( ( R Rÿ Q s
c C sÅ t | d ƒ i } t | d ƒ i } t i | | ƒ } t | ƒ } t | ƒ } | d j o | d j o t
| ƒ Sn t | ƒ } t | ƒ } t d t i d | d | ƒ d d ƒ} t
| | ƒ S( s: outerproduct(a,b) = {a[i]*b[j]}, has shape (len(a),len(b))i i RE N( RI R8 R R/ RW R0 R* R R… R< Rz R{ R Rp R@ RL R= ( R8 RW Rz R… R{ R= R/ R0 ( ( R R ] s *c C s. t t | d ƒ t i t | d ƒ d d ƒ ƒ S( sØ dot(a,b) returns matrix-multiplication between a and b. The product-sum
is over the last dimension of a and the second-to-last dimension of b.
Masked values are replaced by zeros. See also innerproduct.
i iÿÿÿÿiþÿÿÿN( Rÿ RI R8 R* t swapaxesRW ( R8 RW ( ( R R k s iÿÿÿÿc C si t | d ƒ } t | ƒ } | d j o t i | | | ƒ } n t i | t | ƒ | ƒ } t | | ƒ S( sq Select those parts of x for which condition is true.
Masked values in condition are considered false.
i N( RI R> RŒ R< R- R= R R* Ró t dimensionR… Rp ( R> R- R\ RŒ R… R= ( ( R Ró r s
t _minimum_operationc B s8 t Z d Z d „ Z e d „ Z d d „ Z d „ Z RS( s Object to calculate minimac C s d S( s` minimum(a, b) or minimum(a)
In one argument case returns the scalar minimum.
N( ( R ( ( R R s c C sŸ | d j oq t | ƒ } | d j o t t | ƒ i ƒ } | Sn | i ƒ } t | ƒ d j o t Sq› t | i
ƒ ƒ Sn t t | | ƒ | | ƒ d Sd S( s Execute the call behavior.i .N( RW R R< R8 R= t minRI R R… Rô t acR¬ Rr R‹ R~ R^ ( R R8 RW R… R= R_ ( ( R R_ … s
i c C s• t | ƒ } | d j o) t | ƒ } t t i i | | ƒ ƒ SnP t i i t | t | ƒ ƒ | ƒ } t i i | | ƒ } t | | t
| ƒ ƒ Sd S( s# Reduce target along the given axis.N( R< R€ R= R RI R} Rp R* t minimumR R‚ R5 Rƒ R: ( R R€ R‚ R= R} ( ( R R ” s
$c C s” t | ƒ } t | ƒ } | d j o | d j o
d } n+ t | ƒ } t | ƒ } t i | | ƒ } t
i i t | ƒ t | ƒ ƒ } t | | ƒ S( s< Return the function applied to the outer product of a and b.N( R< R8 Rz RW R{ R R= R@ RP R„ R* R` RI R… Rp ( R R8 RW Rz R… R{ R= ( ( R R„ Ÿ s
!( R R
R R R R_ R R„ ( ( ( R R] } s
t _maximum_operationc B s8 t Z d Z d „ Z e d „ Z d d „ Z d „ Z RS( s Object to calculate maximac C s d S( s` maximum(a, b) or maximum(a)
In one argument case returns the scalar maximum.
N( ( R ( ( R R ° s c C sŸ | d j oq t | ƒ } | d j o t t | ƒ i ƒ } | Sn | i ƒ } t | ƒ d j o t Sq› t | i
ƒ ƒ Sn t t | | ƒ | | ƒ d Sd S( s Execute the call behavior.i .N( RW R R< R8 R= t maxRI R R… Rô R_ R¬ Rr R‹ R~ R] ( R R8 RW R… R= R_ ( ( R R_ ¶ s
i c C s• t | ƒ } | d j o) t | ƒ } t t i i | | ƒ ƒ SnP t i i t | t | ƒ ƒ | ƒ } t i i | | ƒ } t | | t
| ƒ ƒ Sd S( s# Reduce target along the given axis.N( R< R€ R= R RI R} Rp R* t maximumR R‚ R6 Rƒ R: ( R R€ R‚ R= R} ( ( R R Å s
$c C s” t | ƒ } t | ƒ } | d j o | d j o
d } n+ t | ƒ } t | ƒ } t i | | ƒ } t
i i t | ƒ t | ƒ ƒ } t | | ƒ S( s< Return the function applied to the outer product of a and b.N( R< R8 Rz RW R{ R R= R@ RP R„ R* Rc RI R… Rp ( R R8 RW Rz R… R{ R= ( ( R R„ Ð s
!( R R
R R R R_ R R„ ( ( ( R Ra ® s
c C sr | d j o t | ƒ } n t | | ƒ } t i | | ƒ } t
| ƒ d j o t | ƒ Sn t | | d d ƒS( s± If x does not have a mask, return a masked array formed from the
result of Numeric.sort(x, axis).
Otherwise, fill x with fill_value. Sort it.
Set a mask where the result is equal to fill_value.
Note that this may have unintended consequences if the data contains the
fill value at a non-masked site.
If fill_value is not given the default fill value for x's type will be
used.
RE i N(
R: R R1 R- RI R… R* t sortR‚ R R< Rp R! ( R- R‚ R: R… R ( ( R Rd ß s
c C sk t i t | ƒ | | | ƒ } t | ƒ } | d j o t | | ƒ Sn t | t i | | | | ƒ ƒ Sd S( s9 diagonal(a,k=0,axis1=0, axis2=1) = the k'th diagonal of aN( R* t diagonalRI R8 t kt axis1t axis2R… R< R= R Rp ( R8 Rf Rg Rh R= R… ( ( R Re ò s
c C s t | | ƒ } t i | | ƒ S( sÊ Treating masked values as if they have the value fill_value,
return sort indices for sorting along given axis.
if fill_value is None, use fill_value(x)
Returns a Numeric array.
N( RI R- R: R… R* t argsortR‚ ( R- R‚ R: R… ( ( R Ri û s c C s t | | ƒ } t i | | ƒ S( s Treating masked values as if they have the value fill_value,
return indices for minimum values along given axis.
if fill_value is None, use fill_value(x).
Returns a Numeric array if x has more than one dimension.
Otherwise, returns a scalar index.
N( RI R- R: R… R* t argminR‚ ( R- R‚ R: R… ( ( R Rj s c C sQ | d j o% t | ƒ } y | } Wq2 q2 Xn t | | ƒ } t i | | ƒ S( s$ Treating masked values as if they have the value fill_value,
return sort indices for maximum along given axis.
if fill_value is None, use -fill_value(x) if it exists.
Returns a Numeric array if x has more than one dimension.
Otherwise, returns a scalar index.
N( R: R R1 R- RI R… R* t argmaxR‚ ( R- R‚ R: R… ( ( R Rk s
c C s t t i | | ƒ ƒ S( s+ apply f to s to create array as in Numeric.N( Rp R* t fromfunctionRi R ( Ri R ( ( R Rl s c C sN t | t ƒ o( | d j p | | i ƒ j o | Sn t | d | d d ƒS( s¥ asarray(data, typecode=None) = array(data, typecode=None, copy=0)
Returns data if typecode if data is a MaskedArray and typecode None
or the same.
R, RE i N( R RŠ R) R, R RF ( RŠ R, ( ( R t asarray" s 0i, c C s» | i } t i i | ƒ } | t j os t | ƒ d j o0 d t i
t t | ƒ d ƒ | i
ƒ | f Sq· t i | t | | | d d ƒ d | t Sn t i | | | | d d ƒ Sd S( s9 Calculate string representation, limiting size of output.i s' array (%s) , type = %s, has %d elementst ,i s
+ %d more elementsN( R8 R? R R* RÓ R t elemst
__MaxElementsR¬ t stringt joint mapR R, t array2stringt max_line_widtht precisiont suppress_small( R8 Ru Rv Rw Ro R ( ( R t limitedArrayRepr/ s
0 c C s” d k } t | ƒ } | a | d j o6 t t _ t t _ | i
t d ƒ | i
t d ƒ n3 t t _ t t _ | i
t d ƒ | i
t d ƒ d S( sB Set the maximum # of elements for printing arrays. <=0 = no limitNi i ( t
multiarrayR± R= R Rp t __original_strR* t array_strt __original_reprt
array_reprt set_string_functionRx ( R= Ry R ( ( R t set_print_limitC s
c C s t S( s3 Get the maximum # of elements for printing arrays. N( Rp ( ( ( R t get_print_limitT s (© R R* Rq R R2 t PrecisionR t piR6 t Int0RB Ru t ExceptionR R
R§ RF RH t Float32R! t Complex32R( R& t UnsignedInt8R$ R0 R1 R5 R6 R9 R< R@ RC RL R> t create_maskRQ R RI R: RY RZ R` Rc Rf Rg Rs Rv R Rá t logt log10t expt conjugatet sinRb t tant arcsint arccost arctant sinht cosht tanhR/ t fabsRÄ t nonzerot aroundt floort ceilRJ R t logical_notRÇ RÐ R RÓ RÕ RØ RÛ RÉ t fmodt hypott arctan2Ré Rë Re Rð R^ R] Rƒ RP t logical_xorR¼ R¿ RÁ t rankR? R” t objectR) R Rr R7 t isarrayt isMAR R R! R R" R R& RÌ RÎ R— R+ R, R- Rù RM Rø Rþ RÞ Rp Rý t productR= R~ R RI RJ RK RL RM RN RO RR RS RT RK RW RX RY R R Rÿ R R Ró R] R` Ra Rc Rd Re Ri Rj Rk Rl Rm Rp Rx R{ Rz R} R| R R€ (š RL Rˆ Rz R^ R RR R£ R? R€ R0 Ra Rp R RÇ R R6 R, R— R™ Rø RŸ Rv Rq R¡ R" R¿ R R* R> R› R- R¼ R~ Rk R) R“ Rœ Rƒ R R R– Rm Rý RÓ Rr R! RI RW RÞ RK RY Re RO R’ R9 RÉ R+ Rb R1 R R Ré Rs Ri R Re R R§ RÕ RZ R$ R Rg RB RC RF R” Rž Rx R‰ RJ RÐ RÄ Ru R RL R` RS RM RY R‚ RŠ Rd R( R R R— R Rþ Rë Rl Rj R” RØ Rc RK R| R] R RÌ Rc RM RP R` RJ RŒ Rÿ RŽ R/ RÁ R] R• R˜ RN RT Rá R‘ R< RÎ RX R6 RÛ RQ R Rf R7 Ró RI R2 R5 R R R& Rš Rð R& R= R: R
R! R‹ R@ R¢ R ( ( R R sH
',K ÿ ÿ ÿ [
l
/ /