Make all larry methods able to handle arbitrary dimensional (nd) input

The larry methods below currently only work on 2d larrys. Extend to work on nd larrys (default input axis is shown):

push, axis=-1
movingrank, axis=-1
movingsum_forward, axis=-1
quantile, axis=-1
last_rank, axis=-1
last_rank_decay, axis=-1
group_ranking, axis=0

and (not yet a method):

la.afunc.geometric_mean