Aspidites/woma/mathutils.py
# cython: language_level=3, annotation_typing=True, c_string_encoding=utf-8, boundscheck=False, wraparound=False, initializedcheck=False
from __future__ import division
from math import ceil as _ceil, floor as _floor
import bisect
import binascii
def clamp(x, lower=float("-inf"), upper=float("inf")):
"""Limit a value to a given range.
Args:
x (int or float): Number to be clamped.
lower (int or float): Minimum value for x.
upper (int or float): Maximum value for x.
The returned value is guaranteed to be between *lower* and
*upper*. Integers, floats, and other comparable types can be
mixed.
>>> clamp(1.0, 0, 5)
1.0
>>> clamp(-1.0, 0, 5)
0
>>> clamp(101.0, 0, 5)
5
>>> clamp(123, upper=5)
5
Similar to `numpy's clip`_ function.
.. _numpy's clip: http://docs.scipy.org/doc/numpy/reference/generated/numpy.clip.html
"""
if upper < lower:
raise ValueError(
"expected upper bound (%r) >= lower bound (%r)" % (upper, lower)
)
return min(max(x, lower), upper)
def ceil(x, options=None):
"""Return the ceiling of *x*. If *options* is set, return the smallest
integer or float from *options* that is greater than or equal to
*x*.
Args:
x (int or float): Number to be tested.
options (iterable): Optional iterable of arbitrary numbers
(ints or floats).
>>> VALID_CABLE_CSA = [1.5, 2.5, 4, 6, 10, 25, 35, 50]
>>> ceil(3.5, options=VALID_CABLE_CSA)
4
>>> ceil(4, options=VALID_CABLE_CSA)
4
"""
if options is None:
return _ceil(x)
options = sorted(options)
i = bisect.bisect_left(options, x)
if i == len(options):
raise ValueError("no ceil options greater than or equal to: %r" % x)
return options[i]
def floor(x, options=None):
"""Return the floor of *x*. If *options* is set, return the largest
integer or float from *options* that is less than or equal to
*x*.
Args:
x (int or float): Number to be tested.
options (iterable): Optional iterable of arbitrary numbers
(ints or floats).
>>> VALID_CABLE_CSA = [1.5, 2.5, 4, 6, 10, 25, 35, 50]
>>> floor(3.5, options=VALID_CABLE_CSA)
2.5
>>> floor(2.5, options=VALID_CABLE_CSA)
2.5
"""
if options is None:
return _floor(x)
options = sorted(options)
i = bisect.bisect_right(options, x)
if not i:
raise ValueError("no floor options less than or equal to: %r" % x)
return options[i - 1]
try:
_int_types = (int, long)
bytes = str
except NameError:
# py3 has no long
_int_types = (int,)
unicode = str