def interpolate_to_mesh_indices(self, loop):
"""
Return interpolated loop indices to the temperature and density meshes defined for
the atomic data. For use with `~scipy.ndimage.map_coordinates`.
"""
nots_itemperature = splrep(self.temperature.value, np.arange(self.temperature.shape[0]))
nots_idensity = splrep(self.density.value, np.arange(self.density.shape[0]))
itemperature = splev(np.ravel(loop.electron_temperature.value), nots_itemperature)
idensity = splev(np.ravel(loop.density.value), nots_idensity)
return itemperature, idensity
评论列表
文章目录