grid_refinement.h

// ---------------------------------------------------------------------
//
// Copyright (C) 2000 - 2016 by the deal.II authors
//
// This file is part of the deal.II library.
//
// The deal.II library is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE at
// the top level of the deal.II distribution.
//
// ---------------------------------------------------------------------

#ifndef dealii__grid_refinement_h
#define dealii__grid_refinement_h


#include <deal.II/base/config.h>
#include <deal.II/base/exceptions.h>
#include <deal.II/grid/tria.h>

#include <vector>
#include <limits>

DEAL_II_NAMESPACE_OPEN

// forward declarations
template <int dim, int spacedim> class Triangulation;
template <class T> class Vector;


namespace GridRefinement
{
  template <int dim>
  std::pair<double, double>
  adjust_refine_and_coarsen_number_fraction (const unsigned int  current_n_cells,
                                             const unsigned int  max_n_cells,
                                             const double        top_fraction_of_cells,
                                             const double        bottom_fraction_of_cells);

  template <int dim, class VectorType, int spacedim>
  void
  refine_and_coarsen_fixed_number
  (Triangulation<dim,spacedim> &triangulation,
   const VectorType            &criteria,
   const double                top_fraction_of_cells,
   const double                bottom_fraction_of_cells,
   const unsigned int          max_n_cells = std::numeric_limits<unsigned int>::max());

  template <int dim, class VectorType, int spacedim>
  void
  refine_and_coarsen_fixed_fraction
  (Triangulation<dim,spacedim> &tria,
   const VectorType            &criteria,
   const double                top_fraction,
   const double                bottom_fraction,
   const unsigned int          max_n_cells = std::numeric_limits<unsigned int>::max());



  template <int dim, class VectorType, int spacedim>
  void
  refine_and_coarsen_optimize (Triangulation<dim,spacedim> &tria,
                               const VectorType            &criteria,
                               const unsigned int          order=2);

  template <int dim, class VectorType, int spacedim>
  void refine (Triangulation<dim,spacedim> &tria,
               const VectorType            &criteria,
               const double                threshold,
               const unsigned int          max_to_mark = numbers::invalid_unsigned_int);

  template <int dim, class VectorType, int spacedim>
  void coarsen (Triangulation<dim,spacedim> &tria,
                const VectorType            &criteria,
                const double                threshold);

  DeclException0(ExcNegativeCriteria);

  DeclException0 (ExcInvalidParameterValue);
}



DEAL_II_NAMESPACE_CLOSE

#endif //dealii__grid_refinement_h