;ò …Rã?c@sDdZdkZdfd„ƒYZd„Zd„Zdd„ZdS(s`This module provides automatic differentiation for functions with any number of variables. Instances of the class DerivVar represent the values of a function and its partial derivatives with respect to a list of variables. All common mathematical operations and functions are available for these numbers. There is no restriction on the type of the numbers fed into the code; it works for real and complex numbers as well as for any Python type that implements the necessary operations. This module is as far as possible compatible with the n-th order derivatives module Derivatives. If only first-order derivatives are required, this module is faster than the general one. Example: >>>print sin(DerivVar(2)) produces the output >>>(0.909297426826, [-0.416146836547]) The first number is the value of sin(2); the number in the following list is the value of the derivative of sin(x) at x=2, i.e. cos(2). When there is more than one variable, DerivVar must be called with an integer second argument that specifies the number of the variable. Example: >>>x = DerivVar(7., 0) >>>y = DerivVar(42., 1) >>>z = DerivVar(pi, 2) >>>print (sqrt(pow(x,2)+pow(y,2)+pow(z,2))) produces the output >>>(42.6950770511, [0.163953328662, 0.98371997197, 0.0735820818365]) The numbers in the list are the partial derivatives with respect to x, y, and z, respectively. Note: It doesn't make sense to use DerivVar with different values for the same variable index in one calculation, but there is no check for this. I.e. >>>print DerivVar(3, 0)+DerivVar(5, 0) produces >>>(8, [2]) but this result is meaningless. NsDerivVarcBsUtZdZddd„Zd„Zd„Zd„Zd„Zd„Zd „Z d „Z d „Z d „Z d „Z e Zd„Zd„Zd„ZeZd„Zd„Zed„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Z d„Z!d„Z"d „Z#d!„Z$d"„Z%d#„Z&d$„Z'RS(%såVariable with derivatives Constructor: DerivVar(|value|, |index| = 0) Arguments: |value| -- the numerical value of the variable |index| -- the variable index (an integer), which serves to distinguish between variables and as an index for the derivative lists. Each explicitly created instance of DerivVar must have a unique index. Indexing with an integer yields the derivatives of the corresponding order. iicCs~|djo td‚n||_|djo g|_n>t|ƒtgƒjo ||_n|dgdg|_dS(NisOnly first-order derivativesi(sorders ValueErrorsvaluesselfsderivstypesindex(sselfsvaluesindexsorder((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__init__Vs      cCsJ|djp |djo td‚n|djo |iSn|iSdS(NiisIndex out of range(sitems ValueErrorsselfsvaluesderiv(sselfsitem((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys __getitem__as    cCs|i|if SdS(N(sselfsvaluesderiv(sself((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__repr__iscCst|i|ifƒSdS(N(sstrsselfsvaluesderiv(sself((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__str__lscCs2t|ƒo||fSn|t|gƒfSdS(N(s isDerivVarsothersselfsDerivVar(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys __coerce__os cCst|i|iƒSdS(N(scmpsselfsvaluesother(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__cmp__uscCs$t|i td„|iƒƒSdS(NcCs| S(N(sa(sa((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysys(sDerivVarsselfsvaluesmapsderiv(sself((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__neg__xscCs|SdS(N(sself(sself((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__pos__{scCs9t|iƒ}t|t|i|d„|iƒƒSdS(NcCs||S(N(sdsa(sasd((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys€s(sabssselfsvaluesabsvaluesDerivVarsmapsderiv(sselfsabsvalue((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__abs__~scCs|idjSdS(Ni(sselfsvalue(sself((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys __nonzero__‚scCs0t|i|itd„|i|iƒƒSdS(NcCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys‡s(sDerivVarsselfsvaluesothers _mapderivsderiv(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__add__…scCs0t|i|itd„|i|iƒƒSdS(NcCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysŒs(sDerivVarsselfsvaluesothers _mapderivsderiv(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__sub__ŠscCs0t|i|itd„|i|iƒƒSdS(NcCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyss(sDerivVarsothersvaluesselfs _mapderivsderiv(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__rsub__ŽscCsTt|i|itd„t|id„|iƒt|id„|iƒƒƒSdS(NcCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys”scCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys•s(sDerivVarsselfsvaluesothers _mapderivsmapsderiv(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__mul__’s cCs{|i o td‚nd|i}t|i|td„t|d„|iƒt|i||d„|iƒƒƒSdS(NsDerivVar divisionf1.0cCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysžscCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysŸs( sothersvaluesZeroDivisionErrorsinvsDerivVarsselfs _mapderivsmapsderiv(sselfsothersinv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__div__™s     cCs ||SdS(N(sothersself(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__rdiv__¢scCsÛ|tj o td‚nt|i|idƒ}||i}t ||id„|i ƒ}t |ƒot |i ƒdjoHt |ti|iƒd„|i ƒ}t|td„||ƒƒSnt||ƒSdS(Ns'DerivVar does not support ternary pow()icCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysªsicCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys®s(szsNones TypeErrorspowsselfsvaluesothersval1svalsmapsderivsderiv1s isDerivVarslensNumericslogsderiv2sDerivVars _mapderiv(sselfsotherszsvalsderiv2sderiv1sval1((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__pow__¥s   #  cCst||ƒSdS(N(spowsothersself(sselfsother((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys__rpow__²scCs5ti|iƒ}t|t|d„|iƒƒSdS(NcCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys·s(sNumericsexpsselfsvaluesvsDerivVarsmapsderiv(sselfsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysexpµscCsBti|iƒ}d|i}t|t|d„|iƒƒSdS(Nf1.0cCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys¼s( sNumericslogsselfsvaluesvsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyslog¹s cCsOti|iƒ}d|itidƒ}t|t|d„|i ƒƒSdS(Nf1.0i cCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysÁs( sNumericslog10sselfsvaluesvslogsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyslog10¾scCs?ti|iƒ}d|}t|t|d„|iƒƒSdS(Nf0.5cCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysÆs( sNumericssqrtsselfsvaluesvsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyssqrtÃs cCs:|idjo td‚ntti|iƒdƒSdS(Nis"can't differentiate sign() at zero(sselfsvalues ValueErrorsDerivVarsNumericssign(sself((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyssignÈs cCsGti|iƒ}ti|iƒ}t|t|d„|i ƒƒSdS(NcCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysÐs( sNumericssinsselfsvaluesvscossdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyssinÍscCsHti|iƒ}ti|iƒ }t|t|d„|i ƒƒSdS(NcCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysÕs( sNumericscossselfsvaluesvssinsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyscosÒscCsHti|iƒ}dt|dƒ}t|t|d„|i ƒƒSdS(Nf1.0icCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysÚs( sNumericstansselfsvaluesvspowsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pystan×scCsGti|iƒ}ti|iƒ}t|t|d„|i ƒƒSdS(NcCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysßs( sNumericssinhsselfsvaluesvscoshsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyssinhÜscCsGti|iƒ}ti|iƒ}t|t|d„|i ƒƒSdS(NcCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysäs( sNumericscoshsselfsvaluesvssinhsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyscosháscCsTti|iƒ}dtti|iƒdƒ}t|t |d„|i ƒƒSdS(Nf1.0icCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysés( sNumericstanhsselfsvaluesvspowscoshsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pystanhæscCsXti|iƒ}dtidt|idƒƒ}t|t |d„|i ƒƒSdS(Nf1.0icCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysîs( sNumericsarcsinsselfsvaluesvssqrtspowsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysarcsinës#cCsXti|iƒ}dtidt|idƒƒ}t|t |d„|i ƒƒSdS(Nf-1.0f1.0icCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysós( sNumericsarccossselfsvaluesvssqrtspowsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysarccosðs#cCsOti|iƒ}ddt|idƒ}t|t|d„|i ƒƒSdS(Nf1.0icCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysøs( sNumericsarctansselfsvaluesvspowsdsDerivVarsmapsderiv(sselfsdsv((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysarctanõscCsŽ|i|i|i|i}|i|}|i|}tti|i|iƒt d„t |d„|i ƒt |d„|i ƒƒƒSdS(NcCs||S(N(sasb(sasb((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysÿscCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyss( sselfsvaluesothersdensssosDerivVarsNumericsarctan2s _mapderivsmapsderiv(sselfsothersosssden((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysarctan2ús   cCsXdkl}l}||iƒ}|||iƒ}t|t|d„|i ƒƒSdS(N(sgammaspsicCs||S(N(sfsx(sxsf((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pyss( stranscendentalsgammaspsisselfsvaluesvsdsDerivVarsmapsderiv(sselfspsisdsvsgamma((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysgammas((s__name__s __module__s__doc__s__init__s __getitem__s__repr__s__str__s __coerce__s__cmp__s__neg__s__pos__s__abs__s __nonzero__s__add__s__radd__s__sub__s__rsub__s__mul__s__rmul__s__div__s__rdiv__sNones__pow__s__rpow__sexpslogslog10ssqrtssignssinscosstanssinhscoshstanhsarcsinsarccossarctansarctan2sgamma(((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pysDerivVarCsJ                               cCs!t|dƒo t|dƒSdS(s&Returns 1 if |x| is a DerivVar object.svaluesderivN(shasattrsx(sx((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys isDerivVar scCsett|ƒt|ƒƒ}||t|ƒdg}||t|ƒdg}t|||ƒSdS(Ni(smaxslensasbsnvarssmapsfunc(sfuncsasbsnvars((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys _mapderivsicCsdkl}t|ƒot|ƒo t|ƒo||||ƒSn4|t||ƒt||dƒt||dƒƒSdS(sùReturns a vector whose components are DerivVar objects. Arguments: |x|, |y|, |z| -- vector components (numbers) |index| -- the DerivVar index for the x component. The y and z components receive consecutive indices. (sVectoriiN(sScientific.Geometry.VectorsVectors isDerivVarsxsyszsDerivVarsindex(sxsyszsindexsVector((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys DerivVectors  '(s__doc__sNumericsDerivVars isDerivVars _mapderivs DerivVector(s _mapderivs isDerivVarsDerivVars DerivVectorsNumeric((sP/mit/python/lib/python2.3/site-packages/Scientific/Functions/FirstDerivatives.pys?;s  È