SciMLBook icon indicating copy to clipboard operation
SciMLBook copied to clipboard

What is E in Chapter 9 - Adaptive Time Stepping?

Open bradcarman opened this issue 1 year ago • 1 comments

I'm assuming $\text{E}$ from here is the residual function value, i.e. $g(u_i) = \text{E}$. Is this correct? And is it scalarized or in vector form?

Does the relative tolerance $\tau_r$ ever take into account the scaling of $u$? For example, a hydraulic problem that has pressure values in the order of 1e7 and valve displacements in the order of 1e-6. If $\text{E}$ is in vector form then is $q$ computed as norm( E ./ (tau_r*u + tau_a) ) ?

bradcarman avatar Dec 29 '23 02:12 bradcarman

No, E is an error estimate between two different methods. It's method-dependent, but generally the difference between two different stepping techniques at different orders.

If is in vector form then is computed as norm( E ./ (tau_r*u + tau_a) ) ?

Yes https://docs.sciml.ai/DiffEqDocs/stable/extras/timestepping/

I'm going to write some notes for the new course that develop this a lot more.

ChrisRackauckas avatar Dec 29 '23 02:12 ChrisRackauckas