From 41b0edbe29b615284de18d003b1edbadcc6f316f Mon Sep 17 00:00:00 2001
From: Guido Kanschat <dr.guido.kanschat@gmail.com>
Date: Thu, 22 May 2014 20:20:24 +0000
Subject: [PATCH] adjust OutputOperator for AnyData

git-svn-id: https://svn.dealii.org/trunk@32967 0785d39b-7218-0410-832d-ea1e28bc413d
---
 deal.II/include/deal.II/algorithms/operator.h |  2 +-
 .../deal.II/algorithms/operator.templates.h   | 50 +++++++------------
 2 files changed, 18 insertions(+), 34 deletions(-)

diff --git a/deal.II/include/deal.II/algorithms/operator.h b/deal.II/include/deal.II/algorithms/operator.h
index 5dbe64a4b7..bad83a96dd 100644
--- a/deal.II/include/deal.II/algorithms/operator.h
+++ b/deal.II/include/deal.II/algorithms/operator.h
@@ -208,7 +208,7 @@ namespace Algorithms
     /**
      * @deprecated Output all the vectors in NamedData.
      */
-    virtual OutputOperator<VECTOR> &operator<< (const NamedData<VECTOR *> &vectors);
+    OutputOperator<VECTOR> &operator<< (const NamedData<VECTOR *> &vectors);
   protected:
     unsigned int step;
   private:
diff --git a/deal.II/include/deal.II/algorithms/operator.templates.h b/deal.II/include/deal.II/algorithms/operator.templates.h
index eed6e1a984..0a941933ed 100644
--- a/deal.II/include/deal.II/algorithms/operator.templates.h
+++ b/deal.II/include/deal.II/algorithms/operator.templates.h
@@ -128,29 +128,27 @@ namespace Algorithms
   {
     if (os == 0)
       {
-        //TODO: make this possible
-        //deallog << ' ' << step;
-        //for (unsigned int i=0;i<vectors.size();++i)
-        //  vectors(i)->print(deallog);
-        //deallog << std::endl;
+        deallog << "Step " << step << std::endl;
+        for (unsigned int i=0; i<vectors.size(); ++i)
+	  {
+	    const VECTOR* v = vectors.try_read_ptr<VECTOR>(i);
+	    if (v == 0) continue;
+	    for (unsigned int j=0; j<v->size(); ++j)
+	      (*os) << ' ' << (*v)(j);
+	    deallog << std::endl;
+	  }
+        deallog << std::endl;
       }
     else
       {
         (*os) << ' ' << step;
         for (unsigned int i=0; i<vectors.size(); ++i)
 	  {
-	    if (vectors.is_type<VECTOR*>(i))
-	      {
-		const VECTOR& v = *vectors.entry<VECTOR*>(i);
-		for (unsigned int j=0; j<v.size(); ++j)
-		  (*os) << ' ' << v(j);
-	      }
-	    else if (vectors.is_type<const VECTOR*>(i))
-	      {
-		const VECTOR& v = *vectors.entry<const VECTOR*>(i);
-		for (unsigned int j=0; j<v.size(); ++j)
-		  (*os) << ' ' << v(j);
-	      }
+	    const VECTOR* v = vectors.try_read_ptr<VECTOR>(i);
+	    if (v == 0) continue;
+	    for (unsigned int j=0; j<v->size(); ++j)
+	      (*os) << ' ' << (*v)(j);
+	    (*os) << std::endl;
 	  }
         (*os) << std::endl;
       }
@@ -162,22 +160,8 @@ namespace Algorithms
   OutputOperator<VECTOR> &
   OutputOperator<VECTOR>::operator<< (const NamedData<VECTOR *> &vectors)
   {
-    if (os == 0)
-      {
-        //TODO: make this possible
-        //deallog << ' ' << step;
-        //for (unsigned int i=0;i<vectors.size();++i)
-        //  vectors(i)->print(deallog);
-        //deallog << std::endl;
-      }
-    else
-      {
-        (*os) << ' ' << step;
-        for (unsigned int i=0; i<vectors.size(); ++i)
-          for (unsigned int j=0; j<vectors(i)->size(); ++j)
-            (*os) << ' ' << (*vectors(i))(j);
-        (*os) << std::endl;
-      }
+    const AnyData newdata = vectors;
+    (*this) << newdata;
     return *this;
   }
 }
-- 
2.39.5