lateral_averaging.h

Go to the documentation of this file.

/*
  Copyright (C) 2011 - 2023 by the authors of the ASPECT code.

  This file is part of ASPECT.

  ASPECT 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, or (at your option)
  any later version.

  ASPECT 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 ASPECT; see the file LICENSE.  If not see
  <http://www.gnu.org/licenses/>.
*/


#ifndef _aspect_lateral_averaging_h
#define _aspect_lateral_averaging_h

#include <aspect/simulator_access.h>

#include <deal.II/fe/fe_values.h>

namespace aspect
{
  using namespace dealii;

  namespace internal
  {
    template <int dim>
    class FunctorBase
    {
      public:
        virtual
        bool
        need_material_properties() const;

        virtual
        void
        create_additional_material_model_outputs (const unsigned int n_points,
                                                  MaterialModel::MaterialModelOutputs<dim> &outputs) const;

        virtual
        void
        setup(const unsigned int q_points);

        virtual
        void
        operator()(const MaterialModel::MaterialModelInputs<dim> &in,
                   const MaterialModel::MaterialModelOutputs<dim> &out,
                   const FEValues<dim> &fe_values,
                   const LinearAlgebra::BlockVector &solution,
                   std::vector<double> &output) = 0;

        virtual
        ~FunctorBase();
    };
  }

  template <int dim>
  class LateralAveraging : public SimulatorAccess<dim>
  {
    public:
      std::vector<std::vector<double>>
      compute_lateral_averages(const unsigned int n_slices,
                               const std::vector<std::string> &property_names) const;

      std::vector<std::vector<double>>
      compute_lateral_averages(const std::vector<double> &depth_bounds,
                               const std::vector<std::string> &property_names) const;

      std::vector<std::vector<double>>
      compute_lateral_averages(const std::vector<double> &depth_bounds,
                               std::vector<std::unique_ptr<internal::FunctorBase<dim>>> &functors) const;

      void
      get_temperature_averages(std::vector<double> &values) const;

      void
      get_composition_averages(const unsigned int composition_index,
                               std::vector<double> &values) const;

      void
      get_viscosity_averages(std::vector<double> &values) const;

      void
      get_log_viscosity_averages(std::vector<double> &values) const;

      void
      get_velocity_magnitude_averages(std::vector<double> &values) const;

      void
      get_sinking_velocity_averages(std::vector<double> &values) const;

      void
      get_rising_velocity_averages(std::vector<double> &values) const;

      void
      get_Vs_averages(std::vector<double> &values) const;

      void
      get_Vp_averages(std::vector<double> &values) const;

      void
      get_vertical_heat_flux_averages(std::vector<double> &values) const;

      void
      get_vertical_mass_flux_averages(std::vector<double> &values) const;
  };
}


#endif