I think so, yes, this is likely 1d.. Perhaps you can use relationships between 1d and 3d, something along these lines: 1d index=i×m×p+j×p+k and back to 3d i=index/(m×p) j=(index mod(m×p))/p k=index mod p where n, m, p are your 3d map dimensions. This is how we'd work with 3d maps in Fortran/C++ as if was 1d, for efficiency. Pavel On 1/20/25 09:39, [email protected] wrote:
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 _______________________________________________ cctbxbb mailing list -- [email protected] To unsubscribe send an email to [email protected]