Ñò šÙÊKc@sB dZddkZddkZddkiiZddkiiZddkl Z l Z l Z ddkl Z l Z ddklZlZddkiiZddkZhZhZeZedƒZdZdefd „ƒYZd fd „ƒYZed ƒZd Zd dZdZdZdZ d„Z!d„Z"d„Z#d„Z$d„Z%d„Z&d„Z'ddd„Z(d„Z)d„Z*e+d„Z,d„Z-d„Z.dfd „ƒYZ/d!fd"„ƒYZ0d#fd$„ƒYZ1d%fd&„ƒYZ2d'fd(„ƒYZ3d)fd*„ƒYZ4d+fd,„ƒYZ5d-fd.„ƒYZ6e3ei7d/e2d/ƒƒZ7e3ei8d0e1d/ƒƒZ8e3ei9d0e1d/ƒƒZ9e3ei:ƒZ:e3ei;ƒZ;e3ei<ƒZ<e3ei=ƒZ=e3ei>d/e0dƒƒZ>e3ei?d/e/d1d0ƒƒZ?e3ei@d/e/d1d0ƒƒZ@e3eiAƒZAe3eiBƒZBe3eiCd0e2d0ƒƒZCe3eiDd/e/d1d2d0d2ƒƒZDe3eiEƒZEe3eiFƒZFe3eiGƒZGe3eiHƒZHe3eiIƒZIe3eiJƒZJd3„ZKe3eiLƒZMe3eiNƒZNe3eiOƒZOe3eiPƒZPe6eiQƒZQe6eiRƒZRe+eR_Se6eiTd4d4ƒZTe5eiUe4ƒdd4ƒZUe5eiVe4ƒdd4ƒZVe5eiWe4ƒdd4ƒZWe5eiXe4ƒdd4ƒZXe5eiYe4ƒdd4ƒZYe6eiZƒZZe6ei[d/d0ƒZ[e+e[_Se6ei\ƒZ\e+e\_Se6ei]ƒZ]e+e]_Se6ei^ƒZ^e+e^_Se6ei_ƒZ_e+e__Se6ei`ƒZ`e+e`_Se6eiaƒZae+ea_Se6eibƒZbe6eibd4d4ƒiSZce6eidƒZdediSZee6eifƒZfe6eigƒZge6eihƒZhe6eiiƒZid5„Zjd6„Zke+d7„Zld8emfd9„ƒYZnenZod:„ZpepZqepZrd4d;d<d=„Zsd4d>„Ztd;d<d4d?„Zud4d@„Zve+d4e+dA„ZwewZxdB„ZydC„ZzdD„Z{dE„Z|e+dF„Z}dG„Z~e+dH„Ze€dI„Ze€dJ„Z‚e+dK„Zƒe+dL„Z„ee+dM„Z…e+e+dN„Z†e+e+dO„Z‡e+e+ddP„ZˆdQ„Z‰e+dRdS„ZŠd4dT„Z‹d4dU„ZŒd4dV„Zd4dW„ZŽd4dX„Zd4dY„Zd4dZ„Z‘d4d[„Z’d4d\„Z“d]„Z”d^„Z•dd_„Z–d`„Z—e+e+dRda„Z˜e+db„Z™dRdc„Zšdd„Z›de„ZœeœZdf„ZžežZŸdg„Z de+dh„Z¡difdj„ƒYZ¢e¢ƒZ£dkfdl„ƒYZ¤e¤ƒZ¥de+dm„Z¦ddd4dn„Z§ddd4e+e+do„Z¨de+e+dp„Z©de+e+dq„Zªde+e+dr„Z«ds„Z¬e+dt„Z­ddukl®Z®dv„Z¯dw„Z°e¯ecƒeo_±e¯eeƒeo_²e¯e«ƒeo_«e¯eªƒeo_ªe¯e©ƒeo_©e³e¯e°ƒƒeo_´e¯e°ƒeo_µdx„Z¶e¯e¶ƒeo_Š[¶e+dy„Z·e¯e·ƒeo_¸e+e+dz„Z¹e¯e¹ƒeo_¡[¹e¯e;ƒeo_ºeo_;e¯e°ƒeo_»e+e+e+d{„Z¼e¯e¼ƒeo_½e+e+e+d|„Z¾e¯e¾ƒeo_¿e¯e§ƒeo_§e¯e°ƒeo_Àe¯e°ƒeo_Áe¯e°ƒeo_Âe³e¯e°ƒƒeo_Ãe¯e•ƒeo_Äe¯e°ƒeo_Åe+e+d}„ZÆe¯eƃeo_Ç[Æe+e+d~„ZÈe¯eȃeo_É[Èe¯eˆƒeo_Êe³e¯e°ƒƒeo_Ëe¯e°ƒeo_Ìe¯eKƒeo_Ke¯e‡ƒeo_Íe+e+d„ZÎe¯e΃eo_Ïe¯e}ƒeo_Ðe¯e|ƒeo_|e¯e°ƒeo_Ñe¯e°ƒeo_Òe¯e°ƒeo_Óe¯e°ƒeo_¦d€„ZÔe¯eÔƒeo_Õe³e¯e°ƒƒeo_Öe¯e†ƒeo_†d„Z×e¯e׃eo_—e¯e˜ƒeo_Øe¯e°ƒeo_Ùe¯e¨ƒeo_¨e¯e™ƒeo_™e+e+e+d‚„ZÚe+e+e+dƒ„ZÛe¯eÚƒeo_Üe¯eÛƒeo_Ýe¯e°ƒeo_Þe¯eMƒeo_ß[¯[®[°endeàd„d4ƒZádd…„ZÐde+dd†„Zâdd‡„ZØdS(ˆs0MA: a facility for dealing with missing observations MA is generally used as a numpy.array look-alike. by Paul F. Dubois. Copyright 1999, 2000, 2001 Regents of the University of California. Released for unlimited redistribution. Adapted for numpy_core 2005 by Travis Oliphant and (mainly) Paul Dubois. iÿÿÿÿN(tnewaxistndarraytinf(tamaxtamin(tbool_t typecodesigzø·¥•ª8tMAErrorcBs#eZdd„Zd„ZeZRS(cCs*t|tƒp |f}n||_dS(sCreate an exceptionN(t isinstancettupletargs(tselfR ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt__init__s cCst|idƒS(s#Calculate the string representationi(tstrR (R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt__str__&sN(t__name__t __module__tNoneR Rt__repr__(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRs  t_MaskedPrintOptioncBsMeZdZd„Zd„Zd„Zd„Zdd„Zd„ZeZ RS(s<One instance of this class, masked_print_option, is created.cCs|i|ƒd|_dS(s&Create the masked print option object.iN(t set_displayt_enabled(R tdisplay((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR -s cCs|iS(s#Show what prints for masked values.(t_display(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR2scCs ||_dS(s2set_display(s) sets what prints for masked values.N(R(R ts((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR6scCs|iS(s(Is the use of the display value enabled?(R(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytenabled:sicCs ||_dS(sSet the enabling flag to flag.N(R(R tflag((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytenable>scCs t|iƒS(N(R R(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRBs( RRt__doc__R RRRRRR(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR+s      s--g@Œµx¯Dyt-i?Bt?cCs#t|tiƒotSt|tiƒpt|tiƒotSt|tiƒotSt|ti ƒot St|t ƒpt|t ƒoˆ|i i}|tdjotS|tdjotS|tdjot S|tdjotS|tdjotitƒStStSdS(s7Function to calculate default fill value for an object.tFloattIntegertComplext CharactertUnsignedIntegerN(Rttypest FloatTypetdefault_real_fill_valuetIntTypetLongTypetdefault_integer_fill_valuet StringTypetdefault_character_fill_valuet ComplexTypetdefault_complex_fill_valuet MaskedArrayRtdtypetcharRtumathtabsolutetdefault_object_fill_value(tobjtx((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytdefault_fill_valueQs,&  cCsÑt|tiƒotiSt|tiƒpt|tiƒotiSt|t ƒpt|t ƒo[|i i }|t djotiS|t djotiS|t djotiSn td‚dS(sDFunction to calculate default fill value suitable for taking minima.RR R#s(Unsuitable type for calculating minimum.N(RR$R%tnumericRR'R(tsystmaxintR.RR/R0Rt TypeError(R4R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytminimum_fill_valueks&   cCsÌt|tiƒot St|tiƒpt|tiƒo ti St|tƒpt|t ƒoW|i i }|t djot S|t djo ti S|t djodSn t d‚dS(sDFunction to calculate default fill value suitable for taking maxima.RR R#is(Unsuitable type for calculating maximum.N(RR$R%RR'R(R8R9R.RR/R0RR:(R4R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmaximum_fill_value|s&     cCs"t|ƒo|i|ƒndS(s,Set fill value of a if it is a masked array.N(t isMaskedArraytset_fill_value(tat fill_value((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR>s cCs#t|tƒo |iƒStSdS(s‰Mask of values in a; could be nomask. Returns nomask if a is not a masked array. To get an array for sure use getmaskarray.N(RR.traw_masktnomask(R?((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytgetmask’s cCs2t|ƒ}|tjott|ƒƒS|SdS(s¦Mask of values in a; an array of zeros if mask is nomask or not a masked array, and is a byte-sized integer. Do not try to add up entries, for example. N(RCRBtmake_mask_nonetshape(R?tm((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt getmaskarray›s  cCs/y|iitjSWntj otSXdS(s;Is m a legal mask? Does not check contents, only type. N(R/ttypetMaskTypetAttributeErrortFalse(RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytis_mask¦scCs¾|tjotSt|tƒoV|iitjo0|oti|dtd|ƒ}qt|}q|itƒ}nt |t ƒitƒ}|ot i t i |ƒƒ otS|SdS(smake_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 nomask. Does not check that contents must be 0s and 1s. if flag, return nomask if m contains no true elements. R/tcopyN(RBRRR/RHRIR7tarraytastypetfilledtTruet fromnumerictsometruetravel(RFRMRtresult((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt make_mask®s  !cCs"ti|dtƒ}||_|S(s&Return a mask of all zeros of shape s.R/(R7tzerosRIRE(RRU((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRDÇs cCse|tjo t|ƒS|tjo t|ƒS||jot|ƒo|Stti||ƒƒS(s…Logical or of the mask candidates m1 and m2, treating nomask as false. Result may equal m1 or m2 if the other is nomask. (RBRVRLR1t logical_or(tm1tm2((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmask_orÍs  cCsvt|tƒo|i|ƒSt|tƒo|ido|St|tiƒoti|dƒSti|ƒSdS(s•a as a contiguous numeric array with any masked areas replaced by value if value is None or the special element "masked", get_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 CONTIGUOUStON( RR.RPRtflagsR$tDictTypeR7RN(R?tvalue((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRPÖs cCs-t|ƒo|iƒ}n t|ƒ}|S(sb The fill value of a, if it has one; otherwise, the default fill value for that type. (R=R@R6(R?RU((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytget_fill_valueês  cCs.t|ƒ}t|ƒ}||jo|SdS(s:The common fill_value of a and b, if there is one, or NoneN(RaR(R?tbtt1tt2((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytcommon_fill_valueõs    tdomain_check_intervalcBs eZdZd„Zd„ZRS(s9domain_check_interval(a,b)(x) = true where x < a or y > bcCs||_||_dS(s9domain_check_interval(a,b)(x) = true where x < a or y > bN(ty1ty2(R RgRh((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR ÿs cCs.titi||iƒti||iƒƒS(sExecute the call behavior.(R1RXtgreaterRhtlessRg(R R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt__call__s(RRRR Rk(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRfýs t domain_tancBs eZdZd„Zd„ZRS(s/domain_tan(eps) = true where abs(cos(x)) < eps)cCs ||_dS(s/domain_tan(eps) = true where abs(cos(x)) < eps)N(teps(R Rm((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR scCs%tititi|ƒƒ|iƒS(sExecute the call behavior.(R1RjR2tcosRm(R R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRks(RRRR Rk(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRl s tdomain_greatercBs eZdZd„Zd„ZRS(s(domain_greater(v)(x) = true where x <= vcCs ||_dS(s(domain_greater(v)(x) = true where x <= vN(tcritical_value(R Rp((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR scCsti||iƒS(sExecute the call behavior.(R1t less_equalRp(R R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRks(RRRR Rk(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRos tdomain_greater_equalcBs eZdZd„Zd„ZRS(s-domain_greater_equal(v)(x) = true where x < vcCs ||_dS(s-domain_greater_equal(v)(x) = true where x < vN(Rp(R Rp((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR scCsti||iƒS(sExecute the call behavior.(R1RjRp(R R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRk$s(RRRR Rk(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRrs tmasked_unary_operationcBs)eZddd„Zd„Zd„ZRS(icCsl||_||_||_t|dt|ƒƒ|_t|dt|ƒƒ|_|t|<|ft|d4„Z?d5„Z@d6„ZAd7„ZBd^d8„ZCd9„ZDd:„ZEd^d;„ZFd<„ZGd=„ZHd>„ZId?„ZJd@„ZKdA„ZLdB„ZMdC„ZNeOdDeNdEdFƒZPdG„ZQeOdDeQdEdHƒZRdI„ZSd^dJ„ZTdK„ZUeOeUdEeViWiXiƒZXdL„ZYeOdDeYdEdMƒZZdN„Z[eOdDe[dEdOƒZ\dP„Z]dQ„Z^d^dR„Z_d^dS„Z`dT„ZadU„ZbdV„ZcdW„ZdeOe e dEdXƒZeeOe edEdYƒZfeOeedEdZƒZgeOeedEd[ƒZhehZieOecd^dEd\ƒZjeOedd^dEd]ƒZkRS(_sœArrays with possibly masked values. Masked values of 1 exclude the corresponding element from any computation. Construction: x = array(data, dtype=None, copy=True, order=False, mask = nomask, fill_value=None) If copy=False, every effort is made not to copy the data: If data is a MaskedArray, and argument mask=nomask, then the candidate data is data.data and the mask used is data.mask. If data is a numeric array, it is used as the candidate raw data. If dtype is not None and is != data.dtype.char 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, dtype=dtype.char, copy=copy) If mask is nomask there are no masked values. Otherwise mask must be convertible to an array of booleans with 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. g333333$@c Csù|djo d}nti|ƒ}|}t|tƒor|i} |djo | i}n|| ijo t}n|tjo |i}q,|tj o t}q,nt|t ƒo>|} |djo | i}q,|| ijo t}q,n1t }ti |d|dtd|ƒ} | i}|oN|| ijo(ti | d|dtd|ƒ|_ qŠ| i |ƒ|_ n | |_ |tjot|_d|_n<t|ƒ|_|itjo d|_n|i|j|_t|iƒ} t|i ƒ} | | jo†| djo+ti|i|i iƒ|_d|_qè| djo4ti|i |iiƒ|_ |ii|i _qètd‚nM| djo?t|iƒt|i ƒjo |iƒ|i i|i_n|i|ƒdS(s°array(data, dtype=None, copy=True, order=False, mask=nomask, fill_value=None) If data already a numeric array, its dtype becomes the default value of dtype. R/RMtorderiisMask and data not compatible.N(RR7R/RR.tdataRQRBtmaskRRKRNt_dataROt_maskt _shared_maskRVR™RRtresizeRERt unshare_maskR>( R R›R/RMRšRœR@ttctneed_data_copiedtctnmtnd((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR sb           ! (          , cCsÕ|itj o¢ti|iƒiƒos|djotidƒ|iS|\}}}t i |ƒ}|djot d|‚n|i ||ƒSq²t|_d|_ n|o|ii|ƒS|iSdS(s:Special hook for numeric. Converts to numeric if possible.siCannot automatically convert masked array to numeric because data is masked in one or more locations.s%s not known to maiN(RžRBRRRTtanyRtwarningstwarnRRytgetRRPRŸRO(R R‰tcontexttfuncR titfills((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt __array__]s     c CsI|djot|dtdtƒS|d \}}t|}ttg}|D]}|t|ƒqO~ƒ}|dj o@t||g} |D]}| t|d|ƒq~ Œƒ}n|tj oky |i } Wnt j oq3X|i | jo4ttg} |D]}| t |ƒq ~ ƒ}q3nt|dtd|ƒS(spSpecial hook for ufuncs. Wraps the numpy array and sets the mask according to context. RMRœiRN( RR.RKRBRxRR[RCRwRERJRG( R RNR«R¬R Rvt_[1]R?RFt_[2]REt_[3]((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt__array_wrap__ys   0 0  8cCs |iiS(sReturn the current shape.(RRE(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _get_shape’scCsB||i_|itj o"|iiƒ|_||i_ndS(sSet the array's shape.N(RRERžRBRM(R tnewshape((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _set_shape–s cCsh|itjo&t|iiƒdtd|iƒƒSt|iiƒd|iiƒd|iƒƒSdS(s"Calculate the flat value. RœR@N(RžRBR{RRTR@(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _get_flats cCs|iƒ}||(dS(sx.flat = valueN(RT(R R`ty((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _set_flat¨s cCs\|itjo#t|iidtd|iƒƒSt|iid|id|iƒƒSdS(s%Get the real part of a complex array.RœR@N(RžRBR{RtrealR@(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _get_real­s cCs|i}||džsN(RžRRERÃRBRRtlisttonesRRRTRORÍR™R1Rãtsum( R R‘RFRtlsRÈR‰twtn1tn2((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytcount“s0              cCs t||ƒS(s%s.dot(other) = innerproduct(s, other)(t innerproduct(R RÙ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytdot±scCs|iS(sGet the current fill value.(RÏ(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR@µsc Cs |i}|i}|tjo!|ido|S|iƒSnÉ|djo |i}n|}|tjoti |ƒ}n…y,ti |d|i ddƒ}|||"s cCs |iiS(N(RR‹(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _get_ndim(scCs |iiS(N(RR™(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _get_size,ss Number of elements in the array.cCs |iiS(N(RR/(R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _get_dtype1sstype of the array elements.cGsa|itj oA|ii|Œ}y|dotSWqQtj otSXn|ii|ŒS(s Return Python scalar if possiblei(RžRBRËRt IndexErrorR(R R RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRË5s   cGs"|d}|d }|||R,R7RR‹R-R™R.R/RËR0R1R2RÉR¡R4R6REtflatRºR¿R½R3tT(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR.þsÔ A        /                            + + - '            "                              cCs t|tƒS(s9Is x a masked array, that is, an instance of MaskedArray?(RR.(R5((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR=zsgñh㈵øä>g:Œ0âŽyE>c CsÍtt|ƒt|ƒƒ}t|ƒ}t|ƒ}tt|ddd|ƒ|ƒitƒ}tt|ddd|ƒdƒitƒ} titi|| ƒ||ti| ƒƒ} 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. RMiRœi( R[RCRPRNRORÊR1RqR2RRtalltrueRT( R?RbR@trtoltatolRFR}RˆR5R¸R“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytallcloses   **0cCsÚtt|ƒt|ƒƒ}|tjoAt|ƒ}t|ƒ}ti||ƒ}titi|ƒƒS|obt|ƒ}t|ƒ}ti||ƒ}t |d|ddƒ}titit|dƒƒƒSdSdS(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œRMiiN( R[RCRBRPR1RRRR@RTRN(R?RbR@RFR5R¸R“R((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytallequal’s      cCs¡ti}t||ƒ}t|iitiƒoYti|||ƒ||||ƒƒ}t |ddƒ}t |d|d|d|ƒSt ||d|ƒSdS(sª masked_values(data, value, rtol=1.e-5, atol=1.e-8) Create a masked array; mask is nomask 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) RiRœRMR@N( R1R2RPt issubclassR/RHR7tfloatingRqRVRNt masked_object(R›R`RARBRMtabsR“RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt masked_values¦s *cCsIt||ƒ}tti||ƒddƒ}t|d|d|d|ƒS(s5Create array masked where exactly data equal to valueRiRœRMR@(RPRVR1RRN(R›R`RMR“R((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRG»scCstti||||ƒƒS(soJust like range() except it returns a array whose type can be specified by the keyword argument dtype. (RNR7R"(tstarttstoptstepR/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR"ÁscCstti||ƒƒS(sAConstruct a masked array from a string. Result will have no mask.(R{R7t fromstring(RR‰((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRMÉscCsht|ƒ}|tjo#tit|ƒ|ƒ}t|ƒStit|dƒ|ƒ}t||ƒSdS(sLeft shift n bitsiN(RCRBR1RèRPR{(R?RÈRFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRèÍs    cCsht|ƒ}|tjo#tit|ƒ|ƒ}t|ƒStit|dƒ|ƒ}t||ƒSdS(sRight shift n bitsiN(RCRBR1RêRPR{(R?RÈRFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRê×s    cCsjt|ƒ}|tj oti||ƒ}nttit|ƒ|ƒd|ƒ}|it|ƒƒ|S(swresize(a, new_shape) returns a new array with the specified shape. The original array's total size can be any size.Rœ(RCRBRRR RNRPR>Ra(R?t new_shapeRFRU((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR ás   $cCsÊt|ƒ}t|tiƒo>|djo |i}n|i|}t|g|ƒ}nt|ƒ}|t j ot i |||ƒ}nt i |||ƒ}t ||ƒ}|i t|ƒƒ|S(srepeat 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(RPRR$R'RR™RER RCRBRRtrepeatR{R>Ra(R?trepeatsR‘taftnumRFR“RU((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt new_repeatës      cCstti|ƒƒS(s<identity(n) returns the identity matrix of shape n x n. (RNR7tidentity(RÈ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRTscCstti||ƒƒS(sindices(dimensions,dtype=None) returns an array representing a grid of indices with row-only, and column-only variation. (RNR7tindices(t dimensionsR/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRUscCstti||ƒƒS(sPzeros(n, dtype=float) = an array of all zeros of the given length or shape.(RNR7RW(RER/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRW scCstti||ƒƒS(sNones(n, dtype=float) = an array of all ones of the given length or shape.(RNR7R(RER/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRscCst|ƒ}|i|ƒS(s?Count of the non-masked elements in a, or along a certain axis.(R{R(R?R‘((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRs c Cs|dj o td‚nt|ƒ}t|ƒ}t||ƒ}t|dƒ}t|dƒ}|iitdjott i ||ƒ|ƒSt t i |dƒddƒ}t||ƒ}|t jott i ||ƒƒSti|d|ƒ}tt i ||ƒ|ƒSdS(sa**bs3-argument power not supported.iR iRN(RRRCR[RPR/R0RR{R1RúRVRjRBR7R…( R?RbRûR†R‡RFtfatfbtmd((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRús     cCst|d|ddd|ƒS(s]masked_array(a, mask=nomask) = array(a, mask=mask, copy=0, fill_value=fill_value) RœRMiR@(RN(R?RœR@((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR{/scCs6|djot|ƒ}d}nti|||ƒS(Ni(RRTRãR(RR‘R/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR5s   cCs6|djot|ƒ}d}nti|||ƒS(Ni(RRTRïR(RR‘R/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytproduct;s   c Cst|ƒ}|i}|i}|d jo d}n|d joD|tjo‘|d jo4ti|iƒiƒƒ}td„|dƒ}q‚t |dƒiƒ}t ii|iƒiƒ|ƒ}t ii|ƒ}~qQ|d joFti|iƒƒ}t i |dƒiƒ}t ii|ƒ}~qQt t |dƒtd|ƒiƒ}ti|iƒ|ƒ}ti|ƒ}~nÌ|tjoY|d jo-||d}t ii|iƒ|ƒ}qQt |dƒ}|i} | djo d} n| |joGti |tddƒ}ti|||ƒ}ti||ƒ}~qQ| ||fjo{tgt|ƒ} td d dƒ| | value(RrRi(R5R`RM((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmasked_greaterØscCstt||ƒ||ƒS(s:masked_greater_equal(x, value) = x masked where x >= value(RrR (R5R`RM((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmasked_greater_equalÜscCstt||ƒ||ƒS(s0masked_less(x, value) = x masked where x < value(RrRj(R5R`RM((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt masked_lessàscCstt||ƒ||ƒS(s7masked_less_equal(x, value) = x masked where x <= value(RrRq(R5R`RM((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmasked_less_equaläscCsLt|dƒ}ti||ƒ}t|t|ƒƒ}t|d|d|ƒS(s6masked_not_equal(x, value) = x masked where x != valueiRœRM(RPR1RR[RCRN(R5R`RMR“R¤RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmasked_not_equalèscCsLt|dƒ}ti||ƒ}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. iRœRM(RPR1RR[RCRN(R5R`RMR“R¤RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt masked_equalïscCs‡||jo|}|}|}nt|dƒ}titi||ƒti||ƒƒ}t|t|ƒƒ}t|d|d|ƒS(sjx with mask of all values of x that are inside [v1,v2] v1 and v2 can be given in either order. iRœRM(RPR1RŽRqR R[RCRN(R5tv1tv2RMR‰R“R¤RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt masked_insideøs  *cCs‡||jo|}|}|}nt|dƒ}titi||ƒti||ƒƒ}t|t|ƒƒ}t|d|d|ƒS(skx with mask of all values of x that are outside [v1,v2] v1 and v2 can be given in either order. iRœRM(RPR1RXRjRiR[RCRN(R5RyRzRMR‰R“R¤RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytmasked_outsides  *cGsYt|ƒ}t|ƒi|Œ}|tjo t|ƒSt|dti||ŒƒSdS(sCopy of a with a new shape.RœN(RCRPRŒRBR{R7(R?RµRFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRŒs    cCsVt|ƒ}tit|ƒƒ}|tjo t|ƒSt|dti|ƒƒSdS(s-a as one-dimensional, may share data and maskRœN(RCRRRTRPRBR{R7(R?RFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRTs    c Cs»g}x!|D]}|it|ƒƒq Wti||ƒ}x0|D]}t|ƒtj oPqCqCWt|ƒSg}x!|D]}|it|ƒƒq|Wti||ƒ}t|d|ƒS(s+Concatenate the arrays along the given axisRœ(tappendRPR7t concatenateRCRBR{RG(tarraysR‘R“R5R((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR~$s   cCszt|ƒ}t|ƒi}|tjotdti|||ƒƒStdti|||ƒdti|||ƒƒSdS(NR›Rœ(RCR{R›RBR7tswapaxes(R?taxis1taxis2RFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR€4s   cCstt|ƒ}t|ƒi}|tjotti|||ƒƒStti|||ƒdti|||ƒƒSdS(s"returns selection of items from a.RœN(RCR{R›RBR7ttake(R?RUR‘RnRoRFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytnew_take>s   cCsht|ƒ}t|ƒ}|tjotti||ƒƒStti||ƒdti||ƒƒSdS(s!reorder dimensions per tuple axesRœN(RCRPRBR{R7R5(R?taxesRFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR5Is    cCs{|iƒ}t|ƒ}t|ƒ}ti|||ƒt|ƒ}|tj o'|iƒti|iƒ|dƒndS(snsets storage-indexed locations to corresponding values. Values and indices are filled if necessary. iN(R(RPR7R&RCRBR¡RA(R?RUR#RoR“R%R+RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR&Ts      cCsp|tjodSti|iƒ||ƒt|ƒ}|tjodS|iƒti|iƒ|dƒdS(s3putmask(a, mask, values) sets a where mask is true.Ni(RBR7R'R(RCR¡RA(R?RœR#RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR'cs    cCszt|dƒ}t|dƒ}t|iƒdjo d|_nt|iƒdjo d|_ntti||ƒƒS(s÷inner(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. ii(i(i(RPRÃRER{R7tinner(R?RbRWRX((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR†ms  cCsÈt|dƒiƒ}t|dƒiƒ}ti||ƒ}t|ƒ}t|ƒ}|tjo|tjo t|ƒSt|ƒ}t|ƒ}tdtid|d|ƒddƒ}t||ƒS(s3outer(a,b) = {a[i]*b[j]}, has shape (len(a),len(b))iiRM( RPRTR7R’RCRBR{RGRV(R?RbRWRXR“R†R‡RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR’{s     *cCs.tt|dƒtit|dƒddƒƒ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. iiÿÿÿÿiþÿÿÿ(RRPR7R€(R?Rb((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR‹scCsit|dƒ}t|ƒ}|tj oti|||ƒ}nti|t|ƒ|ƒ}t||ƒS(sqSelect those parts of x for which condition is true. Masked values in condition are considered false. i(RPRCRBR7RR{(ReR5t dimensionRnR¤RFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR’s   t_minimum_operationcBs8eZdZd„Zdd„Zdd„Zd„ZRS(sObject to calculate minimacCsdS(s`minimum(a, b) or minimum(a) In one argument case returns the scalar minimum. N((R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR ŸscCs˜|djont|ƒ}|tjott|ƒiƒƒ}|S|iƒ}t|ƒdjotSt|i ƒƒSnt t ||ƒ||ƒSdS(sExecute the call behavior.iN( RRCRBRRPRTRRÃRR(R…Rj(R R?RbRFR“tac((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRk¥s    icCs’t|ƒ}|tjo&t|ƒ}ttii||ƒƒStiit|t|ƒƒ|ƒ}tii||ƒ}t||t |ƒƒSdS(s#Reduce target along the given axis.N( RCRBRPR{R1tminimumRR;RŽRa(R RR‘RFR‰((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR´s   $cCs”t|ƒ}t|ƒ}|tjo|tjo t}n+t|ƒ}t|ƒ}ti||ƒ}tiit|ƒt|ƒƒ}t||ƒS(s<Return the function applied to the outer product of a and b.( RCRBRGRXR’R1RŠRPR{(R R?RbR†R‡RFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR’¿s     !N(RRRR RRkRR’(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRˆs    t_maximum_operationcBs8eZdZd„Zdd„Zdd„Zd„ZRS(sObject to calculate maximacCsdS(s`maximum(a, b) or maximum(a) In one argument case returns the scalar maximum. N((R ((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR ÐscCs˜|djont|ƒ}|tjott|ƒiƒƒ}|S|iƒ}t|ƒdjotSt|i ƒƒSnt t ||ƒ||ƒSdS(sExecute the call behavior.iN( RRCRBRRPRTRRÃRR(R…Ri(R R?RbRFR“R‰((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRkÖs    icCs’t|ƒ}|tjo&t|ƒ}ttii||ƒƒStiit|t|ƒƒ|ƒ}tii||ƒ}t||t |ƒƒSdS(s#Reduce target along the given axis.N( RCRBRPR{R1tmaximumRR<RŽRa(R RR‘RFR‰((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRås   $cCs”t|ƒ}t|ƒ}|tjo|tjo t}n+t|ƒ}t|ƒ}ti||ƒ}tiit|ƒt|ƒƒ}t||ƒS(s<Return the function applied to the outer product of a and b.( RCRBRGRXR’R1RŒRPR{(R R?RbR†R‡RFR“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR’ðs     !N(RRRR RRkRR’(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR‹Îs    cCso|djot|ƒ}nt||ƒ}ti||ƒ}t|ƒtjo t|ƒSt||ddƒ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. RMiN( RR6RPRRtsortRCRBR{RI(R5R‘R@R“R((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRÿs  cCshtit|ƒ|||ƒ}t|ƒ}|tjot||ƒSt|ti||||ƒƒSdS(s9diagonal(a,k=0,axis1=0, axis2=1) = the k'th diagonal of aN(RRtdiagonalRPRCRBR{(R?tkRR‚R“RF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRŽs   cCst||||ƒid|ƒS(s~trace(a,offset=0, axis1=0, axis2=1) returns the sum along diagonals (defined by the last two dimenions) of the array. R/(RŽR(R?toffsetRR‚R/Rn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyttracescCst||ƒ}ti||ƒ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 get_fill_value(x) Returns a numpy array. (RPRRtargsort(R5R‘RnR@R“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR’!scCst||ƒ}ti||ƒS(sTreating masked values as if they have the value fill_value, return indices for minimum values along given axis. if fill_value is None, use get_fill_value(x). Returns a numpy array if x has more than one dimension. Otherwise, returns a scalar index. (RPRRtargmin(R5R‘RnR@R“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR“*scCsQ|djo%t|ƒ}y | }Wq2q2Xnt||ƒ}ti||ƒ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 -get_fill_value(x) if it exists. Returns a numpy array if x has more than one dimension. Otherwise, returns a scalar index. N(RR6RPRRtargmax(R5R‘RnR@R“((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR”4s   cCstti||ƒƒS(s)apply f to s to create array as in umath.(R{R7t fromfunction(RtR((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR•DscCsHt|tƒo"|djp||ijo|St|d|ddƒS(s6asarray(data, dtype) = array(data, dtype, copy=0) R/RMiN(RR.RR/RN(R›R/((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyR•Hs(t MethodTypecCst|dtƒS(N(R–RRN(Rt((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_mTscOs td‚dS(Ns'not yet implemented for numpy.ma arrays(tNotImplementedError(R tkwds((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytnot_implementedVscOs t||ƒS(N(R(R R R™((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_choose`sc CsUtd|iit|ƒit|ƒiƒdt|itt|ƒt|ƒƒƒƒS(NR›Rœ(R.R›tclipR•R[RœRC(R ta_minta_maxRn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_clipes cCst|||ƒS(N(R(R tcondR‘Rn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _compresslscCsZ|i}|tj otii|i|ƒ}ntd|idƒi||ƒd|ƒS(NR›iRœ(RœRBR1RXR”R.RPtcumprod(R R‘R/RnRF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_cumprodts  cCsZ|i}|tj otii|i|ƒ}ntd|idƒi||ƒd|ƒS(NR›iRœ(RœRBR1RXR”R.RPtcumsum(R R‘R/RnRF((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_cumsum{s  cCsI|dj otdƒ‚n|djo t|ƒSti||ƒSdS(Ns+Output arrays Unsupported for masked arrays(RR:RŒR(R?R‘Rn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_maxŠs    cCsI|dj otdƒ‚n|djo t|ƒSti||ƒSdS(Ns+Output arrays Unsupported for masked arrays(RR:RŠR(R?R‘Rn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_min“s    cCs|i||ƒ|i|ƒS(N(tmaxtmin(R?R‘Rn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_ptp¢scCsTy+td|iiƒd|iiƒƒ}Wn"tj ot|dƒ}nX|S(NR›Rœtsqueeze(R.R›R«RœRJt_wrapit(R RU((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_squeeze¬s cCs1td|ii||ƒd|ii||ƒƒS(NR›Rœ(R.R›R€Rœ(R RR‚((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt _swapaxes·scCsˆ|djoti|iƒƒiƒS|i|dƒ}||iddƒ}||9}||iddƒ:}|id|ƒi|ƒS(NiR‘( RR7R•RtvarR€tmeanRR(R R‘R/RnR?((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_varÀs  cCs|i||ƒdS(Ngà?(R¯(R R‘R/Rn((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt_stdÈsRœcCst|||ƒS(N(RS(R?RPR‘((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyROÔscCst||||ƒS(N(Rd(R?R‘R_R`((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pytaverage×scCst|||ƒS(N(R„(R?RUR‘((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyRƒÚs(ãRR$R8tnumpy.core.umathtcoreR1tnumpy.core.fromnumericRRtnumpy.core.numericRRRRRtnumpy.core.numerictypesRRR7R¨RxRyRIRBR€t ExceptionRRRÁR&R-R+R)R3R6R;R<R>RCRGRLRVRDR[RRPRaReRfRlRoRrRsR~RRŠRýtlogtlog10texpt conjugatetsinRnttantarcsintarccostarctantarcsinhtarccoshtarctanhtsinhtcoshttanhR2tfabsRàR–tround_taroundtfloortceilt logical_notRãRìRRïRñRôR÷Råtfmodthypottarctan2RRRqR RjRiRŽR@RXRSt logical_xorRØRÛRÝR—RER™R˜R.RNR=tisarraytisMARCRDRIRGR"t arrayrangeRMRèRêR RSRTRURÊRWRRRúR{RRZRdR…RRrRsRtRuRvRwRxR{R|RŒRTR~R€R„R5R&R'R†RR’R!RRRˆRŠR‹RŒRRŽR‘R’R“R”R•R•R–R—RštallR§R=tbasetbyteswapR›RŸRœR¡tconjRMR£R¢R¥R¤tdumptdumpsRuR^tflattentgetfieldR¦R¨R§R©R°tnbytest newbyteordertprodRªtptpROt searchsortedtsetfieldtsetflagsR­R«tstridesR®RƒttofileR±R²R¯tstdtviewtroundRÍRR³(((s9/usr/lib64/python2.6/site-packages/numpy/oldnumeric/ma.pyt sÞ                 #E&             ÿÿÿ}           k              / /