IR-e;\ddlZddlZddlmZmZmZmZmZm Z m Z ddl m Z m Z ddlmZmZmZgdZGdde ZGd d eZGd d eZGd deZGdde ZGddeZGddeZGddeZee_ee_ee_ee jZde_ee jZ de _ee j!Z"de j!j#e"_ee j$Z%de j$j#e%_dS)N) CompositeUnitUnitUnitConversionError UnitsError UnitTypeErrordimensionless_unscaled photometric)FunctionQuantityFunctionUnitBase)dBdexmag) LogUnitMagUnitDexUnit DecibelUnit LogQuantity MagnitudeDecibelDexSTmagABmagM_bolm_bolcneZdZdZedZedZdZdZdZ dZ dZ d Z d Z d Zd S) ra!Logarithmic unit containing a physical one. Usually, logarithmic units are instantiated via specific subclasses such `~astropy.units.MagUnit`, `~astropy.units.DecibelUnit`, and `~astropy.units.DexUnit`. Parameters ---------- physical_unit : `~astropy.units.Unit` or `string` Unit that is encapsulated within the logarithmic function unit. If not given, dimensionless. function_unit : `~astropy.units.Unit` or `string` By default, the same as the logarithmic unit set by the subclass. ctSNrselfs Blib/python3.11/site-packages/astropy/units/function/logarithmic.py_default_function_unitzLogUnit._default_function_unit6 ctSr)rr s r"_quantity_classzLogUnit._quantity_class:sr%cZtj|jtj|S)zjTransformation from value in physical to value in logarithmic units. Used in equivalency. )rto_function_unitnplog10r!xs r" from_physicalzLogUnit.from_physical>s!vd)28A;;777r%cHd|jt|zS)zjTransformation from value in logarithmic to value in physical units. Used in equivalency. )r*r)rr-s r" to_physicalzLogUnit.to_physicalDs"T(++C3333r%c@ t|d||jn/#t$r tcYSt $rt dwxYwt|dt }td|j|g||g}| |S)uAdd/subtract LogUnit to/from another unit, and adjust physical unit. self and other are multiplied by sign_self and sign_other, resp. We wish to do: ±lu_1 + ±lu_2 -> lu_f (lu=logarithmic unit) and pu_1^(±1) * pu_2^(±1) -> pu_f (pu=physical unit) Raises ------ UnitsError If function units are not equivalent. function_unitz;Can only add/subtract logarithmic units of compatible type. physical_unitr ) getattr_tor*AttributeErrorNotImplementedrrr_physical_unit_copy)r!other sign_self sign_otherother_physical_unitr5s r"_add_and_adjust_physical_unitz%LogUnit._add_and_adjust_physical_unitNs   E?E 2 2 6 6t7J K K K K " " "! ! ! !   M   &e_>TUU% #%89Iz;R  zz-(((s),AAc<||jdzS)N)r;r5r s r"__neg__zLogUnit.__neg__oszz$,4555r%c0||ddSNr r@r!r<s r"__add__zLogUnit.__add__r11%R@@@r%c0||ddSrErFrGs r"__radd__zLogUnit.__radd__w11%R@@@r%c0||ddS)Nr rBrFrGs r"__sub__zLogUnit.__sub__zrLr%c0||ddS)NrBr rFrGs r"__rsub__zLogUnit.__rsub__}rIr%N)__name__ __module__ __qualname____doc__propertyr#r'r/r2r@rCrHrKrNrPr%r"rr#s$XX888 444)))B666AAA AAAAAAAAAAAr%rc>eZdZdZedZedZdS)raLogarithmic physical units expressed in magnitudes. Parameters ---------- physical_unit : `~astropy.units.Unit` or `string` Unit that is encapsulated within the magnitude function unit. If not given, dimensionless. function_unit : `~astropy.units.Unit` or `string` By default, this is ``mag``, but this allows one to use an equivalent unit such as ``2 mag``. ctSr)rr s r"r#zMagUnit._default_function_unitr$r%ctSr)rr s r"r'zMagUnit._quantity_classsr%NrQrRrSrTrUr#r'rVr%r"rrsT  XXr%rcPeZdZdZedZedZdfd ZxZS)raLogarithmic physical units expressed in magnitudes. Parameters ---------- physical_unit : `~astropy.units.Unit` or `string` Unit that is encapsulated within the magnitude function unit. If not given, dimensionless. function_unit : `~astropy.units.Unit` or `string` By default, this is ``dex``, but this allows one to use an equivalent unit such as ``0.5 dex``. ctSrrr s r"r#zDexUnit._default_function_unitr$r%ctSr)rr s r"r'zDexUnit._quantity_classr$r%genericc|dkr1|jtkrdSd|j|dStS)Ncdsz[-][)format])r5r to_stringsuper)r!rb __class__s r"rdzDexUnit.to_strings[ U??!%;;;uI4-77v7FFIIII77$$&& &r%)r^) rQrRrSrTrUr#r'rd __classcell__rfs@r"rrs}  XX''''''''''r%rc>eZdZdZedZedZdS)raLogarithmic physical units expressed in dB. Parameters ---------- physical_unit : `~astropy.units.Unit` or `string` Unit that is encapsulated within the decibel function unit. If not given, dimensionless. function_unit : `~astropy.units.Unit` or `string` By default, this is ``dB``, but this allows one to use an equivalent unit such as ``2 dB``. ctSr)r r s r"r#z"DecibelUnit._default_function_units r%ctSr)rr s r"r'zDecibelUnit._quantity_classsr%NrZrVr%r"rrsT  XXr%rceZdZdZeZdZdZdZdZ dZ dZ fdZ d Z fd Zfd Zfd Zd ZdZddZddZddZddZddZejddDzZxZS)raA representation of a (scaled) logarithm of a number with a unit. Parameters ---------- value : number, `~astropy.units.Quantity`, `~astropy.units.LogQuantity`, or sequence of quantity-like. The numerical value of the logarithmic quantity. If a number or a `~astropy.units.Quantity` with a logarithmic unit, it will be converted to ``unit`` and the physical unit will be inferred from ``unit``. If a `~astropy.units.Quantity` with just a physical unit, it will converted to the logarithmic unit, after, if necessary, converting it to the physical unit inferred from ``unit``. unit : str, `~astropy.units.UnitBase`, or `~astropy.units.FunctionUnitBase`, optional For an `~astropy.units.FunctionUnitBase` instance, the physical unit will be taken from it; for other input, it will be inferred from ``value``. By default, ``unit`` is set by the subclass. dtype : `~numpy.dtype`, optional The ``dtype`` of the resulting Numpy array or scalar that will hold the value. If not provided, is is determined automatically from the input value. copy : bool, optional If `True` (default), then the value is copied. Otherwise, a copy will only be made if ``__array__`` returns a copy, if value is a nested sequence, or if a copy is needed to satisfy an explicitly given ``dtype``. (The `False` option is intended mostly for internal use, to speed up initialization where a copy is known to have been made. Use with care.) Examples -------- Typically, use is made of an `~astropy.units.FunctionQuantity` subclass, as in:: >>> import astropy.units as u >>> u.Magnitude(-2.5) >>> u.Magnitude(10.*u.count/u.second) >>> u.Decibel(1.*u.W, u.DecibelUnit(u.mW)) # doctest: +FLOAT_CMP c|jt|dtz}|jt|d|z}|||SNunit_function_viewror6rrp _new_viewr!r<new_unitresults r"rHzLogQuantity.__add__sK9wuf6LMMM$wu6F'N'NN~~fh///r%c,||Sr)rHrGs r"rKzLogQuantity.__radd__ ||E"""r%c|jt|dtz}|j}|t|d|z }|||Srnror6rrp _set_unitr!r<rt function_views r"__iadd__zLogQuantity.__iadd__Q9wuf6LMMM+ (8%@@@  x    r%c|jt|dtz }|jt|d|z }|||Srnrqrss r"rNzLogQuantity.__sub__sI9wuf6LMMM$wu6F'N'NN~~fh///r%c|jt|dt}|jt|d|}||j}|||Srn)rorPr6rrpr)r4rrrss r"rPzLogQuantity.__rsub__sq9%%geV=S&T&TUU$--ge=Mu.U.UVV8122~~fh///r%c|jt|dtz }|j}|t|d|z}|||Srnryr{s r"__isub__zLogQuantity.__isub__&r~r%c8t|tjr_|jj|z}|t j|z}|||j |St |Sr) isinstancenumbersNumberror5viewr+ndarrayrrr;re__mul__r!r<new_physical_unitrurfs r"rzLogQuantity.__mul__.sz eW^ , , *!% 7 > YYrz**U2F>>&$)//:K*L*LMM M77??5)) )r%c,||Sr)rrGs r"__rmul__zLogQuantity.__rmul__9rwr%ct|tjrJ|jj|z}|j}||z}||j||St |Sr) rrrror5rprzr;re__imul__r!r<rr|rfs r"rzLogQuantity.__imul__<sw eW^ , , + $ 7 >  /M U "M NN49??+<== > > >K77##E** *r%c>t|tjrb|jjd|z z}|t j|z }|||j |St |SrE) rrrror5rr+rrrr;re __truediv__rs r"rzLogQuantity.__truediv__Fs eW^ , , .!% 7AI F YYrz**U2F>>&$)//:K*L*LMM M77&&u-- -r%ct|tjrM|jjd|z z}|j}||z}||j||St |SrE) rrrror5rprzr;re __itruediv__rs r"rzLogQuantity.__itruediv__Qs{ eW^ , , / $ 7AI F  /M U "M NN49??+<== > > >K77''.. .r%c t|}n#t$r tcYSwxYw|j|z}|t j|z}|||Sr)float TypeErrorr9rorr+rrr)r!r<rt new_values r"__pow__zLogQuantity.__pow__[ss "%LLEE " " "! ! ! ! "9e#IIbj))U2 ~~i222s &&cR t|}n#t$r tcYSwxYwt||jstS |jj|j}|tj dxx|j |z cc<nb#t$rU | |}n#t$r tcYcYSwxYw||tj d<YnwxYw|||S)N.)rrr9r _unit_classror5r7rr+rr/r _to_valuerz)r!r<factorvalues r" __ilshift__zLogQuantity.__ilshift__esS "KKEE " " "! ! ! ! "%!122 "! ! JY,001DEEF IIbj ! !# & & &$)*A*A&*I*I I & & & &# / / / &u--& & & &%%%%%% &*/DIIbj ! !# & & & / u sB&&$B00 D;CDC'"D&C''%DDNrcj|jjdz}|tj|||||S)Noutddofro)ror4_wrap_functionr+varr!axisdtyperrros r"rzLogQuantity.vars5y&)""264CdQU"VVVr%c|jt}|tj|||||S)Nr)ror;rrr+stdrs r"rzLogQuantity.stds:y566""264CdQU"VVVr%c|jt}|tj|||S)N)rro)ror;rrr+ptp)r!rrros r"rzLogQuantity.ptps5y566""264St"DDDr%r rBc|jt}|tj|||SN)ro)ror;rrr+diff)r!nrros r"rzLogQuantity.diffs5y566""27At$"???r%c|jt}|tj|||Sr)ror;rrr+ediff1d)r!to_endto_beginros r"rzLogQuantity.ediff1ds5y566""2:vxd"KKKr%c8h|]}tt|SrV)r6r+).0functions r" zLogQuantity.s5DDD"*HDDDr%)rrrrr)NNNr)NN)r rB)rQrRrSrTrrrHrKr}rNrPrrrrrrrrrrrrrr _supported_functionsrgrhs@r"rrs++\K000###000000 * * * * *###+++++ . . . . ./////3338WWWWWWWWEEEE@@@@LLLL,@DD.VDDDr%rceZdZeZdS)rN)rQrRrSrrrVr%r"rrKKKr%rceZdZeZdS)rN)rQrRrSrrrVr%r"rrsKKKr%rceZdZeZdS)rN)rQrRrSrrrVr%r"rrrr%rz6ST magnitude: STmag=-21.1 corresponds to 1 erg/s/cm2/Az7AB magnitude: ABmag=-48.6 corresponds to 1 erg/s/cm2/Hzz=Absolute bolometric magnitude: M_bol=0 corresponds to L_bol0=z=Apparent bolometric magnitude: m_bol=0 corresponds to f_bol0=)&rnumpyr+ astropy.unitsrrrrrrr corer r unitsr rr__all__rrrrrrrr_function_unit_classSTfluxrrTABfluxrBolrsibolrrVr%r"rs54444444    ]A]A]A]A]A]A]A]A@g.'''''g'''@'.FFFFF"FFFR+k #%"  "##H  "##I    XKODVXX     XKODVXX r%