High-order compact operator splitting method for three-dimensional fractional equation with subdiffusion