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

Energy-conserving EQP #177

Draft
wants to merge 39 commits into
base: rom-eqp
Choose a base branch
from
Draft

Energy-conserving EQP #177

wants to merge 39 commits into from

Conversation

cval26
Copy link
Collaborator

@cval26 cval26 commented Jun 28, 2023

I'm opening this draft PR so that you can stay up to date with my changes in the code, and so that I can ask you questions about specific parts of the code that I'm unsure of. I'll add comments to the parts of the code I have questions on.

I'll push new changes as I implement them so that we can discuss them. Once we have a correct implementation I can go back and edit the commit history so that it's easier to review.

The indentation in some parts of the code seems to be off in the diff files. I'm not sure why; it doesn't look like that in my editor.

Changes made so far:

  • formed a new, combined NNLS problem for both velocity and energy;
  • formed the new RHS force vectors used in the online timestepping phase;
  • added the energy identity vector to each time window's reduced energy basis;

@cval26 cval26 changed the title Implementation of energy-conserving EQP Energy-conserving EQP Jun 28, 2023
Vector & rv, Vector & re)
{
const int nqp = ir->GetNPoints();
const int dim = trial_fe.GetDim(); // TODO: shouldn't it be the dimension of test_fe?
Copy link
Collaborator Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Shouldn't dim be the dimension of the velocity space test_fe, instead of the energy space trial_fe? In general, my understanding is that test_fe has dimension 3, while trial_fe has dimension 1.

Variable dim is used to set the size of several vectors holding velocity coefficient values, so I'm not sure why the dimension of trial_fe is used.

rom/laghos_rom.cpp Outdated Show resolved Hide resolved
rom/laghos_solver.cpp Outdated Show resolved Hide resolved
rom/laghos.cpp Outdated Show resolved Hide resolved
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

3 participants