Add CPO serialization.

Author: MFraters

include/aspect/postprocess/crystal_preferred_orientation.h

@@ -77,6 +77,23 @@ namespace aspect
         std::list<std::string>
         required_other_postprocessors () const override;
 
+        /**
+               * Save the state of this object.
+               */
+        void save (std::map<std::string, std::string> &status_strings) const override;
+
+        /**
+         * Restore the state of the object.
+         */
+        void load (const std::map<std::string, std::string> &status_strings) override;
+
+        /**
+         * Serialize the contents of this class as far as they are not read
+         * from input parameter files.
+         */
+        template <class Archive>
+        void serialize (Archive &ar, const unsigned int version);
+
         /**
          * Declare the parameters this class takes through input files.
          */

source/postprocess/crystal_preferred_orientation.cc

@@ -602,7 +602,72 @@ namespace aspect
         return Output::EulerAngles;
       return Output::not_found;
     }
+    template <int dim>
+    template <class Archive>
+    void CrystalPreferredOrientation<dim>::serialize (Archive &ar, const unsigned int)
+    {
+      ar &last_output_time
+      & output_file_number
+      & times_and_pvtu_file_names
+      & output_file_names_by_timestep
+      & xdmf_entries
+      ;
+    }
+
+
+    template <int dim>
+    void
+    CrystalPreferredOrientation<dim>::save (std::map<std::string, std::string> &status_strings) const
+    {
+      for (unsigned int particle_manager = 0; particle_manager < this->n_particle_managers(); ++particle_manager)
+        {
+          std::string particles_output_base_name = "Particles";
+          if (particle_manager > 0)
+            {
+              particles_output_base_name += "-" + Utilities::int_to_string(particle_manager+1);
+            }
+
+          // Serialize into a stringstream. Put the following into a code
+          // block of its own to ensure the destruction of the 'oa'
+          // archive triggers a flush() on the stringstream so we can
+          // query the completed string below.
+          std::ostringstream os;
+          {
+            aspect::oarchive oa (os);
+
+            this->get_particle_manager(particle_manager).save(os);
+            oa << (*this);
+          }
 
+          status_strings[particles_output_base_name] = os.str();
+        }
+    }
+
+
+    template <int dim>
+    void
+    CrystalPreferredOrientation<dim>::load (const std::map<std::string, std::string> &status_strings)
+    {
+      for (unsigned int particle_manager = 0; particle_manager < this->n_particle_managers(); ++particle_manager)
+        {
+          std::string particles_output_base_name = "Particles";
+          if (particle_manager > 0)
+            {
+              particles_output_base_name += "-" + Utilities::int_to_string(particle_manager+1);
+            }
+          // see if something was saved
+          if (status_strings.find(particles_output_base_name) != status_strings.end())
+            {
+              std::istringstream is (status_strings.find(particles_output_base_name)->second);
+              aspect::iarchive ia (is);
+
+              // Load the particle manager
+              this->get_particle_manager(particle_manager).load(is);
+
+              ia >> (*this);
+            }
+        }
+    }
 
 
     template <int dim>