rectilineargridwriter.hh

Go to the documentation of this file.

#pragma once
 
#include <array>
#include <iosfwd>
#include <map>
 
#include <dune/vtk/filewriter.hh>
#include <dune/vtk/types.hh>
#include <dune/vtk/vtkwriterbase.hh>
#include <dune/vtk/datacollectors/structureddatacollector.hh>
#include <dune/vtk/utility/concepts.hh>
 
namespace Dune::Vtk
{
 
  template <class GridView, class DataCollector = StructuredDataCollector<GridView>>
  class RectilinearGridWriter
      : public VtkWriterBase<GridView, DataCollector>
  {
    using Super = VtkWriterBase<GridView, DataCollector>;
    using pos_type = typename Super::pos_type;
 
  public:
    using Super::Super;
 
  private:
    virtual void writeSerialFile (std::ofstream& out) const override;
 
    virtual void writeParallelFile (std::ofstream& out, std::string const& pfilename, int size) const override;
 
    void writeCoordinates (std::ofstream& out, std::vector<pos_type>& offsets,
                           std::optional<std::size_t> timestep = {}) const;
 
    template <class T>
    std::array<std::uint64_t, 3> writeCoordinatesAppended (std::ofstream& out) const;
 
    virtual std::string fileExtension () const override
    {
      return "vtr";
    }
 
    virtual void writeGridAppended (std::ofstream& out, std::vector<std::uint64_t>& blocks) const override;
 
  private:
    using Super::dataCollector_;
    using Super::format_;
    using Super::datatype_;
    using Super::headertype_;
 
    // attached data
    using Super::pointData_;
    using Super::cellData_;
  };
 
  // deduction guides
  template <class GridView, class... Args,
    IsGridView<GridView> = true>
  RectilinearGridWriter(GridView, Args...)
    -> RectilinearGridWriter<GridView, StructuredDataCollector<GridView>>;
 
  template <class DataCollector, class... Args,
    IsDataCollector<DataCollector> = true>
  RectilinearGridWriter(DataCollector&, Args...)
    -> RectilinearGridWriter<typename DataCollector::GridView, DataCollector>;
 
  template <class DataCollector, class... Args,
    IsDataCollector<DataCollector> = true>
  RectilinearGridWriter(DataCollector&&, Args...)
    -> RectilinearGridWriter<typename DataCollector::GridView, DataCollector>;
 
  template <class DataCollector, class... Args,
    IsDataCollector<DataCollector> = true>
  RectilinearGridWriter(std::shared_ptr<DataCollector>, Args...)
    -> RectilinearGridWriter<typename DataCollector::GridView, DataCollector>;
 
} // end namespace Dune::Vtk
 
#include "rectilineargridwriter.impl.hh"