diff --git a/dist/metrolopy-0.5.1.tar.gz b/dist/metrolopy-0.5.1.tar.gz
deleted file mode 100644
index db190a3..0000000
Binary files a/dist/metrolopy-0.5.1.tar.gz and /dev/null differ
diff --git a/dist/metrolopy-0.5.1-py3-none-any.whl b/dist/metrolopy-0.5.2-py3-none-any.whl
similarity index 69%
rename from dist/metrolopy-0.5.1-py3-none-any.whl
rename to dist/metrolopy-0.5.2-py3-none-any.whl
index 4356f80..c3fbe9d 100644
Binary files a/dist/metrolopy-0.5.1-py3-none-any.whl and b/dist/metrolopy-0.5.2-py3-none-any.whl differ
diff --git a/dist/metrolopy-0.5.2.tar.gz b/dist/metrolopy-0.5.2.tar.gz
new file mode 100644
index 0000000..fa13738
Binary files /dev/null and b/dist/metrolopy-0.5.2.tar.gz differ
diff --git a/dist/metrolopy-0.5.1.zip b/dist/metrolopy-0.5.2.zip
similarity index 66%
rename from dist/metrolopy-0.5.1.zip
rename to dist/metrolopy-0.5.2.zip
index 2698c1d..c892511 100644
Binary files a/dist/metrolopy-0.5.1.zip and b/dist/metrolopy-0.5.2.zip differ
diff --git a/docs/_build/doctrees/environment.pickle b/docs/_build/doctrees/environment.pickle
index 0886660..8005389 100644
Binary files a/docs/_build/doctrees/environment.pickle and b/docs/_build/doctrees/environment.pickle differ
diff --git a/docs/_build/doctrees/hand_made_doc.doctree b/docs/_build/doctrees/hand_made_doc.doctree
index 7835bca..33acf27 100644
Binary files a/docs/_build/doctrees/hand_made_doc.doctree and b/docs/_build/doctrees/hand_made_doc.doctree differ
diff --git a/docs/_build/doctrees/index.doctree b/docs/_build/doctrees/index.doctree
index 802131a..0dac062 100644
Binary files a/docs/_build/doctrees/index.doctree and b/docs/_build/doctrees/index.doctree differ
diff --git a/docs/_build/doctrees/metrolopy.doctree b/docs/_build/doctrees/metrolopy.doctree
index 9801f32..5fb0264 100644
Binary files a/docs/_build/doctrees/metrolopy.doctree and b/docs/_build/doctrees/metrolopy.doctree differ
diff --git a/docs/_build/doctrees/metrolopy.tests.doctree b/docs/_build/doctrees/metrolopy.tests.doctree
index 3e59272..8f4f152 100644
Binary files a/docs/_build/doctrees/metrolopy.tests.doctree and b/docs/_build/doctrees/metrolopy.tests.doctree differ
diff --git a/docs/_build/html/_sources/hand_made_doc.rst.txt b/docs/_build/html/_sources/hand_made_doc.rst.txt
index fca3d2a..2ed1099 100644
--- a/docs/_build/html/_sources/hand_made_doc.rst.txt
+++ b/docs/_build/html/_sources/hand_made_doc.rst.txt
@@ -12,7 +12,7 @@ class gummy
 ===========
 
 class metrolopy.\ **gummy(x, u=0, unit=one, dof=float('inf'), k=1, p=None, 
-p_method=None, uunit=None, utype=None, name=None)**
+uunit=None, utype=None, name=None)**
 
 A gummy object represents a numerical value with an uncertainty and (or)
 a unit. They can be used in place of float values in Python expressions
@@ -97,21 +97,6 @@ gummy parameters
        0.22
        >>> g.u
        0.11224696052342387
-
-.. _p_method:
-
--  **p_method**: ({'loc', 'cp', 'gauss',
-   'ccp', 'chebyshev', None}) If the *p* parameter is
-   specified, *p_method* sets the method that is used to calculate the
-   coverage factor. If p_method is omitted or set to 'loc', then the
-   uncertainty is assumed to be represented by a normal probability
-   distribution if *dof* = float('inf') and shifted and scaled Student's
-   t distribution otherwise. If *p_method* = 'gauss' or 'cp' then the
-   Guass inequality is used, and if *p_method* = 'chebyshev' or 'ccp'
-   then the Chebyshev inequality is used. For *p* = 0.95 and *dof* =
-   float('inf'), *p_method* = 'loc' gives *k* = 2.0, while *p_method*
-   = 'gauss' gives *k* = 3.0 and *p_method* = 'chebyshev' gives *k* =
-   4.5.
    
 .. _parameter-uunit:
 
@@ -196,9 +181,20 @@ basic gummy properties
    the relation between the value of this property and the property *k*
    is defined by the *p_method* property
    
-.. _property-p-method:
+.. _p_method:
 
--  **p_method**: see the p_method_ parameter
+-  **p_method**: ({'loc', 'cp', 'gauss',
+   'ccp', 'chebyshev', None}) If the *p* parameter is
+   specified, *p_method* sets the method that is used to calculate the
+   coverage factor. If p_method is omitted or set to 'loc', then the
+   uncertainty is assumed to be represented by a normal probability
+   distribution if *dof* = float('inf') and shifted and scaled Student's
+   t distribution otherwise. If *p_method* = 'gauss' or 'cp' then the
+   Guass inequality is used, and if *p_method* = 'chebyshev' or 'ccp'
+   then the Chebyshev inequality is used. For *p* = 0.95 and *dof* =
+   float('inf'), *p_method* = 'loc' gives *k* = 2.0, while *p_method*
+   = 'gauss' gives *k* = 3.0 and *p_method* = 'chebyshev' gives *k* =
+   4.5.  This property may only be set at the class level.
 
 .. _name:
 
@@ -270,6 +266,9 @@ basic gummy properties
    can give some idea of the magnitude of the floating point errors, but is 
    not a substitute for a full numerical error analysis.  The default value
    is ``False``.
+   
+-  **cmp_k**, **cmp_p**: k or p values for comparisons, e.g. a > b is True
+   if c.x > k*c.u where c = b - a.  Set only at the class level.
 
 basic gummy methods
 -------------------
diff --git a/docs/_build/html/_sources/index.rst.txt b/docs/_build/html/_sources/index.rst.txt
index 8a2106c..6be011f 100644
--- a/docs/_build/html/_sources/index.rst.txt
+++ b/docs/_build/html/_sources/index.rst.txt
@@ -121,6 +121,9 @@ version history
 * Version 0.5.1, built 2 April 2019, fixed a major bug that generated negative 
   uncertainties in some cases and fixed some other minor bugs.  Improved support 
   for fraction.Fraction and mpmath.mpf values.
+* Version 0.5.2, built 5 April 2019, fixed a major bug that propagated
+  uncertainty incorrectly if a gummy was created with an uncertainty set with
+  an integer data type.  Fixed several other minor bugs.
 
 
 author
diff --git a/docs/_build/html/genindex.html b/docs/_build/html/genindex.html
index 343d57c..54aec6b 100644
--- a/docs/_build/html/genindex.html
+++ b/docs/_build/html/genindex.html
@@ -7,7 +7,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>Index &#8212; metrolopy 0.5.1 documentation</title>
+    <title>Index &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -15,7 +15,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -39,7 +39,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -913,12 +913,8 @@ <h2 id="P">P</h2>
         <li><a href="metrolopy.html#metrolopy.fit.SinFit.p0">(metrolopy.fit.SinFit attribute)</a>
 </li>
       </ul></li>
-      <li><a href="metrolopy.html#metrolopy.gummy.gummy.p_method">p_method (metrolopy.gummy.gummy attribute)</a>
-
-      <ul>
-        <li><a href="metrolopy.html#metrolopy.gummy.MetaGummy.p_method">(metrolopy.gummy.MetaGummy attribute)</a>
+      <li><a href="metrolopy.html#metrolopy.gummy.MetaGummy.p_method">p_method (metrolopy.gummy.MetaGummy attribute)</a>
 </li>
-      </ul></li>
       <li><a href="metrolopy.html#metrolopy.fit.DoubleExpFit.pf">pf (metrolopy.fit.DoubleExpFit attribute)</a>
 
       <ul>
@@ -1194,6 +1190,8 @@ <h2 id="T">T</h2>
       <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_abs">test_abs() (in module metrolopy.tests.test_create)</a>
 </li>
       <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_add">test_add() (in module metrolopy.tests.test_create)</a>
+</li>
+      <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_addxmul">test_addxmul() (in module metrolopy.tests.test_create)</a>
 </li>
       <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_ap1sincos">test_ap1sincos() (in module metrolopy.tests.test_create)</a>
 </li>
@@ -1208,6 +1206,10 @@ <h2 id="T">T</h2>
       <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_mod">test_mod() (in module metrolopy.tests.test_create)</a>
 </li>
       <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_mul">test_mul() (in module metrolopy.tests.test_create)</a>
+</li>
+      <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_mulxnpow">test_mulxnpow() (in module metrolopy.tests.test_create)</a>
+</li>
+      <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_mulxpow">test_mulxpow() (in module metrolopy.tests.test_create)</a>
 </li>
       <li><a href="metrolopy.tests.html#metrolopy.tests.test_create.test_neg">test_neg() (in module metrolopy.tests.test_create)</a>
 </li>
@@ -1225,14 +1227,14 @@ <h2 id="T">T</h2>
 </li>
       <li><a href="metrolopy.html#metrolopy.ummy.ummy.thousand_spaces">thousand_spaces (metrolopy.ummy.ummy attribute)</a>
 </li>
+  </ul></td>
+  <td style="width: 33%; vertical-align: top;"><ul>
       <li><a href="metrolopy.html#metrolopy.nonlinearunit.ReciprocalConversion.to">to() (metrolopy.nonlinearunit.ReciprocalConversion method)</a>
 
       <ul>
         <li><a href="metrolopy.html#metrolopy.unit.Conversion.to">(metrolopy.unit.Conversion method)</a>
 </li>
       </ul></li>
-  </ul></td>
-  <td style="width: 33%; vertical-align: top;"><ul>
       <li><a href="metrolopy.html#metrolopy.nonlinearunit.NonlinearUnit.to_uunit">to_uunit() (metrolopy.nonlinearunit.NonlinearUnit method)</a>
 </li>
       <li><a href="metrolopy.html#metrolopy.fit.Fit.toascii">toascii() (metrolopy.fit.Fit method)</a>
@@ -1619,7 +1621,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/hand_made_doc.html b/docs/_build/html/hand_made_doc.html
index 08391c1..3f840fa 100644
--- a/docs/_build/html/hand_made_doc.html
+++ b/docs/_build/html/hand_made_doc.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>MetroloPy &#8212; metrolopy 0.5.1 documentation</title>
+    <title>MetroloPy &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -38,7 +38,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -55,7 +55,7 @@ <h3>Navigation</h3>
 <div class="section" id="class-gummy">
 <span id="gummy"></span><h1>class gummy<a class="headerlink" href="#class-gummy" title="Permalink to this headline">¶</a></h1>
 <p>class metrolopy.<strong>gummy(x, u=0, unit=one, dof=float(&#8216;inf&#8217;), k=1, p=None,
-p_method=None, uunit=None, utype=None, name=None)</strong></p>
+uunit=None, utype=None, name=None)</strong></p>
 <p>A gummy object represents a numerical value with an uncertainty and (or)
 a unit. They can be used in place of float values in Python expressions
 and the uncertainty can be propagated with with both first-order and
@@ -142,20 +142,6 @@ <h2>gummy parameters<a class="headerlink" href="#gummy-parameters" title="Permal
 </div>
 </li>
 </ul>
-<ul class="simple" id="p-method">
-<li><strong>p_method</strong>: ({&#8216;loc&#8217;, &#8216;cp&#8217;, &#8216;gauss&#8217;,
-&#8216;ccp&#8217;, &#8216;chebyshev&#8217;, None}) If the <em>p</em> parameter is
-specified, <em>p_method</em> sets the method that is used to calculate the
-coverage factor. If p_method is omitted or set to &#8216;loc&#8217;, then the
-uncertainty is assumed to be represented by a normal probability
-distribution if <em>dof</em> = float(&#8216;inf&#8217;) and shifted and scaled Student&#8217;s
-t distribution otherwise. If <em>p_method</em> = &#8216;gauss&#8217; or &#8216;cp&#8217; then the
-Guass inequality is used, and if <em>p_method</em> = &#8216;chebyshev&#8217; or &#8216;ccp&#8217;
-then the Chebyshev inequality is used. For <em>p</em> = 0.95 and <em>dof</em> =
-float(&#8216;inf&#8217;), <em>p_method</em> = &#8216;loc&#8217; gives <em>k</em> = 2.0, while <em>p_method</em>
-= &#8216;gauss&#8217; gives <em>k</em> = 3.0 and <em>p_method</em> = &#8216;chebyshev&#8217; gives <em>k</em> =
-4.5.</li>
-</ul>
 <ul class="simple" id="parameter-uunit">
 <li><strong>uunit</strong>: (str or <a class="reference internal" href="#class-unit"><span class="std std-ref">Unit</span></a> ) This represents the units of
 <em>u</em>. It may be a unit with the same dimension as the <em>unit</em>
@@ -228,8 +214,19 @@ <h2>basic gummy properties<a class="headerlink" href="#basic-gummy-properties" t
 the relation between the value of this property and the property <em>k</em>
 is defined by the <em>p_method</em> property</li>
 </ul>
-<ul class="simple" id="property-p-method">
-<li><strong>p_method</strong>: see the <a class="reference internal" href="#p-method">p_method</a> parameter</li>
+<ul class="simple" id="p-method">
+<li><strong>p_method</strong>: ({&#8216;loc&#8217;, &#8216;cp&#8217;, &#8216;gauss&#8217;,
+&#8216;ccp&#8217;, &#8216;chebyshev&#8217;, None}) If the <em>p</em> parameter is
+specified, <em>p_method</em> sets the method that is used to calculate the
+coverage factor. If p_method is omitted or set to &#8216;loc&#8217;, then the
+uncertainty is assumed to be represented by a normal probability
+distribution if <em>dof</em> = float(&#8216;inf&#8217;) and shifted and scaled Student&#8217;s
+t distribution otherwise. If <em>p_method</em> = &#8216;gauss&#8217; or &#8216;cp&#8217; then the
+Guass inequality is used, and if <em>p_method</em> = &#8216;chebyshev&#8217; or &#8216;ccp&#8217;
+then the Chebyshev inequality is used. For <em>p</em> = 0.95 and <em>dof</em> =
+float(&#8216;inf&#8217;), <em>p_method</em> = &#8216;loc&#8217; gives <em>k</em> = 2.0, while <em>p_method</em>
+= &#8216;gauss&#8217; gives <em>k</em> = 3.0 and <em>p_method</em> = &#8216;chebyshev&#8217; gives <em>k</em> =
+4.5.  This property may only be set at the class level.</li>
 </ul>
 <ul id="name">
 <li><p class="first"><strong>name</strong>: An arbitrary string naming the gummy. The name is used when
@@ -302,6 +299,9 @@ <h2>basic gummy properties<a class="headerlink" href="#basic-gummy-properties" t
 not a substitute for a full numerical error analysis.  The default value
 is <code class="docutils literal"><span class="pre">False</span></code>.</p>
 </li>
+<li><p class="first"><strong>cmp_k</strong>, <strong>cmp_p</strong>: k or p values for comparisons, e.g. a &gt; b is True
+if c.x &gt; k*c.u where c = b - a.  Set only at the class level.</p>
+</li>
 </ul>
 </div>
 <div class="section" id="basic-gummy-methods">
@@ -2157,7 +2157,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/index.html b/docs/_build/html/index.html
index 0279ebb..5abdca3 100644
--- a/docs/_build/html/index.html
+++ b/docs/_build/html/index.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>MetroloPy, the docs &#8212; metrolopy 0.5.1 documentation</title>
+    <title>MetroloPy, the docs &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -38,7 +38,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="#">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="#">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -132,6 +132,9 @@ <h2>version history<a class="headerlink" href="#version-history" title="Permalin
 <li>Version 0.5.1, built 2 April 2019, fixed a major bug that generated negative
 uncertainties in some cases and fixed some other minor bugs.  Improved support
 for fraction.Fraction and mpmath.mpf values.</li>
+<li>Version 0.5.2, built 5 April 2019, fixed a major bug that propagated
+uncertainty incorrectly if a gummy was created with an uncertainty set with
+an integer data type.  Fixed several other minor bugs.</li>
 </ul>
 </div>
 <div class="section" id="author">
@@ -200,7 +203,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="#">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="#">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/metrolopy.html b/docs/_build/html/metrolopy.html
index 633d3f0..62c7d4d 100644
--- a/docs/_build/html/metrolopy.html
+++ b/docs/_build/html/metrolopy.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>metrolopy package &#8212; metrolopy 0.5.1 documentation</title>
+    <title>metrolopy package &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -38,7 +38,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -3017,7 +3017,7 @@ <h2>Submodules<a class="headerlink" href="#submodules" title="Permalink to this
 
 <dl class="class">
 <dt id="metrolopy.gummy.gummy">
-<em class="property">class </em><code class="descclassname">metrolopy.gummy.</code><code class="descname">gummy</code><span class="sig-paren">(</span><em>x</em>, <em>u=0</em>, <em>unit=</em>, <em>dof=inf</em>, <em>k=1</em>, <em>p=None</em>, <em>p_method=None</em>, <em>uunit=None</em>, <em>utype=None</em>, <em>name=None</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.gummy.gummy" title="Permalink to this definition">¶</a></dt>
+<em class="property">class </em><code class="descclassname">metrolopy.gummy.</code><code class="descname">gummy</code><span class="sig-paren">(</span><em>x</em>, <em>u=0</em>, <em>unit=</em>, <em>dof=inf</em>, <em>k=1</em>, <em>p=None</em>, <em>uunit=None</em>, <em>utype=None</em>, <em>name=None</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.gummy.gummy" title="Permalink to this definition">¶</a></dt>
 <dd><p>Bases: <a class="reference internal" href="#metrolopy.printing.PrettyPrinter" title="metrolopy.printing.PrettyPrinter"><code class="xref py py-class docutils literal"><span class="pre">metrolopy.printing.PrettyPrinter</span></code></a>, <a class="reference internal" href="#metrolopy.nummy.nummy" title="metrolopy.nummy.nummy"><code class="xref py py-class docutils literal"><span class="pre">metrolopy.nummy.nummy</span></code></a></p>
 <p>A gummy object represents a numerical value with an uncertainty and (or) a
 unit.  They can be used in place of float values in Python expressions and
@@ -3629,22 +3629,6 @@ <h2>Submodules<a class="headerlink" href="#submodules" title="Permalink to this
 </div>
 </dd></dl>
 
-<dl class="attribute">
-<dt id="metrolopy.gummy.gummy.p_method">
-<code class="descname">p_method</code><a class="headerlink" href="#metrolopy.gummy.gummy.p_method" title="Permalink to this definition">¶</a></dt>
-<dd><p><cite>str</cite> in {&#8216;loc&#8217;, &#8216;cp&#8217;, &#8216;gauss&#8217;, &#8216;ccp&#8217;, &#8216;chebyshev&#8217;}</p>
-<p>This sets the default <cite>p_method</cite> attribute for newly created gummys which
-determines how the coverage factor is calculated from a given level of
-confidence or coverage probability <cite>p</cite>.</p>
-<p>If <cite>p_method</cite> = &#8216;loc&#8217;, then the uncertainty is assumed to be represented
-by a normal probability distribution if <cite>dof</cite> = <cite>float(&#8216;inf&#8217;)</cite> and shifted
-and scaled Student&#8217;s t distribution otherwise.  If <cite>p_method</cite> = &#8216;gauss&#8217; or
-&#8216;cp&#8217; then the Gauss inequality is used, and if <cite>p_method</cite> = &#8216;chebyshev&#8217; or
-&#8216;ccp&#8217; then the Chebyshev inequality is used.  For <cite>p</cite> = 0.95 and
-<cite>dof</cite> = <cite>float(&#8216;inf&#8217;)</cite>, <cite>p_method</cite> = &#8216;loc&#8217; gives <cite>k</cite> = 2.0, while
-<cite>p_method</cite> = &#8216;gauss&#8217; gives <cite>k</cite> = 3.0 and p_method = &#8216;chebyshev&#8217; gives <cite>k</cite> = 4.5.</p>
-</dd></dl>
-
 <dl class="method">
 <dt id="metrolopy.gummy.gummy.reduce_unit">
 <code class="descname">reduce_unit</code><span class="sig-paren">(</span><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.gummy.gummy.reduce_unit" title="Permalink to this definition">¶</a></dt>
@@ -4948,7 +4932,7 @@ <h2>Submodules<a class="headerlink" href="#submodules" title="Permalink to this
 
 <dl class="attribute">
 <dt id="metrolopy.ummy.GummyTag.tags">
-<code class="descname">tags</code><em class="property"> = &lt;WeakValueDictionary at 0xfae3358&gt;</em><a class="headerlink" href="#metrolopy.ummy.GummyTag.tags" title="Permalink to this definition">¶</a></dt>
+<code class="descname">tags</code><em class="property"> = &lt;WeakValueDictionary at 0x98f3278&gt;</em><a class="headerlink" href="#metrolopy.ummy.GummyTag.tags" title="Permalink to this definition">¶</a></dt>
 <dd></dd></dl>
 
 </dd></dl>
@@ -5641,7 +5625,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/metrolopy.tests.html b/docs/_build/html/metrolopy.tests.html
index 4ef5bf7..f92f7e0 100644
--- a/docs/_build/html/metrolopy.tests.html
+++ b/docs/_build/html/metrolopy.tests.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>metrolopy.tests package &#8212; metrolopy 0.5.1 documentation</title>
+    <title>metrolopy.tests package &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -38,7 +38,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -84,6 +84,11 @@ <h2>Submodules<a class="headerlink" href="#submodules" title="Permalink to this
 <code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_add</code><span class="sig-paren">(</span><em>n=1000</em>, <em>sim=False</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_add" title="Permalink to this definition">¶</a></dt>
 <dd></dd></dl>
 
+<dl class="function">
+<dt id="metrolopy.tests.test_create.test_addxmul">
+<code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_addxmul</code><span class="sig-paren">(</span><em>n=1000</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_addxmul" title="Permalink to this definition">¶</a></dt>
+<dd></dd></dl>
+
 <dl class="function">
 <dt id="metrolopy.tests.test_create.test_ap1sincos">
 <code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_ap1sincos</code><span class="sig-paren">(</span><em>n=1000</em>, <em>sim=False</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_ap1sincos" title="Permalink to this definition">¶</a></dt>
@@ -119,6 +124,16 @@ <h2>Submodules<a class="headerlink" href="#submodules" title="Permalink to this
 <code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_mul</code><span class="sig-paren">(</span><em>n=1000</em>, <em>sim=False</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_mul" title="Permalink to this definition">¶</a></dt>
 <dd></dd></dl>
 
+<dl class="function">
+<dt id="metrolopy.tests.test_create.test_mulxnpow">
+<code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_mulxnpow</code><span class="sig-paren">(</span><em>n=1000</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_mulxnpow" title="Permalink to this definition">¶</a></dt>
+<dd></dd></dl>
+
+<dl class="function">
+<dt id="metrolopy.tests.test_create.test_mulxpow">
+<code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_mulxpow</code><span class="sig-paren">(</span><em>n=1000</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_mulxpow" title="Permalink to this definition">¶</a></dt>
+<dd></dd></dl>
+
 <dl class="function">
 <dt id="metrolopy.tests.test_create.test_neg">
 <code class="descclassname">metrolopy.tests.test_create.</code><code class="descname">test_neg</code><span class="sig-paren">(</span><em>n=1000</em>, <em>sim=False</em><span class="sig-paren">)</span><a class="headerlink" href="#metrolopy.tests.test_create.test_neg" title="Permalink to this definition">¶</a></dt>
@@ -206,7 +221,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/modules.html b/docs/_build/html/modules.html
index eb816d2..33647b2 100644
--- a/docs/_build/html/modules.html
+++ b/docs/_build/html/modules.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>metrolopy &#8212; metrolopy 0.5.1 documentation</title>
+    <title>metrolopy &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -38,7 +38,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -127,7 +127,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/objects.inv b/docs/_build/html/objects.inv
index 9c56189..c8a4083 100644
Binary files a/docs/_build/html/objects.inv and b/docs/_build/html/objects.inv differ
diff --git a/docs/_build/html/py-modindex.html b/docs/_build/html/py-modindex.html
index 270e0a0..badbb99 100644
--- a/docs/_build/html/py-modindex.html
+++ b/docs/_build/html/py-modindex.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>Python Module Index &#8212; metrolopy 0.5.1 documentation</title>
+    <title>Python Module Index &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -41,7 +41,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="#" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -218,7 +218,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="#" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/search.html b/docs/_build/html/search.html
index 6280b6a..be902c4 100644
--- a/docs/_build/html/search.html
+++ b/docs/_build/html/search.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>Search &#8212; metrolopy 0.5.1 documentation</title>
+    <title>Search &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -46,7 +46,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -97,7 +97,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/_build/html/searchindex.js b/docs/_build/html/searchindex.js
index 35bb685..a4718d3 100644
--- a/docs/_build/html/searchindex.js
+++ b/docs/_build/html/searchindex.js
@@ -1 +1 @@
-Search.setIndex({docnames:["index"],envversion:51,filenames:["index.rst"],objects:{},objnames:{},objtypes:{},terms:{"00e":[],"01e":[],"0xfae3358":[],"13e":[],"5px":[],"8ex":[],"9th":0,"\u03bc":[],"abstract":[],"case":0,"class":[],"const":[],"default":[],"float":[],"function":[],"import":0,"int":[],"new":[],"public":0,"return":[],"short":[],"static":[],"super":[],"switch":[],"throw":[],"true":[],"while":[],But:[],For:[],Not:[],One:[],The:0,Then:[],These:0,Use:[],With:[],__init__:[],__str__:[],_add:[],_appli:[],_compositeunit:[],_cov:[],_fit:[],_intervalunit:[],_nappli:[],_one:[],_radd:[],_simul:[],_sub:[],abov:0,abs:[],absolut:[],absoult:[],access:[],account:[],actual:[],add:[],add_symbol:[],added:[],addit:[],additional_nam:[],additional_short_nam:[],advoc:[],affect:[],after:[],again:[],air:[],algebera:[],algorithm:[],alia:[],alias:[],align:[],all:[],allow:[],along:[],alpha:[],alreadi:[],also:0,alter:[],altern:0,alwai:[],amount:[],analysi:[],analyz:[],angl:[],ani:[],anoth:[],any_but_latex:[],api:0,appear:[],append:[],appli:[],appropri:[],approxim:[],april:0,araound:[],arbitrari:[],arcco:[],arccosh:[],arcsin:[],arcsindist:[],arcsinh:[],arctan2:[],arctan:[],arctanh:[],arg1:[],arg2:[],arg:[],argument:[],arithmet:[],around:[],arrai:[],array_lik:[],ascii:[],ascii_symbol:[],assign:[],associ:[],assum:[],assumpt:[],asterisk:[],asterix:[],astronom:0,attach:[],attempt:[],atto:[],attribut:[],auto:[],autoconvert:[],autocorrel:[],automat:0,avail:[],averag:[],avoid:[],axes:[],axhlin:[],axi:[],axvlin:[],back:[],background:[],band:[],bar:[],base:[],base_descript:[],base_prefix:[],bayesian:[],bayesian_default:[],becaus:[],been:[],befor:[],begin:[],behavior:[],below:0,best:0,between:[],bin:[],binari:[],binary_func:[],binaryprefixedunit:[],binomi:[],binomialdist:[],bipm:[],boadcast:[],bold:[],bool:[],border:[],both:[],bottom:[],bound:[],bracket:[],break_on_printing_error:[],broadcast:[],brochur:0,budget:0,budget_first_col_header_cel:[],budget_first_col_x_cel:[],budget_header_cel:[],budget_header_row:[],budget_row:[],budget_t:[],budget_x_cel:[],budgetwarn:[],bug:0,built:0,builtin:[],calcul:[],calibr:0,call:[],can:0,canada:0,cancel:[],cannot:[],capabl:[],care:[],carlo:0,carri:[],cartesian:[],caus:[],cbrt:[],ccp:[],ccp_from_k:[],ceil:[],cell:[],celsiu:[],center:0,centi:[],certain:[],chain:[],chang:[],charact:[],chebyshev:[],chi:[],chi_correct:[],chisquar:[],chisquareddist:[],choic:[],choos:[],ci_mark:[],ci_marker_opt:[],ciformat:[],cik:[],cimethod:[],cioption:[],cip:[],circ:[],circular:[],circularunitconversionerror:[],cisim:[],cisym:[],classmethod:[],clear:[],clformat:[],clk:[],cloption:[],clp:[],cls:[],cm2:0,cmp_k:[],cmp_p:[],codata:0,code:0,coeffici:[],color:[],column:[],column_nam:[],combin:[],command:[],commonli:[],compar:[],comparison:[],compat:[],complet:[],complex:0,compon:[],composit:[],comprehens:0,comput:[],concis:[],confid:[],confus:0,conj:[],conjug:[],conservative_coverage_prob:[],consol:[],constant:0,constitu:[],conta:[],contact:0,contain:0,content:[],contribut:[],control:[],convei:0,convers:0,convert:[],convolut:[],coordin:[],copi:[],core:[],correct:[],correctli:[],correl:0,correlation_matrix:[],correlation_matrix_sim:[],correlation_sim:[],correspond:[],cos:[],cosh:[],cosin:[],council:0,count:[],cov:[],covari:[],covariance_matrix:[],covariance_matrix_sim:[],covariance_sim:[],coverag:[],coverage_factor:[],coverage_prob:[],covplot:[],covsim:[],covsim_matrix:[],cp_from_k:[],creat:[],creation:[],cross:[],css:[],cube:[],cubic:[],cumpord:[],cumprod:[],cumsum:[],curv:0,curvlinear:[],curvlineartrapdist:[],custom:[],custom_head:[],custom_math_mod:[],customari:0,data:[],data_format:[],data_opt:[],datafram:[],date:[],deal:0,deca:[],deci:[],decibel:[],decid:[],decreas:[],def:[],default_column:[],default_column_nam:[],default_css:[],defin:0,definit:0,defualt:[],deg:[],degf:[],degr:[],degre:0,delta:[],delta_diff:[],delta_diff_mean:[],delta_sum:[],delta_sum_mean:[],denomin:[],densiti:[],depend:[],deriv:[],describ:[],descript:[],description_math_mod:[],desir:[],detail:[],determin:[],devat:[],develop:0,deviat:[],df1:[],df2:[],df_html:[],df_latex:[],df_str:[],dfunc:[],dhypot:[],dict:[],dictionari:[],diff:[],differ:0,digit:[],dimens:[],dimension:[],dimensionless:[],direchletdist:[],directli:[],dirichlet:[],disabl:[],displai:[],distribut:0,divid:[],divis:[],divmod:[],document:0,doe:[],dof:[],doffrom:[],don:[],dot:[],doubl:[],doubleexpfit:[],down:[],download:0,dp1:[],dp2:[],draft:0,draw:[],drawn:[],drift:[],dure:[],dx1:[],dx2:[],each:[],ediff1d:[],edit:0,effect:[],either:[],element:[],els:[],enclos:[],encod:[],end:[],entri:[],enumer:[],epsilon:[],equal:[],equat:[],equiv:[],equival:[],eras:[],erg:[],error:[],error_bar:[],error_bar_k:[],errorbar:[],estim:[],even:[],evenli:[],event:[],exa:[],exact:[],exactli:[],exampl:[],exbi:[],except:[],exception_on_fmt_error:[],exception_on_warn:[],exclud:[],execut:[],exist:[],exit:[],exp2:[],exp:[],expan:[],expand:0,expfit:[],explicitli:[],expm1:[],expon:[],exponenti:[],exponentialdist:[],express:0,extens:[],extra:[],fab:[],facilit:[],factor:[],fahrenheit:[],fail:[],fals:[],feet:[],femto:[],few:[],fhtml:[],figur:[],fill:[],find:[],finfo:[],finit:[],fionli:[],first:0,fit:0,fit_format:[],fit_opt:[],fit_output:[],fitwarn:[],fix:0,flag:[],flase:[],flatex:[],floor:[],floor_divid:[],fmt:[],follow:[],font:[],foot:[],form:[],format:[],format_latex:[],forward:[],found:[],four:[],fraction:0,fragment:[],free:0,freedom:0,frm:[],from:[],from_unit:[],from_uunit:[],full:[],func:[],funicod:[],gamma:[],gammadist:[],gauss:[],gconst:[],gener:[],get:[],get_composit:[],get_p0:[],get_puint:[],get_punit:[],get_valu:[],gibi:[],giga:[],github:0,give:[],given:[],gnu:0,good:[],gradient:[],graft:[],greater:[],guass:[],guess:[],guid:0,gum:0,gummi:[],gummytag:[],gummywarn:[],half_width:0,handl:0,handmad:0,harold:0,has:[],have:[],head:[],header:[],heavisid:[],heavysid:[],hecto:[],held:[],help:[],here:[],hist:[],histogram:[],hold:[],holomorph:[],horizont:[],how:[],howev:0,html:[],html_symbol:[],hyperbol:[],hypot:[],iau:0,idea:[],ident:[],ignor:[],ignore_nan:[],imag:[],imaginari:[],immedi:[],implement:0,impos:0,improv:0,inch:[],includ:0,incompat:[],incompatibleunitserror:[],incorpor:[],independ:[],index:0,indic:[],inequ:[],inf:[],infinit:[],inform:[],inherit:[],init:[],initi:[],inlin:[],input:[],insensit:[],insent:[],insert:[],instal:0,instanc:[],instead:[],institut:0,integ:[],integr:[],intend:[],interact:[],intern:0,interpol:[],interpret:[],interv:[],invers:[],involv:[],ipython:0,is_dimensionless:[],isindepend:[],iso:0,issu:0,iter:[],its:[],itself:[],jac:[],jacobian:[],jgcm:0,joint:[],jummi:[],juptyer:[],jupyt:0,just:[],keep:[],kei:[],keyword:[],keywork:[],kibi:[],kilo:[],kilogram:[],known:0,kpa:0,ksim:[],kwd:[],lab:0,label:[],lam:[],laplac:[],laplacedist:[],largest:[],last:[],later:[],latex:[],latex_math:[],latex_math_plot:[],latex_norm:[],latex_norm_plot:[],latex_symbol:[],least:[],leastsq:[],left:[],left_width:[],len:[],length:[],less:[],level:[],librari:[],library_nam:[],lie:[],lies:[],light:[],like:[],limit:[],limit_half_rang:[],line:[],linear:[],linewidth:[],list:0,load:[],loc:[],loc_from_k:[],locat:[],log10:[],log1p:[],log2:[],log:[],log_bas:[],log_func:[],logaddexp2:[],logaddexp:[],logarithm:[],logconvers:[],loge:[],lognormaldist:[],logrithm:[],logunit:[],look:0,lot:[],lower:[],lower_limit:[],lowest:[],machin:[],machineri:[],made:[],magnitud:[],mai:[],major:0,make:[],make_gummi:[],make_numb:[],mani:0,march:0,markdown:[],marker:[],mask:[],match:[],math:[],mathemat:[],mathrm:[],matrix:[],max_digit:[],max_dof:[],maximum:[],maxit:[],mcisim:[],mcisym:[],mcpropag:[],mean:[],mean_datetim:[],mean_mark:[],mean_marker_opt:[],measur:0,measurand:[],mebi:[],mega:[],metagummi:[],metanummi:[],metaprettyprint:[],meter:[],method:[],metr:[],metrolog:0,metrologia:[],mfraction:[],micro:[],microinch:[],microsecond:[],milli:[],millisecond:[],min:[],minim:[],minor:0,minu:[],miscellan:[],miss:[],mix:[],mod:[],mode:[],modf:[],modifi:0,modul:0,modulu:[],mone:[],mont:0,more:0,most:[],mpf:0,mpmath:0,much:0,mulit:[],mulsep:[],multi:[],multinormaldist:[],multipl:[],multipli:[],multitdist:[],multitel:[],multivari:[],multivariatedistribut:[],multivariateel:[],multvariatedistribut:[],must:0,n_eff:[],n_sigma:[],n_sigma_trim:[],name:[],nan:[],nano:[],nappli:[],nation:0,natur:[],ndarrai:[],nearest:[],necessari:[],need:[],neg:0,negat:[],neither:[],neper:[],never:[],newli:[],nist:0,nls:[],non:[],none:[],nonlinear:[],nonlinearconvers:[],nonlinearunit:[],nor:[],norm:[],normal:[],normaldist:[],nosimulateddataerror:[],notat:[],note:0,notebook:0,notimplementederror:[],nounitconversionfounderror:[],now:[],nparam:[],nprop:[],nrc:0,nsig:[],number:[],numer:[],nummi:[],numpi:0,object:0,observ:[],odd:[],odr:[],off:[],offset:[],offsetconvers:[],offsetunit:[],often:[],older:[],ols:[],omit:[],ommit:[],onc:[],one:[],oneovertfit:[],onli:[],oper:[],operand:[],opertaion:[],optim:[],option:[],order:[],origin:[],other:0,otherwis:[],out:[],outlier:[],outlin:0,output:[],outsid:[],over:[],over_plot:[],overrid:[],overridden:[],own:[],p_method:[],packag:0,pad:[],page:0,pair:[],panda:0,param:[],paramet:[],paramt:[],parent:[],park:0,parksh:0,pars:[],part:[],partial:[],particular:[],pass:[],pebi:[],perform:[],perhap:0,peta:[],phi:[],physic:0,pick:[],pico:[],pip:0,place:[],pleas:0,plot:[],plot_opt:[],plot_point:[],plu:[],plusmn:0,pmethod:[],pmi:[],pmsim:[],pmsimi:[],point:[],poisson:[],poissondist:[],polar:[],polyfit:[],polynomi:[],popul:[],posit:[],possibl:[],pow:[],power:[],ppm:[],practic:0,predict:[],prefer:[],prefix:[],prefix_definit:[],prefixedunit:[],pressur:[],pretti:[],prettyprint:[],previou:[],previous:[],primarili:[],print:[],print_html:[],print_latex:[],print_markdown:[],printer:[],prnt:[],probabl:[],problem:[],procedur:[],prod:[],propag:0,properti:[],proport:[],provid:[],ptostr:[],publish:[],punit:[],pure:0,put:[],pyplot:[],python:0,quadrant:[],quaitii:[],quanit:[],quantiti:0,rais:[],random:[],random_st:[],randomst:[],rare:[],rate:[],rather:[],ratiounit:[],raw:[],rchain:[],read:[],readi:[],real:[],reciproc:[],reciprocalconvers:[],recogn:[],recommend:0,red:[],reduc:[],reduce_unit:[],refer:0,referenc:0,rel:[],relat:[],relationship:[],releas:0,relunit:[],remaind:[],remov:[],reorder:[],repres:0,represent:[],requir:0,res:[],research:0,resect:[],resid:[],residu:[],respect:[],respons:[],restrict:0,result:[],retriev:[],right:[],right_width:[],rint:[],root:[],round:[],rounding_u:[],row:[],rss:[],rule:[],run:[],s276:[],same:0,sampl:[],satterthwait:[],scalar:[],scale:[],scatter:[],sci_not:[],sci_notation_high:[],sci_notation_low:[],scientif:[],scipi:0,search:0,search_unit:[],second:[],section:0,see:0,seed:[],select:[],self:[],sensit:[],separ:[],seri:[],serv:[],session:[],set:0,set_print:[],set_se:[],set_tag:[],sever:[],shadow:[],shadow_unit:[],shadowed_alias:[],shadowed_unit:[],shape:[],share:[],shift:[],short_nam:[],shortest:[],should:[],show:[],show_:[],show_al:[],show_c:[],show_d:[],show_data:[],show_dof:[],show_expanded_u:[],show_fit:[],show_k:[],show_nam:[],show_p:[],show_subtot:[],shown:[],side:[],sigma:[],sigma_is_known:[],sigma_trim:[],sign:[],signal:[],signatur:[],signific:[],sim:[],simdata:[],similar:[],simpli:[],simsort:[],simul:[],simultan:[],sin:[],sinc:[],sine:[],sinfit:[],singl:[],sinh:[],siunit:[],size:[],slash:[],slashaxi:[],slightli:0,slow:[],smallest:[],softwar:0,solid:[],solidu:[],solver:[],some:0,someth:[],sometim:[],sometimes_smal:[],sort:[],sound:[],sourc:0,space:[],special:0,specifi:[],spl:[],spline:[],sqrt:[],squar:[],stabl:[],stack:0,standard:0,start:[],statement:[],staticmethod:[],stdev:[],str:[],straight:[],string:0,strip:[],student:[],style:[],sub:[],subclass:[],submodul:[],subpackag:[],substitut:[],subtot:[],subtract:[],subtrat:[],success:[],sum:[],suppli:[],support:0,symbol:[],symmetr:[],system:0,tab:[],tabl:[],tack:[],tag:[],tag_nam:[],take:[],taken:[],tan:[],tangent:[],tanh:[],tdist:[],tebi:[],teh:[],temperatur:[],tera:[],term:0,test:[],test_ab:[],test_add:[],test_ap1sinco:[],test_ap2sinco:[],test_apnsinco:[],test_creat:[],test_div:[],test_gummy_init:[],test_mod:[],test_mul:[],test_neg:[],test_npsinco:[],test_po:[],test_pow:[],test_sinco:[],test_sub:[],test_ubreakdown:[],text:[],than:0,thei:0,them:[],themont:[],thi:0,thick:[],thin:[],those:[],though:[],thousand_spac:[],ths:[],thu:[],time:[],titl:[],title_styl:[],to_unit:[],to_uunit:[],toascii:[],tofloat:[],tohtml:[],tolatex:[],tool:0,top_to_base_ratio:[],tostr:[],total:[],tounicod:[],toward:[],track:[],tradit:[],tradition:[],transcendent:[],trapezoid:[],trapezoidaldist:[],treat:[],tri:[],trial:[],triangular:[],triangulardist:[],true_divid:[],trunc:[],try_fconvert:[],tupl:[],turn:[],tutori:0,two:[],txt:[],type:[],u00b0f:[],ubreakdown:[],ueq:[],uexp:[],ufrom:[],ufsim:[],ummi:[],uncertainti:0,uncertianti:[],uncertiantyprecisionwarn:[],uncorrel:[],under:0,unicod:[],uniform:[],uniformdist:0,unit:0,uniterror:[],unitliberror:[],unitlibnotfounderror:[],unitlibrari:[],unitnotfounderror:[],units_on_valu:[],unitutil:[],unitwarn:[],unless:[],unpredict:[],unshadow:[],until:[],upa:[],upon:[],upper:[],upper_limit:[],use:[],use_n_eff:[],used:0,useful:[],user:[],uses:[],usim:[],using:0,usual:[],usunit:[],usym:[],utyp:[],uunit:0,uunit_is_rel:[],valid:[],valu:0,valueerror:[],vari:[],variabl:[],varianc:[],variou:[],verbatim:0,veri:[],verison:[],version:[],versu:[],vertic:[],view:[],wai:[],want:0,warn:[],weakvaluedictionari:[],weibul:[],weibulldist:[],weight:[],weird:[],welch:[],well:[],were:0,when:[],whenev:[],where:[],whether:[],which:0,who:[],width:[],willink:[],wish:0,without:[],wmean:[],word:[],work:[],would:[],write:[],xdim:[],xfsim:[],xlabel:[],xlist:[],xmax:[],xmin:[],xname:[],xsim:[],xsym:[],xunit:[],yard:[],ydim:[],year:[],ylabel:[],yname:[],yobi:[],yocto:[],yotta:[],you:0,ypred:[],ypredf:[],yunit:[],zepto:[],zero:[],zetta:[],zhang:[],zibi:[]},titles:["MetroloPy, the docs"],titleterms:{"abstract":[],"class":0,"function":[],"static":[],appli:[],arithmet:[],attribut:[],author:0,auto:0,basic:[],budget:[],built:[],carlo:[],come:0,common:[],content:[],convers:[],curv:[],custom:[],dfunc:[],did:0,displai:[],distribut:[],doc:0,except:[],fit:[],format:[],from:0,further:0,gener:0,get:0,gummi:0,histori:0,indic:0,involv:[],jummi:[],licens:0,logunit:[],mathemat:[],mean:[],method:[],metrolopi:0,modul:[],mont:[],multi:[],name:0,nonlinearunit:[],numer:[],nummi:[],offsetunit:[],oper:[],packag:[],paramet:[],pmethod:[],prefixedunit:[],print:[],properti:[],read:0,relat:[],relunit:[],search:[],simul:[],siunit:[],some:[],start:0,sub:[],submodul:[],subpackag:[],tabl:0,test:[],test_creat:[],ummi:[],unit:[],unitutil:[],usunit:[],variat:[],version:0,where:0}})
\ No newline at end of file
+Search.setIndex({docnames:["hand_made_doc","index","metrolopy","metrolopy.tests","modules","todo"],envversion:51,filenames:["hand_made_doc.rst","index.rst","metrolopy.rst","metrolopy.tests.rst","modules.rst","todo.rst"],objects:{"":{metrolopy:[2,0,0,"-"]},"metrolopy.budget":{Budget:[2,1,1,""]},"metrolopy.budget.Budget":{column_names:[2,2,1,""],columns:[2,2,1,""],default_column_names:[2,2,1,""],default_columns:[2,2,1,""],default_css:[2,2,1,""],df:[2,2,1,""],df_html:[2,2,1,""],df_latex:[2,2,1,""],df_str:[2,2,1,""],k:[2,2,1,""],p:[2,2,1,""],show_c:[2,2,1,""],show_d:[2,2,1,""],show_s:[2,2,1,""],tostring:[2,3,1,""],units_on_values:[2,2,1,""]},"metrolopy.dfunc":{Dfunc:[2,1,1,""]},"metrolopy.dfunc.Dfunc":{apply:[2,4,1,""],napply:[2,4,1,""],tofloat:[2,3,1,""]},"metrolopy.distributions":{ArcSinDist:[2,1,1,""],BinomialDist:[2,1,1,""],Convolution:[2,1,1,""],CurvlinearTrapDist:[2,1,1,""],Distribution:[2,1,1,""],ExponentialDist:[2,1,1,""],GammaDist:[2,1,1,""],LaplaceDist:[2,1,1,""],LogNormalDist:[2,1,1,""],MultiNormalDist:[2,1,1,""],MultiTDist:[2,1,1,""],MultiTElement:[2,1,1,""],MultivariateDistribution:[2,1,1,""],MultivariateElement:[2,1,1,""],NormalDist:[2,1,1,""],PoissonDist:[2,1,1,""],TDist:[2,1,1,""],TrapezoidalDist:[2,1,1,""],TriangularDist:[2,1,1,""],UniformDist:[2,1,1,""],WeibullDist:[2,1,1,""]},"metrolopy.distributions.ArcSinDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.BinomialDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.Convolution":{isindependent:[2,2,1,""],random:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.CurvlinearTrapDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.Distribution":{apply:[2,5,1,""],ci:[2,3,1,""],cisym:[2,3,1,""],clear:[2,3,1,""],covplot:[2,5,1,""],covsim:[2,3,1,""],covsim_matrix:[2,5,1,""],hist:[2,3,1,""],isindependent:[2,2,1,""],mean:[2,2,1,""],random:[2,3,1,""],random_state:[2,5,1,""],set_seed:[2,5,1,""],simdata:[2,2,1,""],simsorted:[2,2,1,""],simulate:[2,5,1,""],stdev:[2,2,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.ExponentialDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.GammaDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.LaplaceDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.LogNormalDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.MultiNormalDist":{cov:[2,2,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.MultiTDist":{cov:[2,2,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.MultivariateDistribution":{clear:[2,3,1,""],cov:[2,2,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.MultivariateElement":{clear:[2,3,1,""],cov:[2,3,1,""],random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.NormalDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.PoissonDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.TDist":{bayesian_default:[2,2,1,""],random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.TrapezoidalDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.TriangularDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.UniformDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.distributions.WeibullDist":{random:[2,3,1,""],u:[2,3,1,""],x:[2,3,1,""]},"metrolopy.exceptions":{BudgetWarning:[2,6,1,""],CircularUnitConversionError:[2,6,1,""],FitWarning:[2,6,1,""],GummyWarning:[2,6,1,""],IncompatibleUnitsError:[2,6,1,""],NoSimulatedDataError:[2,6,1,""],NoUnitConversionFoundError:[2,6,1,""],UncertiantyPrecisionWarning:[2,6,1,""],UnitError:[2,6,1,""],UnitLibError:[2,6,1,""],UnitLibNotFoundError:[2,6,1,""],UnitNotFoundError:[2,6,1,""],UnitWarning:[2,6,1,""]},"metrolopy.fit":{DoubleExpFit:[2,1,1,""],ExpFit:[2,1,1,""],Fit:[2,1,1,""],OneOverTFit:[2,1,1,""],PolyFit:[2,1,1,""],SinFit:[2,1,1,""]},"metrolopy.fit.DoubleExpFit":{count:[2,2,1,""],cov:[2,2,1,""],f:[2,3,1,""],fhtml:[2,3,1,""],fit_output:[2,2,1,""],flatex:[2,3,1,""],funicode:[2,3,1,""],get_p0:[2,3,1,""],get_punits:[2,3,1,""],jac:[2,3,1,""],nparam:[2,2,1,""],p0:[2,2,1,""],p:[2,2,1,""],pf:[2,2,1,""],punits:[2,2,1,""],res:[2,2,1,""],s:[2,2,1,""],solver:[2,2,1,""],ux:[2,2,1,""],uy:[2,2,1,""],x:[2,2,1,""],xdim:[2,2,1,""],xf:[2,2,1,""],y:[2,2,1,""],ydim:[2,2,1,""],yf:[2,2,1,""]},"metrolopy.fit.ExpFit":{count:[2,2,1,""],cov:[2,2,1,""],f:[2,3,1,""],fhtml:[2,3,1,""],fit_output:[2,2,1,""],flatex:[2,3,1,""],funicode:[2,3,1,""],get_p0:[2,3,1,""],get_punits:[2,3,1,""],jac:[2,3,1,""],nparam:[2,2,1,""],p0:[2,2,1,""],p:[2,2,1,""],pf:[2,2,1,""],punits:[2,2,1,""],res:[2,2,1,""],s:[2,2,1,""],solver:[2,2,1,""],ux:[2,2,1,""],uy:[2,2,1,""],x:[2,2,1,""],xdim:[2,2,1,""],xf:[2,2,1,""],y:[2,2,1,""],ydim:[2,2,1,""],yf:[2,2,1,""]},"metrolopy.fit.Fit":{f:[2,3,1,""],fhtml:[2,3,1,""],flatex:[2,3,1,""],funicode:[2,3,1,""],get_p0:[2,3,1,""],get_punits:[2,3,1,""],html:[2,3,1,""],jac:[2,3,1,""],latex:[2,3,1,""],latex_math:[2,2,1,""],ptostring:[2,3,1,""],toascii:[2,3,1,""],tohtml:[2,3,1,""],tolatex:[2,3,1,""],tostring:[2,3,1,""],ypred:[2,3,1,""],ypredf:[2,3,1,""]},"metrolopy.fit.OneOverTFit":{count:[2,2,1,""],cov:[2,2,1,""],f:[2,3,1,""],fhtml:[2,3,1,""],fit_output:[2,2,1,""],flatex:[2,3,1,""],funicode:[2,3,1,""],get_p0:[2,3,1,""],get_punits:[2,3,1,""],jac:[2,3,1,""],nparam:[2,2,1,""],p0:[2,2,1,""],p:[2,2,1,""],pf:[2,2,1,""],punits:[2,2,1,""],res:[2,2,1,""],s:[2,2,1,""],solver:[2,2,1,""],ux:[2,2,1,""],uy:[2,2,1,""],x:[2,2,1,""],xdim:[2,2,1,""],xf:[2,2,1,""],y:[2,2,1,""],ydim:[2,2,1,""],yf:[2,2,1,""]},"metrolopy.fit.PolyFit":{f:[2,3,1,""],fhtml:[2,3,1,""],flatex:[2,3,1,""],funicode:[2,3,1,""],get_p0:[2,3,1,""],get_punits:[2,3,1,""],jac:[2,3,1,""]},"metrolopy.fit.SinFit":{count:[2,2,1,""],cov:[2,2,1,""],f:[2,3,1,""],fhtml:[2,3,1,""],fit_output:[2,2,1,""],flatex:[2,3,1,""],funicode:[2,3,1,""],get_p0:[2,3,1,""],get_punits:[2,3,1,""],jac:[2,3,1,""],nparam:[2,2,1,""],p0:[2,2,1,""],p:[2,2,1,""],pf:[2,2,1,""],punits:[2,2,1,""],res:[2,2,1,""],s:[2,2,1,""],solver:[2,2,1,""],ux:[2,2,1,""],uy:[2,2,1,""],x:[2,2,1,""],xdim:[2,2,1,""],xf:[2,2,1,""],y:[2,2,1,""],ydim:[2,2,1,""],yf:[2,2,1,""]},"metrolopy.functions":{absolute:[2,7,1,""],add:[2,7,1,""],angle:[2,7,1,""],arccos:[2,7,1,""],arccosh:[2,7,1,""],arcsin:[2,7,1,""],arcsinh:[2,7,1,""],arctan2:[2,7,1,""],arctan:[2,7,1,""],arctanh:[2,7,1,""],around:[2,7,1,""],cbrt:[2,7,1,""],ceil:[2,7,1,""],conj:[2,7,1,""],cos:[2,7,1,""],cosh:[2,7,1,""],cross:[2,7,1,""],cumprod:[2,7,1,""],cumsum:[2,7,1,""],diff:[2,7,1,""],divide:[2,7,1,""],divmod:[2,7,1,""],ediff1d:[2,7,1,""],exp2:[2,7,1,""],exp:[2,7,1,""],expm1:[2,7,1,""],fix:[2,7,1,""],floor:[2,7,1,""],floor_divide:[2,7,1,""],gradient:[2,7,1,""],heaviside:[2,7,1,""],imag:[2,7,1,""],log10:[2,7,1,""],log1p:[2,7,1,""],log2:[2,7,1,""],log:[2,7,1,""],logaddexp2:[2,7,1,""],logaddexp:[2,7,1,""],mod:[2,7,1,""],modf:[2,7,1,""],multiply:[2,7,1,""],negative:[2,7,1,""],power:[2,7,1,""],prod:[2,7,1,""],real:[2,7,1,""],reciprocal:[2,7,1,""],remainder:[2,7,1,""],rint:[2,7,1,""],sign:[2,7,1,""],sin:[2,7,1,""],sinh:[2,7,1,""],sqrt:[2,7,1,""],square:[2,7,1,""],subtract:[2,7,1,""],sum:[2,7,1,""],tan:[2,7,1,""],tanh:[2,7,1,""],true_divide:[2,7,1,""],trunc:[2,7,1,""]},"metrolopy.gummy":{MetaGummy:[2,1,1,""],gummy:[2,1,1,""],jummy:[2,1,1,""]},"metrolopy.gummy.MetaGummy":{cmp_k:[2,2,1,""],cmp_p:[2,2,1,""],p_method:[2,2,1,""],style:[2,2,1,""]},"metrolopy.gummy.gummy":{U:[2,2,1,""],Usim:[2,2,1,""],ascii:[2,3,1,""],autoconvert:[2,2,1,""],budget:[2,3,1,""],c:[2,2,1,""],cisim:[2,2,1,""],convert:[2,3,1,""],covplot:[2,4,1,""],create:[2,4,1,""],exception_on_fmt_error:[2,2,1,""],graft:[2,3,1,""],hist:[2,3,1,""],html:[2,3,1,""],imag:[2,2,1,""],k:[2,2,1,""],latex:[2,3,1,""],mulsep:[2,2,1,""],p:[2,2,1,""],reduce_unit:[2,3,1,""],show_dof:[2,2,1,""],show_k:[2,2,1,""],show_name:[2,2,1,""],show_p:[2,2,1,""],sim:[2,3,1,""],simulate:[2,5,1,""],slashaxis:[2,2,1,""],solidus:[2,2,1,""],style:[2,2,1,""],toascii:[2,3,1,""],tohtml:[2,3,1,""],tolatex:[2,3,1,""],tostring:[2,3,1,""],u:[2,2,1,""],ubreakdown:[2,2,1,""],unicode:[2,3,1,""],unit:[2,2,1,""],usim:[2,2,1,""],uunit:[2,2,1,""],uunit_is_rel:[2,2,1,""],x:[2,2,1,""],xsim:[2,2,1,""]},"metrolopy.gummy.jummy":{angle:[2,3,1,""],conjugate:[2,3,1,""],copy:[2,3,1,""],cov:[2,2,1,""],imag:[2,2,1,""],real:[2,2,1,""],tostring:[2,3,1,""],unit:[2,2,1,""],x:[2,2,1,""]},"metrolopy.logunit":{LogConversion:[2,1,1,""],LogUnit:[2,1,1,""]},"metrolopy.logunit.LogConversion":{copy:[2,3,1,""]},"metrolopy.logunit.LogUnit":{get_composite:[2,3,1,""],zero:[2,3,1,""]},"metrolopy.mean":{autocorrelation:[2,7,1,""],delta_diff:[2,7,1,""],delta_diff_mean:[2,7,1,""],delta_sum:[2,7,1,""],delta_sum_mean:[2,7,1,""],mean:[2,7,1,""],mean_datetime:[2,7,1,""],n_eff:[2,7,1,""],sigma_trim:[2,7,1,""],wmean:[2,7,1,""]},"metrolopy.nonlinearunit":{NonlinearConversion:[2,1,1,""],NonlinearUnit:[2,1,1,""],ReciprocalConversion:[2,1,1,""]},"metrolopy.nonlinearunit.NonlinearConversion":{chain:[2,3,1,""],copy:[2,3,1,""],linear:[2,2,1,""],pow:[2,3,1,""],rchain:[2,3,1,""]},"metrolopy.nonlinearunit.NonlinearUnit":{from_uunit:[2,3,1,""],get_composite:[2,3,1,""],linear:[2,2,1,""],to_uunit:[2,3,1,""],zero:[2,3,1,""]},"metrolopy.nonlinearunit.ReciprocalConversion":{chain:[2,3,1,""],copy:[2,3,1,""],frm:[2,3,1,""],to:[2,3,1,""]},"metrolopy.nummy":{MetaNummy:[2,1,1,""],nummy:[2,1,1,""]},"metrolopy.nummy.MetaNummy":{bayesian:[2,2,1,""],cimethod:[2,2,1,""],mcpropagate:[2,2,1,""]},"metrolopy.nummy.nummy":{Usim:[2,2,1,""],bayesian:[2,2,1,""],cimethod:[2,2,1,""],cisim:[2,2,1,""],correlation_matrix_sim:[2,5,1,""],correlation_sim:[2,3,1,""],covariance_matrix_sim:[2,5,1,""],covariance_sim:[2,3,1,""],covplot:[2,5,1,""],create:[2,4,1,""],distribution:[2,2,1,""],dof:[2,2,1,""],hist:[2,3,1,""],independent:[2,2,1,""],ksim:[2,2,1,""],set_seed:[2,5,1,""],simdata:[2,2,1,""],simsorted:[2,2,1,""],simulate:[2,5,1,""],ufrom:[2,3,1,""],usim:[2,2,1,""],xsim:[2,2,1,""]},"metrolopy.offsetunit":{OffsetConversion:[2,1,1,""],OffsetUnit:[2,1,1,""]},"metrolopy.offsetunit.OffsetConversion":{copy:[2,3,1,""]},"metrolopy.offsetunit.OffsetUnit":{zero:[2,3,1,""]},"metrolopy.pmethod":{ccp_from_k:[2,7,1,""],conservative_coverage_probability:[2,7,1,""],coverage_factor:[2,7,1,""],coverage_probability:[2,7,1,""],cp_from_k:[2,7,1,""],loc_from_k:[2,7,1,""]},"metrolopy.prefixedunit":{BinaryPrefixedUnit:[2,1,1,""],PrefixedUnit:[2,1,1,""]},"metrolopy.prefixedunit.BinaryPrefixedUnit":{prefix_definitions:[2,2,1,""]},"metrolopy.prefixedunit.PrefixedUnit":{prefix_definitions:[2,2,1,""]},"metrolopy.printing":{MetaPrettyPrinter:[2,1,1,""],PrettyPrinter:[2,1,1,""],print_html:[2,7,1,""],print_latex:[2,7,1,""],print_markdown:[2,7,1,""],set_printer:[2,7,1,""]},"metrolopy.printing.MetaPrettyPrinter":{printer:[2,2,1,""]},"metrolopy.printing.PrettyPrinter":{ascii:[2,3,1,""],break_on_printing_error:[2,2,1,""],html:[2,3,1,""],latex:[2,3,1,""],latex_math:[2,3,1,""],latex_math_plot:[2,3,1,""],latex_norm:[2,3,1,""],latex_norm_plot:[2,3,1,""],printer:[2,2,1,""],toascii:[2,3,1,""],tohtml:[2,3,1,""],tolatex:[2,3,1,""],tounicode:[2,3,1,""],unicode:[2,3,1,""]},"metrolopy.relunits":{RatioUnit:[2,1,1,""]},"metrolopy.tests":{test_create:[3,0,0,"-"]},"metrolopy.tests.test_create":{binary_func:[3,7,1,""],display:[3,7,1,""],make_gummy:[3,7,1,""],make_number:[3,7,1,""],test_abs:[3,7,1,""],test_add:[3,7,1,""],test_addxmul:[3,7,1,""],test_ap1sincos:[3,7,1,""],test_ap2sincos:[3,7,1,""],test_apnsincos:[3,7,1,""],test_div:[3,7,1,""],test_gummy_init:[3,7,1,""],test_mod:[3,7,1,""],test_mul:[3,7,1,""],test_mulxnpow:[3,7,1,""],test_mulxpow:[3,7,1,""],test_neg:[3,7,1,""],test_npsincos:[3,7,1,""],test_pos:[3,7,1,""],test_pow:[3,7,1,""],test_sincos:[3,7,1,""],test_sub:[3,7,1,""],test_ubreakdown:[3,7,1,""]},"metrolopy.ummy":{GummyTag:[2,1,1,""],MFraction:[2,1,1,""],ummy:[2,1,1,""]},"metrolopy.ummy.GummyTag":{get_values:[2,3,1,""],set_tag:[2,5,1,""],tags:[2,2,1,""]},"metrolopy.ummy.ummy":{"const":[2,2,1,""],angle:[2,3,1,""],conjugate:[2,3,1,""],copy:[2,3,1,""],correlation:[2,3,1,""],correlation_matrix:[2,5,1,""],covariance:[2,3,1,""],covariance_matrix:[2,5,1,""],create:[2,4,1,""],dof:[2,2,1,""],doffrom:[2,3,1,""],finfo:[2,2,1,""],imag:[2,2,1,""],max_digits:[2,2,1,""],max_dof:[2,2,1,""],nsig:[2,2,1,""],real:[2,2,1,""],rounding_u:[2,2,1,""],sci_notation:[2,2,1,""],sci_notation_high:[2,2,1,""],sci_notation_low:[2,2,1,""],thousand_spaces:[2,2,1,""],tofloat:[2,3,1,""],u:[2,2,1,""],ufrom:[2,3,1,""],utype:[2,2,1,""],x:[2,2,1,""]},"metrolopy.unit":{Conversion:[2,1,1,""],Unit:[2,1,1,""]},"metrolopy.unit.Conversion":{chain:[2,3,1,""],copy:[2,3,1,""],frm:[2,3,1,""],linear:[2,2,1,""],pow:[2,3,1,""],to:[2,3,1,""],unit:[2,2,1,""]},"metrolopy.unit.Unit":{alias:[2,5,1,""],aliases:[2,2,1,""],conversion:[2,2,1,""],convert:[2,3,1,""],format_latex:[2,5,1,""],is_dimensionless:[2,2,1,""],linear:[2,2,1,""],load:[2,5,1,""],reorder:[2,5,1,""],shadowed_aliases:[2,2,1,""],tostring:[2,3,1,""],unit:[2,5,1,""],units:[2,2,1,""]},"metrolopy.unitutils":{convert:[2,7,1,""],search_units:[2,7,1,""],shadowed_units:[2,7,1,""]},metrolopy:{budget:[2,0,0,"-"],dfunc:[2,0,0,"-"],distributions:[2,0,0,"-"],exceptions:[2,0,0,"-"],fit:[2,0,0,"-"],functions:[2,0,0,"-"],gummy:[2,0,0,"-"],logunit:[2,0,0,"-"],mean:[2,0,0,"-"],nonlinearunit:[2,0,0,"-"],nummy:[2,0,0,"-"],offsetunit:[2,0,0,"-"],pmethod:[2,0,0,"-"],prefixedunit:[2,0,0,"-"],printing:[2,0,0,"-"],relunits:[2,0,0,"-"],siunits:[2,0,0,"-"],tests:[3,0,0,"-"],ummy:[2,0,0,"-"],unit:[2,0,0,"-"],unitutils:[2,0,0,"-"],usunits:[2,0,0,"-"],version:[2,0,0,"-"]}},objnames:{"0":["py","module","Python module"],"1":["py","class","Python class"],"2":["py","attribute","Python attribute"],"3":["py","method","Python method"],"4":["py","classmethod","Python class method"],"5":["py","staticmethod","Python static method"],"6":["py","exception","Python exception"],"7":["py","function","Python function"]},objtypes:{"0":"py:module","1":"py:class","2":"py:attribute","3":"py:method","4":"py:classmethod","5":"py:staticmethod","6":"py:exception","7":"py:function"},terms:{"00e":[0,2],"01e":[0,2],"0x98f3278":2,"13e":[0,2],"5px":2,"8ex":2,"9th":[1,2],"\u03bc":2,"abstract":2,"case":[0,1,2],"class":2,"const":[0,2],"default":[0,2],"float":[0,2],"function":4,"import":[0,1,2],"int":[0,2],"new":[0,2],"public":[1,2],"return":[0,2],"short":[0,2],"static":2,"super":[0,2],"switch":2,"throw":[0,2],"true":[0,2,3],"while":[0,2],But:[0,2],For:[0,2],Not:2,One:[0,2],The:[0,1,2],Then:0,These:[0,1,2],Use:2,With:[0,2],__init__:[0,2],__str__:2,_add:2,_appli:2,_compositeunit:0,_cov:0,_fit:2,_intervalunit:0,_nappli:2,_one:0,_radd:2,_simul:0,_sub:2,abov:[0,1,2],abs:[0,2],absolut:[0,2],absoult:2,access:[0,2],account:[0,2],actual:[0,2],add:[0,2],add_symbol:[0,2],added:0,addit:[0,2],additional_nam:[0,2],additional_short_nam:2,advoc:[0,2],affect:[0,2],after:0,again:[0,2],air:0,algebera:2,algorithm:0,alia:[0,2],alias:[0,2],align:2,all:[0,2],allow:[0,2],along:[0,2],alpha:0,alreadi:[0,2],also:[0,1,2],alter:2,altern:[0,1,2],alwai:2,amount:2,analysi:[0,2],analyz:0,angl:[0,2],ani:[0,2],anoth:[0,2],any_but_latex:[0,2],api:1,appear:0,append:0,appli:2,appropri:2,approxim:2,april:1,araound:0,arbitrari:[0,2],arcco:[0,2],arccosh:[0,2],arcsin:[0,2],arcsindist:[0,2],arcsinh:[0,2],arctan2:[0,2],arctan:[0,2],arctanh:[0,2],arg1:[0,2],arg2:[0,2],arg:[0,2],argument:[0,2],arithmet:2,around:2,arrai:[0,2],array_lik:2,ascii:[0,2],ascii_symbol:[0,2],assign:[0,2],associ:[0,2],assum:[0,2],assumpt:[0,2],asterisk:0,asterix:0,astronom:1,attach:[0,2],attempt:[0,2],atto:2,attribut:2,auto:[0,2],autoconvert:2,autocorrel:2,automat:[0,1,2],avail:[0,2],averag:2,avoid:[0,2],axes:0,axhlin:[0,2],axi:[0,2],axvlin:[0,2],back:[0,2],background:2,band:0,bar:0,base:[0,2],base_descript:[0,2],base_prefix:2,bayesian:[0,2],bayesian_default:2,becaus:[0,2],been:[0,2],befor:[0,2],begin:[0,2],behavior:0,below:[0,1,2],best:[0,1,2],between:[0,2],bin:[0,2],binari:[0,2],binary_func:3,binaryprefixedunit:[0,2],binomi:0,binomialdist:[0,2],bipm:2,boadcast:0,bold:2,bool:[0,2],border:2,both:[0,2],bottom:2,bound:0,bracket:[0,2],break_on_printing_error:2,broadcast:[0,2],brochur:[1,2],budget:[0,1,3,4],budget_first_col_header_cel:2,budget_first_col_x_cel:2,budget_header_cel:2,budget_header_row:2,budget_row:2,budget_t:2,budget_x_cel:2,budgetwarn:2,bug:[1,5],built:[1,2],builtin:2,calcul:[0,2],calibr:1,call:[0,2],can:[0,1,2],canada:1,cancel:[0,2],cannot:2,capabl:0,care:2,carlo:[1,2],carri:[0,2],cartesian:0,caus:2,cbrt:[0,2],ccp:[0,2],ccp_from_k:2,ceil:[0,2],cell:2,celsiu:[0,2],center:[0,1,2],centi:2,certain:0,chain:[0,2],chang:[0,2],charact:[0,2],chebyshev:[0,2],chi:[0,2],chi_correct:2,chisquar:[0,2],chisquareddist:[0,2],choic:0,choos:2,ci_mark:[0,2],ci_marker_opt:[0,2],ciformat:0,cik:0,cimethod:[0,2],cioption:0,cip:0,circ:0,circular:[0,2],circularunitconversionerror:[0,2],cisim:[0,2],cisym:2,classmethod:[0,2],clear:2,clformat:0,clk:0,cloption:0,clp:0,cls:2,cm2:1,cmp_k:[0,2],cmp_p:[0,2],codata:1,code:[0,1,2],coeffici:[0,2],color:[0,2],column:[0,2],column_nam:[0,2],combin:[0,2],command:0,commonli:[0,2],compar:2,comparison:[0,2],compat:0,complet:[0,2],complex:[0,1,2],compon:[0,2],composit:[0,2],comprehens:1,comput:[0,2],concis:[0,2],confid:[0,2],confus:1,conj:[0,2],conjug:[0,2],conservative_coverage_prob:2,consol:[0,2],constant:[0,1,2],constitu:[0,2],conta:0,contact:1,contain:[0,1,2],content:4,contribut:[0,2],control:[0,2],convei:1,convers:[1,2],convert:[0,2],convolut:[0,2],coordin:[0,2],copi:[0,2],core:2,correct:0,correctli:2,correl:[0,1,2],correlation_matrix:[0,2],correlation_matrix_sim:2,correlation_sim:2,correspond:[0,2],cos:[0,2],cosh:[0,2],cosin:[0,2],council:1,count:[0,2],cov:[0,2],covari:[0,2],covariance_matrix:[0,2],covariance_matrix_sim:2,covariance_sim:2,coverag:[0,2],coverage_factor:2,coverage_prob:2,covplot:[0,2],covsim:2,covsim_matrix:2,cp_from_k:2,creat:[0,1,2],creation:2,cross:[0,2],css:[0,2],cube:[0,2],cubic:2,cumpord:0,cumprod:[0,2],cumsum:[0,2],curv:[1,2],curvlinear:0,curvlineartrapdist:[0,2],custom:2,custom_head:[0,2],custom_math_mod:[0,2],customari:1,data:[0,1,2],data_format:0,data_opt:0,datafram:[0,2],date:[0,2],deal:[0,1],deca:2,deci:2,decibel:[0,2],decid:2,decreas:[0,2],def:[0,2],default_column:2,default_column_nam:2,default_css:[0,2],defin:[0,1,2],definit:[0,1],defualt:0,deg:[0,2],degf:0,degr:0,degre:[0,1,2],delta:2,delta_diff:2,delta_diff_mean:2,delta_sum:2,delta_sum_mean:2,denomin:0,densiti:[0,2],depend:[0,2],deriv:[0,2],describ:[0,2],descript:[0,2],description_math_mod:[0,2],desir:[0,2],detail:[0,2],determin:[0,2],devat:2,develop:1,deviat:[0,2],df1:[0,2],df2:[0,2],df_html:[0,2],df_latex:[0,2],df_str:[0,2],dfunc:[0,4],dhypot:[0,2],dict:[0,2],dictionari:[0,2],diff:[0,2],differ:[0,1,2],digit:[0,2],dimens:[0,2],dimension:[0,2],dimensionless:[0,2],direchletdist:0,directli:[0,2],dirichlet:0,disabl:0,displai:[2,3],distribut:[1,4],divid:[0,2],divis:[0,2],divmod:[0,2],document:[0,1,2],doe:[0,2],dof:[0,2,3],doffrom:[0,2],don:2,dot:[0,2],doubl:0,doubleexpfit:[0,2],down:2,download:1,dp1:[0,2],dp2:[0,2],draft:[1,2],draw:[0,2],drawn:[0,2],drift:2,dure:[0,2],dx1:[0,2],dx2:[0,2],each:[0,2],ediff1d:[0,2],edit:[0,1,2],effect:[0,2],either:[0,2],element:[0,2],els:0,enclos:[0,2],encod:2,end:0,entri:[0,2],enumer:0,epsilon:0,equal:[0,2],equat:[0,2],equiv:2,equival:[0,2],eras:[0,2],erg:0,error:[0,2],error_bar:0,error_bar_k:0,errorbar:0,estim:2,even:2,evenli:2,event:2,exa:2,exact:2,exactli:[0,2],exampl:[0,2],exbi:2,except:[0,4],exception_on_fmt_error:2,exception_on_warn:3,exclud:2,execut:[0,2],exist:[0,2],exit:[0,2],exp2:[0,2],exp:[0,2,3],expan:0,expand:[0,1,2],expfit:[0,2],explicitli:0,expm1:[0,2],expon:[0,2],exponenti:0,exponentialdist:[0,2],express:[0,1,2],extens:0,extra:[0,2],fab:0,facilit:2,factor:[0,2],fahrenheit:[0,2],fail:0,fals:[0,2,3],feet:2,femto:2,few:[0,2],fhtml:[0,2],figur:0,fill:2,find:[0,2,5],finfo:2,finit:[0,2],fionli:3,first:[0,1,2],fit:[1,4],fit_format:0,fit_opt:0,fit_output:[0,2],fitwarn:2,fix:[0,1,2],flag:[0,2],flase:0,flatex:[0,2],floor:[0,2],floor_divid:[0,2],fmt:[0,2],follow:[0,2],font:2,foot:[0,2],form:[0,2],format:2,format_latex:2,forward:[0,2],found:[0,2],four:[0,2],fraction:[0,1,2],fragment:[0,2],free:[0,1,2],freedom:[0,1,2],frm:2,from:[0,2],from_unit:2,from_uunit:2,full:0,func:[0,2],funicod:[0,2],gamma:0,gammadist:[0,2],gauss:[0,2],gconst:3,gener:[0,2],get:[0,2],get_composit:2,get_p0:[0,2],get_puint:0,get_punit:[0,2],get_valu:2,gibi:[0,2],giga:[0,2],github:1,give:[0,2],given:[0,2],gnu:1,good:0,gradient:[0,2],graft:[0,2],greater:[0,2],guass:0,guess:[0,2],guid:1,gum:[0,1,2],gummi:4,gummytag:2,gummywarn:2,half_width:[0,1,2],handl:[1,2],handmad:1,harold:1,has:[0,2],have:[0,2],head:[0,2],header:[0,2],heavisid:[0,2],heavysid:[0,2],hecto:2,held:[0,2],help:5,here:[0,2,5],hist:[0,2],histogram:[0,2],hold:[0,2],holomorph:[0,2],horizont:[0,2],how:[0,2],howev:[0,1,2],html:[0,2],html_symbol:[0,2],hyperbol:[0,2],hypot:[0,2],iau:1,idea:0,ident:[0,2],ignor:[0,2],ignore_nan:2,imag:[0,2],imaginari:[0,2],immedi:0,implement:[0,1,2],impos:1,improv:1,inch:[0,2],includ:[0,1,2],incompat:2,incompatibleunitserror:2,incorpor:[0,2],incorrectli:1,independ:[0,2],index:[1,2],indic:[0,2],inequ:[0,2],inf:[0,2],infinit:2,inform:[0,2],inherit:[0,2],init:[0,2],initi:[0,2],inlin:[0,2],input:[0,2],insensit:0,insent:2,insert:[0,2],instal:[0,1],instanc:[0,2],instead:2,institut:1,integ:[0,1,2],integr:2,intend:[0,2],interact:0,intern:[1,2],interpol:2,interpret:2,interv:[0,2],invers:[0,2],involv:2,ipython:[0,1,2],is_dimensionless:[0,2],isindepend:2,iso:1,issu:1,iter:[0,2],its:[0,2],itself:0,jac:[0,2],jacobian:[0,2],jgcm:1,joint:[0,2],jummi:2,juptyer:2,jupyt:[0,1,2],just:0,keep:0,kei:[0,2],keyword:[0,2],keywork:0,kibi:[0,2],kilo:[0,2],kilogram:[0,2],known:[0,1],kpa:1,ksim:[0,2],kwd:[0,2],lab:1,label:[0,2],lam:[0,2],laplac:0,laplacedist:[0,2],largest:0,last:[0,2],later:0,latex:[0,2],latex_math:2,latex_math_plot:2,latex_norm:2,latex_norm_plot:2,latex_symbol:[0,2],least:0,leastsq:[0,2],left:[0,2],left_width:[0,2],len:0,length:[0,2],less:[0,2],level:[0,2],librari:[0,2],library_nam:2,lie:[0,2],lies:0,light:[0,2],like:[0,2],limit:[0,2],limit_half_rang:[0,2],line:[0,2],linear:[0,2],linewidth:[0,2],list:[0,1,2],load:[0,2],loc:[0,2],loc_from_k:2,locat:0,log10:[0,2],log1p:[0,2],log2:[0,2],log:[0,2],log_bas:[0,2],log_func:[0,2],logaddexp2:[0,2],logaddexp:[0,2],logarithm:[0,2],logconvers:[0,2],loge:0,lognormaldist:[0,2],logrithm:[0,2],logunit:[0,4],look:[0,1,2],lot:5,lower:[0,2],lower_limit:[0,2],lowest:0,machin:0,machineri:2,made:[0,2],magnitud:0,mai:[0,2],major:1,make:[0,2],make_gummi:3,make_numb:3,mani:[0,1],march:1,markdown:[0,2],marker:[0,2],mask:2,match:[0,2],math:[0,2],mathemat:2,mathrm:2,matrix:[0,2],max_digit:[0,2],max_dof:2,maximum:[0,2],maxit:[0,2],mcisim:[0,2],mcisym:[0,2],mcpropag:2,mean:[0,4],mean_datetim:2,mean_mark:[0,2],mean_marker_opt:[0,2],measur:[0,1,2],measurand:[0,2],mebi:[0,2],mega:[0,2],metagummi:2,metanummi:2,metaprettyprint:2,meter:0,method:2,metr:[0,2],metrolog:1,metrologia:2,mfraction:2,micro:[0,2],microinch:0,microsecond:[0,2],milli:2,millisecond:[0,2],min:2,minim:2,minor:1,minu:2,miscellan:2,miss:[0,2],mix:[0,2],mod:[0,2],mode:[0,2],modf:[0,2],modifi:[0,1,2],modul:[0,1,4],modulu:0,mone:2,mont:[1,2],more:[0,1,2],most:[0,2],mpf:[0,1,2,3],mpmath:[0,1,2],much:1,mulit:0,mulsep:[0,2],multi:2,multinormaldist:[0,2],multipl:[0,2],multipli:[0,2],multitdist:[0,2],multitel:2,multivari:[0,2],multivariatedistribut:[0,2],multivariateel:2,multvariatedistribut:0,must:[0,1,2],n_eff:2,n_sigma:2,n_sigma_trim:2,name:[0,2],nan:2,nano:2,nappli:[0,2],nation:1,natur:[0,2],ndarrai:[0,2],nearest:[0,2],necessari:0,need:[0,2,5],neg:[0,1,2],negat:0,neither:[0,2],neper:[0,2],never:0,newli:[0,2],nist:[1,2],nls:[0,2],non:[0,2],none:[0,2,3],nonlinear:[0,2],nonlinearconvers:2,nonlinearunit:[0,4],nor:2,norm:2,normal:[0,2],normaldist:[0,2],nosimulateddataerror:[0,2],notat:[0,2],note:[0,1,2],notebook:[0,1,2],notimplementederror:[0,2],nounitconversionfounderror:[0,2],now:[0,2],nparam:[0,2],nprop:[0,2],nrc:1,nsig:[0,2],number:[0,2],numer:2,nummi:4,numpi:[0,1,2],object:[0,1,2],observ:[0,2],odd:2,odr:[0,2],off:2,offset:[0,2],offsetconvers:[0,2],offsetunit:[0,4],often:[0,2],older:[0,2],ols:0,omit:[0,2],ommit:0,onc:[0,2],one:[0,2],oneovertfit:[0,2],onli:[0,2],oper:2,operand:0,opertaion:2,optim:[0,2],option:[0,2],order:[0,2],origin:[0,2],other:[0,1,2],otherwis:[0,2],out:[0,2],outlier:2,outlin:1,output:[0,2],outsid:2,over:[0,2],over_plot:0,overrid:[0,2],overridden:[0,2],own:2,p_method:[0,2],packag:[0,1,4],pad:2,page:1,pair:0,panda:[0,1,2],param:2,paramet:2,paramt:0,parent:2,park:1,parksh:1,pars:2,part:[0,2],partial:[0,2],particular:[0,2],pass:[0,2],pebi:2,perform:[0,2],perhap:1,peta:2,phi:[0,2],physic:[0,1,2],pick:[0,2],pico:2,pip:1,place:[0,2],pleas:1,plot:[0,2,3],plot_opt:[0,2],plot_point:0,plu:2,plusmn:1,pmethod:4,pmi:[0,2],pmsim:[0,2],pmsimi:[0,2],point:[0,2],poisson:0,poissondist:[0,2],polar:0,polyfit:[0,2],polynomi:0,popul:[0,2],posit:2,possibl:[0,2],pow:2,power:[0,2],ppm:[0,2],practic:1,predict:[0,2],prefer:[0,2],prefix:[0,2],prefix_definit:2,prefixedunit:[0,4],pressur:0,pretti:[0,2],prettyprint:2,previou:0,previous:[0,2],primarili:2,print:[0,4],print_html:2,print_latex:2,print_markdown:2,printer:[0,2],prnt:[0,2,3],probabl:[0,2],problem:2,procedur:2,prod:[0,2],propag:[0,1,2],properti:2,proport:0,provid:2,ptostr:2,publish:2,punit:[0,2],pure:1,put:[0,2,5],pyplot:[0,2],python:[0,1,2],quadrant:[0,2],quaitii:0,quanit:[0,2],quantiti:[0,1,2],rais:[0,2],random:[0,2],random_st:[0,2],randomst:[0,2],rare:0,rate:[0,2],rather:[0,2],ratiounit:2,raw:[0,2],rchain:2,read:[0,2],readi:2,real:[0,2],reciproc:[0,2],reciprocalconvers:2,recogn:0,recommend:1,red:[0,2],reduc:[0,2],reduce_unit:[0,2],refer:[0,1,2],referenc:[0,1,2],rel:[0,2],relat:2,relationship:0,releas:1,relunit:4,remaind:[0,2],remov:[0,2],reorder:[0,2],repres:[0,1,2],represent:[0,2],requir:[0,1,2],res:[0,2],research:1,resect:2,resid:2,residu:[0,2],respect:[0,2],respons:[0,2],restrict:1,result:[0,2],retriev:[0,2],right:2,right_width:[0,2],rint:[0,2],root:[0,2],round:[0,2],rounding_u:[0,2],row:[0,2],rss:0,rule:[0,2],run:[0,2],s276:2,same:[0,1,2],sampl:[0,2],satterthwait:2,scalar:[0,2],scale:[0,2],scatter:[0,2],sci_not:[0,2],sci_notation_high:[0,2],sci_notation_low:[0,2],scientif:[0,2],scipi:[0,1,2],search:[1,2],search_unit:[0,2],second:[0,2],section:[0,1],see:[0,1,2],seed:[0,2],select:[0,2],self:[0,2],sensit:[0,2],separ:[0,2],seri:2,serv:[0,2],session:[0,2],set:[0,1,2],set_print:2,set_se:[0,2],set_tag:2,sever:[0,1,2],shadow:[0,2],shadow_unit:2,shadowed_alias:[0,2],shadowed_unit:2,shape:[0,2],share:[0,2],shift:[0,2],short_nam:[0,2],shortest:[0,2],should:[0,2],show:[0,2],show_:[0,2],show_al:[0,2],show_c:[0,2],show_d:[0,2],show_data:0,show_dof:[0,2],show_expanded_u:[0,2],show_fit:0,show_k:[0,2],show_nam:[0,2],show_p:[0,2],show_subtot:[0,2],shown:[0,2],side:2,sigma:[0,2],sigma_is_known:[0,2],sigma_trim:2,sign:[0,2,3],signal:2,signatur:[0,2],signific:[0,2],sim:[0,2,3],simdata:[0,2],similar:[0,2],simpli:2,simsort:[0,2],simul:2,simultan:[0,2],sin:[0,2],sinc:0,sine:[0,2],sinfit:[0,2],singl:[0,2],sinh:[0,2],siunit:[0,4],size:[0,2],slash:[0,2],slashaxi:[0,2],slightli:1,slow:2,smallest:0,softwar:1,solid:2,solidu:[0,2],solver:[0,2],some:[1,2],someth:5,sometim:[0,2],sometimes_smal:3,sort:[0,2],sound:0,sourc:1,space:[0,2],special:[1,2],specifi:[0,2],spl:0,spline:2,sqrt:[0,2],squar:[0,2],stabl:[0,2],stack:1,standard:[0,1,2],start:[0,5],statement:2,staticmethod:[0,2],stdev:2,str:[0,2],straight:0,string:[0,1,2],strip:2,student:[0,2],style:[0,2],sub:2,subclass:[0,2],submodul:4,subpackag:4,substitut:0,subtot:[0,2],subtract:[0,2],subtrat:0,success:0,sum:[0,2],suppli:[0,2],support:[0,1,2],symbol:[0,2],symmetr:[0,2],system:[1,2],tab:0,tabl:[0,2],tack:2,tag:2,tag_nam:2,take:[0,2],taken:[0,2],tan:[0,2],tangent:[0,2],tanh:[0,2],tdist:[0,2],tebi:2,teh:2,temperatur:0,tera:2,term:[0,1,2],test:[2,4],test_ab:3,test_add:3,test_addxmul:3,test_ap1sinco:3,test_ap2sinco:3,test_apnsinco:3,test_creat:[2,4],test_div:3,test_gummy_init:3,test_mod:3,test_mul:3,test_mulxnpow:3,test_mulxpow:3,test_neg:3,test_npsinco:3,test_po:3,test_pow:3,test_sinco:3,test_sub:3,test_ubreakdown:3,text:[0,2],than:[0,1,2],thei:[0,1,2],them:[0,2],themont:0,thi:[0,1,2],thick:[0,2],thin:2,those:[0,2],though:[0,2],thousand_spac:2,ths:2,thu:[0,2],time:0,titl:[0,2],title_styl:[0,2],to_unit:2,to_uunit:2,toascii:[0,2],tofloat:[0,2],tohtml:[0,2],tolatex:[0,2],tool:1,top_to_base_ratio:[0,2],tostr:[0,2],total:[0,2],tounicod:[0,2],toward:[0,2],track:[0,2],tradit:[0,2],tradition:[0,2],transcendent:0,trapezoid:0,trapezoidaldist:[0,2],treat:[0,2],tri:0,trial:0,triangular:0,triangulardist:[0,2],true_divid:[0,2],trunc:[0,2],try_fconvert:0,tupl:[0,2],turn:2,tutori:1,two:[0,2],txt:[0,2],type:[0,1,2],u00b0f:0,ubreakdown:[0,2],ueq:[0,2],uexp:3,ufrom:[0,2],ufsim:[0,2],ummi:[0,4],uncertainti:[0,1,2],uncertianti:0,uncertiantyprecisionwarn:2,uncorrel:[0,2],under:[0,1,2],unicod:[0,2],uniform:0,uniformdist:[0,1,2],unit:[1,3,4],uniterror:2,unitliberror:[0,2],unitlibnotfounderror:2,unitlibrari:2,unitnotfounderror:[0,2],units_on_valu:[0,2],unitutil:4,unitwarn:2,unless:[0,2],unpredict:2,unshadow:[0,2],until:0,upa:0,upon:[0,2],upper:[0,2],upper_limit:[0,2],use:[0,2],use_n_eff:2,used:[0,1,2],useful:[0,2],user:[0,2],uses:[0,2],usim:[0,2],using:[0,1,2],usual:[0,2],usunit:[0,4],usym:2,utyp:[0,2],uunit:[0,1,2],uunit_is_rel:[0,2],valid:[0,2],valu:[0,1,2],valueerror:2,vari:[0,2],variabl:[0,2],varianc:[0,2],variou:0,verbatim:1,veri:[0,2],verison:2,version:[0,4],versu:2,vertic:[0,2],view:[0,2],wai:[0,2],want:[0,1,2],warn:[0,2],weakvaluedictionari:2,weibul:0,weibulldist:[0,2],weight:[0,2],weird:0,welch:2,well:2,were:[0,1],when:[0,2],whenev:0,where:[0,2],whether:[0,2],which:[0,1,2],who:[0,2],width:2,willink:2,wish:[0,1,2],without:[0,2],wmean:2,word:[0,2],work:[0,2],would:[0,2],write:[0,2],xdim:[0,2],xfsim:[0,2],xlabel:[0,2],xlist:[0,2],xmax:0,xmin:0,xname:[0,2],xsim:[0,2],xsym:2,xunit:[0,2],yard:[0,2],ydim:[0,2],year:[0,2],ylabel:[0,2],yname:[0,2],yobi:2,yocto:2,yotta:2,you:[0,1,2],ypred:[0,2],ypredf:[0,2],yunit:[0,2],zepto:2,zero:[0,2],zetta:2,zhang:2,zibi:2},titles:["MetroloPy","MetroloPy, the docs","metrolopy package","metrolopy.tests package","metrolopy","to do"],titleterms:{"abstract":0,"class":[0,1],"function":[0,2],"static":0,appli:0,arithmet:0,attribut:0,author:1,auto:1,basic:0,budget:2,built:0,carlo:0,come:1,common:0,content:[2,3],convers:0,curv:0,custom:0,dfunc:2,did:1,displai:0,distribut:[0,2],doc:1,except:2,fit:[0,2],format:0,from:1,further:1,gener:1,get:1,gummi:[0,1,2],histori:1,indic:1,involv:0,jummi:0,licens:1,logunit:2,mathemat:0,mean:2,method:0,metrolopi:[0,1,2,3,4],modul:[2,3],mont:0,multi:0,name:1,nonlinearunit:2,numer:0,nummi:2,offsetunit:2,oper:0,packag:[2,3],paramet:0,pmethod:2,prefixedunit:2,print:2,properti:0,read:1,relat:0,relunit:2,search:0,simul:0,siunit:2,some:0,start:1,sub:0,submodul:[2,3],subpackag:2,tabl:1,test:3,test_creat:3,ummi:2,unit:[0,2],unitutil:2,usunit:2,variat:0,version:[1,2],where:1}})
\ No newline at end of file
diff --git a/docs/_build/html/todo.html b/docs/_build/html/todo.html
index f808bd7..301a769 100644
--- a/docs/_build/html/todo.html
+++ b/docs/_build/html/todo.html
@@ -6,7 +6,7 @@
   <head>
     <meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
     
-    <title>to do &#8212; metrolopy 0.5.1 documentation</title>
+    <title>to do &#8212; metrolopy 0.5.2 documentation</title>
     
     <link rel="stylesheet" href="_static/nature.css" type="text/css" />
     <link rel="stylesheet" href="_static/pygments.css" type="text/css" />
@@ -14,7 +14,7 @@
     <script type="text/javascript">
       var DOCUMENTATION_OPTIONS = {
         URL_ROOT:    './',
-        VERSION:     '0.5.1',
+        VERSION:     '0.5.2',
         COLLAPSE_INDEX: false,
         FILE_SUFFIX: '.html',
         HAS_SOURCE:  true,
@@ -38,7 +38,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>  
 
@@ -91,7 +91,7 @@ <h3>Navigation</h3>
         <li class="right" >
           <a href="py-modindex.html" title="Python Module Index"
              >modules</a> |</li>
-        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.1 documentation</a> &#187;</li> 
+        <li class="nav-item nav-item-0"><a href="index.html">metrolopy 0.5.2 documentation</a> &#187;</li> 
       </ul>
     </div>
     <div class="footer" role="contentinfo">
diff --git a/docs/hand_made_doc.rst b/docs/hand_made_doc.rst
index fca3d2a..2ed1099 100644
--- a/docs/hand_made_doc.rst
+++ b/docs/hand_made_doc.rst
@@ -12,7 +12,7 @@ class gummy
 ===========
 
 class metrolopy.\ **gummy(x, u=0, unit=one, dof=float('inf'), k=1, p=None, 
-p_method=None, uunit=None, utype=None, name=None)**
+uunit=None, utype=None, name=None)**
 
 A gummy object represents a numerical value with an uncertainty and (or)
 a unit. They can be used in place of float values in Python expressions
@@ -97,21 +97,6 @@ gummy parameters
        0.22
        >>> g.u
        0.11224696052342387
-
-.. _p_method:
-
--  **p_method**: ({'loc', 'cp', 'gauss',
-   'ccp', 'chebyshev', None}) If the *p* parameter is
-   specified, *p_method* sets the method that is used to calculate the
-   coverage factor. If p_method is omitted or set to 'loc', then the
-   uncertainty is assumed to be represented by a normal probability
-   distribution if *dof* = float('inf') and shifted and scaled Student's
-   t distribution otherwise. If *p_method* = 'gauss' or 'cp' then the
-   Guass inequality is used, and if *p_method* = 'chebyshev' or 'ccp'
-   then the Chebyshev inequality is used. For *p* = 0.95 and *dof* =
-   float('inf'), *p_method* = 'loc' gives *k* = 2.0, while *p_method*
-   = 'gauss' gives *k* = 3.0 and *p_method* = 'chebyshev' gives *k* =
-   4.5.
    
 .. _parameter-uunit:
 
@@ -196,9 +181,20 @@ basic gummy properties
    the relation between the value of this property and the property *k*
    is defined by the *p_method* property
    
-.. _property-p-method:
+.. _p_method:
 
--  **p_method**: see the p_method_ parameter
+-  **p_method**: ({'loc', 'cp', 'gauss',
+   'ccp', 'chebyshev', None}) If the *p* parameter is
+   specified, *p_method* sets the method that is used to calculate the
+   coverage factor. If p_method is omitted or set to 'loc', then the
+   uncertainty is assumed to be represented by a normal probability
+   distribution if *dof* = float('inf') and shifted and scaled Student's
+   t distribution otherwise. If *p_method* = 'gauss' or 'cp' then the
+   Guass inequality is used, and if *p_method* = 'chebyshev' or 'ccp'
+   then the Chebyshev inequality is used. For *p* = 0.95 and *dof* =
+   float('inf'), *p_method* = 'loc' gives *k* = 2.0, while *p_method*
+   = 'gauss' gives *k* = 3.0 and *p_method* = 'chebyshev' gives *k* =
+   4.5.  This property may only be set at the class level.
 
 .. _name:
 
@@ -270,6 +266,9 @@ basic gummy properties
    can give some idea of the magnitude of the floating point errors, but is 
    not a substitute for a full numerical error analysis.  The default value
    is ``False``.
+   
+-  **cmp_k**, **cmp_p**: k or p values for comparisons, e.g. a > b is True
+   if c.x > k*c.u where c = b - a.  Set only at the class level.
 
 basic gummy methods
 -------------------
diff --git a/docs/index.rst b/docs/index.rst
index 8a2106c..6be011f 100644
--- a/docs/index.rst
+++ b/docs/index.rst
@@ -121,6 +121,9 @@ version history
 * Version 0.5.1, built 2 April 2019, fixed a major bug that generated negative 
   uncertainties in some cases and fixed some other minor bugs.  Improved support 
   for fraction.Fraction and mpmath.mpf values.
+* Version 0.5.2, built 5 April 2019, fixed a major bug that propagated
+  uncertainty incorrectly if a gummy was created with an uncertainty set with
+  an integer data type.  Fixed several other minor bugs.
 
 
 author
diff --git a/metrolopy/gummy.py b/metrolopy/gummy.py
index 35f82a8..0d1e9f7 100644
--- a/metrolopy/gummy.py
+++ b/metrolopy/gummy.py
@@ -154,7 +154,7 @@ def cmp_k(cls,v):
         if v <= 0:
             raise ValueError('cmp_k <= 0')
         gummy._cmp_k = v
-        gummy._cmp_p = None # k and p may not both be defined
+        gummy._cmp_p = None
         
     # The default level of confidence for new ummies
     @property
@@ -169,7 +169,7 @@ def cmp_p(cls,v):
         if v <= 0 or v >= 1:
             raise ValueError('cmp_p is not in the interval (0,1)')
         gummy._cmp_p = v
-        gummy._cmp_k = None # k and p may not both be defined
+        gummy._cmp_k = None
         
     @property
     def p_method(cls):
@@ -187,9 +187,11 @@ def p_method(cls):
         `dof` = `float('inf')`, `p_method` = 'loc' gives `k` = 2.0, while
         `p_method` = 'gauss' gives `k` = 3.0 and p_method = 'chebyshev' gives `k` = 4.5.
         """
-        return gummy._Pmthd.method
+        return gummy._p_method.method
     @p_method.setter
     def p_method(cls,v):
+        if v is None:
+            v = 'loc'
         gummy._p_method = _Pmthd(v)
         
     
@@ -306,8 +308,8 @@ class gummy(PrettyPrinter,nummy,metaclass=MetaGummy):
                                     # print the will cause _ret_something() to 
                                     # be called.  Can be set to True for debugging
     
-    def __init__(self,x,u=0,unit=one,dof=float('inf'),k=1,p=None,p_method=None,
-                 uunit=None,utype=None,name=None):
+    def __init__(self,x,u=0,unit=one,dof=float('inf'),k=1,p=None,uunit=None,
+                 utype=None,name=None):
         if isinstance(x,ummy):
             self._copy(x,self,formatting=False)
             return
@@ -328,9 +330,6 @@ def __init__(self,x,u=0,unit=one,dof=float('inf'),k=1,p=None,p_method=None,
             else:
                 uunit = None
             
-        if p_method is not None:
-            self._p_method = _Pmthd(p_method)
-            
         if p is not None:
             p = float(p)
             self._k = self._p_method.fptok(p,dof,gummy.bayesian)
@@ -661,29 +660,6 @@ def uunit_is_rel(self):
         except:
             return False
         
-    @property
-    def p_method(self):
-        """`str` in {'loc', 'cp', 'gauss', 'ccp', 'chebyshev'}
-
-        This sets the default `p_method` attribute for newly created gummys which
-        determines how the coverage factor is calculated from a given level of
-        confidence or coverage probability `p`.
-
-        If `p_method` = 'loc', then the uncertainty is assumed to be represented
-        by a normal probability distribution if `dof` = `float('inf')` and shifted
-        and scaled Student's t distribution otherwise.  If `p_method` = 'gauss' or
-        'cp' then the Gauss inequality is used, and if `p_method` = 'chebyshev' or
-        'ccp' then the Chebyshev inequality is used.  For `p` = 0.95 and
-        `dof` = `float('inf')`, `p_method` = 'loc' gives `k` = 2.0, while
-        `p_method` = 'gauss' gives `k` = 3.0 and p_method = 'chebyshev' gives `k` = 4.5.
-        """
-        
-        return self._p_method.method
-    @p_method.setter
-    def p_method(self,v):
-        self._p_method = _Pmthd(v)
-        self.p = self._p
-        
     @property
     def k(self):
         """
@@ -1053,8 +1029,6 @@ def graft(self,unit,uunit=None):
             g._k = self._k
         g._pm = self._pm
         g._p = self._p
-        if self._p_method is not gummy._p_method:
-            g._p_method = self._p_method
         if self.thousand_spaces != gummy.thousand_spaces:
             g.thousand_spaces = self.thousand_spaces
         if self.mulsep != gummy.mulsep:
@@ -1094,8 +1068,6 @@ def _copy(s,r,formatting=True,tofloat=False):
                     r._k = s._k
                 r._p = s._p
                 r._pm = s._pm
-                if s._p_method is not gummy._p_method:
-                    r._p_method = s._p_method
                 r._set_k = s._set_k
                 if tofloat:
                     r._set_U(None,None)
@@ -2767,20 +2739,8 @@ def __eq__(self, v):
                     return False
             else:
                 s = self
-            
-        if v is s:
-            return True
-        if isinstance(v,gummy):
-            if s._u == 0 and v._u == 0:
-                return self._x == v._x
-            if s._ref is v._ref and self._refs == v._refs:
-                if s._x == v._x and self._u == v._u:
-                    return True
-            return False
-        elif s._u == 0:
-            return s._x == v
-        else:
-            return False
+
+        return super(gummy,s).__eq__(v)
     
     def __ne__(self, v):
         if isinstance(v,gummy):
@@ -2799,91 +2759,63 @@ def __ne__(self, v):
             
         if s.__eq__(v):
             return False
-        if s.__lt__(v) or s.__gt__(v):
-            return True
-        return False
+        return s.__lt__(v) or s.__gt__(v)
         
     def __lt__(self, v):
         if isinstance(v,gummy):
             try:
                 s = self.convert(v._unit)
             except NoUnitConversionFoundError:
-                raise IncompatibleUnitsError('values with incompatible units cannot be compared')
+                raise IncompatibleUnitsError('values with incompatible units cannot be compared with > or <')
         else:
             if self._unit is not one:
                 try:
                     s = self.convert(one)
                 except NoUnitConversionFoundError:
-                    raise IncompatibleUnitsError('values with incompatible units cannot be compared')
+                    raise IncompatibleUnitsError('values with incompatible units cannot be compared with > or <')
             else:
                 s = self
             
         df = s - v
-        if s._cmp_k is None:
-            k = self._p_method.fptok(gummy._cmp_p,df.dof,df.bayesian)
-        if df._x < -k*s._u:
-            return True
-        return False
-        
-    def __le__(self, v):
-        if isinstance(v,gummy):
-            try:
-                s = self.convert(v._unit)
-            except NoUnitConversionFoundError:
-                raise IncompatibleUnitsError('values with incompatible units cannot be compared')
+        if self._cmp_k is None:
+            k = self._p_method.fptok(self._cmp_p,df.dof,df.bayesian)
         else:
-            if self._unit is not one:
-                try:
-                    s = self.convert(one)
-                except NoUnitConversionFoundError:
-                    raise IncompatibleUnitsError('values with incompatible units cannot be compared')
-            else:
-                s = self
+            k = self._cmp_k
             
-        if s.__eq__(v):
+        return (df._x < -k*df._u)
+        
+    def __le__(self, v):
+        if self.__eq__(v):
             return True
-        return s.__lt__(v)
+        return self.__lt__(v)
         
     def __gt__(self, v):
         if isinstance(v,gummy):
             try:
                 s = self.convert(v._unit)
             except NoUnitConversionFoundError:
-                raise IncompatibleUnitsError('values with incompatible units cannot be compared')
+                raise IncompatibleUnitsError('values with incompatible units cannot be compared with > or <')
         else:
             if self._unit is not one:
                 try:
                     s = self.convert(one)
                 except NoUnitConversionFoundError:
-                    raise IncompatibleUnitsError('values with incompatible units cannot be compared')
+                    raise IncompatibleUnitsError('values with incompatible units cannot be compared with > or <')
             else:
                 s = self
             
-        df = self - v
-        if s._cmp_k is None:
-            k = self._p_method.fptok(gummy._cmp_p,df.dof,df.bayesian)
-        if df._x > k*s._u:
-            return True
-        return False
-        
-    def __ge__(self, v):
-        if isinstance(v,gummy):
-            try:
-                s = self.convert(v._unit)
-            except NoUnitConversionFoundError:
-                raise IncompatibleUnitsError('values with incompatible units cannot be compared')
+        df = s - v
+        if self._cmp_k is None:
+            k = self._p_method.fptok(self._cmp_p,df.dof,df.bayesian)
         else:
-            if self._unit is not one:
-                try:
-                    s = self.convert(one)
-                except NoUnitConversionFoundError:
-                    raise IncompatibleUnitsError('values with incompatible units cannot be compared')
-            else:
-                s = self
+            k = self._cmp_k
             
-        if s.__eq__(v):
+        return (df._x > k*df._u)
+        
+    def __ge__(self, v):
+        if self.__eq__(v):
             return True
-        return s.__gt__(v)
+        return self.__gt__(v)
     
     @property
     def imag(self):
diff --git a/metrolopy/tests/test_create.py b/metrolopy/tests/test_create.py
index d7f7514..1c9d2fa 100644
--- a/metrolopy/tests/test_create.py
+++ b/metrolopy/tests/test_create.py
@@ -13,6 +13,7 @@
 rand = np.random.RandomState()
 
 def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
+    print('uc.gummy.p_method may be modified\nset uc.gummy.p_method = \'loc\' to recover the default value')
     from scipy.stats import norm, t
     
     if n is None:
@@ -73,9 +74,9 @@ def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
             
             if rand.randint(2):
                 pmethods = ['loc','level of confidence','cp','coverage probability','gauss','ccp','ccp','conservative coverage probability','chebyshev']
-                p_method = pmethods[rand.randint(len(pmethods))]
+                uc.gummy.p_method = pmethods[rand.randint(len(pmethods))]
             else:
-                p_method = None
+                uc.gummy.p_method = None
             
             if unit == 'dB(SPL)':
                 if rand.randint(10) == 0:
@@ -116,6 +117,8 @@ def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
                 else:
                     if rand.randint(20) == 0:
                         U = 1e6*abs(x)*rand.rand()
+                    elif rand.randint(20) == 0:
+                        U = rand.randint(20)+1
                     else:
                         U = 1.1*abs(x)*rand.rand()
                 u = U
@@ -147,9 +150,11 @@ def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
             uc.gummy.bayesian = bayesian
             if uunit is not None and rand.randint(2):
                 uu = uc.gummy(U,unit=uunit)
-                g = uc.gummy(x,u=uu,unit=unit,dof=dof,k=k,p=p,p_method=p_method,name=name,utype=utype)
+                g = uc.gummy(x,u=uu,unit=unit,dof=dof,k=k,p=p,name=name,
+                             utype=utype)
             else:
-                g = uc.gummy(x,u=U,unit=unit,dof=dof,k=k,p=p,p_method=p_method,uunit=uunit,name=name,utype=utype)
+                g = uc.gummy(x,u=U,unit=unit,dof=dof,k=k,p=p,uunit=uunit,
+                             name=name,utype=utype)
                 
             if not const:
                 u = u/g.k
@@ -245,7 +250,6 @@ def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
                 reinit = True
                 uunit = None
                 name = None
-                p_method = 'loc'
             else:
                 reinit = False
     
@@ -304,7 +308,7 @@ def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
             else:
                 assert g.p == p
             
-            if p_method in ['loc','level of confidence',None]:
+            if uc.gummy.p_method in ['loc','level of confidence',None]:
                 if dof == float('inf'):
                     assert abs(g.k - norm.ppf(0.5 + g.p/2)) < 1e-6
                 else:
@@ -314,7 +318,7 @@ def test_gummy_init(n=None,exception_on_warning=True,prnt=False,plot=False):
                     else:
                         assert abs(g.k - t.ppf(0.5 + g.p/2,dof)) < 1e-6
                         
-            elif p_method in ['cp','coverage probability','gauss']:
+            elif uc.gummy.p_method in ['cp','coverage probability','gauss']:
                 if dof == float('inf') or bayesian:
                     if g.p is not None:
                         assert abs(g.k - 2/(3*(np.sqrt(1 - g.p)))) < 1e-6
@@ -393,20 +397,26 @@ def display(g):
         g.ascii()
                 
 def make_gummy(sign=None,exp=None,uexp=-6,sometimes_small=True,dof=None,
-               unit=None,units=None,mpf=False):
+               unit=None,units=None,mpf=False,allowzero=True,allowlargeu=True):
     
     utypes = ['A','B','D',None]
     if units is None:
         units = ['m','lb','m**2 s**3/kg**4','m**3','s**-2','1']
         
-    x = rand.rand()*0.9 + 0.1
-    
-    if exp is None:
-        exp = rand.randint(1,6)
-        if rand.randint(2):
-            exp = -exp
-            
-    x = x*10**exp
+    if rand.randint(10) == 0:
+        if exp is None:
+            x = rand.randint(6000)
+        else:
+            x = rand.randint(int(10**exp))
+    else:
+        x = rand.rand()*0.9 + 0.1
+        
+        if exp is None:
+            exp = rand.randint(1,6)
+            if rand.randint(2):
+                exp = -exp
+                
+        x = x*10**exp
     
     if sign is None:
         if rand.randint(2):
@@ -414,10 +424,16 @@ def make_gummy(sign=None,exp=None,uexp=-6,sometimes_small=True,dof=None,
     else:
         x *= sign
         
-    u = rand.rand()*0.5 + 0.5
-    if sometimes_small and not rand.randint(10):
-        uexp -= 4
-    u = abs(x)*u*10**uexp
+    if not allowzero and x == 0:
+        x = x + 1
+        
+    if allowlargeu and rand.randint(10) == 0:
+        u = rand.randint(20) + 1
+    else:
+        u = rand.rand()*0.5 + 0.5
+        if sometimes_small and not rand.randint(10):
+            uexp -= 4
+        u = abs(x)*u*10**uexp
     
     if mpf and mp is not None and not rand.randint(4):
         x = mp.mpf(x)
@@ -439,7 +455,7 @@ def make_gummy(sign=None,exp=None,uexp=-6,sometimes_small=True,dof=None,
     
     return (g,x,u,dof,unit,utype)
 
-def make_number(sign=None,gconst=True,exp=None,fionly=True):
+def make_number(sign=None,gconst=True,exp=None,fionly=True,allowzero=True):
     q = rand.randint(4)
     if fionly and q > 1:
         q = 1
@@ -479,17 +495,23 @@ def make_number(sign=None,gconst=True,exp=None,fionly=True):
     
     if gconst and rand.randint(2):
         x = uc.gummy(x)
+    
+    if not allowzero and x == 0:
+        x = x + 1
         
     return x
 
-def binary_func(f,df,sim=False,exp=None,fionly=False,uexp=-6):
+def binary_func(f,df,sim=False,exp=None,fionly=False,uexp=-6,allowazero=True,
+                allowbzero=True,allowlargeu=True):
     if sim:
         dof = dof=float('inf')
     else:
         dof = None
-    a,ax,au,adof,_,_ = make_gummy(unit=1,dof=dof,exp=exp,mpf=(not fionly),uexp=uexp)
+    a,ax,au,adof,_,_ = make_gummy(unit=1,dof=dof,exp=exp,mpf=(not fionly),
+                                  uexp=uexp,allowzero=allowazero,
+                                  allowlargeu=allowlargeu)
     if rand.randint(2):
-        b = make_number(exp=exp,fionly=fionly)
+        b = make_number(exp=exp,fionly=fionly,allowzero=allowbzero)
         if isinstance(b,uc.gummy):
             bx = b.x
         else:
@@ -497,7 +519,9 @@ def binary_func(f,df,sim=False,exp=None,fionly=False,uexp=-6):
         bu = 0
         bdof = float('inf')
     else:
-        b,bx,bu,bdof,_,_ = make_gummy(unit=1,dof=dof,exp=exp,mpf=(not fionly),uexp=uexp)
+        b,bx,bu,bdof,_,_ = make_gummy(unit=1,dof=dof,exp=exp,mpf=(not fionly),
+                                      uexp=uexp,allowzero=allowbzero,
+                                      allowlargeu=allowlargeu)
         
     ax = float(ax)
     bx = float(bx)
@@ -527,7 +551,7 @@ def binary_func(f,df,sim=False,exp=None,fionly=False,uexp=-6):
     if x == 0:
         assert abs(g.x) < 1e-15
     else:
-        assert abs((g.x - x)/x) < 1e-15
+        assert abs((g.x - x)/x) < 1e-14
         
     if u == 0:
         assert abs(u) < 1e-6
@@ -537,7 +561,7 @@ def binary_func(f,df,sim=False,exp=None,fionly=False,uexp=-6):
     if dof is not None and not np.isinf(dof):
         assert abs((g.dof - dof)/dof) < 0.01
     
-    if sim and abs(u/x) > 1e-10:
+    if sim and x != 0 and abs(u/x) > 1e-10:
         uc.gummy.simulate([a,b,g])
         
         assert abs((g.xsim - x)/(g.u/np.sqrt(1e5))) < 5
@@ -578,8 +602,12 @@ def f(a,b):
         return a/b
     def df(ax,bx):
         return (1/bx,-ax/bx**2)
+    if sim:
+        lu = False
+    else:
+        lu = True
     for i in range(n):
-        binary_func(f,df,sim=sim)
+        binary_func(f,df,sim=sim,allowbzero=False,allowlargeu=lu)
         
 def test_pow(n=1000,sim=False):
     def f(a,b):
@@ -587,8 +615,12 @@ def f(a,b):
     def df(ax,bx):
         ax = abs(ax)
         return (bx*ax**(bx-1),np.log(ax)*ax**bx)
+    if sim:
+        lu = False
+    else:
+        lu = True
     for i in range(n):
-        binary_func(f,df,sim=sim,exp=0)
+        binary_func(f,df,sim=sim,exp=0,allowazero=False,allowlargeu=lu)
         
 def test_mod(n=1000,sim=False):
     def f(a,b):
@@ -603,8 +635,13 @@ def f(a,b):
         return abs(a)*abs(b)
     def df(ax,bx):
         return (abs(bx),abs(ax))
+    if sim:
+        lu = False
+    else:
+        lu = True
     for i in range(n):
-        binary_func(f,df,sim=sim)
+        binary_func(f,df,sim=sim,allowlargeu=lu,allowazero=lu,
+                    allowbzero=lu)
         
 def test_neg(n=1000,sim=False):
     def f(a,b):
@@ -644,7 +681,7 @@ def f(a,b):
     def df(ax,bx):
         return (np.cos(ax),-np.sin(bx))
     for i in range(n):
-        binary_func(f,df,sim=sim,fionly=True,exp=0)
+        binary_func(f,df,sim=sim,fionly=True,exp=0,allowlargeu=False)
         
 def test_ap2sincos(n=1000,sim=False):
     def ff(a,b):
@@ -666,8 +703,95 @@ def f(a,b):
         return uc.gummy.napply(ff,a,b)
 
     for i in range(n):
-        binary_func(f,df,sim=sim,fionly=True,exp=0)
+        binary_func(f,df,sim=sim,fionly=True,exp=0,allowlargeu=False)
+        
+def test_addxmul(n=1000):
+    def f(*x):
+        r = x[0]
+        for i in range(len(x)-1):
+            r += x[i+1]
+        return r
+        
+    for j in range(n):
+        k = rand.randint(10) + 2
+        x = make_gummy(unit=1)[0]
+        a = f(*(k*[x]))
+        b = k*x
+        c = uc.gummy.apply(f,lambda *x: len(x)*[1],*(k*[x]))
+        d = uc.gummy.napply(f,*(k*[x]))
+        
+        assert abs((a.x - b.x)/b.x) < 1e-14
+        assert (a.u - b.u)/b.u < 1e-3
+        assert a.correlation(b) == 1
+        assert a.correlation(x) == 1
+        
+        assert abs((c.x - b.x)/b.x) < 1e-14
+        assert (c.u - b.u)/b.u < 1e-3
+        assert c.correlation(b) == 1
+        
+        assert abs((d.x - b.x)/b.x) < 1e-14
+        assert (d.u - b.u)/b.u < 1e-3
+        assert d.correlation(b) == 1
+        
+def test_mulxpow(n=1000):
+    def f(*x):
+        r = x[0]
+        for i in range(len(x)-1):
+            r *= x[i+1]
+        return r
+        
+    for j in range(n):
+        k = rand.randint(10) + 2
+        x = make_gummy(unit=1)[0]
+        a = f(*(k*[x]))
+        b = x**k
+        xx = x.x
+        c = uc.gummy.apply(f,lambda *x: len(x)*[xx**(k-1)],*(k*[x]))
+        d = uc.gummy.napply(f,*(k*[x]))
+        
+        assert abs((a.x - b.x)/b.x) < 1e-14
+        assert abs((a.u - b.u)/b.u) < 1e-3
+        assert a.correlation(b) == 1
+        
+        if k%2 == 0 and x.x < 0:
+            assert a.correlation(x) == -1
+        else:
+            assert a.correlation(x) == 1
+            
+        assert abs((c.x - b.x)/b.x) < 1e-14
+        assert abs((c.u - b.u)/b.u) < 1e-3
+        assert c.correlation(b) == 1
+        
+        assert abs((d.x - b.x)/b.x) < 1e-14
+        assert abs((d.u - b.u)/b.u) < 1e-3
+        assert d.correlation(b) == 1
+        
+def test_mulxnpow(n=1000):
+    def f(*x):
+        r = 1/x[0]
+        for i in range(len(x)-1):
+            r *= 1/x[i+1]
+        return r
+        
+    for j in range(n):
+        k = rand.randint(10) + 2
+        x = make_gummy(unit=1,allowzero=False,allowlargeu=False)[0]
+        a = f(*(k*[x]))
+        b = x**-k
+        d = uc.gummy.napply(f,*(k*[x]))
+        
+        assert abs((a.x - b.x)/b.x) < 1e-14
+        assert abs((a.u - b.u)/b.u) < 1e-3
+        assert a.correlation(b) == 1
+        
+        if k%2 == 0 and x.x < 0:
+            assert a.correlation(x) == 1
+        else:
+            assert a.correlation(x) == -1
         
+        assert abs((d.x - b.x)/b.x) < 1e-14
+        assert abs((d.u - b.u)/b.u) < 1e-3
+        assert d.correlation(b) == 1
 
 def test_ubreakdown(n=None,prnt=False,budget=False):
     
@@ -769,7 +893,7 @@ def d(k,*p):
                 
             gn = uc.gummy.napply(f,*g)
             assert abs(gn.x - xr)/abs(xr) < 1e-10
-            assert abs(gn.u - ur)/ur < 1e-5
+            assert abs(gn.u - ur)/ur < 1e3
             if dofr > 100:
                 assert gn.dof > 100
             else:
diff --git a/metrolopy/ummy.py b/metrolopy/ummy.py
index 1f29923..c02b6a3 100644
--- a/metrolopy/ummy.py
+++ b/metrolopy/ummy.py
@@ -104,12 +104,7 @@ def _combu(a,b,c):
         # maybe possibly get here due to a floating point rounding error somewhere
         return 0
     
-    x = np.sqrt(x)
-    try:
-        x = type(a)(x)
-    except:
-        pass
-    return abs(a)*x
+    return abs(a)*x**0.5
 
 def _isscalar(x):
     try:
@@ -176,6 +171,8 @@ def __init__(self,x,u=0,dof=float('inf'),utype=None):
             
         if isinstance(x,Fraction):
             x = MFraction(x)
+        if isinstance(u,Fraction):
+            u = MFraction(u)
             
         self._x = x
         self._u = u
@@ -358,11 +355,7 @@ def _set_covariance_matrix(gummys, matrix):
             else:
                 if gummys[i]._ref is None:
                     gummys[i]._ref = _GummyRef()  
-                e = np.sqrt(m[i][i])
-                try:
-                    e = type(m[i][i])(e)
-                except:
-                    e = float(e)
+                e = m[i][i]**0.5
                 gummys[i]._u = e
         
         # Set the correlations
@@ -606,29 +599,41 @@ def _format_mantissa(self,fmt,x,sig,parenth=False):
         return s + ellipses
         
     @classmethod
-    def _get_dof(cls,dof1, dof2, u, u1, u2, d1, d2, c):
+    def _get_dof(cls,dof1, dof2, u, du1, du2,c):
         # Use the Welch-Satterhwaite approximation to combine dof1 and dof2.
         # See [R. Willink, Metrologia, 44, 340 (2007)] concerning the use
         # of the Welch-Satterwaite approximation with correlated values.
         # Note that we allow the return value to be less than one,
         # but when the dof is retrieved using the dof property get method, 
         # values less than one will be rounded up to one.
+        
         if (isinf(dof1) and isinf(dof2)) or u == 0:
             return float('inf')
-        if isinf(dof1) or isinf(dof2):
-            xt = 0
-        else:
-            xt = (d1*d2*c*(u1/u)*(u2/u))**2
-        d = 0
-        if isfinite(dof1):
-            d += ((d1*u1/u)**4 + xt)/dof1
-        if isfinite(dof2):
-            d += ((d2*u2/u)**4 + xt)/dof2
+            
+        du1 = du1/u
+        du2 = du2/u
+            
+        # Willink's formula, [R. Willink, Metrologia, 44, 340 (2007)], for
+        # the W-S approximations for correlated uncertainties is:
+        # d = (du1**2 + c*du1*du2)**2/dof1 + (du2**2 + c*du1*du2)**2/dof2
+        # Willink's formula is used in the _apply method but not here.
+        
+        d = du1**4/dof1 + du2**4/dof2
+        d += 2*c**2*(((du1 + du2)**4 - du1**4 - du2**4)/
+                     ((1 - 2*c + 2*c**2)*(dof1+dof2)))
+
         if d == 0:
             return float('inf')
+        
         r = 1/d
         if r > cls.max_dof:
             return float('inf')
+        
+        if r < -1e-6:
+            raise ValueError('dof is negative')
+        if r < 0:
+            r == 1
+        
         return r
     
     @classmethod
@@ -674,11 +679,7 @@ def _apply(cls,function,derivative,*args,fxdx=None):
                     raise FloatingPointError('u is sqrt of ' + str(mu**2*u))
                 u = 0
                 
-            u = np.sqrt(u)
-            try:
-                u = type(fx)(u)
-            except:
-                u = float(u)
+            u = u**0.5
             u = mu*u
             
         r = cls(fx,u)
@@ -686,10 +687,13 @@ def _apply(cls,function,derivative,*args,fxdx=None):
             return r
             
         dm = 0
+
         for i,a in enumerate(args):
             if isinstance(a,ummy):
                 if not isinf(a.dof):
-                    dm += np.sum((c.dot(du)[i]/u)**4)/a.dof
+                    # See [R. Willink, Metrologia, 44, 340 (2007)] for this
+                    # extension of the W-S approximation to correlated values.
+                    dm += np.sum(((c.dot(du)[i]/u)*(du[i]/u))**2)/a.dof
                 if a._ref is not None:
                     a._ref.combl(r,a._refs*c[i].dot(du)/u,a._refs*du[i]/u,args)
         if r._ref is not None:
@@ -745,7 +749,7 @@ def _add(self,b):
         self._ref.comb(r,self._refs*(b._u*c + self._u)/r._u,self._refs*self._u/r._u)
         b._ref.comb(r,b._refs*(self._u*c + b._u)/r._u,b._refs*b._u/r._u,self._ref)
 
-        r._dof = ummy._get_dof(self._dof, b._dof, r._u, self._u, b._u, 1, 1, c)
+        r._dof = self._get_dof(self._dof, b._dof, r._u, self._u, b._u,c)
         
         return r
     
@@ -786,7 +790,7 @@ def _sub(self,b):
         self._ref.comb(r,self._refs*(-b._u*c + self._u)/r._u,self._refs*self._u/r._u)
         b._ref.comb(r,b._refs*(self._u*c - b._u)/r._u,-b._refs*b._u/r._u,self._ref)
 
-        r._dof = ummy._get_dof(self._dof, b._dof, r._u, self._u, b._u, 1, -1, c)
+        r._dof = self._get_dof(self._dof, b._dof, r._u, self._u, -b._u, c)
         return r
     
     def _rsub(self,b):
@@ -828,7 +832,7 @@ def _mul(self,b):
         b._ref.comb(r,b._refs*(b._x*self._u*c + self._x*b._u)/r._u,
                         b._refs*self._x*b._u/r._u,self._ref)
 
-        r._dof = ummy._get_dof(self._dof,b._dof,r._u,self._u,b._u,b._x,self._x,c)
+        r._dof = self._get_dof(self._dof,b._dof,r._u,b._x*self._u,self._x*b._u,c)
         return r
     
     def _rmul(self,b):
@@ -863,7 +867,7 @@ def _truediv(self,b):
             else:
                 x = self._x/b._x
             
-        u = _combu(self._u/b._x,self._x*b._u/b._x**2,c)
+        u = _combu(self._u/b._x,-self._x*b._u/b._x**2,c)
 
         r = type(b)(x,u)
         
@@ -885,8 +889,8 @@ def _truediv(self,b):
         b._ref.comb(r,b._refs*(self._u*c/b._x - self._x*b._u/b._x**2)/r._u,
                         b._refs*(-self._x*b._u/(b._x**2*r._u)),self._ref)
 
-        r._dof = ummy._get_dof(self._dof, b._dof, r._u, self._u, b._u, 1/b._x,
-                               -self._x/b._x**2, c)
+        r._dof = self._get_dof(self._dof, b._dof, r._u, self._u/b._x, 
+                               -b._u*self._x/b._x**2, c)
         return r
     
     def _rtruediv(self,b):
@@ -974,7 +978,7 @@ def _pow(self,b):
         b._ref.comb(r,b._refs*(da*self._u*c + db*b._u)/r._u,b._refs*db*b._u/r._u,
                         self._ref)
 
-        r._dof = ummy._get_dof(self._dof, b._dof, r._u, self._u, b._u, da, db, c)
+        r._dof = self._get_dof(self._dof, b._dof, r._u, da*self._u, db*b._u, c)
         return r
     
     def _rpow(self,b):
@@ -1031,6 +1035,21 @@ def __abs__(self):
             if self._ref is not None:
                 r._refs = -self._refs
         return r
+    
+    def __eq__(self,v):
+        if self is v:
+            return True
+        if isinstance(v,ummy):
+            if self._u == 0 and v._u == 0:
+                return self._x == v._x
+            if self._ref is v._ref and self._refs == v._refs:
+                if self._x == v._x and self._u == v._u:
+                    return True
+            return False
+        elif self._u == 0:
+            return self._x == v
+        else:
+            return False
         
     def __float__(self):
         return float(self.x)
@@ -1065,7 +1084,7 @@ def utype(self):
 
         An arbitrary string value labeling the uncertainty type.
         """
-        if self._ref._tag is None:
+        if self._ref is None or self._ref._tag is None:
             return None
         return self._ref._tag.name
         
diff --git a/metrolopy/version.py b/metrolopy/version.py
index f3a5fc7..f766ddd 100644
--- a/metrolopy/version.py
+++ b/metrolopy/version.py
@@ -20,4 +20,4 @@
 # You should have received a copy of the GNU General Public License along with 
 # MetroloPy. If not, see <http://www.gnu.org/licenses/>.
 
-__version__ = '0.5.1'
\ No newline at end of file
+__version__ = '0.5.2'
\ No newline at end of file