richvalues 4.2.3__tar.gz → 4.2.5__tar.gz

{richvalues-4.2.3 → richvalues-4.2.5}/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: richvalues
-Version: 4.2.3
+Version: 4.2.5
 Summary: Python library for working with uncertainties and upper/lower limits
 Home-page: https://github.com/andresmegias/richvalues/
 Author: Andrés Megías Toledano

{richvalues-4.2.3 → richvalues-4.2.5}/pyproject.toml

@@ -4,7 +4,7 @@ build-backend = "setuptools.build_meta"
 
 [project]
 name = "richvalues"
-version = "4.2.3"
+version = "4.2.5"
 description = "Python library for working with uncertainties and upper/lower limits"
 license = {file="LICENSE"}
 authors = [{name="Andrés Megías Toledano"}]

{richvalues-4.2.3 → richvalues-4.2.5}/richvalues/__init__.py

@@ -13,12 +13,12 @@ modification, are permitted provided that the following conditions are
 met:
 
     (1) Redistributions of source code must retain the above copyright
-    notice, this list of conditions and the following disclaimer.
+    notice, this list of conditions and the following disclaimer. 
 
     (2) Redistributions in binary form must reproduce the above copyright
     notice, this list of conditions and the following disclaimer in
     the documentation and/or other materials provided with the
-    distribution.
+    distribution.  
     
     (3) The name of the author may not be used to endorse or promote
     products derived from this software without specific prior written
@@ -37,7 +37,7 @@ IN ANY WAY OUT OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE
 POSSIBILITY OF SUCH DAMAGE.
 """
 
-__version__ = '4.2.2'
+__version__ = '4.2.5'
 __author__ = 'Andrés Megías Toledano'
 
 import copy
@@ -303,18 +303,21 @@ def round_sf_unc(x, dx, n=None, min_exp=None, max_dec=None, extra_sf_lim=None):
                 min_exp = np.inf
                 y, dy = round_sf_unc(x, dx, n, min_exp, max_dec, extra_sf_lim)
         y = y.replace('e+', 'e').replace('e00', 'e0')
-        dy = dy.replace('e+','e').replace('e00', 'e0')
-    d = len(y.split('e')[0].split('.')[-1])
-    if d > max_dec and ')' not in dy:
+        dy = dy.replace('e+','e').replace('e00', 'e0') .replace('(', '').replace(')', '')
+    d = len(y.split('e')[0].split('.')[-1]) if '.' in y else 0
+    if max_dec == 0 or d > max_dec and ')' not in dy:
         if 'e' in y:
             y, exp = y.split('e')
             dy, _ = dy.split('e')
         else:
             exp = None
-        d = len(y.split('.')[-1])
-        d_ = len(dy.split('.')[-1])
+        d = len(y.split('.')[-1]) if '.' in y else 0
+        d_ = len(dy.split('.')[-1]) if '.' in dy else 0
         dy_ = round_sf(float(dy)*10**(d-d_), n, 0, extra_sf_lim)
-        dy = '(' + dy_.split('e')[0] + ')'
+        dy, exp_ = dy_.split('e')
+        if int(exp_) > 0:
+            dy = round_sf(float(dy)*10**int(exp_), n, np.inf, extra_sf_lim)
+        dy = '(' + dy + ')'
         if exp is not None:
             y = '{}e{}'.format(y, exp)
     return y, dy
@@ -632,11 +635,11 @@ class RichValue():
         """Constant value."""
         isconst = self.is_exact and self.domain[0] == self.domain[1]
         return isconst
-    @property
+    @property    
     def center(self):
         """Central value."""
         cent = self.main if self.is_centr else np.nan
-        return cent
+        return cent  
     @property
     def unc_eb(self):
         """Uncertainties with shape (2,1)."""
@@ -663,7 +666,7 @@ class RichValue():
         else:
             s_n = [np.nan]*2
         return s_n
-    @property
+    @property    
     def ampl(self):
         """Amplitudes."""
         m, b = self.main, self.domain
@@ -784,7 +787,7 @@ class RichValue():
                 m = np.nanmedian(distr)
         else:
             m = np.nan
-        return m
+        return m         
     
     def mean(self, num_points=int(1e4), sigmas=8.):
         """Mean of the PDF of the rich value."""
@@ -840,7 +843,7 @@ class RichValue():
             y = self.pdf(x)
             c = np.trapz(y*x, x) if central or standarized else 0.
             s = np.sqrt(np.trapz(y*(x-c)**2)) if standarized else 1.
-            moment = np.trapz((y*(x-c)**n)) / s**n
+            moment = np.trapz((y*(x-c)**n)) / s**n 
         else:
             moment = np.nan
         return moment
@@ -1094,7 +1097,7 @@ class RichValue():
                 else:
                     if is_lolim:
                         symbol = '>'
-                        y = int(np.floor(x)) if is_int else x
+                        y = int(np.floor(x)) if is_int else x               
                     elif is_uplim:
                         symbol = '<'
                         y = int(np.ceil(x)) if is_int else x
@@ -1144,10 +1147,10 @@ class RichValue():
                 else:
                     if is_lolim:
                         symbol = '>'
-                        y = int(np.floor(x)) if is_int else x
+                        y = int(np.floor(x)) if is_int else x  
                     elif is_uplim:
                         symbol = '<'
-                        y = int(np.ceil(x)) if is_int else x
+                        y = int(np.ceil(x)) if is_int else x  
                     y = round_sf(y, n, min_exp, extra_sf_lim)
                     if 'e' in y:
                         y, a = y.split('e')
@@ -1431,7 +1434,7 @@ class RichValue():
         if domain[0] > 0 and other_.sign() != -1:
             domain[0] = 0
         if domain[1] < 0 and other_.sign() != 1:
-            domain[1] = 0
+            domain[1] = 0    
         rvalue = function_with_rich_values('{}%{}', [self, other_],
                                domain=domain, is_domain_cyclic=True)
         return rvalue
@@ -1546,7 +1549,7 @@ class RichValue():
                 x = [x]
             x = np.array(x)
             y = np.zeros(len(x))
-            if self.is_exact:
+            if self.is_exact:    
                 ind = np.argmin(abs(x - main))
                 if hasattr(ind, '__iter__'):
                     ind = ind[0]
@@ -1607,7 +1610,7 @@ class RichValue():
                     if not self.is_lim:
                         distr = general_distribution(main, unc, domain, N)
                     else:
-                        x1, x2 = self.interval()
+                        x1, x2 = self.interval()    
                         distr = loguniform_distribution(x1, x2, N)
                 elif not is_finite_interv and all(np.isinf(unc)):
                     distr = loguniform_distribution(-np.inf, np.inf, N)
@@ -1857,17 +1860,17 @@ class RichArray(np.ndarray):
         return np.array([x.main for x in self.flat]).reshape(self.shape)
     @property
     def uncs(self):
-        return np.array([x.unc for x in self.flat]).reshape([*self.shape,2])
+        return np.array([x.unc for x in self.flat]).reshape([*self.shape,2]) 
     @property
     def are_lolims(self):
         return np.array([x.is_lolim for x in self.flat]).reshape(self.shape)
     @property
     def are_uplims(self):
-        return np.array([x.is_uplim for x in self.flat]).reshape(self.shape)
+        return np.array([x.is_uplim for x in self.flat]).reshape(self.shape) 
     @property
     def are_ranges(self):
         return np.array([x.is_range for x in self.flat]).reshape(self.shape)
-    @property
+    @property 
     def domains(self):
         return np.array([x.domain for x in self.flat]).reshape([*self.shape,2])
     @property
@@ -1908,7 +1911,7 @@ class RichArray(np.ndarray):
         return np.array([x.is_inf for x in self.flat]).reshape(self.shape)
     @property
     def centers(self):
-        return np.array([x.center for x in self.flat]).reshape(self.shape)
+        return np.array([x.center for x in self.flat]).reshape(self.shape) 
     @property
     def rel_uncs(self):
         return (np.array([x.rel_unc for x in self.flat])
@@ -1982,27 +1985,27 @@ class RichArray(np.ndarray):
     
     def medians(self, num_points=None):
         return np.array([x.median(num_points)
-                         for x in self.flat]).reshape(self.shape)
+                         for x in self.flat]).reshape(self.shape) 
 
     def means(self, num_points=int(1e4), sigmas=8.):
         return np.array([x.mean(num_points, sigmas)
-                         for x in self.flat]).reshape(self.shape)
+                         for x in self.flat]).reshape(self.shape) 
 
     def modes(self, num_points=int(1e4), sigmas=8.):
         return np.array([x.mode(num_points, sigmas)
-                         for x in self.flat]).reshape(self.shape)
+                         for x in self.flat]).reshape(self.shape) 
 
     def variances(self, num_points=int(1e4), sigmas=8.):
         return np.array([x.var(num_points, sigmas)
-                         for x in self.flat]).reshape(self.shape)
+                         for x in self.flat]).reshape(self.shape) 
 
     def stds(self, num_points=int(1e4), sigmas=8.):
         return np.array([x.std(num_points, sigmas)
-                         for x in self.flat]).reshape(self.shape)
+                         for x in self.flat]).reshape(self.shape) 
     
     def moments(self, n, central=True, standarized=False):
         return np.array([x.moments(n, central, standarized)
-                                   for x in self.flat]).reshape(self.shape)
+                                   for x in self.flat]).reshape(self.shape) 
     
     def set_params(self, params):
         """Set the rich value parameters of each entry of the rich array."""
@@ -2011,7 +2014,7 @@ class RichArray(np.ndarray):
                          'minimum exponent for scientific notation': 'min_exp',
                          'maximum number of decimals': 'max_dec',
                          'limit for extra significant figure': 'extra_sf_lim'}
-        attributes = ['domain'] + list(abbreviations.values())
+        attributes = ['domain'] + list(abbreviations.values()) 
         for entry in abbreviations:
             name = abbreviations[entry]
             if entry in params:
@@ -2163,7 +2166,7 @@ class RichDataFrame(pd.DataFrame):
         exec(code, {**{'self': self}, **globals()}, output)
         return output['df']
 
-    @property
+    @property    
     def mains(self): return self._property('mains')
     @property
     def uncs(self): return self._property2('uncs')
@@ -2595,7 +2598,7 @@ class ComplexRichValue():
         Parameters
         ----------
         real : rich value
-            Real part of the complex rich value.
+            Real part of the complex rich value. 
         imag : rich value
             Imaginary part of the complex rich value
         """
@@ -2676,14 +2679,14 @@ class ComplexRichValue():
         self.imag.num_sf = x
         
     @property
-    def min_exp(self): return round(np.mean([self.real.min_exp, self.imag.min_exp]))
+    def min_exp(self): return round(np.mean([self.real.min_exp, self.imag.min_exp]))  
     @min_exp.setter
     def min_exp(self, x):
         self.real.min_exp = x
         self.imag.min_exp = x
         
     @property
-    def max_dec(self): return min(self.real.max_dec, self.imag.max_dec)
+    def max_dec(self): return min(self.real.max_dec, self.imag.max_dec)  
     @max_dec.setter
     def max_dec(self, x):
         self.real.max_dec = x
@@ -2723,7 +2726,7 @@ class ComplexRichValue():
     def main(self):
         """Main value."""
         x = self.real.main + 1j*self.imag.main
-        return x
+        return x 
 
     @property
     def unc(self):
@@ -2943,7 +2946,7 @@ class ComplexRichValue():
         """Apply a function to the rich value"""
         return function_with_rich_values(function, self, **kwargs)
     
-    # Instance variable acronyms.
+    # Instance variable acronyms.  
 
     @property
     def real_part(self): return self.real
@@ -3287,7 +3290,7 @@ def rich_value(text=None, domain=None, is_int=None, pdf=None,
         domain = read_domain(text)
         if domain is not None:
             text = text.split('[')[0][:-1]
-        if not '--' in text:
+        if not '--' in text:    
             if text.startswith('+'):
                 text = text[1:]
             if 'e' in text:
@@ -3397,11 +3400,11 @@ def rich_value(text=None, domain=None, is_int=None, pdf=None,
                 if eval(dx1) == eval(dx2) == 0:
                     extra_sf_lim = 1 - 1e-8
                 else:
-                    extra_sf_lim = default_extra_sf_lim
+                    extra_sf_lim = default_extra_sf_lim 
                     base = float('{:e}'.format(eval(val)).split('e')[0])
                     if base <= default_extra_sf_lim:
                         if num_sf < default_num_sf + 1:
-                            extra_sf_lim = base - 1e-8
+                            extra_sf_lim = base - 1e-8        
             else:
                 extra_sf_lim = default_extra_sf_lim
             x = x.replace('e0','')
@@ -3409,7 +3412,7 @@ def rich_value(text=None, domain=None, is_int=None, pdf=None,
             unc = [eval(dx1), eval(dx2)]
             is_range = False
             if domain is None and np.isfinite(main) and unc[0] == 0. == unc[1]:
-                domain = [main]*2
+                domain = [main]*2          
         else:
             text = text.replace(' --','--').replace('-- ','--')
             text1, text2 = text.split('--')
@@ -3456,7 +3459,7 @@ def rich_value(text=None, domain=None, is_int=None, pdf=None,
                 imag = rich_value(text_imag, *args)
             else:
                 text = text.replace(' + ', '+').replace(' - ', '-')
-                val = complex(text)
+                val = complex(text) 
                 real = val.real
                 imag = val.imag
             rvalue = ComplexRichValue(real, imag, domain, is_int)
@@ -4173,7 +4176,7 @@ def loguniform_distribution(low=-1, high=1, size=1,
             log_x2 = _log10(abs(x2))
         if log_x1 < zero_log:
             log_x1 = zero_log
-        if log_x2 > inf_log:
+        if log_x2 > inf_log: 
             log_x2 = inf_log
         if x1 < 0:
             if x2 <= 0:
@@ -4315,7 +4318,7 @@ def add_zero_infs(interval, zero_log, inf_log):
     elif x2 > 0 and x2 > 10**inf_log:
         x2 = np.inf
     new_interval = [x1, x2]
-    return new_interval
+    return new_interval 
 def remove_zero_infs(interval, zero_log, inf_log):
     """Replace 0 and infinity for the given values in the input interval."""
     x1, x2 = interval
@@ -5427,7 +5430,7 @@ def curve_fit(x, y, function, guess, num_samples=3000,
         lim1, lim2 = 0.2, 1.2
         if disp_coef <= lim1:
             frac1 = 1.
-        elif disp_coef < lim2:
+        elif disp_coef < lim2: 
             frac1 = 1. - disp_coef / (lim2 - lim1)
         else:
             frac1 = 0.
@@ -5440,7 +5443,7 @@ def curve_fit(x, y, function, guess, num_samples=3000,
     result = {'parameters': params_fit, 'dispersion': dispersion, 'loss': loss,
               'parameters samples': samples, 'dispersion sample': dispersions,
               'loss sample': losses, 'number of fails': num_fails}
-    return result
+    return result   
 
 def point_fit(y, function, guess, num_samples=3000,
               loss=lambda a,b: (a-b)**2, lim_loss_factor=4.,
@@ -5517,7 +5520,7 @@ def point_fit(y, function, guess, num_samples=3000,
         lim1, lim2 = 0.2, 1.2
         if disp_coef <= lim1:
             frac1 = 1.
-        elif disp_coef < lim2:
+        elif disp_coef < lim2: 
             frac1 = 1. - disp_coef / (lim2 - lim1)
         else:
             frac1 = 0.
@@ -5760,4 +5763,4 @@ evaluate_distribution = evaluate_distr
 center_and_uncertainties = center_and_uncs
 is_not_a_number = is_nan = isnan
 is_infinite = is_inf = isinf
-is_finite = is_finite = isfinite
+is_finite = is_finite = isfinite

{richvalues-4.2.3 → richvalues-4.2.5}/richvalues.egg-info/PKG-INFO

@@ -1,6 +1,6 @@
 Metadata-Version: 2.1
 Name: richvalues
-Version: 4.2.3
+Version: 4.2.5
 Summary: Python library for working with uncertainties and upper/lower limits
 Home-page: https://github.com/andresmegias/richvalues/
 Author: Andrés Megías Toledano

{richvalues-4.2.3 → richvalues-4.2.5}/setup.py

@@ -5,7 +5,7 @@ with open('README.md', 'r') as file:
 
 setuptools.setup(
     name = 'richvalues',
-    version = '4.2.3',
+    version = '4.2.5',
     license = 'BSD-3-Clause',
     author = 'Andrés Megías Toledano',
     description = 'Python library for working with uncertainties and upper/lower limits',