fe_output.hh

Go to the documentation of this file.

#pragma once
 
// standard library
#include <thread>
 
// external libraries
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/vector_memory.h>
#include <deal.II/numerics/data_out.h>
 
// DiFfRG
#include <DiFfRG/common/utils.hh>
 
namespace DiFfRG
{
  using namespace dealii;
 
  template <uint dim, typename VectorType> class FEOutput
  {
  public:
    FEOutput(std::string top_folder, std::string output_name, std::string output_folder, const JSONValue &json);
 
    FEOutput();
    ~FEOutput();
 
    void attach(const DoFHandler<dim> &dof_handler, const VectorType &solution, const std::string &name);
 
    void attach(const DoFHandler<dim> &dof_handler, const VectorType &solution, const std::vector<std::string> &names);
 
    void flush(double time);
 
  private:
    const std::string top_folder;
    const std::string output_name;
    const std::string output_folder;
    constexpr static uint safe_dim = dim == 0 ? 1 : dim;
    DataOut<safe_dim> data_out;
 
    uint series_number, subdivisions;
    std::vector<std::pair<double, std::string>> time_series;
 
    std::thread output_thread;
    std::vector<typename VectorMemory<VectorType>::Pointer> attached_solutions;
 
    GrowingVectorMemory<VectorType> mem;
  };
 
  template <> class FEOutput<0, Vector<double>>
  {
  public:
    void attach(const DoFHandler<0> &dof_handler, const Vector<double> &solution, const std::string &name);
    void attach(const DoFHandler<0> &dof_handler, const Vector<double> &solution,
                const std::vector<std::string> &names);
 
    void flush(double time);
  };
 
} // namespace DiFfRG