nhmc / barak Goto Github PK

-with open('README') as fh:
+with open('README.rst') as fh:

import barak
from barak import convolve

low = 7040.
high = 7050.
subpixkms = 1.3 # km/s
instrument_profile = 2.5 # km/s

wave = barak.sed.make_constant_dv_wa_scale(low, high, subpixkms)
flux = np.ones_like(wave)

# This works as expected:
convolved_flux = convolve.convolve_constant_dv(wave, flux, wa_dv=subpixkms, vfwhm=instrument_profile)

#This one returns an error: 
convolved_flux_fast = convolve.convolve_constant_dv(wave, flux, wa_dv=subpixkms, npix=4.0)

#This is the error: 
---------------------------------------------------------------------------
TypeError                                 Traceback (most recent call last)
<ipython-input-1-343170092fbe> in <module>()
     11 
     12 convolved_flux = convolve.convolve_constant_dv(wave, flux, wa_dv=subpixkms, vfwhm=instrument_profile)
---> 13 convolved_flux_fast = convolve.convolve_constant_dv(wave, flux, wa_dv=subpixkms, npix=4.0)

/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/barak/convolve.pyc in convolve_constant_dv(wa, fl, wa_dv, npix, vfwhm)
    137         wa_dv = make_constant_dv_wa_scale(wa[0], wa[-1], float(vfwhm)/npix)
    138     fl_dv = np.interp(wa_dv, wa, fl)
--> 139     fl_dv_smoothed = convolve_psf(fl_dv, npix)
    140     fl_out = np.interp(wa, wa_dv, fl_dv_smoothed)
    141     return fl_out

/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/barak/convolve.pyc in convolve_psf(a, fwhm, edge, replace_nan, debug)
     44     gauss = np.exp(-0.5 * (x / sigma) ** 2 )
     45 
---> 46     return convolve_window(a, gauss, edge=edge)
     47 
     48 def convolve_window(a, window, edge='invert'):

/Library/Frameworks/EPD64.framework/Versions/7.2/lib/python2.7/site-packages/barak/convolve.pyc in convolve_window(a, window, edge)
     80     # Add edges to either end of the array to reduce edge effects in
     81     # the convolution.
---> 82     if len(a) < 2*n:
     83         raise ValueError('Window is too big for the array! %i %i' % (len(a), n))
     84     if edge == 'invert':

TypeError: object of type 'numpy.float64' has no len()

In [53]: starburst_Calzetti00(3669)
---------------------------------------------------------------------------
ValueError                                Traceback (most recent call last)
<ipython-input-53-660921431fc5> in <module>()
----> 1 starburst_Calzetti00(3669)

/usr/local/lib/python2.7/site-packages/barak/extinction.pyc in starburst_Calzetti00(wa, Rv, EBmV)
    218         if c1.any():
    219             AlamAv[c0] = MW_Cardelli89(wa[c0], 3.1).AlamAv
--> 220             AlamAv[c1] = interp_Akima(wa[c1], wa[~c1], AlamAv[~c1])
    221
    222     if len(wa) == 1:

/usr/local/lib/python2.7/site-packages/barak/interp.pyc in interp_Akima(x_new, x, y)
    344     >>> plt.show()
    345     """
--> 346     interpolator = AkimaSpline(x, y)
    347     return interpolator(x_new)
    348

/usr/local/lib/python2.7/site-packages/barak/interp.pyc in __init__(self, xvals, yvals)
    173         n = len(x)
    174         if n < 3:
--> 175             raise ValueError("Array too small")
    176         if n != len(y):
    177             raise ValueError("Size of x-array must match data shape")

ValueError: Array too small