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

[dealii.demon at gmail.com]

Revision 2368

Rewrite some code to use introspection instead of hard coded numbers.

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


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

Diff:
Modified: trunk/aspect/source/postprocess/composition_statistics.cc
===================================================================
--- trunk/aspect/source/postprocess/composition_statistics.cc	2014-03-26 20:26:20 UTC (rev 2367)
+++ trunk/aspect/source/postprocess/composition_statistics.cc	2014-03-26 20:33:43 UTC (rev 2368)
@@ -94,11 +94,11 @@
 
       for (unsigned int c=0; c<this->n_compositional_fields(); ++c)
         {
-          IndexSet range = this->get_solution().block(3+c).locally_owned_elements();
+          IndexSet range = this->get_solution().block(this->introspection().block_indices.compositional_fields[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(3+c)(idx);
+              const double val =  this->get_solution().block(this->introspection().block_indices.compositional_fields[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);