[cig-commits] commit 2404 by bangerth to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2404

Add a missing factor of 2.

U   trunk/aspect/source/postprocess/dynamic_topography.cc
U   trunk/aspect/source/postprocess/visualization/dynamic_topography.cc


http://www.dealii.org/websvn/revision.php?repname=Aspect+Repository&path=%2F&rev=2404&peg=2404

Diff:
Modified: trunk/aspect/source/postprocess/dynamic_topography.cc
===================================================================
--- trunk/aspect/source/postprocess/dynamic_topography.cc	2014-04-03 12:57:22 UTC (rev 2403)
+++ trunk/aspect/source/postprocess/dynamic_topography.cc	2014-04-04 03:28:01 UTC (rev 2404)
@@ -99,7 +99,7 @@
                   const double viscosity = out.viscosities[q];
                   const double density   = out.densities[q];
 
-                  const SymmetricTensor<2,dim> stress = viscosity * in.strain_rate[q];
+                  const SymmetricTensor<2,dim> stress = 2 * viscosity * in.strain_rate[q];
 
                   const Tensor<1,dim> gravity = this->get_gravity_model().gravity_vector(location);
                   const Tensor<1,dim> gravity_direction = gravity/gravity.norm();
@@ -185,7 +185,7 @@
                                   "that sit along the top surface, we evaluate the stress and "
                                   "evaluate the component of it in the direction in which "
                                   "gravity acts. In other words, we compute "
-                                  "$\sigma_{rr}={\hat g}^T(\eta \varepsilon(\mathbf u))\hat g$ "
+                                  "$\sigma_{rr}={\hat g}^T(2 \eta \varepsilon(\mathbf u))\hat g$ "
                                   "where $\hat g = \mathbf g/\|\mathbf g\|$ is the direction of "
                                   "the gravity vector $\mathbf g$. From this, the dynamic "
                                   "topography is computed using the formula "

Modified: trunk/aspect/source/postprocess/visualization/dynamic_topography.cc
===================================================================
--- trunk/aspect/source/postprocess/visualization/dynamic_topography.cc	2014-04-03 12:57:22 UTC (rev 2403)
+++ trunk/aspect/source/postprocess/visualization/dynamic_topography.cc	2014-04-04 03:28:01 UTC (rev 2404)
@@ -91,7 +91,7 @@
             const double viscosity = out.viscosities[q];
             const double density   = out.densities[q];
 
-            const SymmetricTensor<2,dim> stress = viscosity * in.strain_rate[q];
+            const SymmetricTensor<2,dim> stress = 2 * viscosity * in.strain_rate[q];
 
             const Tensor<1,dim> gravity = this->get_gravity_model().gravity_vector(location);
             const Tensor<1,dim> gravity_direction = gravity/gravity.norm();
@@ -121,7 +121,7 @@
                                                   "dynamic topography requires us to compute the stress tensor and "
                                                   "evaluate the component of it in the direction in which "
                                                   "gravity acts. In other words, we compute "
-                                                  "$\sigma_{rr}={\hat g}^T(\eta \varepsilon(\mathbf u))\hat g$ "
+                                                  "$\sigma_{rr}={\hat g}^T(2 * \eta \varepsilon(\mathbf u))\hat g$ "
                                                   "where $\hat g = \mathbf g/\|\mathbf g\|$ is the direction of "
                                                   "the gravity vector $\mathbf g$. From this, the dynamic "
                                                   "topography is computed using the formula "