From d3388a5387d9fb172085238082940bde5fabf270 Mon Sep 17 00:00:00 2001
From: bangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Date: Tue, 4 Oct 2011 18:49:45 +0000
Subject: [PATCH] Clean up code somewhat.

git-svn-id: https://svn.dealii.org/trunk@24529 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/examples/step-32/step-32.cc | 32 +++++++++++++++++++++--------
 1 file changed, 24 insertions(+), 8 deletions(-)

diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc
index 6d5fc9cbca..14cc93c82f 100644
--- a/deal.II/examples/step-32/step-32.cc
+++ b/deal.II/examples/step-32/step-32.cc
@@ -3275,7 +3275,7 @@ namespace Step32
 	const double old_Ts
 	  = (use_bdf2_scheme ?
 	     (scratch.old_temperature_values[q] *
-	      (time_step + old_time_step) / old_time_step
+	      (1 + time_step/old_time_step)
 	      -
 	      scratch.old_old_temperature_values[q] *
 	      (time_step * time_step) /
@@ -3283,30 +3283,46 @@ namespace Step32
 	     :
 	     scratch.old_temperature_values[q]);
 
+	const double ext_T
+	  = (use_bdf2_scheme ?
+	     (scratch.old_temperature_values[q] *
+	      (1 + time_step/old_time_step)
+	      -
+	      scratch.old_old_temperature_values[q] *
+	      time_step/old_time_step)
+	     :
+	     scratch.old_temperature_values[q]);
+
 	const Tensor<1,dim> ext_grad_T
 	  = (use_bdf2_scheme ?
 	     (scratch.old_temperature_grads[q] *
-	      (1+time_step/old_time_step)
+	      (1 + time_step/old_time_step)
 	      -
 	      scratch.old_old_temperature_grads[q] *
-	      time_step / old_time_step)
+	      time_step/old_time_step)
 	     :
 	     scratch.old_temperature_grads[q]);
 
 	const Tensor<1,dim> extrapolated_u
 	  = (use_bdf2_scheme ?
-	     (scratch.old_velocity_values[q] * (1+time_step/old_time_step) -
-	      scratch.old_old_velocity_values[q] * time_step/old_time_step)
+	     (scratch.old_velocity_values[q] *
+	      (1 + time_step/old_time_step)
+	      -
+	      scratch.old_old_velocity_values[q] *
+	      time_step/old_time_step)
 	     :
 	     scratch.old_velocity_values[q]);
+
 	const SymmetricTensor<2,dim> extrapolated_strain_rate
 	  = (use_bdf2_scheme ?
-	     (scratch.old_strain_rates[q] * (1+time_step/old_time_step) -
-	      scratch.old_old_strain_rates[q] * time_step/old_time_step)
+	     (scratch.old_strain_rates[q] *
+	      (1 + time_step/old_time_step)
+	      -
+	      scratch.old_old_strain_rates[q] *
+	      time_step/old_time_step)
 	     :
 	     scratch.old_strain_rates[q]);
 
-// @todo ?????? why old_Ts?
 	const double gamma
 	  = ((EquationData::radiogenic_heating * EquationData::density(old_Ts)
 	      +
-- 
2.39.5