1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
|
// Copyright (C) 2018-2024 Chris N. Richardson and Garth N. Wells
//
// This file is part of DOLFINx (https://www.fenicsproject.org)
//
// SPDX-License-Identifier: LGPL-3.0-or-later
#include "MPICommWrapper.h"
#include "array.h"
#include "mesh.h"
#include <concepts>
#include <cstdint>
#include <dolfinx/mesh/Mesh.h>
#include <dolfinx/mesh/MeshTags.h>
#include <dolfinx/refinement/interval.h>
#include <dolfinx/refinement/option.h>
#include <dolfinx/refinement/refine.h>
#include <dolfinx/refinement/utils.h>
#include <functional>
#include <nanobind/nanobind.h>
#include <nanobind/ndarray.h>
#include <nanobind/stl/function.h>
#include <nanobind/stl/optional.h>
#include <nanobind/stl/shared_ptr.h>
#include <nanobind/stl/tuple.h>
#include <nanobind/stl/vector.h>
#include <optional>
#include <span>
namespace nb = nanobind;
namespace
{
template <std::floating_point T>
void export_refinement(nb::module_& m)
{
m.def(
"refine",
[](const dolfinx::mesh::Mesh<T>& mesh,
std::optional<
nb::ndarray<const std::int32_t, nb::ndim<1>, nb::c_contig>>
edges,
std::optional<
dolfinx_wrappers::part::impl::PythonCellPartitionFunction>
partitioner,
dolfinx::refinement::Option option)
{
std::optional<std::span<const std::int32_t>> cpp_edges(std::nullopt);
if (edges.has_value())
{
cpp_edges.emplace(
std::span(edges.value().data(), edges.value().size()));
}
dolfinx_wrappers::part::impl::CppCellPartitionFunction cpp_partitioner
= partitioner.has_value()
? dolfinx_wrappers::part::impl::create_cell_partitioner_cpp(
partitioner.value())
: nullptr;
auto [mesh1, cell, facet] = dolfinx::refinement::refine(
mesh, cpp_edges, cpp_partitioner, option);
std::optional<nb::ndarray<std::int32_t, nb::numpy>> python_cell(
std::nullopt);
if (cell.has_value())
{
python_cell.emplace(
dolfinx_wrappers::as_nbarray(std::move(cell.value())));
}
std::optional<nb::ndarray<std::int8_t, nb::numpy>> python_facet(
std::nullopt);
if (facet.has_value())
{
python_facet.emplace(
dolfinx_wrappers::as_nbarray(std::move(facet.value())));
}
return std::tuple{std::move(mesh1), std::move(python_cell),
std::move(python_facet)};
},
nb::arg("mesh"), nb::arg("edges").none(), nb::arg("partitioner").none(),
nb::arg("option"));
}
} // namespace
namespace dolfinx_wrappers
{
void refinement(nb::module_& m)
{
export_refinement<float>(m);
export_refinement<double>(m);
nb::enum_<dolfinx::refinement::Option>(m, "RefinementOption")
.value("none", dolfinx::refinement::Option::none)
.value("parent_facet", dolfinx::refinement::Option::parent_facet)
.value("parent_cell", dolfinx::refinement::Option::parent_cell)
.value("parent_cell_and_facet",
dolfinx::refinement::Option::parent_cell_and_facet);
m.def(
"transfer_facet_meshtag",
[](const dolfinx::mesh::MeshTags<std::int32_t>& parent_meshtag,
std::shared_ptr<const dolfinx::mesh::Topology> topology1,
nb::ndarray<const std::int32_t, nb::ndim<1>, nb::c_contig> parent_cell,
nb::ndarray<const std::int8_t, nb::ndim<1>, nb::c_contig> parent_facet)
{
int tdim = parent_meshtag.topology()->dim();
if (parent_meshtag.dim() != tdim - 1)
throw std::runtime_error("Input meshtag is not facet-based");
auto [entities, values] = dolfinx::refinement::transfer_facet_meshtag(
parent_meshtag, *topology1,
std::span(parent_cell.data(), parent_cell.size()),
std::span(parent_facet.data(), parent_facet.size()));
return dolfinx::mesh::MeshTags<std::int32_t>(
topology1, tdim - 1, std::move(entities), std::move(values));
},
nb::arg("parent_meshtag"), nb::arg("refined_mesh"),
nb::arg("parent_cell"), nb::arg("parent_facet"));
m.def(
"transfer_cell_meshtag",
[](const dolfinx::mesh::MeshTags<std::int32_t>& parent_meshtag,
std::shared_ptr<const dolfinx::mesh::Topology> topology1,
nb::ndarray<const std::int32_t, nb::ndim<1>, nb::c_contig> parent_cell)
{
int tdim = parent_meshtag.topology()->dim();
if (parent_meshtag.dim() != tdim)
throw std::runtime_error("Input meshtag is not cell-based");
if (parent_meshtag.topology()->index_map(tdim)->num_ghosts() > 0)
throw std::runtime_error("Ghosted meshes are not supported");
auto [entities, values] = dolfinx::refinement::transfer_cell_meshtag(
parent_meshtag, *topology1,
std::span(parent_cell.data(), parent_cell.size()));
return dolfinx::mesh::MeshTags<std::int32_t>(
topology1, tdim, std::move(entities), std::move(values));
},
nb::arg("parent_meshtag"), nb::arg("refined_mesh"),
nb::arg("parent_cell"));
}
} // namespace dolfinx_wrappers
|
