TreePiece::adjust(): dEtauDot also limits adiabatic heating. (#282)

N-BodyShop · Dec 13, 2021 · bdf58f3 · bdf58f3
1 parent ee9fd39
commit bdf58f3
Showing 1 changed file with 2 additions and 1 deletion.
diff --git a/TreePiece.cpp b/TreePiece.cpp
@@ -1819,7 +1819,8 @@ void TreePiece::adjust(int iKickRung, int bEpsAccStep, int bGravStep,
           }
 #endif
 
-	  if (dEtauDot > 0.0 && p->PdV() < 0.0) { /* Prevent rapid adiabatic cooling */
+	  if (dEtauDot > 0.0) { /* Prevent rapid adiabatic cooling
+					  or heating */
 	      assert(p->PoverRho2() > 0.0);
 	      // Use P/rho as internal energy estimate since "u" may
 	      // be driven to 0 with cooling.