³ò éeJc @sƒdZddkZddkZddkZddkZddkZddkTd„Zd„Z d dd d d „Z e d „Z dS(s= lib/simfit.py Simultaneous Larentzian Fit Lev, 8 May 2009 iÿÿÿÿN(t*cCsAtii||ti| |ƒ||dti|||ƒS(sw Return spectrum of a Lorentzian described by (f0, A, T2, phi) delayed by 'delay' at given frequencies 'f' i(tnmrfitt Lorentziantltmathtexptpi(tf0tAtT2tphitdelaytf((sc:\py\lib\simfit.pytdelayed_lorentzian sc Csw|\}}}}g}xX|iƒD]J}||}ti|t|it||||||iƒƒgƒ}q%W|S(s) Calculate residuals for the fit (tkeystnumpyt concatenatetabstwR R ( tparamtspectraRRR R tresR tspectrum((sc:\py\lib\simfit.pytsim_lorentzian_residualss  >iigL Ag@$Agˆ³@c CsÌtii|i|ƒid|d|d|ƒƒ}|i|}|i|}h} xRtd|ƒD]A} |i|| |ƒid|d|d|ƒ| || |s0