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
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