
20 Jan
2025
20 Jan
'25
5:39 p.m.
It's close but as far as I can tell from having tried using this function sel = maptbx.grid_indices_around_sites( unit_cell = unit_cell, fft_n_real = map_1.focus(), fft_m_real = map_1.all(), sites_cart = flex.vec3_double([site_cart]), site_radii = flex.double(1, radius)) the sel array is still a 1d array corresponding to 3d grid indices. Or am I wrong? Rob