Using ligo.skymap.moc.rasterize ends up with Aborted (core dumped)
Hi,
I am a relatively new user of ligo.skymaps and I am trying to take transform a moc skyampap into a nested probability map of nside 64 (order 6). I try to use rectorize function for this in the script below. If I use order more than 2 the it crashes with: munmap_chunk(): invalid pointer Aborted (core dumped)
import numpy as np
from ligo.skymap.moc import rasterize
from astropy.table import Table
skymap = Table.read('https://gracedb.ligo.org/api/superevents/S231226av/files/Bilby.multiorder.fits')
moc_data = np.array([np.array(skymap['UNIQ']),np.array(skymap['PROBDENSITY'])], dtype=[('UNIQ', int), ('PROBDENSITY', float)])
P = rasterize(moc_data, order=2)
Edited by Leo P. Singer