demo for polarization ordering

extensions/base_problems/BaseProblem.hh

@@ -95,10 +95,11 @@ Flag BaseProblem<ProblemType>::initDataFromFile(AdaptInfo *adaptInfo)
     if (!file_exists(filename))
       throw(std::runtime_error("The file '" + filename + "' does not exist!"));
     
-    std::vector<DOFVector<double>*> solutions;
-    for (int i = 0; i < prob->getNumComponents(); i++)
-      solutions.push_back(prob->getSolution()->getDOFVector(i));
-    ArhReader::read(filename, prob->getMesh(), solutions);
+//     std::vector<DOFVector<double>*> solutions;
+//     for (int i = 0; i < prob->getNumComponents(); i++)
+//       solutions.push_back(prob->getSolution()->getDOFVector(i));
+//     ArhReader::read(filename, prob->getMesh(), solutions);
+    readFile(filename, prob->getSolution());
   }
   else if (readFormat == "dat") 
   {