add trilinear readout in cic; rename cic.cic

cic.cic -> cic.paint.
add cic.readout.

We also need a k-space cic kernel.
add a gather function in domain.

gather can gather data in the exchanged(shuffled) layout.
three modes:
mean, sum and any.