potential issue with interpolation?

[caseylam]

Hi,

I've found that when I use "get_eep", it fails to return an EEP where I think it should have no problem: for example, mist.get_eep(0.15, 10, 0, accurate=True) gives me RuntimeError: EEP minimization not successful: (0.15, 10, 0).

As a test, I ran get_eep on a grid of masses, logages, and Fe/H (mass_array = 10**np.linspace(np.log10(0.1), np.log10(20), 500), logage_array = np.arange(5, 10.2, 0.1), feh_array = np.arange(-3.5, 0.31, 0.1)).

I then plotted the resultant EEP of that for each value of Fe/H in a grid of mass vs age. If there is an EEP runtime error returned, that piece of the grid is white. The white region in the upper right makes sense (those stars are "dead" so they're not on the isochrones), but a lot of stars which I would expect to be on the MS are not there, especially at [Fe/H] > 0. There also looks to be some slightly weird behavior at low [Fe/H] around a solar mass (there is a "notch" in the distribution.)

Is this expected?