smoothness_estimator.h

Go to the documentation of this file.

// ------------------------------------------------------------------------
//
// SPDX-License-Identifier: LGPL-2.1-or-later
// Copyright (C) 2020 - 2023 by the deal.II authors
//
// This file is part of the deal.II library.
//
// Part of the source code is dual licensed under Apache-2.0 WITH
// LLVM-exception OR LGPL-2.1-or-later. Detailed license information
// governing the source code and code contributions can be found in
// LICENSE.md and CONTRIBUTING.md at the top level directory of deal.II.
//
// ------------------------------------------------------------------------
 
#ifndef dealii_smoothness_estimator_h
#define dealii_smoothness_estimator_h
 
 
#include <deal.II/base/config.h>
 
#include <deal.II/fe/component_mask.h>
 
#include <deal.II/numerics/vector_tools_common.h>
 
#include <functional>
#include <vector>
 
 
DEAL_II_NAMESPACE_OPEN
 
 
// forward declarations
#ifndef DOXYGEN
template <typename Number>
class Vector;
 
template <int dim, int spacedim>
DEAL_II_CXX20_REQUIRES((concepts::is_valid_dim_spacedim<dim, spacedim>))
class DoFHandler;
 
namespace FESeries
{
  template <int dim, int spacedim>
  class Fourier;
  template <int dim, int spacedim>
  class Legendre;
} // namespace FESeries
 
namespace hp
{
  template <int dim, int spacedim>
  class FECollection;
} // namespace hp
#endif
 
 
namespace SmoothnessEstimator
{
  namespace Legendre
  {
    template <int dim, int spacedim, typename VectorType>
    void
    coefficient_decay(FESeries::Legendre<dim, spacedim> &fe_legendre,
                      const DoFHandler<dim, spacedim>   &dof_handler,
                      const VectorType                  &solution,
                      Vector<float>                     &smoothness_indicators,
                      const VectorTools::NormType        regression_strategy =
                        VectorTools::Linfty_norm,
                      const double smallest_abs_coefficient = 1e-10,
                      const bool   only_flagged_cells       = false);
 
    template <int dim, int spacedim, typename VectorType>
    void
    coefficient_decay_per_direction(
      FESeries::Legendre<dim, spacedim> &fe_legendre,
      const DoFHandler<dim, spacedim>   &dof_handler,
      const VectorType                  &solution,
      Vector<float>                     &smoothness_indicators,
      const ComponentMask               &coefficients_predicate   = {},
      const double                       smallest_abs_coefficient = 1e-10,
      const bool                         only_flagged_cells       = false);
 
    template <int dim, int spacedim>
    FESeries::Legendre<dim, spacedim>
    default_fe_series(
      const hp::FECollection<dim, spacedim> &fe_collection,
      const unsigned int component = numbers::invalid_unsigned_int);
  } // namespace Legendre
 
 
 
  namespace Fourier
  {
    template <int dim, int spacedim, typename VectorType>
    void
    coefficient_decay(FESeries::Fourier<dim, spacedim> &fe_fourier,
                      const DoFHandler<dim, spacedim>  &dof_handler,
                      const VectorType                 &solution,
                      Vector<float>                    &smoothness_indicators,
                      const VectorTools::NormType       regression_strategy =
                        VectorTools::Linfty_norm,
                      const double smallest_abs_coefficient = 1e-10,
                      const bool   only_flagged_cells       = false);
 
    template <int dim, int spacedim, typename VectorType>
    void
    coefficient_decay_per_direction(
      FESeries::Fourier<dim, spacedim> &fe_fourier,
      const DoFHandler<dim, spacedim>  &dof_handler,
      const VectorType                 &solution,
      Vector<float>                    &smoothness_indicators,
      const ComponentMask              &coefficients_predicate   = {},
      const double                      smallest_abs_coefficient = 1e-10,
      const bool                        only_flagged_cells       = false);
 
    template <int dim, int spacedim>
    FESeries::Fourier<dim, spacedim>
    default_fe_series(
      const hp::FECollection<dim, spacedim> &fe_collection,
      const unsigned int component = numbers::invalid_unsigned_int);
  } // namespace Fourier
} // namespace SmoothnessEstimator
 
 
DEAL_II_NAMESPACE_CLOSE
 
#endif