diff --git a/laghos.cpp b/laghos.cpp
index 68515e55..2c56452e 100644
--- a/laghos.cpp
+++ b/laghos.cpp
@@ -321,8 +321,8 @@ int main(int argc, char *argv[])
    // Configure the device from the command line options
    Device backend;
    backend.Configure(device, dev);
-   if (Mpi::Root()) { backend.Print(); }
    backend.SetGPUAwareMPI(gpu_aware_mpi);
+   if (Mpi::Root()) { backend.Print(); }
 
 #ifdef LAGHOS_USE_DEVICE_UMPIRE
    // Warm up the device pool on the configured GPU to avoid first-use latency.
diff --git a/laghos_solver.cpp b/laghos_solver.cpp
index caa0ba8e..a6f2ded0 100644
--- a/laghos_solver.cpp
+++ b/laghos_solver.cpp
@@ -222,7 +222,7 @@ LagrangianHydroOperator::LagrangianHydroOperator(const int size,
 
    // Values of rho0DetJ0 and Jac0inv at all quadrature points.
    // Initial local mesh size (assumes all mesh elements are the same).
-   int Ne, ne = NE;
+   HYPRE_BigInt Ne, ne = NE;
    double Volume, vol = 0.0;
    if (dim > 1 && p_assembly)
    {
@@ -249,7 +249,7 @@ LagrangianHydroOperator::LagrangianHydroOperator(const int size,
       for (int e = 0; e < NE; e++) { vol += pmesh->GetElementVolume(e); }
    }
    MPI_Allreduce(&vol, &Volume, 1, MPI_DOUBLE, MPI_SUM, pmesh->GetComm());
-   MPI_Allreduce(&ne, &Ne, 1, MPI_INT, MPI_SUM, pmesh->GetComm());
+   MPI_Allreduce(&ne, &Ne, 1, HYPRE_MPI_BIG_INT, MPI_SUM, pmesh->GetComm());
    switch (pmesh->GetElementBaseGeometry(0))
    {
       case Geometry::SEGMENT: qdata.h0 = Volume / Ne; break;