jacopo-chevallard / PyP-BEAGLE

4 stars 1 forks source link

bangs_photometry singular matrix #1

Closed eclake closed 9 years ago

eclake commented 9 years ago

@jacopo-chevallard - here's the bug I put in the e-mail (you can get issues to e-mail the other person by referencing with @username)

I'm running into a problem when running bangs_photometry.py:

The first band, model_sed.data.field(0) is an array of zeros

127 xdata = model_sed.data.field(i) / nanoJy

which means that

135 pdf = WeightedKDE.gaussian_kde(xdata, weights=probability)

gives an error because it's a singular matrix.

an example object for reproducing this problem can be found here: https://www.dropbox.com/s/sih6bj0of6do7ii/for_jacopo.zip?dl=0

jacopo-chevallard commented 9 years ago

@eclake are you now able to plot the marginal photometry?