Multi-dimensional High-order Discretizations and Fast Solution Algorithms for Simplified Spherical Harmonic Approximations