Skip to content

Commit 79063bb

Browse files
committed
Added Q_D coefficient update to IMEX sweepers
1 parent cf4bd32 commit 79063bb

File tree

3 files changed

+5
-10
lines changed

3 files changed

+5
-10
lines changed

pySDC/implementations/sweeper_classes/imex_1st_order.py

Lines changed: 2 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -61,16 +61,14 @@ def update_nodes(self):
6161
Returns:
6262
None
6363
"""
64-
65-
# get current level and problem description
6664
L = self.level
6765
P = L.prob
66+
M = self.coll.num_nodes
6867

6968
# only if the level has been touched before
7069
assert L.status.unlocked
7170

72-
# get number of collocation nodes for easier access
73-
M = self.coll.num_nodes
71+
self.updateVariableCoeffs(L.status.sweep)
7472

7573
# gather all terms which are known already (e.g. from the previous iteration)
7674
# this corresponds to u0 + QF(u^k) - QIFI(u^k) - QEFE(u^k) + tau

pySDC/implementations/sweeper_classes/imex_1st_order_MPI.py

Lines changed: 1 addition & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -43,14 +43,13 @@ def update_nodes(self):
4343
Returns:
4444
None
4545
"""
46-
4746
L = self.level
4847
P = L.prob
4948

5049
# only if the level has been touched before
5150
assert L.status.unlocked
5251

53-
# get number of collocation nodes for easier access
52+
self.updateVariableCoeffs(L.status.sweep)
5453

5554
# gather all terms which are known already (e.g. from the previous iteration)
5655
# this corresponds to u0 + QF(u^k) - QdF(u^k) + tau

pySDC/implementations/sweeper_classes/imex_1st_order_mass.py

Lines changed: 2 additions & 4 deletions
Original file line numberDiff line numberDiff line change
@@ -15,16 +15,14 @@ def update_nodes(self):
1515
Returns:
1616
None
1717
"""
18-
19-
# get current level and problem description
2018
L = self.level
2119
P = L.prob
20+
M = self.coll.num_nodes
2221

2322
# only if the level has been touched before
2423
assert L.status.unlocked
2524

26-
# get number of collocation nodes for easier access
27-
M = self.coll.num_nodes
25+
self.updateVariableCoeffs(L.status.sweep)
2826

2927
# gather all terms which are known already (e.g. from the previous iteration)
3028
# this corresponds to u0 + QF(u^k) - QIFI(u^k) - QEFE(u^k) + tau

0 commit comments

Comments
 (0)