Homec4science
Diffusion akantu (bugfixes/bug_fix_in_synchronization_of_mat_vec_mul_product_when_assembled_to_residual) Browse

rAKA/python/py_mesh.cc
b4da2bec742fbugfixes/bug_fix_in_syn…

Last Change
Raw File

py_mesh.cc
View Options

/* -------------------------------------------------------------------------- */
#include "aka_config.hh"
/* -------------------------------------------------------------------------- */
#include "py_aka_array.hh"
/* -------------------------------------------------------------------------- */
#include <mesh.hh>
#include <mesh_utils.hh>
#include <mesh_accessor.hh>
/* -------------------------------------------------------------------------- */
#include <pybind11/pybind11.h>
/* -------------------------------------------------------------------------- */
namespace py = pybind11;
/* -------------------------------------------------------------------------- */
namespace akantu {
/* -------------------------------------------------------------------------- */
template <typename T>
void register_element_type_map_array(py::module & mod,
const std::string & name) {
py::class_<ElementTypeMapArray<T>, std::shared_ptr<ElementTypeMapArray<T>>>(
mod, ("ElementTypeMapArray" + name).c_str())
.def(
"__call__",
[](ElementTypeMapArray<T> & self, ElementType type,
GhostType ghost_type) -> decltype(auto) {
return self(type, ghost_type);
},
py::arg("type"), py::arg("ghost_type") = _not_ghost,
py::return_value_policy::reference)
.def(
"elementTypes",
[](ElementTypeMapArray<T> & self, UInt _dim, GhostType _ghost_type,
ElementKind _kind) -> decltype(auto) {
auto types = self.elementTypes(_dim, _ghost_type, _kind);
std::vector<ElementType> _types;
for (auto && t : types) {
_types.push_back(t);
}
return _types;
},
py::arg("dim") = _all_dimensions, py::arg("ghost_type") = _not_ghost,
py::arg("kind") = _ek_regular);
}
/* -------------------------------------------------------------------------- */
void register_mesh(py::module & mod) {
register_element_type_map_array<Real>(mod, "Real");
register_element_type_map_array<UInt>(mod, "UInt");
py::class_<MeshData>(mod, "MeshData")
.def(
"getElementalDataUInt",
[](MeshData & _this, const ID & name) -> ElementTypeMapArray<UInt> & {
return _this.getElementalData<UInt>(name);
},
py::return_value_policy::reference);
py::class_<Mesh, GroupManager, Dumpable, MeshData>(mod, "Mesh",
py::multiple_inheritance())
.def(py::init<UInt, const ID &, const MemoryID &>(),
py::arg("spatial_dimension"), py::arg("id") = "mesh",
py::arg("memory_id") = 0)
.def("read", &Mesh::read, py::arg("filename"),
py::arg("mesh_io_type") = _miot_auto, "read the mesh from a file")
.def(
"getNodes",
[](Mesh & self) -> decltype(auto) { return self.getNodes(); },
py::return_value_policy::reference)
.def("getNbNodes", &Mesh::getNbNodes)
.def(
"getConnectivity",
[](Mesh & self, ElementType type) -> decltype(auto) {
return self.getConnectivity(type);
},
py::return_value_policy::reference)
.def("addConnectivityType",
[](Mesh & self, ElementType type, GhostType ghost_type) -> void {
self.addConnectivityType(type, ghost_type);
},
py::arg("type"), py::arg("ghost_type") = _not_ghost)
.def("distribute", [](Mesh & self) { self.distribute(); })
.def("makePeriodic",
[](Mesh & self, const SpatialDirection & direction) {
self.makePeriodic(direction);
})
.def(
"getNbElement",
[](Mesh & self, const UInt spatial_dimension,
GhostType ghost_type, ElementKind kind) {
return self.getNbElement(spatial_dimension, ghost_type, kind);
},
py::arg("spatial_dimension") = _all_dimensions,
py::arg("ghost_type") = _not_ghost, py::arg("kind") = _ek_not_defined)
.def(
"getNbElement",
[](Mesh & self, ElementType type,
GhostType ghost_type) {
return self.getNbElement(type, ghost_type);
},
py::arg("type"), py::arg("ghost_type") = _not_ghost)
.def_static("getSpatialDimension", [](ElementType & type) {
return Mesh::getSpatialDimension(type);
})
/* Extra Normal specific stuff
*/
.def("hasExtraNormal",
[](Mesh & self, ElementType type, GhostType ghost_type) -> bool {
return self.hasData<Real>("extra_normal", type, ghost_type);
},
py::arg("type"), py::arg("ghost_type") = _not_ghost)
.def("makeExtraNormal",
[](Mesh & self, ElementType type, UInt dim) -> decltype(auto) {
if(not self.hasData<Real>("extra_normal", type)){
self.getElementalData<Real>("extra_normal")
.initialize(self, _element_kind = _ek_structural,
_nb_component = dim, _with_nb_element = true,
_default_value = 0.);
};
return self.getData<Real>("extra_normal", type);
},
R"(If the "extra_normal" mesh data does not exists create it and return it.)",
py::arg("type"), py::arg("spatial_dimension") = 3,
py::return_value_policy::reference)
.def("getExtraNormal",
[](Mesh & self, ElementType type) -> decltype(auto) {
return self.getData<Real>("extra_normal", type);
},
py::arg("type"),
py::return_value_policy::reference)
;
/* ------------------------------------------------------------------------ */
py::class_<MeshUtils>(mod, "MeshUtils")
.def_static("buildFacets", &MeshUtils::buildFacets);
py::class_<MeshAccessor>(mod, "MeshAccessor")
.def(
py::init<Mesh &>(), py::arg("mesh"))
.def("resizeConnectivity",
[](MeshAccessor & self, UInt new_size, ElementType type, GhostType gt) -> void {
self.resizeConnectivity(new_size, type, gt);},
py::arg("new_size"), py::arg("type"), py::arg("ghost_type") = _not_ghost)
.def("resizeNodes",
[](MeshAccessor& self, UInt new_size) -> void {
self.resizeNodes(new_size);},
py::arg("new_size"))
.def("makeReady", &MeshAccessor::makeReady);
}
} // namespace akantu
c4science · Help