Skip to content

Switch to adapter-based checkpointing in the two-scale heat conduction tutorial #676

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
merged 5 commits into from
Oct 23, 2025
Merged

Switch to adapter-based checkpointing in the two-scale heat conduction tutorial #676

Changes from 2 commits
Commits
Show all changes
5 commits
Select commit Hold shift + click to select a range
cb858fc
Switch to DuMuX adapter checkpointing functionality
IshaanDesai Oct 21, 2025
e9365a1
Remove manual checkpointing of time as the adapter handles this inter…
IshaanDesai Oct 21, 2025
4afaba8
Checkpoint timestepping
IshaanDesai Oct 22, 2025
e05c405
Add changelog entry
IshaanDesai Oct 22, 2025
7d8a5d0
Move initialize() call before time settings
IshaanDesai Oct 22, 2025
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
22 changes: 8 additions & 14 deletions two-scale-heat-conduction/macro-dumux/appl/main.cc
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -173,9 +173,7 @@ int main(int argc, char **argv)
problem->applyInitialSolution(x);
auto xOld = x;

auto xCheckpoint = x;
double timeCheckpoint = 0.0;
int timeStepCheckpoint = 0;
int timeStepCheckpoint = 0;

// initialize the coupling data
std::vector<double> temperatures;
Expand Down Expand Up @@ -226,8 +224,11 @@ int main(int argc, char **argv)
// output every vtkOutputInterval time step
const int vtkOutputInterval = getParam<int>("TimeLoop.OutputInterval");

// initialize preCICE
couplingParticipant.initialize();
// initialize preCICE and the adapter checkpointing
if (runWithCoupling) {
couplingParticipant.initialize();
couplingParticipant.initializeCheckpoint(x, *gridVariables);
}

// time loop parameters
const auto tEnd = getParam<Scalar>("TimeLoop.TEnd");
Expand Down Expand Up @@ -272,9 +273,7 @@ int main(int argc, char **argv)
break;

// write checkpoint
if (couplingParticipant.requiresToWriteCheckpoint()) {
xCheckpoint = x;
timeCheckpoint = timeLoop->time();
if (couplingParticipant.writeCheckpointIfRequired()) {
timeStepCheckpoint = timeLoop->timeStepIndex();
}

Expand Down Expand Up @@ -357,14 +356,9 @@ int main(int argc, char **argv)
couplingParticipant.advance(dt);

// reset to checkpoint if not converged
if (couplingParticipant.requiresToReadCheckpoint()) {
x = xCheckpoint;
xOld = x;
timeLoop->setTime(timeCheckpoint, timeStepCheckpoint);

if (couplingParticipant.readCheckpointIfRequired()) {
// TODO: previousTimeStep might be more appropriate, last one could be small
timeLoop->setTimeStepSize(dt);
gridVariables->update(x);
gridVariables->advanceTimeStep();
continue;
}
Expand Down