<ol>
-<li> New: DoFTools::make_periodicity_constraints implemented which
-inserts algebraic constraints due to periodic boundary conditions
-into a ConstraintMatrix.
-<br>
-(Matthias Maier, 2012/05/22)
-
<li>
New: step-48 demonstrates the solution of a nonlinear wave equation
with an explicit time stepping method. The usage of Gauss-Lobatto
(Katharina Kormann and Martin Kronbichler, 2012/05/05)
<li>
-New: step-44 demonstrates one approach to modelling large
+New: step-44 demonstrates one approach to modeling large
deformations of nearly-incompressible elastic solids. The
elastic response is governed by a non-linear, hyperelastic
free-energy function. The geometrical response is also
<h3>Specific improvements</h3>
<ol>
+<li> New: The function DataOutInterface::write_pvd_record can be used
+to provide Paraview with metadata that describes which time in a
+simulation a particular output file corresponds to.
+<br>
+(Marco Engelhard 2012/05/30)
+
<li> Fixed: Bug in 3d with hanging nodes in GridTools::find_cells_adjacent_to_vertex()
that caused find_active_cell_around_point() to fail in those cases.
<br>
(Timo Heister, Wolfgang Bangerth 2012/05/30)
+<li> New: DoFTools::make_periodicity_constraints implemented which
+inserts algebraic constraints due to periodic boundary conditions
+into a ConstraintMatrix.
+<br>
+(Matthias Maier, 2012/05/22)
+
<li> New: The GridIn::read_unv function can now read meshes generated
by the Salome framework, see http://www.salome-platform.org/ .
<br>
void write_pvtu_record (std::ostream &out,
const std::vector<std::string> &piece_names) const;
+ /**
+ * In ParaView it is possible to visualize time-
+ * dependent data tagged with the current
+ * integration time of a time dependent simulation. To use this
+ * feature you need a <code>.pvd</code>
+ * file that describes which VTU or PVTU file
+ * belongs to which timestep. This function writes a file that
+ * provides this mapping, i.e., it takes a list of pairs each of
+ * which indicates a particular time instant and the corresponding
+ * file that contains the graphical data for this time instant.
+ *
+ * A typical use case, in program that computes a time dependent
+ * solution, would be the following (<code>time</code> and
+ * <code>time_step</code> are member variables of the class with types
+ * <code>double</code> and <code>unsigned int</code>, respectively;
+ * the variable <code>times_and_names</code> is of type
+ * <code>std::vector@<std::pair@<double,std::string@> @></code>):
+ *
+ * @code
+ * template <int dim>
+ * void MyEquation<dim>::output_results () const
+ * {
+ * DataOut<dim> data_out;
+ *
+ * data_out.attach_dof_handler (dof_handler);
+ * data_out.add_data_vector (solution, "U");
+ * data_out.build_patches ();
+ *
+ * const std::string filename = "solution-" +
+ * Utilities::int_to_string (timestep_number, 3) +
+ * ".vtu";
+ * std::ofstream output (filename.c_str());
+ * data_out.write_vtu (output);
+ *
+ * times_and_names.push_back (std::pair<double,std::string> (time, filename));
+ * std::ofstream pvd_output ("solution.pvd");
+ * data_out.write_pvd_record (pvd_output, times_and_names);
+ * }
+ * @endcode
+ *
+ * @note See DataOutBase::write_vtu or
+ * DataOutBase::write_pvtu_record for
+ * writing solutions at each timestep.
+ *
+ * @note The second element of each pair, i.e., the file in which
+ * the graphical data for each time is stored, may itself be again
+ * a file that references other files. For example, it could be
+ * the name for a <code>.pvtu</code> file that references multiple
+ * parts of a parallel computation.
+ *
+ * @author Marco Engelhard, 2012
+ */
+ void write_pvd_record (std::ostream &out,
+ const std::vector<std::pair<double,std::string> > ×_and_names) const;
+
/**
* This function is the exact
* equivalent of the
}
+template <int dim, int spacedim>
+void
+DataOutInterface<dim,spacedim>::
+write_pvd_record (std::ostream &out,
+ const std::vector<std::pair<double,std::string> > ×_and_names) const
+{
+ AssertThrow (out, ExcIO());
+
+ out << "<?xml version=\"1.0\"?>\n";
+
+ std::time_t time1= std::time (0);
+ std::tm *time = std::localtime(&time1);
+ out << "<!--\n";
+ out << "#This file was generated by the deal.II library on "
+ << time->tm_year+1900 << "/"
+ << time->tm_mon+1 << "/"
+ << time->tm_mday << " at "
+ << time->tm_hour << ":"
+ << std::setw(2) << time->tm_min << ":"
+ << std::setw(2) << time->tm_sec
+ << "\n-->\n";
+
+ out << "<VTKFile type=\"Collection\" version=\"0.1\" ByteOrder=\"LittleEndian\">\n";
+ out << " <Collection>\n";
+
+ for(unsigned int i=0; i<times_and_names.size(); ++i)
+ out << " <DataSet timestep=\"" << times_and_names[i].first
+ << "\" group=\"\" part=\"0\" file=\"" << times_and_names[i].second
+ << "\"/>\n";
+
+ out << " </Collection>\n";
+ out << "</VTKFile>\n";
+
+ out.flush();
+
+ AssertThrow (out, ExcIO());
+}
+
template <int dim, int spacedim>
void
--- /dev/null
+//-----------------------------------------------------------------------------
+// $Id$
+// Version: $Name$
+//
+// Copyright (C) 2006, 2007, 2010, 2012 by the deal.II authors
+//
+// This file is subject to QPL and may not be distributed
+// without copyright and license information. Please refer
+// to the file deal.II/doc/license.html for the text and
+// further information on this license.
+//
+//-----------------------------------------------------------------------------
+
+// write the pvd master record for parallel visualization through the
+// vtu file format
+
+#include "../tests.h"
+#include <deal.II/base/data_out_base.h>
+#include <deal.II/base/logstream.h>
+
+#include <vector>
+#include <iomanip>
+#include <fstream>
+#include <string>
+#include <stdio.h>
+
+#include "patches.h"
+
+
+
+std::vector<DataOutBase::Patch<2,2> > patches;
+
+class DataOutX : public DataOutInterface<2,2>
+{
+ virtual
+ const std::vector< ::DataOutBase::Patch<2,2> > &
+ get_patches () const
+ {
+ return patches;
+ }
+
+ virtual
+ std::vector<std::string>
+ get_dataset_names () const
+ {
+ return std::vector<std::string>();
+ }
+};
+
+
+template <int dim, int spacedim>
+void check(std::ostream& out)
+{
+ std::vector<std::pair<double,std::string> > names(5);
+ names[0] = std::make_pair(0,"x1");
+ names[1] = std::make_pair(1,"x2");
+ names[2] = std::make_pair(1e1,"x3");
+ names[3] = std::make_pair(3.141,"d");
+ names[4] = std::make_pair(42e19,"i");
+
+ DataOutX x;
+ x.write_pvd_record (out, names);
+}
+
+
+
+int main()
+{
+ std::ofstream logfile("data_out_base_pvd/output");
+ check<2,2>(logfile);
+}
--- /dev/null
+<?xml version="1.0"?>
+<!--
+#This file was generated
+-->
+<VTKFile type="Collection" version="0.1" ByteOrder="LittleEndian">
+ <Collection>
+ <DataSet timestep="0" group="" part="0" file="x1"/>
+ <DataSet timestep="1" group="" part="0" file="x2"/>
+ <DataSet timestep="10" group="" part="0" file="x3"/>
+ <DataSet timestep="3.141" group="" part="0" file="d"/>
+ <DataSet timestep="4.2e+20" group="" part="0" file="i"/>
+ </Collection>
+</VTKFile>
projects / dealii-svn.git / commitdiff