Skip to content

Add timing info to DG local advection prom example #251

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

Merged
ckendrick merged 4 commits into from
Dec 18, 2023
Merged

Add timing info to DG local advection prom example #251

Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
b6b8e57
Add assemble timer and timing prints
ckendrick Nov 15, 2023
d2a4a70
Fix issue with -online_interp when running in parallel
ckendrick Nov 20, 2023
73ff70a
Update timing print
ckendrick Nov 30, 2023
b975105
Fix minor typo in examples
ckendrick Dec 7, 2023
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
26 changes: 19 additions & 7 deletions examples/prom/dg_advection_local_rom_matrix_interp.cpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -32,7 +32,7 @@
// dg_advection_local_rom_matrix_interp -interp_prep -ff 1.08 -rdim 40
// dg_advection_local_rom_matrix_interp -fom -ff 1.05
// dg_advection_local_rom_matrix_interp -online_interp -ff 1.05 -rdim 40 (interpolate using a linear solve)
// dg_advection_local_rom_matrix_interp -online_interp -ff 1.05 -rdim 40 -im "LP" (interpolate using lagragian polynomials)
// dg_advection_local_rom_matrix_interp -online_interp -ff 1.05 -rdim 40 -im "LP" (interpolate using lagrangian polynomials)
// dg_advection_local_rom_matrix_interp -online_interp -ff 1.05 -rdim 40 -im "IDW" (interpolate using inverse distance weighting)
//
// Sample runs:
Expand Down Expand Up @@ -527,7 +527,9 @@ int main(int argc, char *argv[])
VectorFunctionCoefficient velocity(dim, velocity_function);
FunctionCoefficient inflow(inflow_function);
FunctionCoefficient u0(u0_function);
StopWatch fom_timer, assemble_timer;

assemble_timer.Start();
ParBilinearForm *m = new ParBilinearForm(fes);
ParBilinearForm *k = new ParBilinearForm(fes);
if (pa)
Expand Down Expand Up @@ -574,6 +576,8 @@ int main(int argc, char *argv[])
u->ProjectCoefficient(u0);
HypreParVector *U = u->GetTrueDofs();

assemble_timer.Stop();

{
ostringstream mesh_name, sol_name;
mesh_name << "dg_advection_local_rom_matrix_interp-mesh." << setfill('0') <<
Expand Down Expand Up @@ -681,9 +685,7 @@ int main(int argc, char *argv[])
}
}

StopWatch fom_timer;
double t = 0.0;

int max_num_snapshots = t_final / dt + 1;
bool update_right_SV = false;
bool isIncremental = false;
Expand Down Expand Up @@ -714,6 +716,8 @@ int main(int argc, char *argv[])

if (online)
{
assemble_timer.Start();

if (!online_interp)
{
CAROM::BasisReader reader(basisName);
Expand Down Expand Up @@ -805,14 +809,11 @@ int main(int argc, char *argv[])
std::vector<CAROM::Matrix*> M_hats;
std::vector<CAROM::Vector*> b_hats;
std::vector<CAROM::Vector*> u_init_hats;
std::ofstream fout;
fout.open("frequencies.txt");
for(auto it = frequencies.begin(); it != frequencies.end(); it++)
{
CAROM::Vector* point = new CAROM::Vector(1, false);
point->item(0) = *it;
parameter_points.push_back(point);
fout << *it << std::endl;

std::string parametricBasisName = "basis_" + std::to_string(*it);
CAROM::BasisReader reader(parametricBasisName);
Expand All @@ -839,7 +840,6 @@ int main(int argc, char *argv[])
parametricuinithat->read("u_init_hat_" + std::to_string(*it));
u_init_hats.push_back(parametricuinithat);
}
fout.close();
if (myid == 0) printf("spatial basis dimension is %d x %d\n", numRowRB,
numColumnRB);

Expand Down Expand Up @@ -884,6 +884,8 @@ int main(int argc, char *argv[])

u_in = new Vector(numColumnRB);
*u_in = 0.0;

assemble_timer.Stop();
}

TimeDependentOperator* adv;
Expand Down Expand Up @@ -1037,6 +1039,16 @@ int main(int argc, char *argv[])
osol.close();
}

// 13. print timing info
if (myid == 0)
{
const std::string type = online ? "ROM" : "FOM";
std::cout << "Elapsed time for assembling " << type << ": " << std::scientific
<< std::setprecision(8) << assemble_timer.RealTime() << " second\n";
std::cout << "Elapsed time for solving " << type << ": " << std::scientific <<
std::setprecision(8) << fom_timer.RealTime() << " second\n";
}

// 16. Free the used memory.
delete U;
delete u;
Expand Down