1
// Copyright (C) 2012 Anders Logg
3
// This file is part of DOLFIN.
5
// DOLFIN is free software: you can redistribute it and/or modify
6
// it under the terms of the GNU Lesser General Public License as published by
7
// the Free Software Foundation, either version 3 of the License, or
8
// (at your option) any later version.
10
// DOLFIN is distributed in the hope that it will be useful,
11
// but WITHOUT ANY WARRANTY; without even the implied warranty of
12
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13
// GNU Lesser General Public License for more details.
15
// You should have received a copy of the GNU Lesser General Public License
16
// along with DOLFIN. If not, see <http://www.gnu.org/licenses/>.
18
// Modified by Benjamin Kehlet, 2013
20
// First added: 2012-04-11
21
// Last changed: 2013-04-18
23
#include "CSGGeometry.h"
24
#include <dolfin/common/NoDeleter.h>
26
using namespace dolfin;
28
//-----------------------------------------------------------------------------
29
CSGGeometry::CSGGeometry()
33
//-----------------------------------------------------------------------------
34
CSGGeometry::~CSGGeometry()
38
//-----------------------------------------------------------------------------
39
void CSGGeometry::set_subdomain(std::size_t i, std::shared_ptr<CSGGeometry> s)
41
dolfin_assert(dim() == s->dim());
45
error("Setting reserved CSG subdomain (0)");
48
// Check if i already used
49
std::list<std::pair<std::size_t, std::shared_ptr<const CSGGeometry> > >::iterator it = subdomains.begin();
50
while (it != subdomains.end())
54
warning("Double declaration of CSG subdomain with index %u.", i);
56
// Remove existing declaration
57
it = subdomains.erase(it);
63
subdomains.push_back(std::make_pair(i, s));
65
//-----------------------------------------------------------------------------
66
void CSGGeometry::set_subdomain(std::size_t i, CSGGeometry& s)
68
set_subdomain(i, reference_to_no_delete_pointer(s));
70
//-----------------------------------------------------------------------------
71
bool CSGGeometry::has_subdomains() const
73
return subdomains.size() > 0;