Skip to content

Commit 3bed49c

Browse files
committed
boris: Add reevaluate method.
Signed-off-by: Matthew Emmett <[email protected]>
1 parent 6cf30c0 commit 3bed49c

File tree

1 file changed

+9
-2
lines changed

1 file changed

+9
-2
lines changed

examples/boris/boris_sweeper.hpp

Lines changed: 9 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -380,7 +380,7 @@ namespace pfasst
380380

381381
// building rules for Q_E and Q_I:
382382
// Q_E is striclty lower diagonal matrix with delta nodes of column index
383-
// Q_I is lower diagonal matrix with first row and column all zero and delta nodes of
383+
// Q_I is lower diagonal matrix with first row and column all zero and delta nodes of
384384
// column index minus one
385385
Matrix<time> qe_mat = Matrix<time>(nnodes, nnodes);
386386
qe_mat.fill(time(0.0));
@@ -423,7 +423,7 @@ namespace pfasst
423423
this->energy_evals.front() = this->energy_evals.back();
424424
}
425425

426-
virtual void evaluate(size_t m) override
426+
void evaluate(size_t m)
427427
{
428428
// Vector3d<scalar> vel; vel.fill(scalar(0.0));
429429
// Vector3d<scalar> B; B.fill(scalar(0.0));
@@ -448,6 +448,13 @@ namespace pfasst
448448
this->f_evals++;
449449
}
450450

451+
virtual void reevaluate(bool initial) override
452+
{
453+
for (size_t m = 0; m < this->get_quadrature()->get_num_nodes(); m++) {
454+
this->evaluate(m);
455+
}
456+
}
457+
451458
virtual void predict(bool initial) override
452459
{
453460
UNUSED(initial);

0 commit comments

Comments
 (0)