output_post_processor.hpp

Go to the documentation of this file.

#pragma once
 
#include <deal.II/dofs/dof_handler.h>
 
#include <deal.II/lac/la_parallel_vector.h>
 
#include <deal.II/numerics/data_component_interpretation.h>
#include <deal.II/numerics/data_postprocessor.h>
 
#include <meltpooldg/post_processing/generic_data_out.hpp>
#include <meltpooldg/post_processing/postprocessor.hpp>
 
#include <meltpooldg/utilities/cpp23_functions.h>
 
#include <string>
#include <vector>
 
namespace MeltPoolDG::SpeciesTransport
{
  template <int dim, int n_species, typename number, typename View>
    requires requires { typename View::state_type; }
  class MassFractionPostProcessor : public dealii::DataPostprocessor<dim>
  {
  public:
    MassFractionPostProcessor(const std::function<View(typename View::state_type &)> &view_creator,
                              const std::vector<std::string>                         &species_names)
      : view_creator(view_creator)
      , species_names(species_names)
    {
      Assert(
        species_names.size() == n_species,
        dealii::ExcMessage(
          "The number of provided species names must match the number of species in the simulation."));
      for (std::string &name : this->species_names)
        {
          name += "_mass_fraction";
        }
    }
 
    std::vector<std::string>
    get_names() const override
    {
      return species_names;
    }
 
    std::vector<dealii::DataComponentInterpretation::DataComponentInterpretation>
    get_data_component_interpretation() const override
    {
      return std::vector<dealii::DataComponentInterpretation::DataComponentInterpretation>(
        n_species, dealii::DataComponentInterpretation::component_is_scalar);
    }
 
    dealii::UpdateFlags
    get_needed_update_flags() const override
    {
      return dealii::update_values;
    }
 
    void
    evaluate_vector_field(const dealii::DataPostprocessorInputs::Vector<dim> &inputs,
                          std::vector<dealii::Vector<number>> &computed_quantities) const override
    {
      const unsigned int n_evaluation_points = inputs.solution_values.size();
 
      for (unsigned int p = 0; p < n_evaluation_points; ++p)
        {
          typename View::state_type w(dealii::make_array_view(inputs.solution_values[p]));
          const View                accessor = view_creator(w);
 
          for (int i = 0; i < n_species; ++i)
            {
              computed_quantities[p](i) = accessor.mass_fraction(i);
            }
        }
    }
 
    const std::function<View(typename View::state_type &)> view_creator;
 
    std::vector<std::string> species_names;
  };
 
  template <int dim, int n_species, typename number, typename View>
    requires requires { typename View::state_type; }
  class PartialDensityPostProcessor : public dealii::DataPostprocessor<dim>
  {
  public:
    PartialDensityPostProcessor(
      const std::function<View(typename View::state_type &)> &view_creator,
      const std::vector<std::string>                         &species_names)
      : view_creator(view_creator)
      , species_names(species_names)
    {
      Assert(
        species_names.size() == n_species,
        dealii::ExcMessage(
          "The number of provided species names must match the number of species in the simulation."));
      for (std::string &name : this->species_names)
        {
          name += "_partial_density";
        }
    }
 
    std::vector<std::string>
    get_names() const override
    {
      return species_names;
    }
 
    std::vector<dealii::DataComponentInterpretation::DataComponentInterpretation>
    get_data_component_interpretation() const override
    {
      return std::vector<dealii::DataComponentInterpretation::DataComponentInterpretation>(
        n_species, dealii::DataComponentInterpretation::component_is_scalar);
    }
 
    dealii::UpdateFlags
    get_needed_update_flags() const override
    {
      return dealii::update_values;
    }
 
    void
    evaluate_vector_field(const dealii::DataPostprocessorInputs::Vector<dim> &inputs,
                          std::vector<dealii::Vector<number>> &computed_quantities) const override
    {
      const unsigned int n_evaluation_points = inputs.solution_values.size();
 
      for (unsigned int p = 0; p < n_evaluation_points; ++p)
        {
          typename View::state_type w(dealii::make_array_view(inputs.solution_values[p]));
          const View                accessor = view_creator(w);
 
          for (int i = 0; i < n_species; ++i)
            {
              computed_quantities[p](i) = accessor.partial_density(i);
            }
        }
    }
 
    const std::function<View(typename View::state_type &)> view_creator;
 
    std::vector<std::string> species_names;
  };
} // namespace MeltPoolDG::SpeciesTransport