Homec4science
Diffusion akantu (bugfixes/60-hybrid-phase-field-formulation-unstable-in-dynamics) Browse

rAKA/python/py_aka_array.hh
d0573c64ba77bugfixes/60-hybrid…

Last Change
Raw File

py_aka_array.hh
View Options

/* -------------------------------------------------------------------------- */
#include "aka_array.hh"
/* -------------------------------------------------------------------------- */
#include <pybind11/numpy.h>
#include <pybind11/pybind11.h>
/* -------------------------------------------------------------------------- */
namespace py = pybind11;
namespace _aka = akantu;
namespace akantu {
template <typename VecType> class Proxy : public VecType {
protected:
using T = typename VecType::value_type;
// deallocate the memory
void deallocate() override final {}
// allocate the memory
void allocate(__attribute__((unused)) UInt size,
__attribute__((unused)) UInt nb_component) override final {}
// allocate and initialize the memory
void allocate(__attribute__((unused)) UInt size,
__attribute__((unused)) UInt nb_component,
__attribute__((unused)) const T & value) override final {}
public:
Proxy(T * data, UInt size, UInt nb_component) {
this->values = data;
this->size_ = size;
this->nb_component = nb_component;
}
Proxy(const Array<T> & src) {
this->values = src.storage();
this->size_ = src.size();
this->nb_component = src.getNbComponent();
}
~Proxy() { this->values = nullptr; }
void resize(UInt /*size*/, const T & /*val */) override final {
AKANTU_EXCEPTION("cannot resize a temporary array");
}
void resize(UInt /*new_size*/) override final {
AKANTU_EXCEPTION("cannot resize a temporary array");
}
void reserve(UInt /*size*/, UInt /*new_size*/) override final {
AKANTU_EXCEPTION("cannot resize a temporary array");
}
};
template <typename T> using vec_proxy = Vector<T>;
template <typename T> using mat_proxy = Matrix<T>;
template <typename T> using array_proxy = Proxy<Array<T>>;
template <typename array> struct ProxyType { using type = Proxy<array>; };
template <typename T> struct ProxyType<Vector<T>> { using type = Vector<T>; };
template <typename T> struct ProxyType<Matrix<T>> { using type = Matrix<T>; };
template <typename array> using ProxyType_t = typename ProxyType<array>::type;
} // namespace akantu
namespace pybind11 {
namespace detail {
template <typename U>
using array_type = array_t<U, array::c_style | array::forcecast>;
template <typename T>
void create_proxy(std::unique_ptr<_aka::vec_proxy<T>> & proxy,
array_type<T> ref) {
proxy =
std::make_unique<_aka::vec_proxy<T>>(ref.mutable_data(), ref.shape(0));
}
template <typename T>
void create_proxy(std::unique_ptr<_aka::mat_proxy<T>> & proxy,
array_type<T> ref) {
proxy = std::make_unique<_aka::mat_proxy<T>>(ref.mutable_data(),
ref.shape(0), ref.shape(1));
}
template <typename T>
void create_proxy(std::unique_ptr<_aka::array_proxy<T>> & proxy,
array_type<T> ref) {
proxy = std::make_unique<_aka::array_proxy<T>>(ref.mutable_data(),
ref.shape(0), ref.shape(1));
}
/* ------------------------------------------------------------------------ */
template <typename T>
py::handle aka_array_cast(const _aka::Array<T> & src,
py::handle base = handle(), bool writeable = true) {
array a;
a = array_type<T>({src.size(), src.getNbComponent()}, src.storage(), base);
if (not writeable)
array_proxy(a.ptr())->flags &= ~detail::npy_api::NPY_ARRAY_WRITEABLE_;
return a.release();
}
template <typename U>
using tensor_type = array_t<U, array::f_style | array::forcecast>;
template <typename T>
py::handle aka_array_cast(const _aka::Vector<T> & src,
py::handle base = handle(), bool writeable = true) {
array a;
a = tensor_type<T>({src.size()}, src.storage(), base);
if (not writeable)
array_proxy(a.ptr())->flags &= ~detail::npy_api::NPY_ARRAY_WRITEABLE_;
return a.release();
}
template <typename T>
py::handle aka_array_cast(const _aka::Matrix<T> & src,
py::handle base = handle(), bool writeable = true) {
array a;
a = tensor_type<T>({src.size(0), src.size(1)}, src.storage(), base);
if (not writeable)
array_proxy(a.ptr())->flags &= ~detail::npy_api::NPY_ARRAY_WRITEABLE_;
return a.release();
}
/* ------------------------------------------------------------------------ */
template <typename VecType>
class [[gnu::visibility("default")]] my_type_caster {
protected:
using T = typename VecType::value_type;
using type = VecType;
using proxy_type = _aka::ProxyType_t<VecType>;
type value;
public:
static PYBIND11_DESCR name() { return type_descr(_("Toto")); };
/**
* Conversion part 1 (Python->C++)
*/
bool load(handle src, bool convert) {
bool need_copy = not isinstance<array_type<T>>(src);
auto && fits = [&](auto && aref) {
auto && dims = aref.ndim();
if (dims < 1 || dims > 2)
return false;
return true;
};
if (not need_copy) {
// We don't need a converting copy, but we also need to check whether
// the strides are compatible with the Ref's stride requirements
auto aref = py::cast<array_type<T>>(src);
if (not fits(aref)) {
return false;
}
copy_or_ref = std::move(aref);
} else {
if (not convert) {
return false;
}
auto copy = array_type<T>::ensure(src);
if (not copy) {
return false;
}
if (not fits(copy)) {
return false;
}
copy_or_ref = std::move(array_type<T>::ensure(src));
loader_life_support::add_patient(copy_or_ref);
}
create_proxy(array_proxy, copy_or_ref);
return true;
}
operator type *() { return array_proxy.get(); }
operator type &() { return *array_proxy; }
template <typename _T>
using cast_op_type = pybind11::detail::cast_op_type<_T>;
/**
* Conversion part 2 (C++ -> Python)
*/
static handle cast(const type & src, return_value_policy policy,
handle parent) {
switch (policy) {
case return_value_policy::copy:
return aka_array_cast<T>(src);
case return_value_policy::reference_internal:
return aka_array_cast<T>(src, parent);
case return_value_policy::reference:
case return_value_policy::automatic:
case return_value_policy::automatic_reference:
return aka_array_cast<T>(src, none());
default:
pybind11_fail("Invalid return_value_policy for ArrayProxy type");
}
}
protected:
std::unique_ptr<proxy_type> array_proxy;
array_type<T> copy_or_ref;
};
/* ------------------------------------------------------------------------ */
// specializations
/* ------------------------------------------------------------------------ */
template <typename T>
struct type_caster<_aka::Array<T>> : public my_type_caster<_aka::Array<T>> {};
template <typename T>
struct type_caster<_aka::Vector<T>> : public my_type_caster<_aka::Vector<T>> {
};
template <typename T>
struct type_caster<_aka::Matrix<T>> : public my_type_caster<_aka::Matrix<T>> {
};
} // namespace detail
} // namespace pybind11
c4science · Help