Efficient expectation propagation for posterior approximation in high-dimensional probit models