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.

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

Add option of eigenmodes oversampling #303

Merged
merged 1 commit into from
Oct 29, 2024
Merged

Add option of eigenmodes oversampling #303

Changes from all commits
Commits
Show all changes
1 commit
Select commit
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
35 changes: 21 additions & 14 deletions examples/prom/elliptic_eigenproblem_global_rom.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -93,6 +93,7 @@ int main(int argc, char *argv[])
bool dirichlet = true;
int order = 1;
int nev = 4;
int nev_os = 0;
int seed = 75;
bool prescribe_init = false;
int lobpcg_niter = 200;
Expand Down Expand Up @@ -133,6 +134,8 @@ int main(int argc, char *argv[])
"Order (degree) of the finite elements.");
args.AddOption(&nev, "-n", "--num-eigs",
"Number of desired eigenmodes.");
args.AddOption(&nev_os, "-nos", "--num-eigs-os",
"Number of oversampled eigenmodes.");
args.AddOption(&seed, "-s", "--seed",
"Random seed used to initialize LOBPCG.");
args.AddOption(&amplitude, "-a", "--amplitude",
Expand Down Expand Up @@ -172,7 +175,7 @@ int main(int argc, char *argv[])
"Number of iterations for the LOBPCG solver.");
args.AddOption(&lobpcg_tol, "-tol", "--fom-tol",
"Tolerance for the LOBPCG solver.");
args.AddOption(&eig_tol, "-tol", "--fom-tol",
args.AddOption(&eig_tol, "-eig-tol", "--eig-tol",
"Tolerance for eigenvalues to be considered equal.");
#ifdef MFEM_USE_SUPERLU
args.AddOption(&slu_solver, "-slu", "--superlu", "-no-slu",
Expand Down Expand Up @@ -298,7 +301,7 @@ int main(int argc, char *argv[])
// 8. Set BasisGenerator if offline
if (offline)
{
options = new CAROM::Options(fespace->GetTrueVSize(), nev, nev,
options = new CAROM::Options(fespace->GetTrueVSize(), nev,
update_right_SV);
std::string snapshot_basename = baseName + "par" + std::to_string(id);
generator = new CAROM::BasisGenerator(*options, isIncremental,
Expand All @@ -309,7 +312,7 @@ int main(int argc, char *argv[])
if (merge)
{
mergeTimer.Start();
options = new CAROM::Options(fespace->GetTrueVSize(), max_num_snapshots, nev,
options = new CAROM::Options(fespace->GetTrueVSize(), max_num_snapshots,
dylan-copeland marked this conversation as resolved.
Show resolved Hide resolved
update_right_SV);
generator = new CAROM::BasisGenerator(*options, isIncremental, basis_filename);
for (int paramID=0; paramID<nsets; ++paramID)
Expand Down Expand Up @@ -457,7 +460,7 @@ int main(int argc, char *argv[])
}

lobpcg = new HypreLOBPCG(MPI_COMM_WORLD);
lobpcg->SetNumModes(nev);
lobpcg->SetNumModes(nev + nev_os);
lobpcg->SetRandomSeed(seed);
lobpcg->SetPreconditioner(*precond);
lobpcg->SetMaxIter(lobpcg_niter);
Expand All @@ -469,8 +472,8 @@ int main(int argc, char *argv[])

if (prescribe_init && (fom || (offline && id > 0)))
{
HypreParVector** snapshot_vecs = new HypreParVector*[nev];
for (int i = 0; i < nev; i++)
HypreParVector** snapshot_vecs = new HypreParVector*[nev + nev_os];
for (int i = 0; i < nev + nev_os; i++)
{
std::string snapshot_filename = baseName + "ref_snapshot_" + std::to_string(i);
std::ifstream snapshot_infile(snapshot_filename + "." + std::to_string(myid));
Expand All @@ -479,7 +482,7 @@ int main(int argc, char *argv[])
snapshot_vecs[i]->Read(MPI_COMM_WORLD, snapshot_filename.c_str());
if (myid == 0) std::cout << "Loaded " << snapshot_filename << std::endl;
}
lobpcg->SetInitialVectors(nev, snapshot_vecs);
lobpcg->SetInitialVectors(nev + nev_os, snapshot_vecs);
if (myid == 0) std::cout << "LOBPCG initial vectors set" << std::endl;
}

Expand All @@ -496,25 +499,29 @@ int main(int argc, char *argv[])
{
if (myid == 0)
{
std::cout << " Eigenvalue " << i << ": " << eigenvalues[i] << "\n";
std::cout << "Eigenvalue " << i << ": " << eigenvalues[i] << "\n";
}
if (offline)
{
eigenfunction_i = lobpcg->GetEigenvector(i);
eigenfunction_i /= sqrt(InnerProduct(eigenfunction_i, eigenfunction_i));
generator->takeSample(eigenfunction_i.GetData());
if (prescribe_init && id == 0)
}
}

if (offline)
{
if (prescribe_init && id == 0)
{
for (int i = 0; i < nev + nev_os; i++)
{
std::string snapshot_filename = baseName + "ref_snapshot_" + std::to_string(i);
const HypreParVector snapshot_vec = lobpcg->GetEigenvector(i);
snapshot_vec.Print(snapshot_filename.c_str());
if (myid == 0) std::cout << "Saved " << snapshot_filename << std::endl;
if (myid == 0) std::cout << "Saved " << snapshot_filename <<
" for LOBPCG initialization" << std::endl;
}
}
}

if (offline)
{
generator->writeSnapshot();
delete generator;
delete options;
Expand Down
Loading