hpx spatial model integration
@@ -387,6 +387,11 @@ class PointSpatialModel(SpatialModel):
Predicted flux map
"""
geom_image = geom.to_image()
+ if geom.is_hpx:
+ idx, weights = geom_image.interp_weights({"skycoord": self.position})
+ data = np.zeros(geom_image.data_shape)
+ data[tuple(idx)] = weights
+ else:
x, y = geom_image.get_pix()
x0, y0 = self.position.to_pixel(geom.wcs)
data = self._grid_weights(x, y, x0, y0)
