Skip to content

Commit 01b98db

Browse files
committed
advection_diffusion: Save errors in a map: step, iter -> error.
1 parent 53c60b2 commit 01b98db

File tree

1 file changed

+7
-2
lines changed

1 file changed

+7
-2
lines changed

examples/advection_diffusion/advection_diffusion_sweeper.hpp

Lines changed: 7 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -5,10 +5,11 @@
55
#ifndef _ADVECTION_DIFFUSION_SWEEPER_HPP_
66
#define _ADVECTION_DIFFUSION_SWEEPER_HPP_
77

8-
#include <complex>
9-
#include <vector>
108
#include <cassert>
9+
#include <complex>
10+
#include <map>
1111
#include <ostream>
12+
#include <vector>
1213

1314
#include <pfasst/encap/imex_sweeper.hpp>
1415

@@ -29,6 +30,7 @@ class AdvectionDiffusionSweeper
2930
FFT fft;
3031

3132
vector<complex<double>> ddx, lap;
33+
map<pair<size_t, size_t>, double> errors;
3234

3335
double v = 1.0;
3436
time t0 = 1.0;
@@ -95,6 +97,9 @@ class AdvectionDiffusionSweeper
9597
cout << "err: " << n << " " << k << " " << scientific << max
9698
<< " (" << qend->size() << ", " << predict << ")"
9799
<< endl;
100+
101+
errors.insert(pair<pair<size_t, size_t>, double>
102+
(pair<size_t, size_t>(n, k), max));
98103
}
99104

100105
void predict(bool initial)

0 commit comments

Comments
 (0)