[cig-commits] commit 2358 by heister to /var/svn/dealii/aspect

[dealii.demon at gmail.com]

Revision 2358

fix bug in compositional statistics

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


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

Diff:
Modified: trunk/aspect/source/postprocess/composition_statistics.cc
===================================================================
--- trunk/aspect/source/postprocess/composition_statistics.cc	2014-03-26 14:58:29 UTC (rev 2357)
+++ trunk/aspect/source/postprocess/composition_statistics.cc	2014-03-26 15:19:20 UTC (rev 2358)
@@ -94,13 +94,11 @@
 
       for (unsigned int c=0; c<this->n_compositional_fields(); ++c)
         {
-          double local_min_temperature = std::numeric_limits<double>::max();
-          double local_max_temperature = -std::numeric_limits<double>::max();
           IndexSet range = this->get_solution().block(3+c).locally_owned_elements();
           for (unsigned int i=0;i<range.n_elements();++i)
             {
               const unsigned int idx = range.nth_index_in_set(i);
-              const double val =  this->get_solution().block(2)(idx);
+              const double val =  this->get_solution().block(3+c)(idx);
 
               local_min_compositions[c] = std::min<double> (local_min_compositions[c], val);
               local_max_compositions[c] = std::max<double> (local_max_compositions[c], val);