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

[dealii.demon at gmail.com]

Revision 1975

Augment a comment.

U   trunk/aspect/source/postprocess/heat_flux_statistics.cc


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

Diff:
Modified: trunk/aspect/source/postprocess/heat_flux_statistics.cc
===================================================================
--- trunk/aspect/source/postprocess/heat_flux_statistics.cc	2013-10-18 15:15:59 UTC (rev 1974)
+++ trunk/aspect/source/postprocess/heat_flux_statistics.cc	2013-10-18 15:28:29 UTC (rev 1975)
@@ -87,7 +87,13 @@
                       composition_values[c]);
 
                 in.position = fe_face_values.get_quadrature_points();
-                in.strain_rate.resize(0);// we are not reading the viscosity
+
+		// since we are not reading the viscosity and the viscosity
+		// is the only coefficient that depends on the strain rate,
+		// we need not compute the strain rate. set the corresponding
+		// array to empty, to prevent accidental use
+                in.strain_rate.resize(0);
+		
                 for (unsigned int i=0; i<fe_face_values.n_quadrature_points; ++i)
                   {
                     for (unsigned int c=0; c<this->n_compositional_fields(); ++c)