;ò …Rã?c@s†dkZdkZdklZdkZdfd„ƒYZee_defd„ƒYZee_d„Zd„Z e djod kTe d d d ƒZ e e ƒZee feƒZed ƒGe d ƒGHeiƒZedƒGde dƒGHeiƒZedƒGdddGHeiƒGHee dd…efƒee ƒZee e feƒZeddƒGedƒedƒGHndS(N(sindex_expressionsInterpolatingFunctioncBsþtZdZed„Zd„Zd„Zd„Zd„Zd„Z dd„Z dd „Z dd „Z dd „Z d „Zd „Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„Zd„ZRS(sjFunction defined by values on a grid using interpolation An interpolating function of n variables with m-dimensional values is defined by an (n+m)-dimensional array of values and n one-dimensional arrays that define the variables values corresponding to the grid points. The grid does not have to be equidistant. Constructor: InterpolatingFunction(|axes|, |values|, |default|=None) Arguments: |axes| -- a sequence of one-dimensional arrays, one for each variable, specifying the values of the variables at the grid points |values| -- an array containing the function values on the grid |default| -- the value of the function outside the grid. A value of 'None' means that the function is undefined outside the grid and that any attempt to evaluate it there yields an exception. Evaluation: 'function(x1, x2, ...)' yields the function value obtained by linear interpolation. Indexing: all array indexing operations except for the NexAxis operator are supported. cCs~t|ƒt|iƒjo td‚nt|ƒ|_||_||_f|_x$|iD]}|i|i|_q]WdS(NsInconsistent arguments( slensaxessvaluessshapes ValueErrorslistsselfsdefaultsaxis(sselfsaxessvaluessdefaultsaxis((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys__init__/s     cGst|ƒt|iƒjo td‚nytt||iƒ}Wn9tj o-}|i t j o |i Sq~t|‚nXf}x1|D])}|t|d|ddd…}q‹W|i|}x4|D],}d|d|d|d|d}qÌW|SdS(NsWrong number of argumentsiif1.0i(slenspointssselfsaxess TypeErrorsmaps_lookups neighbourss ValueErrorstextsdefaultsNonesslicessitemsindex_expressionsvalues(sselfspointssslicesstextsitems neighbourssvalues((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys__call__9s"  ' *cCst|idƒSdS(Ni(slensselfsaxes(sself((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys__len__KscCsit|ƒ}|tdƒjoXt|iƒdjo |id||i|fSqe|i|id|i|ƒSnï|ttt ƒƒjo:|id|g|id}|i||i|ƒSnœ|tfƒjog}|i}xG|D]?}t|ƒtdƒjo|i |d|ƒn|d=qôW||}|i||i|ƒSn t d‚dS(Niisillegal index type(stypesistislensselfsaxessvaluess _constructorsslicesNonesrestsitemsappends TypeError(sselfsisaxessitemsreststi((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys __getitem__Ns&  "   cCs@|id||!g|id}|i||i||!ƒSdS(Nii(sselfsaxessisjs _constructorsvalues(sselfsisjsaxes((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys __getslice__ds"cCs|djof|i}y |i}Wntj onX|i}y |i}WnnX|i|i||ƒSn¨|djo‘y|ii}Wn tj od|i}nX|i}y|ii}Wn(yd|i}Wq÷t }q÷XnX|i|i||ƒSn t |‚dS(Nsrealsimagi( sattrsselfsvaluessreals ValueErrorsdefaults _constructorsaxessimagsNonesAttributeError(sselfsattrsdefaultsvalues((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys __getattr__hs6       icCs•|i|}titi||ƒti||ƒƒ}|i| ti ||ƒg|i|d}ti ||i |ƒ}|i||tƒSdS(s«Returns a new InterpolatingFunction whose grid is restricted to the interval from |first| to |last| along the variable whose number is |variable|. iN(sselfsaxessvariablesxsNumerics logical_ands greater_equalsfirsts less_equalslastscscompresssi_axessvaluessi_valuess _constructorsNone(sselfsfirstslastsvariablescsi_valuessxsi_axes((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pysselectInterval†s /c Cs?|i|}|tddd…tddd…td}|tddd…tddd…td}tddd…t|ii ƒ|dtt i }|i||i||d|d |}d|d|d }|i| |g|i|d}t}|itj o d}n|i|||ƒSdS(swReturns a new InterpolatingFunction describing the derivative with respect to |variable| (an integer). Ni.iÿÿÿÿf0.5f0.0(sselfsaxessvariablesdiffaxissindex_expressionsuislislensvaluessshapesNumericsNewAxissaisd_valuessd_axessNones d_defaultsdefaults _constructor( sselfsvariables d_defaultsaisd_axessdiffaxisslisd_valuessui((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys derivative“s 006,# c Csq|i|}|tddd…tddd…td}|tddd…tddd…td}tddd…t|ii ƒ|dtt i }tddd…t|ii ƒ|dtt i }dt ii|i||i||||||ƒ}t|ii ƒ}d||6s(sScientific.IO.NetCDFs NetCDFFilesfilenamesselfsfilesmaps axesnamessaxess variabless variablenamesvaluessdefaultsshapesaxis(sselfsfilenames axesnamess variablenamesdefaults NetCDFFilesaxis((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys__init__3s    (s__name__s __module__s__doc__sNones__init__(((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pysNetCDFInterpolatingFunctions cCs¼tiiti||ƒƒ}|t|ƒjo8|||djo|d|ddfSqid}n|djo td‚n|d}|||||||}|||fSdS(Niif1.0isPoint outside grid of values( sNumericsaddsreduces less_equalsaxisspointsjslens ValueErrorsisweight(spointsaxissweightsisj((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys_lookupBs    cCsƒt|ƒdjotd„|dƒSnUt|dƒ}g}x7|dD]+}x"|D]}|i|f|ƒqYWqLW|SdS(NicCs|fS(N(sx(sx((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pysQsi( slensaxessmaps _combinationssrestslsxsysappend(saxesslsrestsysx((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys _combinationsOs s__main__(s*if1.1000000000000001f0.10000000000000001f0.22f0.34999999999999998f0.5f0.41999999999999998f1.5f0.23000000000000001f0.77000000000000002(sNumerics PolynomialsScientific.indexingsindex_expressionsoperatorsInterpolatingFunctions _constructorsNetCDFInterpolatingFunctions_lookups _combinationss__name__sarangesaxisssqrtsvaluessss derivativessdsintegralssisdefiniteIntegralssinsNewAxisscosssc(s _combinationss_lookupssisNumericsoperatorsssvaluessindex_expressionssds PolynomialsaxissscsNetCDFInterpolatingFunctionsInterpolatingFunction((sM/mit/python/lib/python2.3/site-packages/Scientific/Functions/Interpolation.pys?s0    ÿ '       &