Computing connolly surfaces