HydPy 6.2.dev1__cp313-cp313-win_amd64.whl

hydpy/cythons/autogen/c_conv_idw.pyx

@@ -0,0 +1,264 @@
+#!python
+# distutils: define_macros=NPY_NO_DEPRECATED_API=NPY_1_7_API_VERSION
+# cython: language_level=3
+# cython: cpow=True
+# cython: boundscheck=False
+# cython: wraparound=False
+# cython: initializedcheck=False
+# cython: cdivision=True
+from typing import Optional
+import numpy
+cimport numpy
+from libc.math cimport exp, fabs, log, sin, cos, tan, tanh, asin, acos, atan, isnan, isinf
+from libc.math cimport NAN as nan
+from libc.math cimport INFINITY as inf
+import cython
+from cpython.mem cimport PyMem_Malloc
+from cpython.mem cimport PyMem_Realloc
+from cpython.mem cimport PyMem_Free
+from hydpy.cythons.autogen cimport configutils
+from hydpy.cythons.autogen cimport interfaceutils
+from hydpy.cythons.autogen cimport interputils
+from hydpy.cythons.autogen import pointerutils
+from hydpy.cythons.autogen cimport pointerutils
+from hydpy.cythons.autogen cimport quadutils
+from hydpy.cythons.autogen cimport rootutils
+from hydpy.cythons.autogen cimport smoothutils
+from hydpy.cythons.autogen cimport masterinterface
+
+
+cdef void do_nothing(Model model)  noexcept nogil:
+    pass
+
+cpdef get_wrapper():
+    cdef CallbackWrapper wrapper = CallbackWrapper()
+    wrapper.callback = do_nothing
+    return wrapper
+
+@cython.final
+cdef class Parameters:
+    pass
+@cython.final
+cdef class ControlParameters:
+    pass
+@cython.final
+cdef class DerivedParameters:
+    pass
+@cython.final
+cdef class Sequences:
+    pass
+@cython.final
+cdef class InletSequences:
+    cpdef inline alloc(self, name, numpy.int64_t length):
+        if name == "inputs":
+            self._inputs_length_0 = length
+            self._inputs_ready = numpy.full(length, 0, dtype=numpy.int64)
+            self.inputs = <double**> PyMem_Malloc(length * sizeof(double*))
+    cpdef inline dealloc(self, name):
+        if name == "inputs":
+            PyMem_Free(self.inputs)
+    cpdef inline set_pointer1d(self, str name, pointerutils.Double value, numpy.int64_t idx):
+        cdef pointerutils.PDouble pointer = pointerutils.PDouble(value)
+        if name == "inputs":
+            self.inputs[idx] = pointer.p_value
+            self._inputs_ready[idx] = 1
+    cpdef get_value(self, str name):
+        cdef numpy.int64_t idx
+        if name == "inputs":
+            values = numpy.empty(self.len_inputs)
+            for idx in range(self.len_inputs):
+                pointerutils.check0(self._inputs_length_0)
+                if self._inputs_ready[idx] == 0:
+                    pointerutils.check1(self._inputs_length_0, idx)
+                    pointerutils.check2(self._inputs_ready, idx)
+                values[idx] = self.inputs[idx][0]
+            return values
+    cpdef set_value(self, str name, value):
+        if name == "inputs":
+            for idx in range(self.len_inputs):
+                pointerutils.check0(self._inputs_length_0)
+                if self._inputs_ready[idx] == 0:
+                    pointerutils.check1(self._inputs_length_0, idx)
+                    pointerutils.check2(self._inputs_ready, idx)
+                self.inputs[idx][0] = value[idx]
+@cython.final
+cdef class FluxSequences:
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t jdx0
+        cdef numpy.int64_t k
+        if self._inputs_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._inputs_length_0):
+                self.inputs[jdx0] = self._inputs_ncarray[k]
+                k += 1
+        elif self._inputs_ramflag:
+            for jdx0 in range(self._inputs_length_0):
+                self.inputs[jdx0] = self._inputs_array[idx, jdx0]
+        if self._outputs_diskflag_reading:
+            k = 0
+            for jdx0 in range(self._outputs_length_0):
+                self.outputs[jdx0] = self._outputs_ncarray[k]
+                k += 1
+        elif self._outputs_ramflag:
+            for jdx0 in range(self._outputs_length_0):
+                self.outputs[jdx0] = self._outputs_array[idx, jdx0]
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil:
+        cdef numpy.int64_t jdx0
+        cdef numpy.int64_t k
+        if self._inputs_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._inputs_length_0):
+                self._inputs_ncarray[k] = self.inputs[jdx0]
+                k += 1
+        if self._inputs_ramflag:
+            for jdx0 in range(self._inputs_length_0):
+                self._inputs_array[idx, jdx0] = self.inputs[jdx0]
+        if self._outputs_diskflag_writing:
+            k = 0
+            for jdx0 in range(self._outputs_length_0):
+                self._outputs_ncarray[k] = self.outputs[jdx0]
+                k += 1
+        if self._outputs_ramflag:
+            for jdx0 in range(self._outputs_length_0):
+                self._outputs_array[idx, jdx0] = self.outputs[jdx0]
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value):
+        pass
+    cpdef inline void update_outputs(self) noexcept nogil:
+        pass
+@cython.final
+cdef class OutletSequences:
+    cpdef inline alloc(self, name, numpy.int64_t length):
+        if name == "outputs":
+            self._outputs_length_0 = length
+            self._outputs_ready = numpy.full(length, 0, dtype=numpy.int64)
+            self.outputs = <double**> PyMem_Malloc(length * sizeof(double*))
+    cpdef inline dealloc(self, name):
+        if name == "outputs":
+            PyMem_Free(self.outputs)
+    cpdef inline set_pointer1d(self, str name, pointerutils.Double value, numpy.int64_t idx):
+        cdef pointerutils.PDouble pointer = pointerutils.PDouble(value)
+        if name == "outputs":
+            self.outputs[idx] = pointer.p_value
+            self._outputs_ready[idx] = 1
+    cpdef get_value(self, str name):
+        cdef numpy.int64_t idx
+        if name == "outputs":
+            values = numpy.empty(self.len_outputs)
+            for idx in range(self.len_outputs):
+                pointerutils.check0(self._outputs_length_0)
+                if self._outputs_ready[idx] == 0:
+                    pointerutils.check1(self._outputs_length_0, idx)
+                    pointerutils.check2(self._outputs_ready, idx)
+                values[idx] = self.outputs[idx][0]
+            return values
+    cpdef set_value(self, str name, value):
+        if name == "outputs":
+            for idx in range(self.len_outputs):
+                pointerutils.check0(self._outputs_length_0)
+                if self._outputs_ready[idx] == 0:
+                    pointerutils.check1(self._outputs_length_0, idx)
+                    pointerutils.check2(self._outputs_ready, idx)
+                self.outputs[idx][0] = value[idx]
+@cython.final
+cdef class Model:
+    cpdef inline void simulate(self, numpy.int64_t idx)  noexcept nogil:
+        self.idx_sim = idx
+        self.update_inlets()
+        self.run()
+        self.update_outlets()
+        self.update_outputs()
+    cpdef void reset_reuseflags(self) noexcept nogil:
+        pass
+    cpdef void save_data(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        self.sequences.fluxes.save_data(idx)
+    cpdef inline void run(self) noexcept nogil:
+        self.calc_outputs_v2()
+    cpdef inline void update_inlets(self) noexcept nogil:
+        self.pick_inputs_v1()
+    cpdef inline void update_outlets(self) noexcept nogil:
+        self.pass_outputs_v1()
+    cpdef inline void update_receivers(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        pass
+    cpdef inline void update_senders(self, numpy.int64_t idx) noexcept nogil:
+        self.idx_sim = idx
+        pass
+    cpdef inline void update_outputs(self) noexcept nogil:
+        pass
+    cpdef inline void pick_inputs_v1(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmbinputs):
+            self.sequences.fluxes.inputs[idx] = self.sequences.inlets.inputs[idx][0]
+    cpdef inline void calc_outputs_v2(self) noexcept nogil:
+        self.interpolate_inversedistance_v1(self.sequences.fluxes.inputs, self.sequences.fluxes.outputs)
+    cpdef inline void interpolate_inversedistance_v1(self, double[:] inputs, double[:] outputs) noexcept nogil:
+        cdef numpy.int64_t idx_in
+        cdef numpy.int64_t idx_try
+        cdef numpy.int64_t counter_inf
+        cdef double d_sumvalues_inf
+        cdef double d_sumvalues
+        cdef double d_sumweights
+        cdef numpy.int64_t idx_out
+        for idx_out in range(self.parameters.derived.nmboutputs):
+            d_sumweights = 0.0
+            d_sumvalues = 0.0
+            d_sumvalues_inf = 0.0
+            counter_inf = 0
+            for idx_try in range(self.parameters.control.maxnmbinputs):
+                idx_in = self.parameters.derived.proximityorder[idx_out, idx_try]
+                if not isnan(inputs[idx_in]):
+                    if isinf(self.parameters.derived.weights[idx_out, idx_try]):
+                        d_sumvalues_inf = d_sumvalues_inf + (inputs[idx_in])
+                        counter_inf = counter_inf + (1)
+                    else:
+                        d_sumweights = d_sumweights + (self.parameters.derived.weights[idx_out, idx_try])
+                        d_sumvalues = d_sumvalues + (self.parameters.derived.weights[idx_out, idx_try] * inputs[idx_in])
+            if counter_inf:
+                outputs[idx_out] = d_sumvalues_inf / counter_inf
+            elif d_sumweights:
+                outputs[idx_out] = d_sumvalues / d_sumweights
+            else:
+                outputs[idx_out] = nan
+    cpdef inline void pass_outputs_v1(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmboutputs):
+            self.sequences.outlets.outputs[idx][0] = self.sequences.fluxes.outputs[idx]
+    cpdef inline void pick_inputs(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmbinputs):
+            self.sequences.fluxes.inputs[idx] = self.sequences.inlets.inputs[idx][0]
+    cpdef inline void calc_outputs(self) noexcept nogil:
+        self.interpolate_inversedistance_v1(self.sequences.fluxes.inputs, self.sequences.fluxes.outputs)
+    cpdef inline void interpolate_inversedistance(self, double[:] inputs, double[:] outputs) noexcept nogil:
+        cdef numpy.int64_t idx_in
+        cdef numpy.int64_t idx_try
+        cdef numpy.int64_t counter_inf
+        cdef double d_sumvalues_inf
+        cdef double d_sumvalues
+        cdef double d_sumweights
+        cdef numpy.int64_t idx_out
+        for idx_out in range(self.parameters.derived.nmboutputs):
+            d_sumweights = 0.0
+            d_sumvalues = 0.0
+            d_sumvalues_inf = 0.0
+            counter_inf = 0
+            for idx_try in range(self.parameters.control.maxnmbinputs):
+                idx_in = self.parameters.derived.proximityorder[idx_out, idx_try]
+                if not isnan(inputs[idx_in]):
+                    if isinf(self.parameters.derived.weights[idx_out, idx_try]):
+                        d_sumvalues_inf = d_sumvalues_inf + (inputs[idx_in])
+                        counter_inf = counter_inf + (1)
+                    else:
+                        d_sumweights = d_sumweights + (self.parameters.derived.weights[idx_out, idx_try])
+                        d_sumvalues = d_sumvalues + (self.parameters.derived.weights[idx_out, idx_try] * inputs[idx_in])
+            if counter_inf:
+                outputs[idx_out] = d_sumvalues_inf / counter_inf
+            elif d_sumweights:
+                outputs[idx_out] = d_sumvalues / d_sumweights
+            else:
+                outputs[idx_out] = nan
+    cpdef inline void pass_outputs(self) noexcept nogil:
+        cdef numpy.int64_t idx
+        for idx in range(self.parameters.derived.nmboutputs):
+            self.sequences.outlets.outputs[idx][0] = self.sequences.fluxes.outputs[idx]

hydpy/cythons/autogen/c_conv_idw_ed.pxd

@@ -0,0 +1,197 @@
+#!python
+# distutils: define_macros=NPY_NO_DEPRECATED_API=NPY_1_7_API_VERSION
+# cython: language_level=3
+# cython: cpow=True
+# cython: boundscheck=False
+# cython: wraparound=False
+# cython: initializedcheck=False
+# cython: cdivision=True
+from typing import Optional
+import numpy
+cimport numpy
+from libc.math cimport exp, fabs, log, sin, cos, tan, tanh, asin, acos, atan, isnan, isinf
+from libc.math cimport NAN as nan
+from libc.math cimport INFINITY as inf
+import cython
+from cpython.mem cimport PyMem_Malloc
+from cpython.mem cimport PyMem_Realloc
+from cpython.mem cimport PyMem_Free
+from hydpy.cythons.autogen cimport configutils
+from hydpy.cythons.autogen cimport interfaceutils
+from hydpy.cythons.autogen cimport interputils
+from hydpy.cythons.autogen import pointerutils
+from hydpy.cythons.autogen cimport pointerutils
+from hydpy.cythons.autogen cimport quadutils
+from hydpy.cythons.autogen cimport rootutils
+from hydpy.cythons.autogen cimport smoothutils
+from hydpy.cythons.autogen cimport masterinterface
+ctypedef void (*CallbackType) (Model)  noexcept nogil
+cdef class CallbackWrapper:
+    cdef CallbackType callback
+@cython.final
+cdef class Parameters:
+    cdef public ControlParameters control
+    cdef public DerivedParameters derived
+@cython.final
+cdef class ControlParameters:
+    cdef public double[:,:] inputcoordinates
+    cdef public double[:,:] outputcoordinates
+    cdef public double[:] inputheights
+    cdef public double[:] outputheights
+    cdef public numpy.int64_t maxnmbinputs
+    cdef public numpy.int64_t minnmbinputs
+    cdef public double defaultconstant
+    cdef public double defaultfactor
+    cdef public double power
+@cython.final
+cdef class DerivedParameters:
+    cdef public numpy.int64_t nmbinputs
+    cdef public numpy.int64_t nmboutputs
+    cdef public double[:,:] distances
+    cdef public numpy.int64_t[:,:] proximityorder
+    cdef public double[:,:] weights
+@cython.final
+cdef class Sequences:
+    cdef public InletSequences inlets
+    cdef public FluxSequences fluxes
+    cdef public OutletSequences outlets
+@cython.final
+cdef class InletSequences:
+    cdef double **inputs
+    cdef public numpy.int64_t len_inputs
+    cdef public numpy.int64_t[:] _inputs_ready
+    cdef public numpy.int64_t _inputs_ndim
+    cdef public numpy.int64_t _inputs_length
+    cdef public numpy.int64_t _inputs_length_0
+    cpdef inline alloc(self, name, numpy.int64_t length)
+    cpdef inline dealloc(self, name)
+    cpdef inline set_pointer1d(self, str name, pointerutils.Double value, numpy.int64_t idx)
+    cpdef get_value(self, str name)
+    cpdef set_value(self, str name, value)
+@cython.final
+cdef class FluxSequences:
+    cdef public double[:] inputs
+    cdef public numpy.int64_t _inputs_ndim
+    cdef public numpy.int64_t _inputs_length
+    cdef public numpy.int64_t _inputs_length_0
+    cdef public bint _inputs_ramflag
+    cdef public double[:,:] _inputs_array
+    cdef public bint _inputs_diskflag_reading
+    cdef public bint _inputs_diskflag_writing
+    cdef public double[:] _inputs_ncarray
+    cdef public double actualconstant
+    cdef public numpy.int64_t _actualconstant_ndim
+    cdef public numpy.int64_t _actualconstant_length
+    cdef public bint _actualconstant_ramflag
+    cdef public double[:] _actualconstant_array
+    cdef public bint _actualconstant_diskflag_reading
+    cdef public bint _actualconstant_diskflag_writing
+    cdef public double[:] _actualconstant_ncarray
+    cdef public bint _actualconstant_outputflag
+    cdef double *_actualconstant_outputpointer
+    cdef public double actualfactor
+    cdef public numpy.int64_t _actualfactor_ndim
+    cdef public numpy.int64_t _actualfactor_length
+    cdef public bint _actualfactor_ramflag
+    cdef public double[:] _actualfactor_array
+    cdef public bint _actualfactor_diskflag_reading
+    cdef public bint _actualfactor_diskflag_writing
+    cdef public double[:] _actualfactor_ncarray
+    cdef public bint _actualfactor_outputflag
+    cdef double *_actualfactor_outputpointer
+    cdef public double[:] inputpredictions
+    cdef public numpy.int64_t _inputpredictions_ndim
+    cdef public numpy.int64_t _inputpredictions_length
+    cdef public numpy.int64_t _inputpredictions_length_0
+    cdef public bint _inputpredictions_ramflag
+    cdef public double[:,:] _inputpredictions_array
+    cdef public bint _inputpredictions_diskflag_reading
+    cdef public bint _inputpredictions_diskflag_writing
+    cdef public double[:] _inputpredictions_ncarray
+    cdef public double[:] outputpredictions
+    cdef public numpy.int64_t _outputpredictions_ndim
+    cdef public numpy.int64_t _outputpredictions_length
+    cdef public numpy.int64_t _outputpredictions_length_0
+    cdef public bint _outputpredictions_ramflag
+    cdef public double[:,:] _outputpredictions_array
+    cdef public bint _outputpredictions_diskflag_reading
+    cdef public bint _outputpredictions_diskflag_writing
+    cdef public double[:] _outputpredictions_ncarray
+    cdef public double[:] inputresiduals
+    cdef public numpy.int64_t _inputresiduals_ndim
+    cdef public numpy.int64_t _inputresiduals_length
+    cdef public numpy.int64_t _inputresiduals_length_0
+    cdef public bint _inputresiduals_ramflag
+    cdef public double[:,:] _inputresiduals_array
+    cdef public bint _inputresiduals_diskflag_reading
+    cdef public bint _inputresiduals_diskflag_writing
+    cdef public double[:] _inputresiduals_ncarray
+    cdef public double[:] outputresiduals
+    cdef public numpy.int64_t _outputresiduals_ndim
+    cdef public numpy.int64_t _outputresiduals_length
+    cdef public numpy.int64_t _outputresiduals_length_0
+    cdef public bint _outputresiduals_ramflag
+    cdef public double[:,:] _outputresiduals_array
+    cdef public bint _outputresiduals_diskflag_reading
+    cdef public bint _outputresiduals_diskflag_writing
+    cdef public double[:] _outputresiduals_ncarray
+    cdef public double[:] outputs
+    cdef public numpy.int64_t _outputs_ndim
+    cdef public numpy.int64_t _outputs_length
+    cdef public numpy.int64_t _outputs_length_0
+    cdef public bint _outputs_ramflag
+    cdef public double[:,:] _outputs_array
+    cdef public bint _outputs_diskflag_reading
+    cdef public bint _outputs_diskflag_writing
+    cdef public double[:] _outputs_ncarray
+    cpdef inline void load_data(self, numpy.int64_t idx)  noexcept nogil
+    cpdef inline void save_data(self, numpy.int64_t idx)  noexcept nogil
+    cpdef inline set_pointeroutput(self, str name, pointerutils.PDouble value)
+    cpdef inline void update_outputs(self) noexcept nogil
+@cython.final
+cdef class OutletSequences:
+    cdef double **outputs
+    cdef public numpy.int64_t len_outputs
+    cdef public numpy.int64_t[:] _outputs_ready
+    cdef public numpy.int64_t _outputs_ndim
+    cdef public numpy.int64_t _outputs_length
+    cdef public numpy.int64_t _outputs_length_0
+    cpdef inline alloc(self, name, numpy.int64_t length)
+    cpdef inline dealloc(self, name)
+    cpdef inline set_pointer1d(self, str name, pointerutils.Double value, numpy.int64_t idx)
+    cpdef get_value(self, str name)
+    cpdef set_value(self, str name, value)
+@cython.final
+cdef class Model:
+    cdef public numpy.int64_t idx_sim
+    cdef public Parameters parameters
+    cdef public Sequences sequences
+    cpdef inline void simulate(self, numpy.int64_t idx)  noexcept nogil
+    cpdef void reset_reuseflags(self) noexcept nogil
+    cpdef void save_data(self, numpy.int64_t idx) noexcept nogil
+    cpdef inline void run(self) noexcept nogil
+    cpdef inline void update_inlets(self) noexcept nogil
+    cpdef inline void update_outlets(self) noexcept nogil
+    cpdef inline void update_receivers(self, numpy.int64_t idx) noexcept nogil
+    cpdef inline void update_senders(self, numpy.int64_t idx) noexcept nogil
+    cpdef inline void update_outputs(self) noexcept nogil
+    cpdef inline void pick_inputs_v1(self) noexcept nogil
+    cpdef inline void calc_actualconstant_actualfactor_v1(self) noexcept nogil
+    cpdef inline void calc_inputpredictions_v1(self) noexcept nogil
+    cpdef inline void calc_outputpredictions_v1(self) noexcept nogil
+    cpdef inline void calc_inputresiduals_v1(self) noexcept nogil
+    cpdef inline void calc_outputresiduals_v1(self) noexcept nogil
+    cpdef inline void calc_outputs_v3(self) noexcept nogil
+    cpdef inline double return_mean_v1(self, double[:] values, double[:] mask, numpy.int64_t number) noexcept nogil
+    cpdef inline void interpolate_inversedistance_v1(self, double[:] inputs, double[:] outputs) noexcept nogil
+    cpdef inline void pass_outputs_v1(self) noexcept nogil
+    cpdef inline void pick_inputs(self) noexcept nogil
+    cpdef inline void calc_actualconstant_actualfactor(self) noexcept nogil
+    cpdef inline void calc_inputpredictions(self) noexcept nogil
+    cpdef inline void calc_outputpredictions(self) noexcept nogil
+    cpdef inline void calc_inputresiduals(self) noexcept nogil
+    cpdef inline void calc_outputresiduals(self) noexcept nogil
+    cpdef inline void calc_outputs(self) noexcept nogil
+    cpdef inline double return_mean(self, double[:] values, double[:] mask, numpy.int64_t number) noexcept nogil
+    cpdef inline void interpolate_inversedistance(self, double[:] inputs, double[:] outputs) noexcept nogil
+    cpdef inline void pass_outputs(self) noexcept nogil