What is the relationship between the computational domain range, computational grid size, and sub-fault size to ensure stable simulation results?