parsed_dirichlet_bcs.h

Go to the documentation of this file.

//-----------------------------------------------------------
//
//    Copyright (C) 2015 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_dirichlet_bcs_h
#define _d2k_parsed_dirichlet_bcs_h

#include <deal2lkit/config.h>
#include <deal2lkit/parameter_acceptor.h>
#include <deal.II/base/exceptions.h>
#include <deal2lkit/parsed_mapped_functions.h>
#include <deal.II/dofs/dof_handler.h>
#include <deal.II/lac/constraint_matrix.h>
#include <deal.II/numerics/vector_tools.h>

using namespace dealii;

D2K_NAMESPACE_OPEN

template <int dim, int spacedim=dim>
class ParsedDirichletBCs : public ParsedMappedFunctions<spacedim>
{
public:
  ParsedDirichletBCs (const std::string &name = "Dirichlet BCs",
                      const unsigned int &n_components=1,
                      const std::string &component_names = "",
                      const std::string &default_id_components = "0=ALL",
                      const std::string &default_id_functions = "",
                      const std::string &default_constants = "");

  virtual void declare_parameters (ParameterHandler &prm);

  virtual void parse_parameters_call_back ();

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

  void interpolate_boundary_values (const Mapping<dim,spacedim> &mapping,
                                    const DoFHandler<dim,spacedim> &dof_handler,
                                    ConstraintMatrix &constraints) const;

  void interpolate_boundary_values (const DoFHandler<dim,spacedim> &dof_handler,
                                    std::map<types::global_dof_index,double> &d_dofs) const;

  void interpolate_boundary_values (const Mapping<dim,spacedim> &mapping,
                                    const DoFHandler<dim,spacedim> &dof_handler,
                                    std::map<types::global_dof_index,double> &d_dofs) const;

  void project_boundary_values (const DoFHandler<dim,spacedim> &dof_handler,
                                const Quadrature<dim-1> &quadrature,
                                ConstraintMatrix &constraints) const;

  void project_boundary_values (const Mapping<dim,spacedim> &mapping,
                                const DoFHandler<dim,spacedim> &dof_handler,
                                const Quadrature<dim-1> &quadrature,
                                ConstraintMatrix &constraints) const;

  void project_boundary_values (const DoFHandler<dim,spacedim> &dof_handler,
                                const Quadrature<dim-1> &quadrature,
                                std::map<types::global_dof_index,double> &projected_bv) const;

  void project_boundary_values (const Mapping<dim,spacedim> &mapping,
                                const DoFHandler<dim,spacedim> &dof_handler,
                                const Quadrature<dim-1> &quadrature,
                                std::map<types::global_dof_index,double> &projected_bv) const;


  void compute_no_normal_flux_constraints (const DoFHandler<dim,spacedim> &dof_handler,
                                           ConstraintMatrix &constraints) const;
  void compute_no_normal_flux_constraints (const DoFHandler<dim,spacedim> &dof_handler,
                                           const Mapping< dim, spacedim > &mapping,
                                           ConstraintMatrix &constraints) const;
  void compute_nonzero_normal_flux_constraints (const DoFHandler<dim,spacedim> &dof_handler,
                                                ConstraintMatrix &constraints) const;

  void compute_nonzero_normal_flux_constraints (const DoFHandler<dim,spacedim> &dof_handler,
                                                const Mapping< dim, spacedim > &mapping,
                                                ConstraintMatrix &constraints) const;
private:
  const unsigned int n_components;

};

D2K_NAMESPACE_CLOSE

#endif