Adapt tests for new numpy 1.25 dtype behaviour

... sometimes resulting in float32 arrays instead of float64; not quite clear why this happens, but it's not our responsibility to do dtype coercion.