cdf_samples = np.random.uniform(0, 1, size=(100, 100))
cdf_samples *= (cdf_stop - cdf_start)[:, (None)]
cdf_samples += cdf_start[:, (None)]
truncnorm_samples = norm.ppf(cdf_samples)
Vectorised code for sampling from truncated normal distributions with different intervals