Hi All,

I am Yuji from Japan.

I tyr to get some benchmarks using GPU on cube_hex24.msh made by Dr. Park. Howenver I get some problems.

When singel calclation is going, the terminal shows how much this calclation is done. like below;

```
$ pyfr run -b cuda -p cube_hex24.pyfrm config.ini
Authorization required, but no authorization protocol specified
Authorization required, but no authorization protocol specified
Authorization required, but no authorization protocol specified
100.0% [===========================================================================>] 3.00/3.00 ela: 00:02:35 rem: 00:00:00
$
```

When more than 2mpi calclation are going, the terminal does not show how much this calclation is done. like below;

```
$mpiexec -n 2 pyfr run -b cuda cube_hex24_p2.pyfr mconfig.ini
Authorization required, but no authorization protocol specified
Authorization required, but no authorization protocol specified
$
```

Do you know why it occurs?

I attach my config.ini below

```
[backend]
precision = double
rank-allocator = linear
[backend-openmp]
;cblas-st = Enter path to local single-threaded BLAS library for OpenMP backend
;cblas-mt = Enter path to local multi-threaded BLAS library for OpenMP backend
[constants]
gamma = 1.4
mu = 0.001
Pr = 0.72
; Convergence paramters
k = 3.14159265358979323846
omega = 3.14159265358979323846
a = 3
[solver]
system = navier-stokes
order = 4
;anti-alias = flux, surf-flux
[solver-time-integrator]
scheme = rk4
controller = none
t0 = 0.0
dt = 2e-4
tstart = 0.0
tend = 3.0
[solver-interfaces]
riemann-solver = rusanov
ldg-beta = 0.5
ldg-tau = 0.1
[solver-interfaces-tri]
flux-pts = williams-shunn
[solver-interfaces-quad]
flux-pts = gauss-legendre
[solver-elements-tet]
soln-pts = shunn-ham
[solver-elements-hex]
soln-pts = gauss-legendre
;[soln-output]
;format = pyfrs-dir
;basedir = .
;basename = run-%(t).2f
;times = range(0.0, 0.1, 2)
[soln-plugin-writer]
dt-out = 0.1
basedir = .
basename = runout-{t:.2f}
[soln-ics]
f = (sin(k*(x + y + z)) + a)
rho = %(f)s
u = 1
v = 1
w = 1
p = (gamma - 1)*(%(f)s*%(f)s - 1.5*%(f)s)
[solver-source-terms]
theta = (k*(x + y + z) - omega*t)
rho = (3.0*k - omega)*cos(%(theta)s)
rhou = (-omega + 3.0*k + 2.0*(gamma-1.0)*a*k - 1.5*k*(gamma-1.0))*cos(%(theta)s) + k*(gamma-1.0)*sin(2.0*%(theta)s)
rhov = (-omega + 3.0*k + 2.0*(gamma-1.0)*a*k - 1.5*k*(gamma-1.0))*cos(%(theta)s) + k*(gamma-1.0)*sin(2.0*%(theta)s)
rhow = (-omega + 3.0*k + 2.0*(gamma-1.0)*a*k - 1.5*k*(gamma-1.0))*cos(%(theta)s) + k*(gamma-1.0)*sin(2.0*%(theta)s)
; Euler prolbem without last term of E
E = (6.0*k*a*gamma - 4.5*(gamma-1.0)*k - 2.0*omega*a)*cos(%(theta)s) + (3.0*gamma*k - omega)*sin(2.0*%(theta)s) + sin(%(theta)s)*(3.0*k*k*mu*gamma/Pr)
```

Best,

Yuji