%PDF- %PDF-
Direktori : /lib/python3/dist-packages/pythran/pythonic/numpy/ |
Current File : //lib/python3/dist-packages/pythran/pythonic/numpy/nanargmax.hpp |
#ifndef PYTHONIC_NUMPY_NANARGMAX_HPP #define PYTHONIC_NUMPY_NANARGMAX_HPP #include "pythonic/include/numpy/nanargmax.hpp" #include "pythonic/utils/functor.hpp" #include "pythonic/types/ndarray.hpp" #include "pythonic/builtins/ValueError.hpp" #include "pythonic/numpy/isnan.hpp" PYTHONIC_NS_BEGIN namespace numpy { namespace { template <class E, class F> void _nanargmax(E begin, E end, F &max, long &index, long &where, utils::int_<1>) { for (; begin != end; ++begin, ++index) { auto curr = *begin; if (!functor::isnan()(curr) && curr > max) { max = curr; where = index; } } } template <class E, class F, size_t N> void _nanargmax(E begin, E end, F &max, long &index, long &where, utils::int_<N>) { for (; begin != end; ++begin) _nanargmax((*begin).begin(), (*begin).end(), max, index, where, utils::int_<N - 1>()); } } template <class E> long nanargmax(E const &expr) { typename E::dtype max = -std::numeric_limits<typename E::dtype>::infinity(); long where = -1; long index = 0; _nanargmax(expr.begin(), expr.end(), max, index, where, utils::int_<E::value>()); if (where >= 0) return where; else throw types::ValueError("empty sequence"); } } PYTHONIC_NS_END #endif