Changeset 2306 in lmdz_wrf for trunk/tools

Timestamp:

Jan 31, 2019, 2:20:35 PM (7 years ago)

Author:

lfita

Message:

Adding:

• `almost_zero': Function to provide if two real values are the same up to certain power

File:

• trunk/tools/generic_tools.py (modified) (2 diffs)

trunk/tools/generic_tools.py

@@ -61 +61 @@
 # all_consecutive_combs: Function to provide all the consecutive possible 
 #   combinations from a series of indices
+# almost_zero: Function to provide if two real values are the same up to certain power
 # angle_DegMinSec: Function to transform an angle to Degrees Minutes Seconds
 # ASCIIfile_stats: Function to provide the statistics of a series of values from an ASCII file
@@ -15078 +15079 @@
     return newarray, newdims
 
+def almost_zero(val1, val2, powv):
+    """ Function to provide if two real values are the same up to certain power
+      val1, val2: values to compare
+      powv: power to use to determine zero (as decimal)
+    >>> almost_zero(1.12345, 1.12346, 3)
+    True
+    >>> almost_zero(1.12345, 1.12346, 6)
+    False
+    """
+    fname = 'almost_zero'
+
+    pot = 10.**powv
+
+    diff = val1 - val2
+    adiff = np.abs(diff)
+    idiff = int(adiff*pot)
+
+    fdiff = idiff/pot
+
+    iszero = np.abs(fdiff - diff) < 1./pot
+    print diff, idiff, fdiff, np.abs(fdiff - diff), ':', 1./pot
+
+    return iszero
+
 #quit()