system_subset_by_subdomain.h

Go to the documentation of this file.

// The libMesh Finite Element Library.
// Copyright (C) 2002-2018 Benjamin S. Kirk, John W. Peterson, Roy H. Stogner

// This library is free software; you can 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.

// This library 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
// Lesser General Public License for more details.

// You should have received a copy of the GNU Lesser General Public
// License along with this library; if not, write to the Free Software
// Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA



#ifndef LIBMESH_SYSTEM_SUBSET_BY_SUBDOMAIN_H
#define LIBMESH_SYSTEM_SUBSET_BY_SUBDOMAIN_H

// Local Includes
#include "libmesh/system_subset.h"
#include "libmesh/id_types.h"
#include "libmesh/parallel_object.h"

// C++ includes
#include <cstddef>
#include <set>

namespace libMesh
{

// Forward Declarations

class SystemSubsetBySubdomain : public SystemSubset,
                                public ParallelObject
{
public:

  class SubdomainSelection : public ReferenceCountedObject<SubdomainSelection>
  {
  public:

    SubdomainSelection ();

    virtual ~SubdomainSelection ();

    virtual bool operator() (const subdomain_id_type & subdomain_id) const = 0;

  private:
    SubdomainSelection (const SubdomainSelection &);

    SubdomainSelection & operator = (const SubdomainSelection &);
  }; // subclass \p SubdomainSelection



  class SubdomainSelectionByList : public SubdomainSelection
  {
  public:
    explicit
    SubdomainSelectionByList (const std::set<subdomain_id_type> & list);

    virtual bool operator() (const subdomain_id_type & subdomain_id) const override;

  protected:
    const std::set<subdomain_id_type> & _list;
  };

  SystemSubsetBySubdomain (const System & system,
                           const SubdomainSelection & subdomain_selection,
                           const std::set<unsigned int> * const var_nums = nullptr);

  SystemSubsetBySubdomain (const System & system,
                           const std::set<subdomain_id_type> & subdomain_ids,
                           const std::set<unsigned int> * const var_nums = nullptr);

  virtual ~SystemSubsetBySubdomain ();

  virtual const std::vector<unsigned int> & dof_ids () const override;

  void init (const SubdomainSelection & subdomain_selection);

  void init (const std::set<subdomain_id_type> & subdomain_ids);

protected:

  void set_var_nums (const std::set<unsigned int> * const var_nums);

  std::set<unsigned int> _var_nums;

  std::vector<unsigned int> _dof_ids;

}; // class SystemSubset

} // namespace libMesh

#endif // LIBMESH_SYSTEM_SUBSET_BY_SUBDOMAIN_H