A Computational Method for Solving the Stochastic Joint Replenishment Problem in High Dimensions