parsed_zero_average_constraints.h

Go to the documentation of this file.

//-----------------------------------------------------------
//
//    Copyright (C) 2016 by the deal2lkit authors
//
//    This file is part of the deal2lkit library.
//
//    The deal2lkit 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 deal2lkit distribution.
//
//-----------------------------------------------------------

#ifndef _d2k_parsed_zero_average_constraints_h
#define _d2k_parsed_zero_average_constraints_h

#include <deal2lkit/config.h>
#include <deal.II/base/exceptions.h>
#include <deal.II/base/parsed_function.h>
#include <deal2lkit/parameter_acceptor.h>
#include <deal.II/fe/component_mask.h>
#include <algorithm>
#include <map>
#include <deal2lkit/utilities.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>


using namespace dealii;


D2K_NAMESPACE_OPEN

template <int dim, int spacedim=dim>
class ParsedZeroAverageConstraints : public ParameterAcceptor
{
public:
  ParsedZeroAverageConstraints  (const std::string &name,
                                 const unsigned int &n_components=1,
                                 const std::string &component_names = "",
                                 const std::string &default_components = "",
                                 const std::string &default_boundary_components ="");


  void apply_zero_average_constraints(const DoFHandler<dim,spacedim> &dof_handler,
                                      ConstraintMatrix &constraints) const;


  ComponentMask get_boundary_mask () const;


  ComponentMask get_mask () const;


  virtual void declare_parameters (ParameterHandler &prm);

  virtual void parse_parameters_call_back ();



  DeclException2(ExcWrongComponent, unsigned int, unsigned int,
                 << "Wrong component number has been used: "
                 <<  arg1 << " is not in the range [0, "
                 << arg2 <<").");

  DeclException2(ExcWrongVariable, std::string, std::vector<std::string>,
                 << "Wrong variabile name has been used: "
                 <<  arg1 << " does not belong to the knwon variables: "
                 << print(unique(arg2)) <<".");

protected:

  void internal_zero_average_constraints(const DoFHandler<dim,spacedim> &dof_handler,
                                         const ComponentMask mask,
                                         const bool at_boundary,
                                         ConstraintMatrix &constraints) const;

  std::string name;
  std::string str_components;
  std::string str_boundary_components;
  std::string str_component_names;
  std::vector<std::string> _component_names;
  std::vector<std::string> boundary_components;
  std::vector<std::string> components;

  std::vector<bool> mask;
  std::vector<bool> boundary_mask;

  const unsigned int n_components;
};



D2K_NAMESPACE_CLOSE

#endif