Sampling from a truncated multivariate normal distribution subject to multiple linear inequality constraints is a recurring problem in many areas in statistics and econometrics, such as the order restricted regressions, censored data models, and shape-restricted nonparametric regressions. However, the sampling problem still appears non-trivial due to the existence of the analytically intractabl...
