ararpy 0.1.15__py3-none-any.whl → 0.1.18__py3-none-any.whl

ararpy/smp/style.py

@@ -58,6 +58,13 @@ def set_plot_style(smp: Sample):
                 setattr(figure, 'title', Plot.Text())
             setattr(getattr(figure, 'title'), 'text', f"{suffix} {getattr(figure, 'name', '')}")
 
+    age_unit = "Undefined"
+    try:
+        age_unit = str(smp.Info.preference['ageUnit']).capitalize()
+    except:
+        print(traceback.format_exc())
+    smp.AgeSpectraPlot.yaxis.title.text = f"Apparent Age ({age_unit})"
+
 
 def reset_plot_scale(smp: Sample, only_figure: str = None):
     """
@@ -184,9 +191,10 @@ def set_table_style(sample: Sample):
     for key, comp in basic.get_components(sample).items():
         if isinstance(comp, Table):
             comp.header = TABLEHEADER(index=int(comp.id))
-            comp.colcount = len(comp.header)
-            comp.coltypes = [{'type': 'numeric'}] * (comp.colcount)
-            textindexs = getattr(comp, 'textindexs', [0]) if hasattr(comp, 'textindexs') else [0]
-            for i in textindexs:
-                comp.coltypes[i] = {'type': 'text'}
+            comp.set_coltypes()
+            # comp.colcount = len(comp.header)
+            # comp.coltypes = [{'type': 'numeric'}] * (comp.colcount)
+            # text_indexes = getattr(comp, 'text_indexes', [0]) if hasattr(comp, 'text_indexes') else [0]
+            # for i in text_indexes:
+            #     comp.coltypes[i] = {'type': 'text'}
 

ararpy/smp/table.py

@@ -9,6 +9,8 @@
 #
 #
 """
+import ast
+import re
 import copy
 from .. import calc
 from . import (sample as samples, basic)
@@ -96,6 +98,12 @@ def update_handsontable(smp: Sample, data: list, id: str):
             'true': True, 'false': False, 'True': True, 'False': False, '1': True, '0': False, 'none': False,
         }
         return [bools_dict.get(str(col).lower(), False) for col in cols]
+
+    def _digitize_data(a):
+        # pattern = r'^[-+]?(\d+(\.\d*)?|\.\d+)([eE][-+]?\d+)?$'
+        # return [[ast.literal_eval(str(cell)) if re.fullmatch(pattern, str(cell)) else np.nan if str(cell) == "" else cell for cell in row] for row in a]
+        return a
+
     try:
         smp.SequenceName = data[0]
     except IndexError:
@@ -104,47 +112,37 @@ def update_handsontable(smp: Sample, data: list, id: str):
     update_all_table = False
     try:
         if data[1] != smp.SequenceValue:
+            update_all_table = True
             smp.SequenceValue = data[1]
     except IndexError:
         pass
-    else:
-        update_all_table = True
 
     if id == '1':  # 样品值
         data = _normalize_data(data, len(samples.SAMPLE_INTERCEPT_HEADERS), 2)
-        smp.SampleIntercept = data
+        smp.SampleIntercept = _digitize_data(data)
     elif id == '2':  # 本底值
         data = _normalize_data(data, len(samples.BLANK_INTERCEPT_HEADERS), 2)
-        smp.BlankIntercept = data
+        smp.BlankIntercept = _digitize_data(data)
     elif id == '3':  # 校正值
         data = _normalize_data(data, len(samples.CORRECTED_HEADERS), 2)
-        smp.CorrectedValues = data
+        smp.CorrectedValues = _digitize_data(data)
     elif id == '4':  # Degas table
         data = _normalize_data(data, len(samples.DEGAS_HEADERS), 2)
-        smp.DegasValues = data
+        smp.DegasValues = _digitize_data(data)
     elif id == '5':  # 发行表
         data = _normalize_data(data, len(samples.PUBLISH_TABLE_HEADERS), 2)
-        smp.PublishValues = data
+        smp.PublishValues = _digitize_data(data)
     elif id == '6':  # 年龄谱
         data = _normalize_data(data, len(samples.SPECTRUM_TABLE_HEADERS), 2)
-        smp.ApparentAgeValues = data
+        smp.ApparentAgeValues = _digitize_data(data)
     elif id == '7':  # 等时线
-        smp.IsochronMark = data[2]
+        smp.IsochronMark = _digitize_data(data)[2]
         data = _normalize_data(data, len(samples.ISOCHRON_TABLE_HEADERS), 3)
-        smp.IsochronValues = data
-        smp.SelectedSequence1 = [
-            i for i in range(len(smp.IsochronMark)) if str(smp.IsochronMark[i]) == "1"]
-        smp.SelectedSequence2 = [
-            i for i in range(len(smp.IsochronMark)) if str(smp.IsochronMark[i]) == "2"]
-        smp.UnselectedSequence = [
-            i for i in range(len(smp.IsochronMark)) if
-            i not in smp.SelectedSequence1 + smp.SelectedSequence2]
-        #
-        smp.Info.results.selection[0]['data'] = smp.SelectedSequence1
-        smp.Info.results.selection[1]['data'] = smp.SelectedSequence2
-        smp.Info.results.selection[2]['data'] = smp.UnselectedSequence
+        smp.IsochronValues = _digitize_data(data)
+        smp.sequence()
     elif id == '8':  # 总参数
         data = _normalize_data(data, len(samples.TOTAL_PARAMS_HEADERS), 2)
+        data = _digitize_data(data)
         data[101: 112] = [_strToBool(i) for i in data[101: 112]]
         smp.TotalParam = data
     else:

ararpy/{argon_diffusion_simulator → thermo}/__init__.py

@@ -9,4 +9,4 @@
 #
 # 
 """
-from . import main
+from . import atomic_level_random_walk as arw, basic

ararpy/thermo/arrhenius.py

@@ -0,0 +1,15 @@
+#!/usr/bin/env python
+# -*- coding: UTF-8 -*-
+"""
+# ==========================================
+# Copyright 2025 Yang 
+# ararpy - arrhenius
+# ==========================================
+#
+#
+# 
+"""
+
+from . import basic
+
+

ararpy/{argon_diffusion_simulator/main.py → thermo/atomic_level_random_walk.py}

@@ -3,7 +3,7 @@
 """
 # ==========================================
 # Copyright 2024 Yang
-# ArgonDiffusionRandomWalk - main
+# argon atomic_level_random_walk model - main
 # ==========================================
 #
 #

ararpy/thermo/basic.py

@@ -0,0 +1,306 @@
+#!/usr/bin/env python
+# -*- coding: UTF-8 -*-
+"""
+# ==========================================
+# Copyright 2025 Yang 
+# ararpy - basic
+# ==========================================
+#
+#
+# 
+"""
+import numpy as np
+from math import log
+
+SEC2YEAR = 3600 * 24 * 365.2425
+CAL2JOULE = 4.184
+CM2UM = 10000
+GAS_CONSTANT = 8.31446261815324  # in J / (mol K)
+
+
+def gammq(a, x):
+    if x < 0 or a <= 0:
+        raise ValueError("ERROR(GAMMQ): x < 0 or a <= 0")
+
+    if x < a + 1:
+        gamser, gln = gser(a, x)
+        return 1.0 - gamser
+    else:
+        gammcf, gln = gcf(a, x)
+        return gammcf
+
+def gser(a, x, itmax=100, eps=3e-7):
+    gln = gammln(a)
+    if x <= 0:
+        if x < 0:
+            raise ValueError("ERROR(GSER): x < 0")
+        return 0.0, gln
+
+    ap = a
+    sum_ = 1.0 / a
+    delta = sum_
+    for n in range(1, itmax + 1):
+        ap += 1
+        delta *= x / ap
+        sum_ += delta
+        if abs(delta) < abs(sum_) * eps:
+            break
+    else:
+        raise RuntimeError("ERROR(GSER): a too large, itmax too small")
+
+    gamser = sum_ * np.exp(-x + a * np.log(x) - gln)
+    return gamser, gln
+
+def gcf(a, x, itmax=100, eps=3e-7):
+    gln = gammln(a)
+    gold = 0.0
+    a0 = 1.0
+    a1 = x
+    b0 = 0.0
+    b1 = 1.0
+    fac = 1.0
+    for n in range(1, itmax + 1):
+        an = float(n)
+        ana = an - a
+        a0 = (a1 + a0 * ana) * fac
+        b0 = (b1 + b0 * ana) * fac
+        anf = an * fac
+        a1 = x * a0 + anf * a1
+        b1 = x * b0 + anf * b1
+        if a1 != 0:
+            fac = 1.0 / a1
+            g = b1 * fac
+            if abs((g - gold) / g) < eps:
+                return g * np.exp(-x + a * np.log(x) - gln), gln
+            gold = g
+    else:
+        raise RuntimeError("ERROR(GCF): a too large, itmax too small")
+
+def gammln(xx):
+    cof = np.array([76.18009173, -86.50532033, 24.01409822,
+                    -1.231739516, 0.00120858003, -0.00000536382])
+    stp = 2.50662827465
+    half = 0.5
+    one = 1.0
+    fpf = 5.5
+
+    x = xx - one
+    tmp = x + fpf
+    tmp = (x + half) * np.log(tmp) - tmp
+    ser = one
+    for j in range(6):
+        x += one
+        ser += cof[j] / x
+
+    return tmp + np.log(stp * ser)
+
+
+def fit(x, y, sigx, sigy, pho=None):
+    """
+    equivalent to York2 regression
+    y = a + bx
+
+    Parameters
+    ----------
+    x
+    y
+    sigx
+    sigy
+    pho
+
+    Returns
+    -------
+    a: intercept
+    b: slope
+    siga: intercept uncertainty
+    sigb: slope uncertainty
+    chi2: chi-square value
+    q: the incomplete gamma function value, should be large for the best fitting.
+
+    """
+
+    print(f"{x = }")
+    print(f"{y = }")
+    print(f"{sigx = }")
+    print(f"{sigy = }")
+    print(f"{pho = }")
+
+    ndata = len(x)
+    if pho is None:
+        pho = np.zeros(ndata)
+    nd = 20
+    imax = 20
+    xerr = 0.001
+    wt = []
+    a = 0
+    b = -1.
+    siga = 0
+    sigb = 0
+    chi2 = 0
+    q = 0
+    iter = 0
+
+    while iter <= imax:
+        sx = 0.
+        sy = 0.
+        st2 = 0.
+        st3 = 0.
+        ss = 0.
+        b0 = b
+
+        for i in range(ndata):
+            wt.append(0)
+            wt[i] = 1. / (sigy[i] ** 2 + b ** 2 * sigx[i] ** 2 - 2 * pho[i] * sigx[i] * sigy[i])
+            ss = ss + wt[i]
+            sx = sx + x[i] * wt[i]
+            sy = sy + y[i] * wt[i]
+
+            # print(f"{x[i] = }, {y[i] = }, {wt[i] = }")
+
+        sxoss = sx / ss
+        syoss = sy / ss
+
+        for i in range(ndata):
+            t1 = (x[i] - sxoss) * sigy[i] ** 2
+            t2 = (y[i] - syoss) * sigx[i] ** 2 * b
+            t3 = sigx[i] * sigy[i] * pho[i]
+            st2 = st2 + wt[i] ** 2 * (y[i] - syoss) * (t1 + t2 - t3 * (y[i] - syoss))
+            st3 = st3 + wt[i] ** 2 * (x[i] - sxoss) * (t1 + t2 - b * t3 * (x[i] - sxoss))
+
+        b = st2 / st3
+        iter = iter + 1
+
+        # print(f"{sxoss = }, {syoss = }, {b = }, {abs(b0 - b) = }")
+
+        if abs(b0 - b) > xerr:
+            continue
+
+        a = (syoss - sxoss * b)
+        # print(f"{a = }, {b = }")
+        sgt1 = 0.
+        sgt2 = 0.
+
+        for i in range(ndata):
+            sgt1 = sgt1 + wt[i] * (x[i] - sxoss) ** 2
+            sgt2 = sgt2 + wt[i] * x[i] ** 2
+
+        sigb = (1. / sgt1) ** 0.5
+        siga = sigb * (sgt2 / ss) ** 0.5
+        chi2 = 0.
+
+        for i in range(ndata):
+            chi2 = chi2 + wt[i] * (y[i] - a - b * x[i]) ** 2
+
+        q = gammq(0.5 * (ndata - 2), 0.5 * chi2)
+
+        if abs(b0 - b) <= xerr:
+            break
+
+    return a, b, siga, sigb, chi2, q
+
+
+def get_tc(da2, sda2, E, sE, cooling_rate=10, temp_in_celsius=True, temp=None, A: float = 27.0, R=None, pho=0):
+    """
+    E / (R Tc) = ln(A t D0/a2)
+    E: diffusion activation energy
+    R: gas constant, 8.314 J/(K*mol)
+    Tc: closure temperature
+    A: geometric constant, A is 55, 27, or 8.7 for volume diffusion from a sphere, cylinder or plane sheet respectively
+    tau: time constant, can be determined by
+        (1) a pre-defined cooling rate following t = -R Tc2 / E * dT / dt, where dT / dt is the average cooling rate of the local geological body
+        (2) the age of dated sample following t = R * age * E * (Tp - Tc), where Tp is the temperature at present that could be 300 K.
+    D0: the diffusion coefficient at the given temperature, which will be the frequency factor (the pre-exponential factor) whe the temperature is infinite high
+    a: the effective diffusion radius, the characteristic dimension of the system, or the domain size, or roughly the grain size
+
+
+    In either way of them, a equation with Tc in both right and left hand will be yielded. Then using iterative calculation, the Tc can be determined.
+
+    Parameters
+    -------
+    da2: D/a2 in 1/a, pear year
+    sda2:
+    E: diffusion activation energy
+    sE:
+    pho: covariance of da2 and E
+    temp: temperature of the given diffusion coefficient, in K
+    A: geometric constant
+    R: gas constant
+    cooling_rate: in degree/m.y., degrees per million year
+    temp_in_celsius: temperature in celsius, True, else in K, False
+
+    Returns
+    -------
+
+    """
+    R = 8.314 if R is None else R # in J/K mol
+    Tc = 600
+    sTc = 0
+    Tm = 99999999 if temp is None else temp
+    Tm = Tm + 273.15 if temp_in_celsius else Tm
+    cooling_rate = cooling_rate / 1000000
+
+    iter_num = 0
+    iter_diff = 100
+    while abs(iter_diff) > 0.01 and iter_num < 100:
+        tau = R * Tc ** 2 / (E * cooling_rate)
+        # new_Tc = (E/R) / log(A * tau * da2)
+        new_Tc = (R / E) * log(A * tau * da2) + 1 / Tm  #
+        d1 = cooling_rate / (A * da2 * Tc ** 2)
+        s1 = d1 ** 2 * sda2 ** 2  # da2
+        d2 = R / E ** 2 * (log(A * tau * da2) + 1)
+        s2 = d2 ** 2 * sE ** 2  # E
+        d3 = 2 * R / (E * Tc)
+        s3 = d3 ** 2 * sTc ** 2  # Tc
+        new_Tc = 1 / new_Tc
+        iter_diff = abs(new_Tc - Tc)
+        Tc = new_Tc
+        sTc = Tc ** 2 * np.sqrt(s1 + s2 + s3 + 2 * d1 * d2 * pho)
+        iter_num += 1
+        # print(f"Get Tc: {iter_num = }, {Tc = }, {tau = }, {Tm = }")
+    return Tc - 273.15 if temp_in_celsius else Tc, sTc
+
+
+if __name__ == "__main__":
+    """
+    tests
+    """
+    # 柴田贤, 钾长石
+    print("# 柴田贤, 钾长石")
+    print(get_tc(da2=5.6 * 3600 * 24 * 365.24, sda2=0, E=28.8 * 1000, R=1.987, sE=0, cooling_rate=5, A=8.7))
+    print(get_tc(da2=5.6 * 3600 * 24 * 365.24, sda2=0, E=28.8 * CAL2JOULE * 1000, sE=0, cooling_rate=5, A=8.7))
+    print(get_tc(da2=5.6 * 3600 * 24 * 365.24, sda2=0, E=28.8 * CAL2JOULE * 1000, sE=0, cooling_rate=30, A=8.7))
+    # 杨静, 斜钾铁矾 K101-1
+    print("# 杨静, 斜钾铁矾 K101-1")
+    print(get_tc(da2=10 ** 15.07 * 3600 * 24 * 365.24, sda2=0, E=76.17 * 1000 * CAL2JOULE, sE=0, cooling_rate=5, A=55))  # 这里跟文章中吻合
+    print(get_tc(da2=10 ** 15.07 * 3600 * 24 * 365.24, sda2=0, E=76.17 * 1000 * CAL2JOULE, sE=0, cooling_rate=10, A=55))  # 10冷却速率不吻合
+    # 杨静, 斜钾铁矾 K101-2
+    print("# 杨静, 斜钾铁矾 K101-2")
+    print(get_tc(da2=10 ** 12.35 * 3600 * 24 * 365.24, sda2=0, E=66.43 * 1000 * CAL2JOULE, sE=0, cooling_rate=5, A=55))  # 这里跟文章中吻合
+    print(get_tc(da2=10 ** 12.35 * 3600 * 24 * 365.24, sda2=0, E=66.43 * 1000 * CAL2JOULE, sE=0, cooling_rate=10, A=55))  # 10冷却速率不吻合
+    # 芦武长
+    print("# 芦武长")
+    print(get_tc(da2=3e-4 * 3600 * 24 * 365.24, sda2=0, E=35 * 1000 * CAL2JOULE, sE=0, cooling_rate=0.1, A=27, temp=600))  # 缺乏样品特征尺寸参数，da2不对，结果不对
+    print(get_tc(da2=3e-4 * 3600 * 24 * 365.24, sda2=0, E=35 * 1000 * CAL2JOULE, sE=0, cooling_rate=1, A=27, temp=600))  # 缺乏样品特征尺寸参数，da2不对，结果不对
+    # Harrison 1981, 角闪石, A = 55
+    print("# Harrison 1981, 角闪石, A = 55")
+    print(get_tc(da2=0.024 * SEC2YEAR / 0.0080 ** 2, sda2=0, E=64.1 * 1000 * CAL2JOULE, sE=0, cooling_rate=500, A=55))
+    print(get_tc(da2=0.024 * SEC2YEAR / 0.0080 ** 2, sda2=0, E=64.1 * 1000 * CAL2JOULE, sE=0, cooling_rate=5, A=55))
+    # Harrison 1985, 黑云母, A = 27
+    print("# Harrison 1985, 黑云母, A = 27")
+    print(get_tc(da2=0.077 * SEC2YEAR / 0.0150 ** 2, sda2=0, E=47.0 * 1000 * CAL2JOULE, sE=0, cooling_rate=100, A=27))
+    print(get_tc(da2=0.077 * SEC2YEAR / 0.0150 ** 2, sda2=0, E=47.0 * 1000 * CAL2JOULE, sE=0, cooling_rate=1, A=27))
+    print("# Blereau 2019, 大隅石，Osumilite, A = 8.7")
+    print(get_tc(da2=8.34e8 * SEC2YEAR / (0.0175 ** 2), sda2=0, E=461 * 1000, sE=0, cooling_rate=10, A=8.7))
+    print(get_tc(da2=6.70e8 * SEC2YEAR / (0.0175 ** 2), sda2=0, E=453 * 1000, sE=0, cooling_rate=10, A=8.7))  # 这两个与原文对不上，且原文描述表格与文字不符
+    print(get_tc(da2=4.49e5 * SEC2YEAR / (0.0175 ** 2), sda2=0, E=369 * 1000, sE=0, cooling_rate=10, A=8.7))  # 这两个与原文对不上，且原文描述表格与文字不符
+    print("# Thern 2020, 电气石, A = 55")  # 大部分都对或者接近，不完全一致的原因大概是他们采用蒙特卡洛方法计算且正负误差不一样
+    print(get_tc(da2=4.86e22 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=678 * 1000, sE=0, cooling_rate=10, A=55))
+    print(get_tc(da2=8.99e23 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=756 * 1000, sE=0, cooling_rate=10, A=55))
+    print(get_tc(da2=3.95e19 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=604 * 1000, sE=0, cooling_rate=10, A=55))
+    print(get_tc(da2=6.95e27 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=815 * 1000, sE=0, cooling_rate=10, A=55))
+    print(get_tc(da2=2.70e14 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=519 * 1000, sE=0, cooling_rate=10, A=55))
+    print(get_tc(da2=1.65e15 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=505 * 1000, sE=0, cooling_rate=10, A=55))
+    print(get_tc(da2=9.09e21 * SEC2YEAR / (0.0100 ** 2), sda2=0, E=656 * 1000, sE=0, cooling_rate=10, A=55))
+    print("# Shi 2020, 黑云母, A = 55")  # 大部分都对或者接近，不完全一致的原因大概是他们采用蒙特卡洛方法计算且正负误差不一样
+    print(get_tc(da2=10 ** 18 * SEC2YEAR, sda2=0, E=87 * 1000 * CAL2JOULE, sE=0, cooling_rate=10, A=6))  # 文中的参数可能是cooling_rate=100, A=8.7，但是10和6的组合与文中报导的结果一致
+    print(get_tc(da2=10 ** 16.6 * SEC2YEAR, sda2=0, E=81.5 * 1000 * CAL2JOULE, sE=0, cooling_rate=10, A=6))

{ararpy-0.1.15.dist-info → ararpy-0.1.18.dist-info}/METADATA

@@ -1,13 +1,13 @@
-Metadata-Version: 2.2
+Metadata-Version: 2.4
 Name: ararpy
-Version: 0.1.15
+Version: 0.1.18
 Summary: A project for Ar-Ar geochronology
 Home-page: https://github.com/wuyangchn/ararpy.git
 Author: Yang Wu
 Author-email: wuycug@hotmail.com
+License: MIT
 Classifier: Development Status :: 3 - Alpha
 Classifier: Programming Language :: Python :: 3
-Classifier: License :: OSI Approved :: MIT License
 Classifier: Operating System :: OS Independent
 Requires-Python: >=3.5
 Description-Content-Type: text/markdown
@@ -18,6 +18,8 @@ Dynamic: classifier
 Dynamic: description
 Dynamic: description-content-type
 Dynamic: home-page
+Dynamic: license
+Dynamic: license-file
 Dynamic: requires-python
 Dynamic: summary
 

{ararpy-0.1.15.dist-info → ararpy-0.1.18.dist-info}/RECORD

@@ -1,19 +1,17 @@
-ararpy/__init__.py,sha256=XKlsgklkzGfvqQkt59y-UST9aVhgCf-eXpEDiaZaynQ,6725
+ararpy/__init__.py,sha256=abmIxFaRq9m7mUI66xHbMGD2wRdzV9v8FAgvlViIlTc,6733
 ararpy/test.py,sha256=4F46-JJ1Ge12HGae0qO44Qc6kiEMHBgn2MsY_5LlHDo,3973
-ararpy/argon_diffusion_simulator/__init__.py,sha256=8VEdUP7sh17lNg6FcRVVK0JynthEsSk_qcTQtwisqU8,229
-ararpy/argon_diffusion_simulator/main.py,sha256=-EN7ndhzVslNctJSaQTn_JwumA6Cnt37yiPxuLtahXQ,25490
 ararpy/calc/__init__.py,sha256=kUjRuLE8TLuKOv3i976RnGJoEMj23QBZDu37LWs81U4,322
-ararpy/calc/age.py,sha256=zaiy_2AgKEe1a-KfZoi6HmrLpyU_aLUuNVmSeJ_X5KQ,5820
+ararpy/calc/age.py,sha256=POBF_UxfXM_GCfXR_sTE__-fOvYZleNXXSuR_SIrFrI,5782
 ararpy/calc/arr.py,sha256=w31bn6MXF8I8qPXYo5kI-TfMKCYspx1rZ5g_UpwwSd8,14939
-ararpy/calc/basic.py,sha256=QN__j0JrEZA-W29iPAnwT7XtR6PCuvgDptLtyTgGp_k,3236
-ararpy/calc/corr.py,sha256=i0ku_L0dV4Cex0SLZxyDHYML8lF_njA0yfIufX8DXXo,17926
+ararpy/calc/basic.py,sha256=uJCCUFaPd9zvfkggrdbFYSGLl2pt7UJ7ENgXanzHy68,4036
+ararpy/calc/corr.py,sha256=HnCbf_84gKD6vMX7PflZym-I1MezRKfNHYfEtkDCCaI,18598
 ararpy/calc/err.py,sha256=63LtprqjemlIb1QGDst4Ggcv5KMSDHdlAIL-nyQs1eA,2691
 ararpy/calc/histogram.py,sha256=0GVbDdsjd91KQ1sa2B7NtZ4KGo0XpRIJapgIrzAwQUo,5777
 ararpy/calc/isochron.py,sha256=ej9G2e68k6yszonWHsLcEubh3TA7eh1upTJP_X0ttAA,5726
 ararpy/calc/jvalue.py,sha256=zHUhJ1iYe5mPrY95mGYxoUPAp7hwu4coUgiHKiruKfM,1138
 ararpy/calc/plot.py,sha256=iWxPYIaT0OmJBiho7X9wFY7m-QrvYexu6nWHp24T_VY,1962
 ararpy/calc/raw_funcs.py,sha256=UC01lvA6GyZ5FJv43jgoUULAFoLnZJMxeSa0BeVFCAM,2637
-ararpy/calc/regression.py,sha256=Xes0rvIX5a9e_gje3n6iTXYbNtRuZZMKRcADp2b-5Dw,40091
+ararpy/calc/regression.py,sha256=w5kni6LGqvISvlvbBZnJ3N2c5eQjgkz3bBbj0PXPyGs,40251
 ararpy/calc/spectra.py,sha256=_Q23eP9necHlaCoHf3_UfW1N3JmVZj5rcWFro8GS-CA,1995
 ararpy/examples/022_VU124-M11a.ahd,sha256=3m0Gd-ZObou3KsnRNFMf77QwzT1Uz3nu3vA33Sqeyng,5414
 ararpy/examples/20WHA0103.age,sha256=cT-a4d7Wt77aotx6v0G47vulY_TZIcZUcaVHB3pqTPM,380416
@@ -37,28 +35,32 @@ ararpy/examples/sample-default.smp,sha256=YNkoQGgPrsL_fXS7ZHxfRtLQWekCDqT9czS6vB
 ararpy/files/__init__.py,sha256=l5B5ZQ01WdtvjjN0aMkyAFNgpwANdM_1I0tQbqnRuEY,69
 ararpy/files/arr_file.py,sha256=KqksGlEA6nmMQofTgi7v45flscQZVtefxaNCKrV3Am4,837
 ararpy/files/basic.py,sha256=nc7Hgo_qLSkdmtKzZmd5SQ8Jy0dhW46ly4gh-oisUDs,2095
-ararpy/files/calc_file.py,sha256=ezKtP2PFSW5W2w6cPY3p4l9VJ4uUPAjc5guHeZxxERA,28303
+ararpy/files/calc_file.py,sha256=Mww0R1RH3BpTCo6EetW9jikd19G6-Ljhgak5PXccRfA,28341
 ararpy/files/new_file.py,sha256=efblARIBROVLWS2w3-98BxLX5VZ8grRpiTkJFtf_rAk,214
 ararpy/files/raw_file.py,sha256=aVHlebiGpxHGj_z9mrTR-Y2JEivxa4k1MJZ1TFXyMnA,22319
 ararpy/files/xls.py,sha256=8ibT4ZRY2grK34ikKm1pDmlWFlVTgDL7vSMO6mphzrY,702
 ararpy/smp/EXPORT_TO_PDF_DATA_PROPERTIES.py,sha256=baDM437tu6hsPv0uYfod0TREXlPd6kvMBFT1S9ZZlkk,3024
 ararpy/smp/__init__.py,sha256=k6_fa27UJsQK7K7oC5GYlwMo6l0Xd8af3QtOrZz2XJk,478
-ararpy/smp/basic.py,sha256=Jufv_4zXO1-XbmwM4UUrvZsp4tRGD05ec-035u8CiEs,21533
+ararpy/smp/basic.py,sha256=1zsBznVojY-SFAjq_1tkHxvpI9-mD_uotcrd9aDd6n4,22923
 ararpy/smp/calculation.py,sha256=LCFJWjLVLEKEQ5b7RFUIxsMahEzgLdodW4kCYXV5Z34,2919
 ararpy/smp/consts.py,sha256=XIdjdz8cYxspG2jMnoItdlUsxr3hKbNFJjMZJh1bpzw,393
-ararpy/smp/corr.py,sha256=ga3OsRSwPeLd4VV7DbOjZ3T4qlOYweUZhFfAM9MkM5s,25484
-ararpy/smp/diffusion_funcs.py,sha256=_8SchTD27CAS8SmQVDoXJscTRB7M-BqtpASEgg_izBM,175748
-ararpy/smp/export.py,sha256=4RtoDvItQ2ysWXY78gS_tDocTgSY9ftP9Vb-YKrOGQI,107646
+ararpy/smp/corr.py,sha256=1MDH41TsBVyuyI45mqaqP2IvEm_jAUMZgR0MG7uNiu4,25598
+ararpy/smp/diffusion_funcs.py,sha256=4-PMMIZWzjk2HOYYWNgSp4GmApygp1MmOxJ2g3xrqWc,175049
+ararpy/smp/export.py,sha256=_kanGNHjCvjzwUS1xNAPTaTk7VdIn88wOZ9mmq2eLPQ,117324
 ararpy/smp/info.py,sha256=iKUELm-BuUduDlJKC1d8tKKNHbwwbNmhUg2pi6bcBvA,489
-ararpy/smp/initial.py,sha256=J30NYDpccRJ1TSP6Fl_nHSZ12J6PDMtIic-iT5fSlC4,15580
+ararpy/smp/initial.py,sha256=eWvrzTh80AigxDI_BQn6dLkaf-oJeam24LFXcz8zB0o,16011
 ararpy/smp/json.py,sha256=BTZCjVN0aj9epc700nwkYEYMKN2lHBYo-pLmtnz5oHY,2300
-ararpy/smp/plots.py,sha256=jYvi4oNEdI1r4371Wh3cB6vkL-1rrzgbyClK87Uw6YM,31617
+ararpy/smp/plots.py,sha256=rb_Hsk2M3Dftbsj4bT6JIdSZNjzxyllnMgXhlVChNjU,31802
 ararpy/smp/raw.py,sha256=51n-rrbW2FqeZHQyevuG7iObPLGvIBzTe414QDVM1FE,6523
-ararpy/smp/sample.py,sha256=2aqyszY81GXzbZWRKip9BEaeYaFV32-gxmhue3gCvfs,55429
-ararpy/smp/style.py,sha256=NZwTeodhtHBBdCyjvupF64nh2B5EDdMHzlT8oE360zg,6746
-ararpy/smp/table.py,sha256=yNr4PLxbH4ezeUXpKBxLjRfq4RuYiW-G3sMWFDtjlVM,6802
-ararpy-0.1.15.dist-info/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
-ararpy-0.1.15.dist-info/METADATA,sha256=_zC9lSFdJP2rU_a6zTZlmrz231SNMTGtLR0iMHnLR1E,24513
-ararpy-0.1.15.dist-info/WHEEL,sha256=jB7zZ3N9hIM9adW7qlTAyycLYW9npaWKLRzaoVcLKcM,91
-ararpy-0.1.15.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
-ararpy-0.1.15.dist-info/RECORD,,
+ararpy/smp/sample.py,sha256=qrsILB0O2v_8HAq_XsWTyF_rEAGs8GYOVVGz0Cb-Qc0,55641
+ararpy/smp/style.py,sha256=Z4exOAdtvw5xwDf6nyIjPxuvww2U5JK6H_Ql02DBSyM,7030
+ararpy/smp/table.py,sha256=9bNAOqAIOc0nSC3LNeqjJKUYSJSM28Ji3o9VimwMU8A,6645
+ararpy/thermo/__init__.py,sha256=6VBuqTRFl403PVqOuMkVrut0nKaQsAosBmfW91X1dMg,263
+ararpy/thermo/arrhenius.py,sha256=Ass1ichHfqIAtpv8eLlgrUc1UOb3Urh1qzr1E3gLB4U,233
+ararpy/thermo/atomic_level_random_walk.py,sha256=Q97zfe2h2RaxADkoBAqd0uEiP16BFOajrTmXHMkL2EQ,25502
+ararpy/thermo/basic.py,sha256=nBGHI9uK7VdJwThwBIOcKAzdnYqPyQseFoY6s4zKizk,11504
+ararpy-0.1.18.dist-info/licenses/LICENSE,sha256=cvG5t_C1qY_zUyJI7sNOa7gCArdngNPaOrfujl2LYuc,1085
+ararpy-0.1.18.dist-info/METADATA,sha256=NJXSnrB9-l5wMxHlvWdB9yA11EYO_G2co1YjZStm4Tk,24516
+ararpy-0.1.18.dist-info/WHEEL,sha256=1tXe9gY0PYatrMPMDd6jXqjfpz_B-Wqm32CPfRC58XU,91
+ararpy-0.1.18.dist-info/top_level.txt,sha256=9iTpsPCYuRYq09yQTk9d2lqB8JtTEOmbN-IcGB-K3vY,7
+ararpy-0.1.18.dist-info/RECORD,,

{ararpy-0.1.15.dist-info → ararpy-0.1.18.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: setuptools (75.8.2)
+Generator: setuptools (77.0.3)
 Root-Is-Purelib: true
 Tag: py3-none-any