diff --git a/src/scheme/FiniteVolumesDiffusion.hpp b/src/scheme/FiniteVolumesDiffusion.hpp
index 548e87c0c0fda773b42ebcf75c12495c3a030cc6..4f66bcb3bc627153a73c78a0bc013e2dae20c7e9 100644
--- a/src/scheme/FiniteVolumesDiffusion.hpp
+++ b/src/scheme/FiniteVolumesDiffusion.hpp
@@ -279,7 +279,7 @@ public:
 	  momentum_fluxes +=  Fl(l)*Cjr(j,R);
 	  energy_fluxes   += (Gl(l), Cjr(j,R));
 	}
-	uj[j] -= std::exp(-t)*dt*Vj(j)*Sj(j)+(dt*inv_mj[j]) * momentum_fluxes;
+	uj[j] -= std::exp(-t)*(dt*inv_mj[j])*Vj(j)*Sj(j) + (dt*inv_mj[j]) * momentum_fluxes;
        	Ej[j] -= (dt*inv_mj[j]) * energy_fluxes;
       });