Changeset 618
- Timestamp:
- 09/16/08 09:28:24 (2 months ago)
- Files:
-
- trunk/pysynphot/lib/spectrum.py (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
- Modified
- Copied
- Moved
trunk/pysynphot/lib/spectrum.py
r617 r618 312 312 hdu.writeto(filename) 313 313 314 315 316 314 317 315 def integrate(self,fluxunits='photlam'): 316 #Extract the flux in the desired units 318 317 u=self.fluxunits 319 318 self.convert(fluxunits) 320 321 319 wave,flux=self.getArrays() 322 320 self.convert(u) 323 321 #then do the integration 324 322 return self.trapezoidIntegration(wave,flux) 325 323 … … 547 545 newasc = False 548 546 549 ## First need to pad the ends of the spectrum with zeros550 tapered = self.taper()551 552 547 ## Use numpy interpolation function 553 if tapered._wavetable[0]<tapered._wavetable[-1]:548 if self._wavetable[0]<self._wavetable[-1]: 554 549 oldasc = True 555 ans = N.interp(newwave, tapered._wavetable,556 tapered._fluxtable)550 ans = N.interp(newwave,self._wavetable, 551 self._fluxtable) 557 552 else: 558 553 oldasc = False 559 rev = N.interp(newwave, tapered._wavetable[::-1],560 tapered._fluxtable[::-1])554 rev = N.interp(newwave,self._wavetable[::-1], 555 self._fluxtable[::-1]) 561 556 ans = rev[::-1] 562 557 … … 1008 1003 newasc = False 1009 1004 1010 #First need to pad the ends with zeros(?)1011 tapered=self.taper()1012 1005 ## Use numpy interpolation function 1013 if tapered._wavetable[0]<tapered._wavetable[-1]:1006 if self._wavetable[0]<self._wavetable[-1]: 1014 1007 oldasc = True 1015 ans = N.interp(newwave, tapered._wavetable,1016 tapered._throughputtable)1008 ans = N.interp(newwave,self._wavetable, 1009 self._throughputtable) 1017 1010 else: 1018 1011 oldasc = False 1019 rev = N.interp(newwave, tapered._wavetable[::-1],1020 tapered._throughputtable[::-1])1012 rev = N.interp(newwave,self._wavetable[::-1], 1013 self._throughputtable[::-1]) 1021 1014 ans = rev[::-1] 1022 1015
