
Fast hydrodynamics on heterogeneous many-core hardware 281
MATLAB
3
is used to solve the linear system at every stage and a compar-
ison is made between single- and unfiltered double-precision calculations. It
is shown in Figure 11.9 a) that without a filter, the single-precision calcu-
lations result in “blow-up” after which the solver fails just before 50 wave
periods of calculation time. However, in Figure 11.9 (b) it is demonstrated
that invoking a smoothening filter, cf. (11.44), stabilizes the accumulation of
round-off errors and the calculations continue and achieve reduced accuracy
compared to the computed double-precision results. Thus, it is confirmed ...