Computation of topological phase diagram of disordered Pb1-x Snx Te using the kernel polynomial method