;ò cöz?c@s# dZdkiZdklZlZlZdkiZ dk i Z dk Z dk Z dkZdkZdkTdkZe iZdZdefd„ƒYZdfd„ƒYZed ƒZeid gƒieƒZeid d gƒieƒZd Z eid gƒie!ƒ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„Z0d„Z1d„Z2dfd„ƒYZ3dfd„ƒYZ4dfd „ƒYZ5d!fd"„ƒYZ6d#fd$„ƒYZ7d%fd&„ƒYZ8d'fd(„ƒYZ9d)fd*„ƒYZ:e7ei;d+e6d+ƒƒZ;e7ei<d,e5d+ƒƒZ<e7ei=d,e5d+ƒƒZ=e7ei>ƒZ>e7ei?ƒZ?e7ei@ƒZ@e7eiAƒZAe7eiBd+e4dƒƒZBe7eiCd+e3d-d,ƒƒZCe7eiDd+e3d-d,ƒƒZDe7eiEƒZEe7eiFƒZFe7eiGƒZGe7eiHƒZHe7eiIƒZIe7eiJƒZJe7eiKƒZKe7eiLƒZLe7eiMƒZMe7eiNƒZNe7eiOƒZOe7eiPƒZPe7eiQd.ƒZQe7eiRƒZRe:eiSƒZSe:eiTƒZTe/eT_Ue:eiVd.d.ƒZVe9eiWe8ƒd d.ƒZWye9eiXe8ƒd d.ƒZXWnnXye9eiYe8ƒd d.ƒZYWnnXe9eiZe8ƒd d.ƒZZe9ei[e8ƒd d.ƒZ[e:ei\ƒZ\e:ei]d+d,ƒZ]e/e]_Ue:ei^ƒZ^e/e^_Ue:ei_ƒZ_e/e__Ue:ei`ƒZ`e/e`_Ue:eiaƒZae/ea_Ue:eibƒZbe/eb_Ue:eicƒZce/ec_Ue:eidƒZde:eieƒZee:eifƒZfe:eigƒZge:eihƒZhe:eiiƒZieijZjeikZkeilZld/e imfd0„ƒYZnenZd1enfd2„ƒYZoeoƒZpd3„ZqeqZreqZsd.d4d5d6„Ztd.d7„Zud4d5d.d d8„Zvd.d d9„Zwe/d.e/d:„ZxexZyd;„Zzd<„Z{d=„Z|d>„Z}d d?„Z~d@„Ze/dA„Z€ed dB„Z‚ed dC„Zƒe/dD„Z„e/dE„Z…e/e/dF„Z†eSiUZ‡eViUZˆd e/d dG„Z‰dH„ZŠdI„Z‹d.dJ„ZŒd.dK„Zd.dL„ZŽd.dM„Zd.dN„Zd.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[„Zd\„Zžd]d^„ZŸd_fd`„ƒYZ e ƒZ¡dafdb„ƒ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ªdjZ«dk„Zke/dl„Zld dm„Z¬dn„Z­e¬e«ƒdS(osþ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(sespisNewAxis(s*f1e-35sMAErrorcBs#tZed„Zd„ZeZRS(NcCs ||_dS(sCreate an exceptionN(sargssself(sselfsargs((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__scCst|iƒSdS(s#Calculate the string representationN(sstrsselfsargs(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__str__s(s__name__s __module__sNones__init__s__str__s__repr__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysMAErrors  s_MaskedPrintOptioncBsGtZdZd„Zd„Zd„Zd„Zdd„Zd„ZRS(s<One instance of this class, masked_print_option, is created.cCs|i|ƒd|_dS(s&Create the masked print option object.iN(sselfs set_displaysdisplays_enabled(sselfsdisplay((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__$s cCs |iSdS(s#Show what prints for masked values.N(sselfs_display(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdisplay)scCs ||_dS(s2set_display(s) sets what prints for masked values.N(sssselfs_display(sselfss((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys set_display-scCs |iSdS(s(Is the use of the display value enabled?N(sselfs_enabled(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysenabled1sicCs ||_dS(sSet the enabling flag to flag.N(sflagsselfs_enabled(sselfsflag((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysenable5scCst|iƒSdS(N(sstrsselfs_display(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__str__9s( s__name__s __module__s__doc__s__init__sdisplays set_displaysenabledsenables__str__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys_MaskedPrintOption"s      s--f1e+20x0.00.0s?icCs+t|tiƒotSn t|tiƒpt|tiƒotSnßt|tiƒot SnÄt|ti ƒot Sn©t|t ƒpt|t iƒo|iƒ}|tdjotSn|tdjotSn|tdjot Sn|tdjot itƒSntSntSdS(s7Function to calculate default fill value for an object.sFloatsIntegersComplexsUnsignedIntegerN(s isinstancesobjstypess FloatTypesdefault_real_fill_valuesIntTypesLongTypesdefault_integer_fill_values StringTypesdefault_character_fill_values ComplexTypesdefault_complex_fill_values MaskedArraysNumerics arraytypestypecodesxs typecodessabsolutesdefault_object_fill_value(sobjsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdefault_fill_valueFs*&# cCsÝt|tiƒotSn¿t|tiƒpt|tiƒo tiSnŽt|t ƒpt|t i ƒoa|i ƒ}|tdjotSn|tdjo tiSn|tdjo tiSqÙn td‚dS(sDFunction to calculate default fill value suitable for taking minima.sFloatsIntegersUnsignedIntegers(Unsuitable type for calculating minimum.N(s isinstancesobjstypess FloatTypesdefault_real_fill_valuesIntTypesLongTypessyssmaxints MaskedArraysNumerics arraytypestypecodesxs typecodess TypeError(sobjsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysminimum_fill_value`s& #  cCsÛt|tiƒo t Sn¼t|tiƒpt|tiƒo t Snt|tƒpt|t i ƒo`|i ƒ}|t djo t Sn|t djo ti Sn|t djodSq×n td‚dS(sDFunction to calculate default fill value suitable for taking maxima.sFloatsIntegersUnsignedIntegeris(Unsuitable type for calculating maximum.N(s isinstancesobjstypess FloatTypesdefault_real_fill_valuesIntTypesLongTypesdefault_integer_fill_values MaskedArraysNumerics arraytypestypecodesxs typecodesssyssmaxints TypeError(sobjsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmaximum_fill_valueqs & #    cCs"t|ƒo|i|ƒndS(s,Set fill value of a if it is a masked array.N(s isMaskedArraysasset_fill_values fill_value(sas fill_value((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysset_fill_value‚s cCs&t|tƒo|iƒSntSdS(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(s isinstancesas MaskedArraysraw_masksNone(sa((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysgetmask‡scCs5t|ƒ}|tjott|ƒƒSn|SdS(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(sgetmasksasmsNonesmake_mask_nonesshape(sasm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys getmaskarrays   cCsC|tjp#t|tiƒo|iƒtjodSndSdS(s;Is m a legal mask? Does not check contents, only type. iiN(smsNones isinstancesNumerics ArrayTypestypecodesMaskType(sm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysis_mask›s3cCsß|tjotSnšt|tiƒoe|iƒtjo?|oti|ddƒ}q‰ti|ddddƒ}q®|i tƒ}n"tit |dƒtddƒ}|oti ti|ƒƒ otSn|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 None. Does not check that contents must be 0s and 1s. if flag, return None if m contains no true elements. s savespaceiscopyiN(smsNones isinstancesNumerics ArrayTypestypecodesMaskTypescopysarraysresultsastypesfilledsflagssometruesravel(smscopysflagsresult((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys make_mask¤s !!cCs#ti|tƒ}||_|SdS(s&Return a mask of all zeros of shape s.N(sNumericszerossssMaskTypesresultsshape(sssresult((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmake_mask_none½s cCsr|tjot|ƒSn|tjot|ƒSn||jo t|ƒo|Sntti||ƒƒSdS(sLogical or of the mask candidates m1 and m2, treating None as false. Result may equal m1 or m2 if the other is None. N(sm1sNones make_masksm2sis_masksNumerics logical_or(sm1sm2((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmask_orÅs  cCst|tƒo|i|ƒSn]t|tiƒo |iƒo|Sn5t|ti ƒoti |dƒSnti |ƒSdS(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. sON( s isinstancesas MaskedArraysfilledsvaluesNumerics ArrayTypes iscontiguousstypessDictTypesarray(sasvalue((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysfilledÎs  cCs1t|ƒo|iƒ}n t|ƒ}|SdS(sb The fill value of a, if it has one; otherwise, the default fill value for that type. N(s isMaskedArraysas fill_valuesresultsdefault_fill_value(sasresult((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys fill_valueâs   cCs5t|ƒ}t|ƒ}||jo|SntSdS(s:The common fill_value of a and b, if there is one, or NoneN(s fill_valuesast1sbst2sNone(sasbst2st1((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyscommon_fill_valueís    sdomain_check_intervalcBs tZdZd„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(sy1sselfsy2(sselfsy1sy2((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__÷s cCs2titi||iƒti||iƒƒSdS(sExecute the call behavior.N(sNumerics logical_orsgreatersxsselfsy2slesssy1(sselfsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__üs(s__name__s __module__s__doc__s__init__s__call__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdomain_check_intervalõs  s domain_tancBs tZdZd„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(sepssself(sselfseps((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__scCs)tititi|ƒƒ|iƒSdS(sExecute the call behavior.N(sNumericslesssabsolutescossxsselfseps(sselfsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__s(s__name__s __module__s__doc__s__init__s__call__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys domain_tans  sdomain_greatercBs tZdZd„Zd„ZRS(s(domain_greater(v)(x) = true where x <= vcCs ||_dS(s(domain_greater(v)(x) = true where x <= vN(scritical_valuesself(sselfscritical_value((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__scCsti||iƒSdS(sExecute the call behavior.N(sNumerics less_equalsxsselfscritical_value(sselfsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__s(s__name__s __module__s__doc__s__init__s__call__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdomain_greater s  sdomain_greater_equalcBs tZdZd„Zd„ZRS(s-domain_greater_equal(v)(x) = true where x < vcCs ||_dS(s-domain_greater_equal(v)(x) = true where x < vN(scritical_valuesself(sselfscritical_value((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__scCsti||iƒSdS(sExecute the call behavior.N(sNumericslesssxsselfscritical_value(sselfsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__s(s__name__s __module__s__doc__s__init__s__call__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdomain_greater_equals  smasked_unary_operationcBs)tZded„Zd„Zd„ZRS(NicCs:||_||_||_t|dt|ƒƒ|_dS(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. s__doc__N(saufuncsselfsfsfillsdomainsgetattrsstrs__doc__(sselfsaufuncsfillsdomain((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__!s    cOst|ƒ}t||iƒ}|itj ot ||i|ƒƒ}n|tjoA|i |||Ž}t|ƒtijot|ƒSq|Sn}t||ƒ}|i t||iƒ||Ž}t|ƒtijot||ƒSn&titi|ƒƒotSn|SdS(sExecute the call behavior.N(sgetmasksasmsfilledsselfsfillsd1sdomainsNonesmask_orsfsargsskwargssresultstypesNumerics ArrayTypes masked_arraysdxssometruesravelsmasked(sselfsasargsskwargssmsresultsdxsd1((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__,s"  !cCsdt|iƒSdS(NsMasked version of (sstrsselfsf(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__str__Cs(s__name__s __module__sNones__init__s__call__s__str__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_unary_operation s sdomain_safe_dividecBstZed„Zd„ZRS(NcCs ||_dS(N(s tolerancesself(sselfs tolerance((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__HscCs'ti|ƒ|iti|ƒjSdS(N(sNumericsabsolutesasselfs tolerancesb(sselfsasb((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__Js(s__name__s __module__sdivide_tolerances__init__s__call__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdomain_safe_divideGs sdomained_binary_operationcBs/tZdZddd„Zd„Zd„ZRS(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. icCsC||_||_||_||_t|dt|ƒƒ|_dS(sjabfunc(fillx, filly) must be defined. abfunc(x, filly) = x for all x to enable reduce. s__doc__N( sabfuncsselfsfsdomainsfillxsfillysgetattrsstrs__doc__(sselfsabfuncsdomainsfillxsfilly((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__Qs     c Cs}t|ƒ}t|ƒ}t||iƒ} t||i ƒ}|i | |ƒ}t it i|ƒƒo(t||i |ƒ}t||ƒ}nt||ƒ}|tjo>|i| |ƒ}t|ƒt ijot|ƒSqé|Sn|i| |ƒ}t|ƒt ijoC|i|ijott|ƒt|ƒƒ}nt||ƒSn&t it i|ƒƒotSn|SdS(sExecute the call behavior.N(sgetmasksasmasbsmbsfilledsselfsfillxsd1sfillysd2sdomainstsNumericssometruesravelswheresmask_orsmsNonesfsresultstypes ArrayTypes masked_arraysshapes getmaskarraysmasked( sselfsasbsmasresultsmbsmstsd2sd1((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__[s.   cCsdt|iƒSdS(NsMasked version of (sstrsselfsf(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__str__vs(s__name__s __module__s__doc__s__init__s__call__s__str__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdomained_binary_operationMs  smasked_binary_operationcBsJtZddd„Zd„Zdd„Zd„Zdd„Zd„ZRS(NicCs:||_||_||_t|dt|ƒƒ|_dS(sjabfunc(fillx, filly) must be defined. abfunc(x, filly) = x for all x to enable reduce. s__doc__N(sabfuncsselfsfsfillxsfillysgetattrsstrs__doc__(sselfsabfuncsfillxsfilly((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__zs    c OsNtt|ƒt|ƒƒ}|tjoht||iƒ}t||i ƒ}|i ||||Ž}t|ƒtijot|ƒSq|Snt||iƒ}t||i ƒ}|i ||||Ž}t|ƒtijoC|i|ijott|ƒt|ƒƒ}nt||ƒSn&titi|ƒƒotSn|SdS(sExecute the call behavior.N(smask_orsgetmasksasbsmsNonesfilledsselfsfillxsd1sfillysd2sfsargsskwargssresultstypesNumerics ArrayTypes masked_arraysshapes getmaskarrayssometruesravelsmasked( sselfsasbsargsskwargssmsresultsd2sd1((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__ƒs& cCs t|ƒ}t||iƒ}|ifjo?df|_|tj o"t |ddƒ}df|_qmn|tjot |i i ||ƒƒSn†t ||ƒ}|i i t||iƒ|ƒ}tii ||ƒ}t|tiƒot ||t|ƒƒSn|otSn|SdS(s6Reduce target along the given axis with this function.iscopyN(sgetmaskstargetsmsfilledsselfsfillystsshapesNones make_masks masked_arraysfsreducesaxissNumerics logical_ands isinstances ArrayTypes fill_valuesmasked(sselfstargetsaxissmst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysreducešs$    !cCs¤t|ƒ}t|ƒ}|tjo |tjo t}n+t|ƒ}t|ƒ}ti ||ƒ}|i i t ||i ƒt ||iƒƒ}t||ƒSdS(s<Return the function applied to the outer product of a and b.N(sgetmasksasmasbsmbsNonesms getmaskarrays logical_orsoutersselfsfsfilledsfillxsfillysds masked_array(sselfsasbsmasdsmbsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysouter°s     -cCs/t||iƒ}t|ii||ƒƒSdS(s=Accumulate target along axis after filling with y fill value.N( sfilledstargetsselfsfillysts masked_arraysfs accumulatesaxis(sselfstargetsaxisst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys accumulate½scCsdt|iƒSdS(NsMasked version of (sstrsselfsf(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__str__Ás(s__name__s __module__s__init__s__call__sreducesouters accumulates__str__(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_binary_operationys    f0.0f1.0f-1.0is MaskedArraycBswtZdZdZdZedeeed„Zed„Zd„Zd„Z d„Z d „Z d „Z d „Z d „Zd „Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„ZeZeZeZd„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)„Z0d*„Z1d+„Z2ed,„Z3d-„Z4d.„Z5d/„Z6d0„Z7d1„Z8d2„Z9d3„Z:d4„Z;d5„Z<d6„Z=d7„Z>d8„Z?d9„Z@d:„ZAed;„ZBd<„ZCd=„ZDed>„ZEd?„ZFd@„ZGdA„ZHdB„ZIdC„ZJdD„ZKdE„ZLdF„ZMdG„ZNdH„ZOdI„ZPedJ„ZQedK„ZRdL„ZOdM„ZSedN„ZTedO„ZUdP„ZVdQ„ZWeXe e dRdSƒZYeXe e dRdTƒZZeXe edRdUƒZ[eXeedRdVƒZ\e\Z]RS(Ws€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. f1.0sMA.MaskedArrayicCs×|} |} |}t|tƒo¾|i ƒ} | iƒ}| t jo |} n| |jo d}n| i ƒ} | t jo | } n| | jo d}nd} |t jo|iƒ}q¦|t j o d}q¦nÇt|tiƒo~|} | iƒ}| t jo |} n| |jo d}n| i ƒ} | t jo | } q¦| | jo d}q¦n6d}| t jo d} nti|| d| ƒ} |oE| |jo"ti| ddd| ƒ|_qû| i| ƒ|_n | |_|t jot |_d|_n©t|tƒ oV|djoti|iitƒ}q…|djoti|iitƒ}q…nt|d| ƒ|_|it jo d|_n|i|j|_t|iƒ}t|iƒ} || jo†|djo+ti |i|iiƒ|_d|_qÆ| djo4ti |i|iiƒ|_|ii|i_qÆt!d‚nM|djot|iƒt|iƒjo |i"ƒ|ii|i_n|i#|ƒdS(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. iis savespacescopysflagsMask and data not compatible.N(%stypecodestcs savespacesssscopysneed_data_copieds isinstancesdatas MaskedArraysraw_datascsctcsNones spacesaverscsssmasksNumerics ArrayTypesarraysselfs_datasastypes_masks _shared_masksonessshapesMaskTypeszeross make_maskssizesnmsndsresizesMAErrors unshare_masksset_fill_values fill_value(sselfsdatastypecodescopys savespacesmasks fill_valuesneed_data_copiedsnmsndstcscssscssssctc((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__/sˆ                         "             , cCsu|itj o?titi|iƒƒo td‚qOt|_d|_n|o|i i |ƒSn|i SdS(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. iN( sselfs_masksNonesNumericssometruesravelsMAErrors _shared_masksts_datasastype(sselfst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys __array__‡s   cCsAdGH|iiƒdGH|itj o|iiƒndGHdS(s)prints out information about masked arraysMaskedArray data info:sMaskedArray mask info:sno maskN(sselfs_datasinfos_masksNone(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysinfo˜s cCs|iiSdS(sReturn the current shape.N(sselfs_datasshape(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys _get_shape¢scCs²|iiƒ o4ti|i|iiƒd|iiƒƒ|_n||i_|i t j oM|i ƒ|i iƒ o"ti|i t ddƒ|_ n||i _ndS(sSet the array's shape.iN( sselfs_datas iscontiguoussNumericsarraystypecodes spacesaversnewshapesshapes_masksNones unshare_masksMaskType(sselfsnewshape((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys _set_shape¦s  "cCsb|itjo&t|iidtd|iƒƒSn)t|iid|iid|iƒƒSdS(s#Calculate the flat value. smasks fill_valueN(sselfs_masksNones masked_arrays_datasflats fill_value(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys _get_flat²s  cCs|i}||(dS(sx.flat = valueN(sselfsflatsysvalue(sselfsvaluesy((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys _set_flat½s cCse|itjo&t|iidtd|iƒƒSn,t|iid|iƒid|iƒƒSdS(s%Get the real part of a complex array.smasks fill_valueN( sselfs_masksNones masked_arrays_datasreals fill_valuesmasksflat(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys _get_realÂs cCs|i}||d¿sN(sselfs_masksms_datasshapessslenslssNonesaxissreducesnsliststsonessNumericsravelsastypesIntswssizesn1sn2saddssum( sselfsaxissmsnssslsswstsn1sn2((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyscount´s2              cCst||ƒSdS(s%s.dot(other) = innerproduct(s, other)N(s innerproductsselfsother(sselfsother((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdotÒscCs |iSdS(sGet the current fill value.N(sselfs _fill_value(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys fill_valueÖscCs |i}|i}|tjo8|iƒo|SqWti|d|i ƒddƒSn|}|tjo |i }n|t jo:t|tƒ oti|ƒ}n|}|i|_n[y8ti|d|i ƒddƒ}ti|||ƒWnti|||fƒ}nX|SdS(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. stypecodescopyiN(sselfs_datasds_masksmsNones iscontiguoussNumericsarraystypecodes fill_valuesvalues _fill_valuesmaskeds isinstances_MaskedPrintOptionsresultsshapesputmasksobjschoose(sselfs fill_valuesresultsdsmsvalue((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysfilledÚs*    #   !cCs t|iƒt|iƒfSdS(s)Return the ids of the data and mask areasN(sidsselfs_datas_mask(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysidsùscCs|iiƒSdS(sIs the data contiguous?N(sselfs_datas iscontiguous(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys iscontiguousýscCs|iiƒSdS(sItem size of each data item.N(sselfs_datasitemsize(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysitemsizescCsG|i}|tjo|Sn&|iƒo|Snti|iƒSdS(s1Return the data mask, or None. Result contiguous.N(sselfs_masksmsNones iscontiguoussNumericsarray(sselfsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasks   cCst||ƒSdS(s's.outer(other) = outerproduct(s, other)N(s outerproductsselfsother(sselfsother((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysouterscCsŽti|iƒƒ}|itjo |}ntid|i|ƒ}t |ƒt|ƒjo t d‚nti |i |t|ƒƒdS(s\Set the non-masked entries of self to filled(values). No change to mask is(x.put(values) incorrect count of values.N(sNumericsarangesselfssizesiotas_masksNonesindscompressslensvaluessMAErrorsputs_datasfilled(sselfsvaluessindsiota((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysputs  cCs–|itj o‚ti|iƒƒ}ti|i|ƒ}t |ƒt|ƒjo t d‚nti |i |t|ƒƒt|_d|_ndS(s\Set the masked entries of self to filled(values). Mask changed to None. s(x.put(values) incorrect count of values.iN(sselfs_masksNonesNumericsarangessizesiotascompresssindslensvaluessMAErrorsputs_datasfilleds _shared_mask(sselfsvaluessindsiota((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysputmask s  cCs |iSdS(s_ The raw data; portions may be meaningless. May be noncontiguous. Expert use only.N(sselfs_data(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysraw_data-scCs |iSdS(sh The raw mask; portions may be meaningless. May be noncontiguous. Expert use only. N(sselfs_mask(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysraw_mask2scCsdSdS(sDeprecated. returns 0.iN((sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys spacesaver8scCsdS(s Deprecated.N((sselfsvalue((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys savespace<scCs0|tjot|iƒƒ}n||_dS(s3Set the fill value to v. Omit v to restore default.N(svsNonesdefault_fill_valuesselfsraw_datas _fill_value(sselfsv((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysset_fill_value@s cCsT|ii}|tjo/t|ƒdjodSqPtd„|ƒSn ||SdS(s5Number of elements in array, or in a particular axis.iicCs||S(N(sxsy(sxsy((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysMsN(sselfs_datasshapesssaxissNoneslensreduce(sselfsaxisss((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyssizeFs  cCs|iiƒSdS(s.spacesaver() queries the spacesaver attribute.N(sselfs_datas spacesaver(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys spacesaverQscCs|iiƒSdS(N(sselfs_datastypecode(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pystypecodeUscCs|i|ƒiƒSdS(sConvert to listN(sselfsfilleds fill_valuestolist(sselfs fill_value((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pystolistXscCs|i|ƒiƒSdS(sConvert to stringN(sselfsfilleds fill_valuestostring(sselfs fill_value((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pystostring\scCsT|itjodSnt|iddƒ}|tjot|_d|_ndS(s%Replace the mask by None if possible.Nsflagii(sselfs_masksNones make_masksms _shared_mask(sselfsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysunmask`s  cCs9|io+t|iddddƒ|_d|_ndS(s'If currently sharing mask, make a copy.scopyisflagiN(sselfs _shared_masks make_masks_mask(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys unshare_maskhs sdocs#tuple giving the shape of the arraysAccess array in flat form.s!Access the real part of the arrays&Access the imaginary part of the array(^s__name__s __module__s__doc__s op_priorityshandler_cache_keysNones__init__s __array__sinfos _get_shapes _set_shapes _get_flats _set_flats _get_reals _set_reals_get_imaginarys_set_imaginarys__str__s__repr__s __float__s__int__s __getitem__s __getslice__s __setitem__s __setslice__s__len__s__and__s__or__s__xor__s__rand__s__ror__s__rxor__s__abs__s__neg__s__pos__s__add__s__radd__s__mod__s__rmod__s __lshift__s __rshift__s__sub__s__rsub__s__mul__s__rmul__s__div__s__rdiv__s __truediv__s __rtruediv__s __floordiv__s __rfloordiv__s__pow__s__sqrt__s__iadd__s__imul__s__isub__s__idiv__s__eq__s__ne__s__lt__s__le__s__gt__s__ge__sastypes byte_swappeds compressedscountsdots fill_valuesfilledsidss iscontiguoussitemsizesmasksoutersputsputmasksraw_datasraw_masks spacesavers savespacesset_fill_valuessizestypecodestoliststostringsunmasks unshare_maskspropertysshapesflatsreals imaginarysimag(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys MaskedArray s¾ X      /                              + + + '                                   s MaskedScalarcBs tZd„ZeeiƒZRS(NcCs8ti|dgddgƒf|i_f|i_dS(Nismaski(s MaskedArrays__init__sselfs_datasshapes_mask(sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__s (s__name__s __module__s__init__spropertys MaskedArrays _get_shapesshape(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys MaskedScalar~s cCst|tƒSdS(s9Is x a masked array, that is, an instance of MaskedArray?N(s isinstancesxs MaskedArray(sx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys isMaskedArrayˆsf1.0000000000000001e-05f1e-08c CsÑtt|ƒt|ƒƒ}t|ƒ} t|ƒ} tt| ddd|ƒ|ƒi t ƒ}tt| ddd|ƒdƒi t ƒ}titi||ƒ||ti|ƒƒ}titi|ƒƒSdS(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. scopyismaskiN(smask_orsgetmasksasbsmsfilledsd1sd2sarrays fill_valuesastypesFloatsxsysNumerics less_equalsabsolutesatolsrtolsdsalltruesravel( sasbs fill_valuesrtolsatolsdsmsysxsd2sd1((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysallcloses  **0cCsàtt|ƒt|ƒƒ}|tjoDt|ƒ}t|ƒ}t i ||ƒ}t i t i |ƒƒSnq|oet|ƒ}t|ƒ}t i ||ƒ}t|d|ddƒ}t i t i t|dƒƒƒSndSdS(s„ True if all entries of a and b are equal, using fill_value as a truth value where either or both are masked. smaskscopyiiN(smask_orsgetmasksasbsmsNonesfilledsxsysNumericsequalsdsalltruesravels fill_valuesarraysdm(sasbs fill_valuesdsmsdmsysx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysallequal s     #c Cs¨ti}t||ƒ}|iƒtdjobti |||ƒ||||ƒƒ}t |ddƒ}t|d|d|d|d|ƒSnt||||ƒSdS( 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) sFloatsflagismasks savespacescopys fill_valueN(sNumericsabsolutesabssfilledsdatasvaluesdstypecodes typecodess less_equalsatolsrtolsms make_masksarrays savespacescopys masked_object( sdatasvaluesrtolsatolscopys savespacesabssmsd((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_values´s  *c CsSt||ƒ}tti||ƒddƒ}t|d|d|d|d|ƒSdS(s5Create array masked where exactly data equal to valuesflagismaskscopys savespaces fill_valueN( sfilledsdatasvaluesds make_masksNumericsequalsdmsarrayscopys savespace(sdatasvaluescopys savespacesdmsd((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_objectÊs cCs tti||||ƒƒSdS(sqJust like range() except it returns a array whose type can be specfied by the keyword argument typecode. N(sarraysNumerics arrayrangesstartsstopsstepstypecode(sstartsstopsstepstypecode((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys arrayrangeÑscCstti||ƒƒSdS(sAConstruct a masked array from a string. Result will have no mask.N(s masked_arraysNumerics fromstringssst(ssst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys fromstringÙscCskt|ƒ}|tjo&tit|ƒ|ƒ}t |ƒSn)tit|dƒ|ƒ}t ||ƒSdS(sLeft shift n bitsiN( sgetmasksasmsNonesNumerics left_shiftsfilledsnsds masked_array(sasnsdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys left_shiftÝs  cCskt|ƒ}|tjo&tit|ƒ|ƒ}t |ƒSn)tit|dƒ|ƒ}t ||ƒSdS(sRight shift n bitsiN( sgetmasksasmsNonesNumerics right_shiftsfilledsnsds masked_array(sasnsdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys right_shiftçs  cCsnt|ƒ}|tj oti||ƒ}nttit|ƒ|ƒd|ƒ}|i t |ƒƒ|SdS(swresize(a, new_shape) returns a new array with the specified shape. The original array's total size can be any size.smaskN( sgetmasksasmsNonesNumericsresizes new_shapesarraysfilledsresultsset_fill_values fill_value(sas new_shapesmsresult((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysresizeñs  $cCs±t|ƒ}t|tiƒo!t|gt|ƒ|ƒ}nt |ƒ}|t j ot i|||ƒ}nt i|||ƒ}t||ƒ}|it|ƒƒ|SdS(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(sfilledsasafs isinstancesrepeatsstypessIntTypestuplesshapesaxissgetmasksmsNonesNumericsrepeatsds masked_arraysresultsset_fill_values fill_value(sasrepeatssaxissdsresultsafsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysrepeatûs !  cCstti|ƒƒSdS(s<identity(n) returns the identity matrix of shape n x n. N(sarraysNumericsidentitysn(sn((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysidentity scCstti||ƒƒSdS(s„indices(dimensions,typecode=None) returns an array representing a grid of indices with row-only, and column-only variation. N(sarraysNumericsindicess dimensionsstypecode(s dimensionsstypecode((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysindicesscCstti|||ƒƒSdS(s_zeros(n, typecode=Int, savespace=0) = an array of all zeros of the given length or shape.N(sarraysNumericszerossshapestypecodes savespace(sshapestypecodes savespace((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyszerosscCstti|||ƒƒSdS(s]ones(n, typecode=Int, savespace=0) = an array of all ones of the given length or shape.N(sarraysNumericsonessshapestypecodes savespace(sshapestypecodes savespace((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysonesscCst|ƒ}|i|ƒSdS(s?Count of the non-masked elements in a, or along a certain axis.N(s masked_arraysascountsaxis(sasaxis((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyscount"s c Cs|tj o td‚nt|ƒ}t|ƒ}t||ƒ}t |dƒ}t |dƒ}|i ƒtdjotti||ƒ|ƒSntti|dƒddƒ}t||ƒ}|tjotti||ƒƒSn/ti|d|ƒ}tti||ƒ|ƒSdS(sa**bs3-argument power not supported.isIntegerisflagN(sthirdsNonesMAErrorsgetmasksasmasbsmbsmask_orsmsfilledsfasfbstypecodes typecodess masked_arraysNumericspowers make_masks less_equalsmdswhere( sasbsthirdsmdsmasmbsmsfasfb((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyspower's      cCsZ|tjo-t|ƒo|iƒ}q:t|ƒ}nt|d|ddd|ƒSdS(s~masked_array(a, mask=None) = array(a, mask=mask, copy=0, fill_value=fill_value) Use fill_value(a) if None. smaskscopyis fill_valueN(s fill_valuesNones isMaskedArraysasdefault_fill_valuesarraysmask(sasmasks fill_value((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_array:s   c Cst|ƒ}|iƒ} |i}|fjo df}n|tjo8| tjoˆ|tjo1ti |i ƒi ƒ} t d„|dƒ}q|t|dƒi } tii|i ƒi | ƒ} tii| ƒ}~ qk|tjoFti |i ƒ} ti| ddfƒi } tii| ƒ}~ qktt|dƒtid| ƒi } ti |i | ƒ} ti | ƒ}~ nì| tjof|tjo-||d}tii|i ƒ|ƒ} qkt|dƒ} | i}|fjo df}n||joGti| tddƒ} ti || |ƒ} ti | |ƒ}~ qk|||fjo…||}tgt|ƒ} tttdƒ| |esf1.0f0.0smaskscopyisw[s] * ones(ash, Float)saverage: weights wrong shape.s(] * masked_array(ones(ash, Float), mask)N(&s masked_arraysasmasksshapesashsaxissNonesweightssaddsreducesraw_datasflatsnsdsfilledswsNumericszreduceschoosesarraysFloatswshsnisNewAxisslensrsslicesevalsreprstuples ValueErrorsdividesresults isinstances MaskedArraysunmasksreturnedsones( sasaxissweightssreturnedsresultsnisashsdswshsmasksnsrsw((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysaverageMs¢         $                   c Cstt|dƒdƒ}|tjod} d}n0t|ƒ} t|ƒ}|t jo d}n|tjod} d}n0t|ƒ} t|ƒ}|t jo d}nt i|| | fƒ}t i|||fƒ}t|ƒ}tt||ƒddddƒ}t||ƒSdS(swhere(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. iiscopysflagN(sfilleds not_equals conditionsfcsxsmaskedsxvsxmsgetmasksNonesysyvsymsNumericschoosesdsmdsms make_masksmask_ors masked_array( s conditionsxsysmdsymsdsxmsmsfcsyvsxv((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyswhere¹s*           !c Cs×d„} d„}t|dƒ}gi}|D]} ||| ƒƒq/~}gi}|D]} || | ƒƒq\~}t i ||ƒ}t i ||ƒ}tt|t|ƒƒddddƒ}t||ƒSdS(s=Returns array shaped like indices with elements chosen from tcCs#|tjodSnt|ƒSdS(Ni(sxsmaskedsfilled(sx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysfmask×s cCs>|tjodSnt|ƒ}|tjodSn|SdS(Nii(sxsmaskedsgetmasksmsNone(sxsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysnmaskÚs    iscopysflagiN(sfmasksnmasksfilledsindicesscsappends_[1]stsxsmaskssasNumericschoosesdsms make_masksmask_orsgetmasks masked_array( sindicesstsascsnmasksdsmsmaskss_[1]sxsfmask((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyschooseÕs  --'cCsDt|dƒ}tt|ƒ|ƒ}tt|ƒd|d|ƒSdS(sjReturn x as an array masked where condition is true. Also masked where x or condition masked. iscopysmaskN( sfilleds conditionscmsmask_orsgetmasksxsmsarrayscopy(s conditionsxscopyscmsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_whereçscCstt||ƒ||ƒSdS(s3masked_greater(x, value) = x masked where x > valueN(s masked_wheresgreatersxsvaluescopy(sxsvaluescopy((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_greaterïscCstt||ƒ||ƒSdS(s:masked_greater_equal(x, value) = x masked where x >= valueN(s masked_wheres greater_equalsxsvaluescopy(sxsvaluescopy((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_greater_equalóscCstt||ƒ||ƒSdS(s0masked_less(x, value) = x masked where x < valueN(s masked_whereslesssxsvaluescopy(sxsvaluescopy((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_less÷scCstt||ƒ||ƒSdS(s7masked_less_equal(x, value) = x masked where x <= valueN(s masked_wheres less_equalsxsvaluescopy(sxsvaluescopy((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_less_equalûscCsPt|dƒ}ti||ƒ}t|t|ƒƒ}t |d|d|ƒSdS(s6masked_not_equal(x, value) = x masked where x != valueismaskscopyN( sfilledsxsdsNumerics not_equalsvaluescsmask_orsgetmasksmsarrayscopy(sxsvaluescopyscsdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_not_equalÿs cCsPt|dƒ}ti||ƒ}t|t|ƒƒ}t |d|d|ƒSdS(s{masked_equal(x, value) = x masked where x == value For floating point consider masked_values(x, value) instead. ismaskscopyN( sfilledsxsdsNumericsequalsvaluescsmask_orsgetmasksmsarrayscopy(sxsvaluescopyscsdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_equals cCs‹||jo|}|}|}nt|dƒ}titi||ƒti ||ƒƒ}t |t |ƒƒ}t|d|d|ƒSdS(sjx with mask of all values of x that are inside [v1,v2] v1 and v2 can be given in either order. ismaskscopyN(sv2sv1stsfilledsxsdsNumerics logical_ands less_equals greater_equalscsmask_orsgetmasksmsarrayscopy(sxsv1sv2scopyscsdsmst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys masked_insides  *cCs‹||jo|}|}|}nt|dƒ}titi||ƒti ||ƒƒ}t |t |ƒƒ}t|d|d|ƒSdS(skx with mask of all values of x that are outside [v1,v2] v1 and v2 can be given in either order. ismaskscopyN(sv2sv1stsfilledsxsdsNumerics logical_orslesssgreaterscsmask_orsgetmasksmsarrayscopy(sxsv1sv2scopyscsdsmst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysmasked_outsides  *cCs_t|ƒ}tit|ƒ|ƒ}|tjot |ƒSnt |dti||ƒƒSdS(sCopy of a with a new shape.smaskN( sgetmasksasmsNumericsreshapesfilledsnewshapesdsNones masked_array(sasnewshapesdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysreshape)s   cCsYt|ƒ}tit|ƒƒ}|tjot|ƒSnt|dti|ƒƒSdS(s-a as one-dimensional, may share data and masksmaskN( sgetmasksasmsNumericsravelsfilledsdsNones masked_array(sasdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysravel2s   cCs¿g}x!|D]}|it|ƒƒq Wti||ƒ}x0|D]}t|ƒt j oPqCqCWt |ƒSg}x!|D]}|it |ƒƒq|Wti||ƒ}t |d|ƒSdS(s+Concatenate the arrays along the given axissmaskN( sdsarrayssxsappendsfilledsNumerics concatenatesaxissgetmasksNones masked_arraysdms getmaskarray(sarrayssaxissdmsdsx((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys concatenate;s"  cCszt|ƒ}t|ƒiƒ}|tjotti|||ƒƒSn/tti|||ƒdti|||ƒƒSdS(s;take(a, indices, axis=0) returns selection of items from a.smaskN( sgetmasksasms masked_arraysraw_datasdsNonesNumericstakesindicessaxis(sasindicessaxissdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pystakeKs  cCskt|ƒ}t|ƒ}|tjotti||ƒƒSn)tti||ƒdti||ƒƒSdS(s9transpose(a, axes=None) reorder dimensions per tuple axessmaskN( sgetmasksasmsfilledsdsNones masked_arraysNumerics transposesaxes(sasaxessdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys transposeUs   cCs{|iƒ}t|ƒ}t|ƒ}ti |||ƒt |ƒ}|t j o'|i ƒti |iƒ|dƒndS(s{put(a, indices, values) sets storage-indexed locations to corresponding values. values and indices are filled if necessary.iN(sasraw_datasdsfilledsindicessindsvaluessvsNumericsputsgetmasksmsNones unshare_masksraw_mask(sasindicessvaluessdsmsindsv((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysput_s      cCsv|tjodSnti|iƒ||ƒt|ƒ}|tjodSn|i ƒti|i ƒ|dƒdS(s0put (a, mask, values) sets a where mask is true.Ni( smasksNonesNumericsputmasksasraw_datasvaluessgetmasksms unshare_masksraw_mask(sasmasksvaluessm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysputmaskjs    cCs„t|dƒ}t|dƒ}t|iƒdjodf|_nt|iƒdjodf|_ntti ||ƒƒSdS(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. iiN( sfilledsasfasbsfbslensshapes masked_arraysNumerics innerproduct(sasbsfasfb((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys innerproducttscCsÉt|dƒi}t|dƒi}ti||ƒ}t |ƒ}t |ƒ}|t jo |t jot |ƒSnt|ƒ}t|ƒ}tdtid|d|ƒddƒ}t ||ƒSdS(s:outerproduct(a,b) = {a[i]*b[j]}, has shape (len(a),len(b))iiscopyN(sfilledsasflatsfasbsfbsNumerics outerproductsdsgetmasksmasmbsNones masked_arrays getmaskarrays make_masksm(sasbsmasdsmbsmsfasfb((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys outerproduct€s    *cCs2tt|dƒtit|dƒddƒƒSdS(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þÿÿÿN(s innerproductsfilledsasNumericsswapaxessb(sasb((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdotŽsiÿÿÿÿcCsmt|dƒ}t|ƒ}|tj oti|||ƒ}nti|t|ƒ|ƒ}t ||ƒSdS(sqSelect those parts of x for which condition is true. Masked values in condition are considered false. iN( sfilleds conditionscsgetmasksxsmsNonesNumericscompresss dimensionsds masked_array(s conditionsxs dimensionscsdsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyscompress•s  s_minimum_operationcBs8tZdZd„Zed„Zdd„Zd„ZRS(sObject to calculate minimacCsdS(s`minimum(a, b) or minimum(a) In one argument case returns the scalar minimum. N((sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__¢scCsŸ|tjoqt|ƒ}|tjott|ƒiƒ}|Sn|i ƒ}t |ƒdjot Sq›t|i ƒƒSntt||ƒ||ƒdSdS(sExecute the call behavior.i.N(sbsNonesgetmasksasmsminsfilledsflatsds compressedsacslensmaskedsraw_dataswheresless(sselfsasbsdsmsac((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__¨s    icCs•t|ƒ}|tjo)t|ƒ}ttii ||ƒƒSnPtii t|t |ƒƒ|ƒ}ti i ||ƒ}t||t |ƒƒSdS(s#Reduce target along the given axis.N(sgetmaskstargetsmsNonesfilledsts masked_arraysNumericsminimumszreducesaxissminimum_fill_values logical_ands fill_value(sselfstargetsaxissmst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysreduce·s   $cCs˜t|ƒ}t|ƒ}|tjo |tjo t}n+t|ƒ}t|ƒ}ti ||ƒ}t i i t |ƒt |ƒƒ}t||ƒSdS(s<Return the function applied to the outer product of a and b.N(sgetmasksasmasbsmbsNonesms getmaskarrays logical_orsoutersNumericsminimumsfilledsds masked_array(sselfsasbsmasdsmbsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysouterÂs     !(s__name__s __module__s__doc__s__init__sNones__call__sreducesouter(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys_minimum_operation s    s_maximum_operationcBs8tZdZd„Zed„Zdd„Zd„ZRS(sObject to calculate maximacCsdS(s`maximum(a, b) or maximum(a) In one argument case returns the scalar maximum. N((sself((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__init__ÓscCsŸ|tjoqt|ƒ}|tjott|ƒiƒ}|Sn|i ƒ}t |ƒdjot Sq›t|i ƒƒSntt||ƒ||ƒdSdS(sExecute the call behavior.i.N(sbsNonesgetmasksasmsmaxsfilledsflatsds compressedsacslensmaskedsraw_dataswheresgreater(sselfsasbsdsmsac((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys__call__Ùs    icCs•t|ƒ}|tjo)t|ƒ}ttii ||ƒƒSnPtii t|t |ƒƒ|ƒ}ti i ||ƒ}t||t |ƒƒSdS(s#Reduce target along the given axis.N(sgetmaskstargetsmsNonesfilledsts masked_arraysNumericsmaximumszreducesaxissmaximum_fill_values logical_ands fill_value(sselfstargetsaxissmst((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysreduceès   $cCs˜t|ƒ}t|ƒ}|tjo |tjo t}n+t|ƒ}t|ƒ}ti ||ƒ}t i i t |ƒt |ƒƒ}t||ƒSdS(s<Return the function applied to the outer product of a and b.N(sgetmasksasmasbsmbsNonesms getmaskarrays logical_orsoutersNumericsmaximumsfilledsds masked_array(sselfsasbsmasdsmbsm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysouterós     !(s__name__s __module__s__doc__s__init__sNones__call__sreducesouter(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys_maximum_operationÑs    cCsv|tjot|ƒ}nt||ƒ}ti||ƒ}t |ƒtjot |ƒSnt ||ddƒSdS(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. scopyiN( s fill_valuesNonesdefault_fill_valuesxsfilledsdsNumericssortsaxissssgetmasks masked_arrays masked_values(sxsaxiss fill_valuesdss((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyssorts  cCsktit|ƒ|||ƒ}t|ƒ}|t jot ||ƒSn t |ti||||ƒƒSdS(s9diagonal(a,k=0,axis1=0, axis2=1) = the k'th diagonal of aN( sNumericsdiagonalsfilledsasksaxis1saxis2sdsgetmasksmsNones masked_array(sasksaxis1saxis2smsd((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysdiagonals   cCs#t||ƒ}ti||ƒSdS(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(sfilledsxs fill_valuesdsNumericsargsortsaxis(sxsaxiss fill_valuesd((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysargsortscCs#t||ƒ}ti||ƒSdS(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 fill_value(x). Returns a Numeric array if x has more than one dimension. Otherwise, returns a scalar index. N(sfilledsxs fill_valuesdsNumericsargminsaxis(sxsaxiss fill_valuesd((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysargmin'scCsU|tjo%t|ƒ}y | }Wq2q2Xnt||ƒ}ti||ƒSdS(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( s fill_valuesNonesdefault_fill_valuesxsfilledsdsNumericsargmaxsaxis(sxsaxiss fill_valuesd((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysargmax1s   cCstti||ƒƒSdS(s+apply f to s to create array as in Numeric.N(s masked_arraysNumerics fromfunctionsfss(sfss((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys fromfunctionAscCsRt|tƒo|tjp||iƒjo|Snt|d|ddƒSdS(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. stypecodescopyiN(s isinstancesdatas MaskedArraystypecodesNonesarray(sdatastypecode((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysasarrayEs0i,cCs7t|tƒo|ii}nti|ƒ}|SdS(s2returns shape tuple for 'a'. MaskedArray version.N(s isinstancesas MaskedArrays_datasshapesssNumeric(sass((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysshapeRs cCs;t|tƒoti|i|ƒSnti||ƒSdS(s2returns shape tuple for 'a'. MaskedArray version.N(s isinstancesas MaskedArraysNumericssizes_datasaxis(sasaxis((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pyssizeZscCsti|ƒdS(sBSet the maximum # of elements for printing arrays. <=0 = no limitN(s arrayprints set_summarysm(sm((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysset_print_limitascCstiƒSdS(s3Get the maximum # of elements for printing arrays. N(s arrayprintsget_summary_threshhold(((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pysget_print_limites(®s__doc__snumarray.numericsnumericsNumericsespisNewAxissnumarray.objectssobjectssobjsnumarray.arrayprints arrayprintsnumarraysstringstypesssyss Precisionswarningss _warningssInt8sMaskTypesdivide_tolerances ExceptionsMAErrors_MaskedPrintOptionsmasked_print_optionsarraysastypesFloat32sdefault_real_fill_values Complex32sdefault_complex_fill_valuesdefault_character_fill_values UnsignedInt8sdefault_integer_fill_valuesdefault_object_fill_valuesdefault_fill_valuesminimum_fill_valuesmaximum_fill_valuesset_fill_valuesgetmasks getmaskarraysis_masks make_masksmake_mask_nones create_masksmask_orsNonesfilleds fill_valuescommon_fill_valuesdomain_check_intervals domain_tansdomain_greatersdomain_greater_equalsmasked_unary_operationsdomain_safe_dividesdomained_binary_operationsmasked_binary_operationssqrtslogslog10sexps conjugatessinscosstansarcsinsarccossarctanssinhscoshstanhsabsolutesfabssnegativesnonzerosaroundsfloorsceilssometruesalltrues logical_notsaddssubtractsreducesmultiplysdivides true_divides floor_divides remaindersfmodshypotsarctan2sequals not_equals less_equals greater_equalslesssgreaters logical_ands logical_ors logical_xors bitwise_ands bitwise_ors bitwise_xorsranksshapessizesUsesOpPrioritys MaskedArrays MaskedScalarsmaskeds isMaskedArraysisarraysisMAsallclosesallequals masked_valuess masked_objects arrayrangesaranges fromstrings left_shifts right_shiftsresizesrepeatsidentitysindicessIntszerossonesscountspowers masked_arrayssumsproductsaverageswhereschooses masked_wheresmasked_greatersmasked_greater_equals masked_lesssmasked_less_equalsmasked_not_equals masked_equals masked_insidesmasked_outsidesreshapesravels concatenatestakes transposesputsputmasks innerproducts outerproductsdotscompresss_minimum_operationsminimums_maximum_operationsmaximumssortsdiagonalsargsortsargminsargmaxs fromfunctionsasarrays __MaxElementssset_print_limitsget_print_limit(œs make_masks create_maskslesssputmasks masked_insidesproductsshapesdefault_object_fill_values_maximum_operations masked_arraysset_print_limitsaddschoosesNewAxissravelsaroundsceilsonessranksdomained_binary_operationsstringsisarrays arrayranges bitwise_orsMAErrorsNumericsmake_mask_nonesfmodsindicess bitwise_andswheresargmaxs MaskedArrayscoshshypots logical_andsarangesallclosesnonzerosasarrayssumsmultiplysmaskedsdefault_real_fill_valuesfilleds concatenatespowersmasked_greater_equals transposesdiagonals masked_equalssinhsset_fill_values remaindersrepeatscossobjsdefault_fill_valuesmasked_binary_operationsputsequalsdomain_safe_dividesargsortsarccoss less_equals masked_objectsmasked_print_optionsdividesdomain_check_intervalsdefault_integer_fill_valuesarctan2smasked_unary_operationsMaskTypesis_masksarrayssizes logical_xorslogssometruessubtractsnegativesdivide_tolerancesalltrues masked_lesssidentitys domain_tansmasked_outsideszerosscommon_fill_valuespislog10s _warningsssortsdefault_complex_fill_valuesarcsins logical_notsresizes MaskedScalarscounts not_equals fromfunctionsargminstanhs true_dividesmaximumsget_print_limitsgreaters outerproducts left_shiftsdomain_greatersmasked_less_equals logical_orsminimumsmasked_greaters conjugates innerproductstansabsolutes bitwise_xors_minimum_operationsfabssnumarraysmasked_not_equalsreshapessqrtsarctansgetmasks right_shiftstakesmaximum_fill_values floor_dividesmask_orssinsdomain_greater_equals isMaskedArrayscompresss masked_wheressyssfloorsminimum_fill_values arrayprintsallequals __MaxElementsstypessdefault_character_fill_valueses greater_equals fromstringsaverages fill_values_MaskedPrintOptions masked_valuessexps getmaskarraysisMAsdot((sG/mit/python/arch/i386_rh9/lib/python2.3/site-packages/numarray/ma/MA.pys? sZ                    ',K           ÿÿÿr           l            / /