Hello everyone.

@Will, thanks, I will try using version `1.12.2`

. Unfortunately, I donâ€™t have web access on the cluster that I am using, therefore I cannot make a git diff between the two versions of the source code. However, adopting one brute force technique with `diff -ru`

of the two directories containing the two versions of the source code, I got a lot of differences (as I donâ€™t find a way to attach a file, I donâ€™t think that it would worth to write all of them here). On the other hand, doing `diff -u`

between the config files that I am used, I obtain the following:

```
--- PyFR_1.12/Vermeire_et_al_2017/sd7003_CompExprt/sd7003.ini 2021-07-02 09:17:50.000000000 +0200
+++ PyFR_1.12.1/Vermeire_et_al_2017/sd7003_new/sd7003.ini 2021-09-23 17:06:02.000000000 +0200
@@ -72,8 +72,32 @@
[soln-plugin-writer]
-dt-out=10.0
+dt-out=1.0
basedir = .
-basename = sd7003-{t:.5f}
+basename = sd7003-{t:.2f}
```

being `---`

the successful test with version `1.12.0`

and `+++`

the second failed test adopting version `1.12.1`

. Therefore, I am basically using the same config file, apart from the `dt`

of the output.

@Gonzalo, I used the same version I am using to run the simulation for genereating the mesh (i.e. `PyFR 1.12.1`

for the latest failed attempts, and `PyFR 1.12.0`

for the previous successful simulation). I did not explicitly set any tolerance for the linear elements but on both the failed and successful simulations, I used the same parameters that where present in the paper supplemental material. In particular, the latest failing config file, looks like that:

```
[backend]
precision = double
[backend-cuda]
device-id = local-rank
[constants]
gamma = 1.4
mu = 3.94405318873308E-6
Pr = 0.72
M = 0.2
[solver-time-integrator]
scheme = rk45
controller = pi
tstart=0.0
tend=260.0
dt = 0.00001
atol = 0.000001
rtol = 0.000001
safety-fact = 0.5
min-fact = 0.3
max-fact = 1.2
[solver]
system = navier-stokes
order = 4
[solver-interfaces]
riemann-solver = rusanov
ldg-beta = 0.5
ldg-tau = 0.1
[solver-interfaces-quad]
flux-pts = gauss-legendre
[solver-elements-hex]
soln-pts = gauss-legendre
[soln-bcs-outlet]
type = char-riem-inv
rho = 1.0
u = 0.2366431913
v = 0.0
w = 0.0
p = 1.0
[soln-bcs-inlet]
type = char-riem-inv
rho = 1.0
u = 0.2366431913
v = 0.0
w = 0.0
p = 1.0
[soln-bcs-wall]
type = no-slp-adia-wall
cpTw = 3.5
[soln-ics]
rho = 1.0
u = 0.2366431913
v = 0.001
w = 0.001*cos(x)*cos(y)
p = 1.0
[soln-plugin-writer]
dt-out=1.0
basedir = .
basename = sd7003-{t:.2f}
```