Commit 3277845a authored by Thibault Hallouin's avatar Thibault Hallouin
Browse files

first version working for 1D arrays only

No related merge requests found
Showing with 287 additions and 0 deletions
+287 -0
.gitignore 0 → 100644
*.egg-info/
.ipynb_checkpoints/
dist/
build/
*.py[cod]
# OS X
.DS_Store
%% Cell type:code id: tags:
``` python
import numpy
import evalhyd
```
%% Cell type:code id: tags:
``` python
sim = numpy.array(
[5.3, 4.2, 5.7, 2.3]
)
obs = numpy.array(
[4.7, 4.3, 5.5, 2.7]
)
print(evalhyd.nse(sim, obs))
```
%% Output
[0.86298077]
%% Cell type:code id: tags:
``` python
sim = numpy.array(
[[5.3, 5.3, 5.3],
[4.2, 4.2, 4.2],
[5.7, 5.7, 5.7],
[2.3, 2.3, 2.3]]
)
obs = numpy.array(
[[4.7],
[4.3],
[5.5],
[2.7]]
)
print(sim.shape, obs.shape)
print(evalhyd.nse(sim, obs, axis=0))
```
%% Output
(4, 3) (4, 1)
[[ 0.86298077 -3.08173077 -8.63221154]]
%% Cell type:code id: tags:
``` python
sim = numpy.array(
[[5.3, 4.2, 5.7, 2.3],
[5.3, 4.2, 5.7, 2.3],
[5.3, 4.2, 5.7, 2.3]]
)
obs = numpy.array(
[[4.7, 4.3, 5.5, 2.7]]
)
print(sim.shape, obs.shape)
print(evalhyd.nse(sim, obs, axis=1))
```
%% Output
(3, 4) (1, 4)
[[ 0.86298077]
[-19.07451923]
[-19.07451923]]
setup.py 0 → 100644
from setuptools import setup, Extension
from setuptools.command.build_ext import build_ext
import sys
import os
import setuptools
__version__ = '0.0.1'
class get_pybind_include(object):
"""Helper class to determine the pybind11 include path
The purpose of this class is to postpone importing pybind11
until it is actually installed, so that the ``get_include()``
method can be invoked. """
def __init__(self, user=False):
self.user = user
def __str__(self):
import pybind11
return pybind11.get_include(self.user)
class get_numpy_include(object):
"""Helper class to determine the numpy include path
The purpose of this class is to postpone importing numpy
until it is actually installed, so that the ``get_include()``
method can be invoked. """
def __init__(self):
pass
def __str__(self):
import numpy as np
return np.get_include()
ext_modules = [
Extension(
'evalhyd',
['src/evalhyd-python.cpp'],
include_dirs=[
# Path to pybind11 headers
get_pybind_include(),
get_pybind_include(user=True),
get_numpy_include(),
os.path.join(os.getcwd(), 'deps', 'evalhyd', 'deps', 'xtl', 'include'),
os.path.join(os.getcwd(), 'deps', 'evalhyd', 'deps', 'xtensor', 'include'),
os.path.join(os.getcwd(), 'deps', 'xtensor-python', 'include'),
os.path.join(os.getcwd(), 'deps', 'evalhyd', 'include'),
os.path.join(sys.prefix, 'include'),
os.path.join(sys.prefix, 'Library', 'include')
],
language='c++'
),
]
def has_flag(compiler, flagname):
"""Return a boolean indicating whether a flag name is supported on
the specified compiler.
"""
import tempfile
with tempfile.NamedTemporaryFile('w', suffix='.cpp') as f:
f.write('int main (int argc, char **argv) { return 0; }')
try:
compiler.compile([f.name], extra_postargs=[flagname])
except setuptools.distutils.errors.CompileError:
return False
return True
def cpp_flag(compiler):
"""Return the -std=c++14 compiler flag and errors when the flag is
no available.
"""
if has_flag(compiler, '-std=c++14'):
return '-std=c++14'
else:
raise RuntimeError('C++14 support is required by xtensor!')
class BuildExt(build_ext):
"""A custom build extension for adding compiler-specific options."""
c_opts = {
'msvc': ['/EHsc'],
'unix': [],
}
if sys.platform == 'darwin':
c_opts['unix'] += ['-stdlib=libc++', '-mmacosx-version-min=10.7']
def build_extensions(self):
ct = self.compiler.compiler_type
opts = self.c_opts.get(ct, [])
if ct == 'unix':
opts.append('-DVERSION_INFO="%s"' % self.distribution.get_version())
opts.append(cpp_flag(self.compiler))
if has_flag(self.compiler, '-fvisibility=hidden'):
opts.append('-fvisibility=hidden')
elif ct == 'msvc':
opts.append('/DVERSION_INFO=\\"%s\\"' % self.distribution.get_version())
for ext in self.extensions:
ext.extra_compile_args = opts
build_ext.build_extensions(self)
setup(
name='evalhyd-python',
version=__version__,
author='Thibault Hallouin',
author_email='thibault.hallouin@inrae.fr',
url='https://gitlab.irstea.fr/thibault.hallouin/evalhyd-python',
description='EvalHyd in Python',
long_description='An evaluator for hydrological simulations/forecasts.',
ext_modules=ext_modules,
install_requires=['pybind11>=2.0.1', 'numpy'],
cmdclass={'build_ext': BuildExt},
zip_safe=False,
)
#include <pybind11/pybind11.h>
#define FORCE_IMPORT_ARRAY
#include <xtensor-python/pyarray.hpp>
#include "nse.hpp"
namespace py = pybind11;
// Python Module and Docstrings
PYBIND11_MODULE(evalhyd, m)
{
xt::import_numpy();
m.doc() = "evaluator for hydrological simulations/forecasts";
m.def("nse", nse<xt::pyarray<double>>, "Return the Nash-Sutcliffe Efficiency (NSE)",
py::arg("sim"), py::arg("obs"), py::arg("axis") = 0);
}
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment