Vectorise everything that can be
Should go through and check that anything which returns a scalar is vectorised in the bindings, so that there are no surprises in the numpy side
Should go through and check that anything which returns a scalar is vectorised in the bindings, so that there are no surprises in the numpy side