[cig-commits] commit: Add generate_pvd

[Mercurial]

changeset:   177:90f25c050c7f
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Dec 23 14:04:19 2009 -0800
files:       tools/generate_pvd.cxx
description:
Add generate_pvd


diff -r 0fe143638a45 -r 90f25c050c7f tools/generate_pvd.cxx
--- /dev/null	Thu Jan 01 00:00:00 1970 +0000
+++ b/tools/generate_pvd.cxx	Wed Dec 23 14:04:19 2009 -0800
@@ -0,0 +1,88 @@
+/* A simple program to generate field and particle pvd files 
+
+   You will need Boost (www.boost.org) to compile this.  Almost any
+   version should do.  On Debian, it is
+
+     g++ generate_pvd.cxx -o generate_pvd -lboost_filesystem
+
+   To generate a static binary, compile it with
+
+     g++ generate_pvd.cxx -o generate_pvd -lboost_system -lboost_filesystem -static
+
+   Usage: generate_pvd START END STEP
+
+   This will output two files, fields.pvd and particles.pvd.  Put them
+   in the same directory as the other fields and particle files.  In
+   paraview, open the pvd files.  You will be able use the movie
+   controls (play, pause, step forward/backward, jump to
+   beginning/end) to examine the time series.
+
+**  Copyright (C) 2009, California Institute of Technology
+
+**  This script is free software; you can redistribute it and/or
+**  modify it under the terms of the GNU General Public License as
+**  published by the Free Software Foundation; either version 2 of the
+**  License, or (at your option) any later version.
+**
+**  This library is distributed in the hope that it will be useful,
+**  but WITHOUT ANY WARRANTY; without even the implied warranty of
+**  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
+**  General Public License for more details.
+**
+**  You should have received a copy of the GNU General Public License
+**  along with this library; if not, write to the Free Software
+**  Foundation, Inc., 51 Franklin St, Fifth Floor, Boston, MA
+**  02110-1301 USA
+*/
+
+
+#include <boost/filesystem.hpp>
+#include <boost/filesystem/fstream.hpp>
+#include <cstdlib>
+#include <iostream>
+
+namespace fs=boost::filesystem;
+using namespace std;
+
+
+int main(int argc, char *argv[])
+{
+  if(argc<3)
+    {
+      cerr << "Not enough arguments\n"
+           << "Usage: generate_pvd [start] [end] [step]\n";
+      exit(1);
+    }
+  int start=atoi(argv[1]);
+  int end=atoi(argv[2]);
+  int step=1;
+  if(argc>3)
+    step=atoi(argv[3]);
+  
+  fs::ofstream particles("particles.pvd"), fields("fields.pvd");
+
+  particles << "<?xml version=\"1.0\"?>\n"
+            << "<VTKFile type=\"Collection\" version=\"0.1\">\n"
+            << "  <Collection>\n";
+  fields << "<?xml version=\"1.0\"?>\n"
+         << "<VTKFile type=\"Collection\" version=\"0.1\">\n"
+         << "  <Collection>\n";
+  for(int i=start; i<=end; i+=step)
+    {
+      particles << "    <DataSet timestep=\"" << i
+                << "\" file=\"particles.";
+      particles.width(5);
+      particles.fill('0');
+      particles << i << ".pvtu\"/>\n";
+
+      fields << "    <DataSet timestep=\"" << i
+              << "\" file=\"fields.";
+      fields.width(5);
+      fields.fill('0');
+      fields << i << ".pvts\"/>\n";
+    }
+  particles << "  </Collection>\n"
+            << "</VTKFile>\n";
+  fields << "  </Collection>\n"
+         << "</VTKFile>\n";
+}