Solving High-dimensional Shape-constrained Convex Regression Problems