Skip to content

Energy-conserving EQP #177

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.

Sign up for GitHub

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

Draft
cval26 wants to merge 40 commits into
base: rom
Choose a base branch
from
Draft

Energy-conserving EQP #177

Show file tree
Hide file tree
Changes from 4 commits
Commits
Show all changes
40 commits
Select commit Hold shift + click to select a range
c7edc36
form combined NNLS problem
cval26 Jun 27, 2023
b78829d
modify read/write of NNLS solution
cval26 Jun 28, 2023
c4280af
form RHS force vectors
cval26 Jun 30, 2023
0da7e93
modify parameter flags
cval26 Jun 30, 2023
90327af
add energy identity to basis
cval26 Jul 6, 2023
0b3907e
review 1 changes
cval26 Jul 7, 2023
9b9b0a7
print energy total & diff
cval26 Jul 11, 2023
3d661d2
review 2 changes
cval26 Jul 11, 2023
36b2232
form inverse reduced mass matrices
cval26 Jul 12, 2023
5adac53
Merge branch 'rom-eqp' into rom-eqp-en
siuwuncheung Jul 12, 2023
684a039
form RHSs for timemarching
cval26 Jul 12, 2023
7cac132
bug fixes & minor changes
cval26 Jul 14, 2023
679c952
fix energy identity mistake
cval26 Jul 25, 2023
c2c7fd5
fix variables offset
cval26 Jul 26, 2023
4ade861
fix restore issue
cval26 Jul 31, 2023
834de51
add NNLS tolerance input option
cval26 Jul 28, 2023
fa591bd
add gitignore file
cval26 Jul 18, 2023
124dd00
use libROM matrix rescaling method
cval26 Jul 31, 2023
3fa45ef
sampling frequency; energy calc every dt
cval26 Aug 1, 2023
4105bae
restore sampling call flags
cval26 Aug 4, 2023
8f5ddc3
minor matrix allocation change
cval26 Aug 4, 2023
cd57d9d
fix online loop timing
cval26 Aug 11, 2023
eda3d08
windowing: online stage changes
cval26 Sep 4, 2023
5e95bae
windowing: offline stage changes
cval26 Sep 4, 2023
85948c9
form ROM basis/operator at time of window change
cval26 Sep 26, 2023
593497e
fix basis formation in restore phase
cval26 Sep 26, 2023
09e7268
add Dylan's EQP bug fixes
cval26 Oct 20, 2023
b0e6d08
fix L2 errors calculation
cval26 Nov 3, 2023
473c0e3
use incremental GS method
cval26 Nov 3, 2023
242155c
use snapshots in EQP derivation
cval26 Nov 3, 2023
cf55f61
snapshots matrix bug fix
cval26 Dec 7, 2023
90b38be
basis orthonormality tests
cval26 Dec 7, 2023
ea2cdd5
Dylan's quadrature data bug fix
cval26 Dec 20, 2023
317ca53
use double-pass GS
cval26 Jan 22, 2024
332d59c
perform LQ decomposition of NNLS matrix
cval26 Feb 6, 2024
42f9b0c
add offset vectors to basis matrices
cval26 Mar 4, 2024
e3a7be9
update libROM sampling func. calls
cval26 Mar 25, 2024
9c8294e
use mass induced Gram-Schmidt
cval26 May 21, 2024
63990ad
get L directly from libROM's LQ
cval26 Jun 25, 2024
298da43
include rom makefile changes
cval26 Jul 24, 2024
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
21 changes: 18 additions & 3 deletions rom/laghos.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -394,8 +394,15 @@ int main(int argc, char *argv[])
romOptions.basisIdentifier = std::string(basisIdentifier);

romOptions.hyperreductionSamplingType = getHyperreductionSamplingType(hyperreductionSamplingType);
romOptions.use_sample_mesh = romOptions.hyperreduce && (romOptions.hyperreductionSamplingType != eqp);
MFEM_VERIFY(!(romOptions.SNS) || (romOptions.hyperreductionSamplingType != eqp), "Using SNS with EQP is prohibited");
romOptions.use_sample_mesh = romOptions.hyperreduce && (romOptions.hyperreductionSamplingType != eqp
|| romOptions.hyperreductionSamplingType != eqp_energy);

MFEM_VERIFY(!(romOptions.SNS) || (romOptions.hyperreductionSamplingType != eqp &&
romOptions.hyperreductionSamplingType != eqp_energy),
"Using SNS with EQP is prohibited");

if (romOptions.hyperreduce && romOptions.hyperreductionSamplingType == eqp_energy)
romOptions.GramSchmidt = true;

romOptions.spaceTimeMethod = getSpaceTimeMethod(spaceTimeMethod);
const bool spaceTime = (romOptions.spaceTimeMethod != no_space_time);
Expand Down Expand Up @@ -1573,7 +1580,15 @@ int main(int argc, char *argv[])
{
romOper[0]->ApplyHyperreduction(romS);
}
double windowEndpoint = 0.0;

// To perform the induced Gram-Schmidt, so that the reduced mass
// matrices are identity.
if (romOptions.hyperreduce && romOptions.hyperreductionSamplingType == eqp_energy)
{
romOper[0]->ApplyHyperreduction(romS);
}

double windowEndpoint = 0.0;
double windowOverlapMidpoint = 0.0;
for (ti = 1; !last_step; ti++)
{
Expand Down
Loading