display state.xn in scientific notation

integration/integrator_setup_strang.H

@@ -210,7 +210,7 @@ void integrator_cleanup (IntegratorT& int_state, BurnT& state,
             std::cout << "temp start = " << std::setprecision(16) << state_save.T_in << std::endl;
             std::cout << "xn start = ";
             for (const double X : state_save.xn_in) {
-                std::cout << std::setprecision(16) << X << " ";
+                std::cout << std::scientific << std::setprecision(16) << X << " ";
             }
             std::cout << std::endl;
             std::cout << "dens current = " << std::setprecision(16) << state.rho << std::endl;