たとえば、98/42
を計算するとき、7/3
ではなく2.3333333
を取得したいのですが、PythonまたはNumpy
?
fractions
モジュールはそれを行うことができます
>>> from fractions import Fraction
>>> Fraction(98, 42)
Fraction(7, 3)
Numpy gcdの here にはレシピがあります。分数を分割するために使用できるもの
>>> def numpy_gcd(a, b):
... a, b = np.broadcast_arrays(a, b)
... a = a.copy()
... b = b.copy()
... pos = np.nonzero(b)[0]
... while len(pos) > 0:
... b2 = b[pos]
... a[pos], b[pos] = b2, a[pos] % b2
... pos = pos[b[pos]!=0]
... return a
...
>>> numpy_gcd(np.array([98]), np.array([42]))
array([14])
>>> 98/14, 42/14
(7, 3)