Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

What is E in Chapter 9 - Adaptive Time Stepping? #137

Open
bradcarman opened this issue Dec 29, 2023 · 1 comment
Open

What is E in Chapter 9 - Adaptive Time Stepping? #137

bradcarman opened this issue Dec 29, 2023 · 1 comment
Labels
enhancement New feature or request

Comments

@bradcarman
Copy link

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 bradcarman added the enhancement New feature or request label Dec 29, 2023
@ChrisRackauckas
Copy link
Member

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.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
enhancement New feature or request
Projects
None yet
Development

No branches or pull requests

2 participants