iobase.hh

Go to the documentation of this file.

/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
/*                                                                           */
/*  This file is part of the library KASKADE 7                               */
/*  https://www.zib.de/research/projects/kaskade7-finite-element-toolbox     */
/*                                                                           */
/*  Copyright (C) 2009-2018 Zuse Institute Berlin                            */
/*                                                                           */
/*  KASKADE 7 is distributed under the terms of the ZIB Academic License.    */
/*    see $KASKADE/academic.txt                                              */
/*                                                                           */
/* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * */
 
#ifndef IOBASE_HH
#define IOBASE_HH
 
#include <string>
 
#include <boost/fusion/include/for_each.hpp>
#include <boost/fusion/include/zip_view.hpp>
#include <boost/fusion/include/vector.hpp>
 
//---------------------------------------------------------------------
//---------------------------------------------------------------------
 
namespace Kaskade
{
  struct UnaryTrue {
    template <class T>
    bool operator()(T const&) const { return true; }
  };
 
 
  namespace IoDetail {
 
 
    template <class GridView,
              class Writer,
              class Filter,
              class RAIter>
    struct AddDataTo {};
 
 
  } // End of namespace IoDetail
  struct IoOptions
  {
    enum OutputType
    {
      ascii,
      binary
    };
    
    enum DataMode
    {
      conforming, 
      nonconforming, 
      inferred
    };
    
    enum Info
    {
      none, summary, detail
    };
    
    IoOptions() = default;
    
    IoOptions& setOrder(int ord) { order = ord; return *this; }
    
    IoOptions& setOutputType(OutputType out) { outputType = out; return *this; }
 
    IoOptions& setPrecision(int pre) { precision = pre; return *this; }
    
    IoOptions& setDataMode(DataMode mode) { dataMode = mode; return *this; }
    
    OutputType outputType = ascii;
    DataMode dataMode     = inferred;
    Info info             = none;
    int precision         = 6;
    int order             = 1;
  };
 
  extern IoOptions ioOptions_default;
  //---------------------------------------------------------------------
  
  std::string paddedString(int n, int places=3);
  
  //---------------------------------------------------------------------
 
  template <class GridView,
            class Writer,
            class VariableSet,
            class Filter,
            class Options>
  void writePartialFile(GridView const& gridView,
                        Writer& writer,
                        VariableSet const& vars,
                        std::string const& filename,
                        Filter const& filter,
                        Options const& options)
  {
    using namespace boost::fusion;
 
    // Bug in boost::fusion? The following line does not work!
    // for_each(zip(description.variables,vars.data),AddDataTo<VTKWriter>(writer));
 
    typedef vector<typename VariableSet::Descriptions::Variables const&,
                   typename VariableSet::Functions const&> ZipVector;
    zip_view<ZipVector> zipView(ZipVector(typename VariableSet::Descriptions::Variables(),vars.data));
    for_each(zipView,IoDetail::AddDataTo<GridView,Writer,Filter,std::string const*>(gridView,writer,filter,&vars.descriptions.names[0]));
 
    writer.write(filename.c_str(),options);
  }
 
} // namespace Kaskade
 
#endif