p = np.polyfit(np.log(simDFNew[xField]), np.log(timeForceSPH), deg=1)