-
Notifications
You must be signed in to change notification settings - Fork 166
Rckirby/timedep rb #4339
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
base: release
Are you sure you want to change the base?
Rckirby/timedep rb #4339
Conversation
demos/time_dependent_rayleigh_benard/timedep-rayleigh-benard.py.rst
Outdated
Show resolved
Hide resolved
demos/time_dependent_rayleigh_benard/timedep-rayleigh-benard.py.rst
Outdated
Show resolved
Hide resolved
demos/time_dependent_rayleigh_benard/timedep-rayleigh-benard.py.rst
Outdated
Show resolved
Hide resolved
demos/time_dependent_rayleigh_benard/timedep-rayleigh-benard.py.rst
Outdated
Show resolved
Hide resolved
mh = MeshHierarchy(base_msh, ref_levels) | ||
msh = mh[-1] | ||
|
||
MC = MeshConstant(msh) |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
This needs a little bit of explanation. I think the only reason to use the R space is to be able to take the adjoint of a simulation for variables that appear inside a Form
and are updated at each time step. This example might be fine with just Constant
.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I thought we had UFL manipulation errors internally with Constant
?
Also, t += dt
works great for MC
but not for `Constant.
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
I'll try this out and see if it passes CI. Worked locally for me.
Description
This is a demo solving the time-dependent Rayleigh-Benard convection problem using Irksome for time stepping and a monolithic multigrid solver for the algebraic systems.