Incorporate compute_expansion()

[seasmith]

Currently compute_expansion() is a placeholder function to remove expansion computation from compute_limits(), and currently commented out.