Skip to content

Commit 9ed98b1

Browse files
authored
Update Time Integration lesson.md
Change step sizes for stability test
1 parent 768ca4f commit 9ed98b1

File tree

1 file changed

+5
-5
lines changed

1 file changed

+5
-5
lines changed

_lessons/time_integrators_sundials/lesson.md

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -236,25 +236,25 @@ mpirun -n 1 ./fcompare plt00001/ reference_solution/
236236
Notice that the computed solution error is rather small (the solution has
237237
magnitude $$\mathcal{O}(1)$$, so we hope for errors well below 0.1).
238238

239-
Now re-run this hands-on code using a larger time step size of 100.0,
239+
Now re-run this hands-on code using a larger time step size of 25.0,
240240

241241
```bash
242-
mpirun -n 1 ./HandsOn1.CUDA.exe inputs-1 fixed_dt=100.0
242+
mpirun -n 1 ./HandsOn1.CUDA.exe inputs-1 fixed_dt=25.0
243243
```
244244

245-
_see how much faster the code ran!_ However, if we check the accuracy of the
245+
_the code now runs 5x faster._ However, if we check the accuracy of the
246246
computed solution,
247247

248248
```bash
249249
mpirun -n 1 ./fcompare plt00001/ reference_solution/
250250
```
251251

252-
we see it has an incredibly large error (mine was $$\mathcal{O}(10^{284})$$).
252+
we see it has an incredibly large error (mine was $$\mathcal{O}(10^{98})$$).
253253

254254
{% include qanda
255255
question='What do you think happened?'
256256
answer='At this mesh size, the explicit algorithm is unstable for
257-
a time step size of 100, but is stable for a time step size of 5.' %}
257+
a time step size of 25, but is stable for a time step size of 5.' %}
258258

259259
Run the code a few more times, trying to identify the largest stable time step
260260
size.

0 commit comments

Comments
 (0)