]> https://gitweb.dealii.org/ - dealii-svn.git/commitdiff
Fix an actual bug: we were taking a completely meaningless temperature
authorbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 4 Oct 2011 18:58:33 +0000 (18:58 +0000)
committerbangerth <bangerth@0785d39b-7218-0410-832d-ea1e28bc413d>
Tue, 4 Oct 2011 18:58:33 +0000 (18:58 +0000)
field to evaluate the density in computing the rhs of the temperature
equation. Use the correctly extrapolated field.

git-svn-id: https://svn.dealii.org/trunk@24532 0785d39b-7218-0410-832d-ea1e28bc413d

deal.II/examples/step-32/step-32.cc

index f1ad6a271858cdb5d2331979a4fdde455b59bf2d..b689c959d6945a19f9bc8edafefa8cfc902180d7 100644 (file)
@@ -3324,10 +3324,10 @@ namespace Step32
             scratch.old_strain_rates[q]);
 
        const double gamma
-         = ((EquationData::radiogenic_heating * EquationData::density(T_term_for_rhs)
+         = ((EquationData::radiogenic_heating * EquationData::density(ext_T)
              +
              2 * EquationData::eta * extrapolated_strain_rate * extrapolated_strain_rate) /
-            (EquationData::density(T_term_for_rhs) * EquationData::specific_heat));
+            (EquationData::density(ext_T) * EquationData::specific_heat));
 
        for (unsigned int i=0; i<dofs_per_cell; ++i)
          {

In the beginning the Universe was created. This has made a lot of people very angry and has been widely regarded as a bad move.

Douglas Adams


Typeset in Trocchi and Trocchi Bold Sans Serif.