PyTransit 2.9.2__py3-none-any.whl → 2.10.0__py3-none-any.whl

pytransit/__init__.py

@@ -41,7 +41,10 @@ Date
 
 """
 
-__version__ = '2.9.2'
+try:
+    from ._version import version as __version__
+except ImportError:
+    __version__ = '0+unknown'
 
 # Generic
 # -------

pytransit/_version.py

@@ -0,0 +1,24 @@
+# file generated by vcs-versioning
+# don't change, don't track in version control
+from __future__ import annotations
+
+__all__ = [
+    "__version__",
+    "__version_tuple__",
+    "version",
+    "version_tuple",
+    "__commit_id__",
+    "commit_id",
+]
+
+version: str
+__version__: str
+__version_tuple__: tuple[int | str, ...]
+version_tuple: tuple[int | str, ...]
+commit_id: str | None
+__commit_id__: str | None
+
+__version__ = version = '2.10.0'
+__version_tuple__ = version_tuple = (2, 10, 0)
+
+__commit_id__ = commit_id = 'gd00e4aa7f'

pytransit/lpf/lpf.py

@@ -505,7 +505,7 @@ class BaseLPF(LogPosteriorFunction):
         self._additional_log_priors.append(T14)
 
     def add_as_prior(self, mean: float, std: float) -> None:
-        """Add a normal prior on the scaled semi-major axis :math:`(a / R_\star)`.
+        r"""Add a normal prior on the scaled semi-major axis :math:`(a / R_\star)`.
 
         Parameters
         ----------

pytransit/models/roadrunner/common.py

@@ -2,7 +2,7 @@ from numba import njit
 from typing import Optional
 
 from numpy import (sqrt, sin, arctan2, pi, nan, zeros, floor, arccos, linspace, ndarray, atleast_1d,
-                   atleast_2d, asarray, full, empty, concatenate, array, broadcast_to)
+                   atleast_2d, asarray, full, empty, concatenate, array, broadcast_to, int64)
 from numpy.polynomial import Polynomial
 from scipy.special import roots_jacobi, roots_legendre
 
@@ -578,6 +578,72 @@ def ldm_table(mu: ndarray, wf: ndarray, t0: float, dt: float, ldp: ndarray, out:
         out[ig] = num / den
 
 
+@njit(fastmath=True)
+def ldm_table_map(mu: ndarray, wf: ndarray, t0: float, dt: float, nmu: int):
+    """The linear map from a tabulated profile to the mean intensity table.
+
+    `ldm_table` is linear in the profile, and its nodes, weights, and interpolation stencils depend
+    only on the planet, so the map can be built once and applied to any number of profiles with
+    `ldm_table_banded`. This is what lets the spectroscopy model share the work between its
+    wavelength bins.
+
+    Each table node only reads the profile across the planet, so its row of the map is nonzero in a
+    narrow band. Returns the ``(ng, nmu)`` map and, for every row, the first and one-past-the-last
+    nonzero column.
+    """
+    ng, nn = mu.shape
+    wt = zeros((ng, nmu))
+    lo = empty(ng, dtype=int64)
+    hi = empty(ng, dtype=int64)
+    for ig in range(ng):
+        den = 0.0
+        for q in range(nn):
+            den += wf[ig, q]
+        for q in range(nn):
+            w = wf[ig, q]
+            if w == 0.0:
+                continue
+            w /= den
+            x = (sqrt(mu[ig, q]) - t0) / dt
+            if x <= 0.0:
+                wt[ig, 0] += w
+            elif x >= nmu - 1:
+                wt[ig, nmu - 1] += w
+            else:
+                i = int(x) - 1
+                if i < 0:
+                    i = 0
+                elif i > nmu - 4:
+                    i = nmu - 4
+                u = x - i
+                wt[ig, i] -= w * (u - 1) * (u - 2) * (u - 3) / 6.0
+                wt[ig, i + 1] += w * u * (u - 2) * (u - 3) / 2.0
+                wt[ig, i + 2] -= w * u * (u - 1) * (u - 3) / 2.0
+                wt[ig, i + 3] += w * u * (u - 1) * (u - 2) / 6.0
+
+        lo[ig] = 0
+        while lo[ig] < nmu - 1 and wt[ig, lo[ig]] == 0.0:
+            lo[ig] += 1
+        hi[ig] = nmu
+        while hi[ig] > lo[ig] + 1 and wt[ig, hi[ig] - 1] == 0.0:
+            hi[ig] -= 1
+    return wt, lo, hi
+
+
+@njit(fastmath=True)
+def ldm_table_banded(wt: ndarray, lo: ndarray, hi: ndarray, ldp: ndarray, scale: float, out: ndarray) -> None:
+    """The mean intensity table of one profile from the banded map of `ldm_table_map`, times `scale`."""
+    for ig in range(wt.shape[0]):
+        # Slicing the band first lets the loop index start from zero, which numba needs to drop its
+        # negative index check and vectorise the sum.
+        w = wt[ig, lo[ig]:hi[ig]]
+        p = ldp[lo[ig]:hi[ig]]
+        s = 0.0
+        for j in range(w.size):
+            s += w[j] * p[j]
+        out[ig] = scale * s
+
+
 def _cubic_matrices() -> ndarray:
     """Maps from four stencil values to the cubic's coefficients in the local interval variable.
 
@@ -608,8 +674,15 @@ def cubic_coefficients(tab: ndarray) -> ndarray:
     Lagrange interpolant on a stencil shifted inwards at the ends, so the interpolation is fourth
     order everywhere including the first and last intervals.
     """
+    c = zeros((tab.size - 1, 4))
+    cubic_coefficients_into(tab, c)
+    return c
+
+
+@njit(fastmath=True)
+def cubic_coefficients_into(tab: ndarray, out: ndarray) -> None:
+    """`cubic_coefficients` written into `out`, which has ``tab.size - 1`` rows."""
     n = tab.size
-    c = zeros((n - 1, 4))
     for i in range(n - 1):
         j = i - 1
         if j < 0:
@@ -618,17 +691,16 @@ def cubic_coefficients(tab: ndarray) -> ndarray:
             j = n - 4
         s = i - j
         for r in range(4):
-            c[i, r] = (CUBIC_MATRICES[s, r, 0] * tab[j] + CUBIC_MATRICES[s, r, 1] * tab[j + 1]
-                       + CUBIC_MATRICES[s, r, 2] * tab[j + 2] + CUBIC_MATRICES[s, r, 3] * tab[j + 3])
-    return c
+            out[i, r] = (CUBIC_MATRICES[s, r, 0] * tab[j] + CUBIC_MATRICES[s, r, 1] * tab[j + 1]
+                         + CUBIC_MATRICES[s, r, 2] * tab[j + 2] + CUBIC_MATRICES[s, r, 3] * tab[j + 3])
 
 
 @njit(fastmath=True)
-def ldm_lookup(g: float, gc: float, n1: int, coef: ndarray) -> float:
-    """The mean intensity under the planet at grazing parameter g from the split cubic table.
+def ldm_interval(g: float, gc: float, n1: int, nc: int):
+    """The interval of the split cubic table holding grazing parameter g, and g's position in it.
 
-    `coef` holds the interval coefficients of both segments back to back: ``n1 - 1`` intervals
-    covering [0, gc] and the rest covering [gc, 1].
+    `nc` is the number of intervals in the table. Returns the interval index and the local
+    position in [0, 1) that `ldm_lookup` evaluates the interval's cubic at.
     """
     if g < gc:
         x = g / gc * (n1 - 1)
@@ -636,14 +708,24 @@ def ldm_lookup(g: float, gc: float, n1: int, coef: ndarray) -> float:
         if i > n1 - 2:
             i = n1 - 2
     else:
-        n2 = coef.shape[0] - (n1 - 1) + 1
+        n2 = nc - (n1 - 1) + 1
         x = (g - gc) / (1.0 - gc) * (n2 - 1)
         i = int(x)
         if i > n2 - 2:
             i = n2 - 2
         x += n1 - 1
         i += n1 - 1
-    a = x - i
+    return i, x - i
+
+
+@njit(fastmath=True)
+def ldm_lookup(g: float, gc: float, n1: int, coef: ndarray) -> float:
+    """The mean intensity under the planet at grazing parameter g from the split cubic table.
+
+    `coef` holds the interval coefficients of both segments back to back: ``n1 - 1`` intervals
+    covering [0, gc] and the rest covering [gc, 1].
+    """
+    i, a = ldm_interval(g, gc, n1, coef.shape[0])
     return coef[i, 0] + a * (coef[i, 1] + a * (coef[i, 2] + a * coef[i, 3]))
 
 
@@ -657,5 +739,5 @@ def split_cubic_coefficients(tab: ndarray, n1: int, out: ndarray) -> None:
     """
     # The split node is stored twice, once as the last node of the first segment and once as
     # the first node of the second, so the second segment starts at index n1.
-    out[:n1 - 1] = cubic_coefficients(tab[:n1])
-    out[n1 - 1:] = cubic_coefficients(tab[n1:])
+    cubic_coefficients_into(tab[:n1], out[:n1 - 1])
+    cubic_coefficients_into(tab[n1:], out[n1 - 1:])

pytransit/models/roadrunner/model_trspec.py

@@ -1,19 +1,68 @@
-from math import fabs, floor
+from math import floor
 
 from meepmeep.backends.numba.point2d import sep_c, solve2d, bounding_box
 from numba import njit, prange, get_num_threads, set_num_threads
-from numpy import zeros, dot, ndarray, isnan, nan, mean, floor, fabs, max, empty
+from numpy import zeros, ndarray, isnan, nan, mean, max, empty, int64
 
-from .common import (g_nodes, ldm_nodes, ldm_table, split_cubic_coefficients, split_point, ldm_lookup,
-                     valid_radius_ratios)
+from .common import (g_nodes, ldm_nodes, ldm_table_map, ldm_table_banded, split_cubic_coefficients, split_point,
+                     ldm_interval, valid_radius_ratios)
 from .common import circle_circle_intersection_area_kite as ccia
 
 
-@njit(parallel=False, fastmath=False)
-def tsmodel_serial(times: ndarray,
-                   k: ndarray, t0: ndarray, p: ndarray, a: ndarray, i: ndarray, e: ndarray, w: ndarray,
-                   nsamples: ndarray, exptimes: ndarray, ldp: ndarray, istar: ndarray,
-                   pt0: float, pdt: float, rules: ndarray, ng: int) -> ndarray:
+@njit(fastmath=False)
+def transit_samples(times: ndarray, t0: float, p: float, xyc: ndarray, k: float, kmax: float,
+                    gc: float, n1: int, nc: int, nsamples: int, exptime: float,
+                    interval: ndarray, ipts: ndarray, basis: ndarray, area: ndarray) -> int:
+    """The geometry of every exposure sample that overlaps the star, shared by all the wavelength bins.
+
+    For each overlapping sample, stores its exposure index, the mean intensity table interval it
+    falls in, the cubic basis ``(a, a**2, a**3)`` at its position in that interval, and the three
+    overlap area terms the bins combine as ``area[0] + area[1] * (k / k_mean)**2 + area[2] * (k - k_mean)``.
+    The area terms are divided by the number of samples, so the bins can sum the samples directly.
+
+    Returns the number of overlapping samples.
+    """
+    bt1, bt4 = bounding_box(k, xyc)
+    bt1 -= 0.0015 + exptime
+    bt4 += 0.0015 + exptime
+
+    n = 0
+    for ipt in range(times.size):
+        epoch = floor((times[ipt] - t0 + 0.5 * p) / p)
+        tc = times[ipt] - (t0 + epoch * p)
+        if not (bt1 <= tc <= bt4):
+            continue
+        for isample in range(1, nsamples + 1):
+            z = sep_c(tc + exptime * ((isample - 0.5) / nsamples - 0.5), xyc)
+            ap0, kappa = ccia(1.0, k, z)
+            dadk = 2.0 * k * kappa
+            if ap0 == 0.0 and dadk == 0.0:
+                continue
+            i, a = ldm_interval(z / (1.0 + k), gc, n1, nc)
+            ipts[n] = ipt
+            interval[n] = i
+            basis[n, 0] = a
+            basis[n, 1] = a * a
+            basis[n, 2] = a * a * a
+
+            # The area is exact while the planet is fully on the disk, and linear in k - k_mean across the limb.
+            if z <= 1.0 - kmax:
+                area[n, 0] = 0.0
+                area[n, 1] = ap0 / nsamples
+                area[n, 2] = 0.0
+            else:
+                area[n, 0] = ap0 / nsamples
+                area[n, 1] = 0.0
+                area[n, 2] = dadk / nsamples
+            n += 1
+    return n
+
+
+def _tsmodel(times: ndarray,
+             k: ndarray, t0: ndarray, p: ndarray, a: ndarray, i: ndarray, e: ndarray, w: ndarray,
+             nsamples: ndarray, exptimes: ndarray, ldp: ndarray, istar: ndarray,
+             pt0: float, pdt: float, rules: ndarray, ng: int) -> ndarray:
+    """The transmission spectroscopy model. Compiled both serially and in parallel over the wavelength bins."""
     if k.ndim != 2:
         raise ValueError(" The radius ratios must be given as a 2D array with shape (npv, npb)")
 
@@ -26,15 +75,22 @@ def tsmodel_serial(times: ndarray,
     npt = times.size
     npv = k.shape[0]
     npb = k.shape[1]
+    nmu = ldp.shape[2]
+    ns = nsamples[0]
+    nc = ng - 2                      # Number of intervals in the split mean intensity table
 
-
-    flux = zeros((npv, npb, npt))  # Model flux
+    flux = empty((npv, npb, npt))  # Model flux
     nq = rules.shape[2]
-    coef = zeros((npb, ng - 2, 4))   # Split cubic coefficients of the mean intensity tables
+    tables = empty((npb, ng))      # Mean intensity tables
+    coef = empty((npb, nc, 4))     # Split cubic coefficients of the mean intensity tables
     mu = empty((ng, 2 * nq))
     wf = zeros((ng, 2 * nq))
-    ldm = empty(ng)
-    xyc = zeros((2, 5))            # Taylor series coefficients for the (x, y) position
+
+    # Geometry of the overlapping exposure samples
+    ipts = empty(npt * ns, dtype=int64)
+    interval = empty(npt * ns, dtype=int64)
+    basis = empty((npt * ns, 3))
+    area = empty((npt * ns, 3))
 
     for ipv in range(npv):
         if isnan(a[ipv]) or (a[ipv] <= 1.0) or (e[ipv] < 0.0) or not valid_radius_ratios(k[ipv]):
@@ -44,141 +100,55 @@ def tsmodel_serial(times: ndarray,
         kmean = mean(k[ipv])
         kmax = max(k[ipv])
         afac = k[ipv] ** 2 / kmean ** 2
+        dk = k[ipv] - kmean
 
         # -----------------------------------#
         # Calculate the limb darkening means #
         # -----------------------------------#
+        # The table is linear in the profile, so the map from a profile to its table is built once
+        # and applied to every bin. Dividing by the stellar intensity here saves a division per bin
+        # and exposure later.
         gs, n1 = g_nodes(kmean, ng)
         gc = split_point(kmean)
         ldm_nodes(kmean, gs, rules, mu, wf)
-        for ipb in range(npb):
-            ldm_table(mu, wf, pt0, pdt, ldp[ipv, ipb, :], ldm)
-            split_cubic_coefficients(ldm, n1, coef[ipb])
-
-        # -----------------------------------------------------#
-        # Calculate the Taylor series expansions for the orbit #
-        # -----------------------------------------------------#
-        xyc[:, :] = solve2d(0.0, p[ipv], a[ipv], i[ipv], e[ipv], w[ipv])
+        wt, lo, hi = ldm_table_map(mu, wf, pt0, pdt, nmu)
+        for ipb in prange(npb):
+            ldm_table_banded(wt, lo, hi, ldp[ipv, ipb], 1.0 / istar[ipv, ipb], tables[ipb])
+            split_cubic_coefficients(tables[ipb], n1, coef[ipb])
+
+        # ------------------------------------------#
+        # Calculate the geometry shared by the bins #
+        # ------------------------------------------#
+        xyc = solve2d(0.0, p[ipv], a[ipv], i[ipv], e[ipv], w[ipv])
+        nsm = transit_samples(times, t0[ipv], p[ipv], xyc, kmean, kmax, gc, n1, nc, ns, exptimes[0],
+                              interval, ipts, basis, area)
 
         # ---------------------------#
-        # Calculate the bounding box #
+        # Calculate the light curves #
         # ---------------------------#
-        bt1, bt4 = bounding_box(kmean, xyc)
-        bt1 -= 0.0015 + exptimes[0]
-        bt4 += 0.0015 + exptimes[0]
-
-        # --------------------------#
-        # Calculate the light curve #
-        # --------------------------#
-        for ipt in range(npt):
-            epoch = floor((times[ipt] - t0[ipv] + 0.5 * p[ipv]) / p[ipv])
-            tc = times[ipt] - (t0[ipv] + epoch * p[ipv])
-            if not (bt1 <= tc <= bt4):
-                flux[ipv, :, ipt] = 1.0
-            else:
-                for isample in range(1, nsamples[0] + 1):
-                    time_offset = exptimes[0] * ((isample - 0.5) / nsamples[0] - 0.5)
-                    z = sep_c(tc + time_offset, xyc)
-                    ap0, kappa = ccia(1.0, kmean, z)
-                    dadk = 2.0*kmean*kappa
-                    if z <= 1.0 - kmax:
-                        for ipb in range(npb):
-                            iplanet = ldm_lookup(z / (1.0 + kmean), gc, n1, coef[ipb])
-                            flux[ipv, ipb, ipt] += (istar[ipv, ipb] - iplanet * ap0 * afac[ipb]) / istar[ipv, ipb]
-                    else:
-                        for ipb in range(npb):
-                            iplanet = ldm_lookup(z / (1.0 + kmean), gc, n1, coef[ipb])
-                            flux[ipv, ipb, ipt] += (istar[ipv, ipb] - iplanet * (ap0 + (k[ipv, ipb]-kmean)*dadk)) / istar[ipv, ipb]
-                flux[ipv, :, ipt] /= nsamples[0]
+        for ipb in prange(npb):
+            f = flux[ipv, ipb]
+            c = coef[ipb]
+            f[:] = 1.0
+            for s in range(nsm):
+                ci = c[interval[s]]
+                iplanet = ci[0] + basis[s, 0] * ci[1] + basis[s, 1] * ci[2] + basis[s, 2] * ci[3]
+                f[ipts[s]] -= iplanet * (area[s, 0] + area[s, 1] * afac[ipb] + area[s, 2] * dk[ipb])
     return flux
 
 
-@njit(parallel=True, fastmath=False)
-def tsmodel_parallel(times: ndarray,
-                   k: ndarray, t0: ndarray, p: ndarray, a: ndarray, i: ndarray, e: ndarray, w: ndarray,
-                   nsamples: ndarray, exptimes: ndarray, ldp: ndarray, istar: ndarray,
-                   pt0: float, pdt: float, rules: ndarray, ng: int,
-                   nthreads: int) -> ndarray:
+tsmodel_serial = njit(parallel=False, fastmath=False)(_tsmodel)
+_tsmodel_parallel = njit(parallel=True, fastmath=False)(_tsmodel)
+
 
+@njit
+def tsmodel_parallel(times: ndarray,
+                     k: ndarray, t0: ndarray, p: ndarray, a: ndarray, i: ndarray, e: ndarray, w: ndarray,
+                     nsamples: ndarray, exptimes: ndarray, ldp: ndarray, istar: ndarray,
+                     pt0: float, pdt: float, rules: ndarray, ng: int,
+                     nthreads: int) -> ndarray:
     nthreads_current = get_num_threads()
     set_num_threads(nthreads)
-
-    if k.ndim != 2:
-        raise ValueError(" The radius ratios must be given as a 2D array with shape (npv, npb)")
-
-    if ldp.ndim != 3:
-        raise ValueError("The limb darkening profiles must be given as a 3D array with shape (npv, npb, nmu)")
-
-    if k.shape[1] != ldp.shape[1]:
-        raise ValueError("The transmission spectrum transit model requires that the number or radius ratios and the number of passbands match.")
-
-    npt = times.size
-    npv = k.shape[0]
-    npb = k.shape[1]
-
-    flux = zeros((npv, npb, npt))  # Model flux
-    nq = rules.shape[2]
-    coef = zeros((npb, ng - 2, 4))   # Split cubic coefficients of the mean intensity tables
-    mu = empty((ng, 2 * nq))
-    wf = zeros((ng, 2 * nq))
-    ldm = empty(ng)
-    xyc = zeros((2, 5))            # Taylor series coefficients for the (x, y) position
-
-    for ipv in range(npv):
-        if isnan(a[ipv]) or (a[ipv] <= 1.0) or (e[ipv] < 0.0) or not valid_radius_ratios(k[ipv]):
-            flux[ipv, :, :] = nan
-            continue
-
-        kmean = mean(k[ipv])
-        kmax = max(k[ipv])
-        afac = k[ipv] ** 2 / kmean ** 2
-
-        # -----------------------------------#
-        # Calculate the limb darkening means #
-        # -----------------------------------#
-        gs, n1 = g_nodes(kmean, ng)
-        gc = split_point(kmean)
-        ldm_nodes(kmean, gs, rules, mu, wf)
-        for ipb in range(npb):
-            ldm_table(mu, wf, pt0, pdt, ldp[ipv, ipb, :], ldm)
-            split_cubic_coefficients(ldm, n1, coef[ipb])
-
-        # -----------------------------------------------------#
-        # Calculate the Taylor series expansions for the orbit #
-        # -----------------------------------------------------#
-        xyc[:, :] = solve2d(0.0, p[ipv], a[ipv], i[ipv], e[ipv], w[ipv])
-
-        # ---------------------------#
-        # Calculate the bounding box #
-        # ---------------------------#
-        bt1, bt4 = bounding_box(kmean, xyc)
-        bt1 -= 0.0015 + exptimes[0]
-        bt4 += 0.0015 + exptimes[0]
-
-        # --------------------------#
-        # Calculate the light curve #
-        # --------------------------#
-        for ipt in prange(npt):
-            epoch = floor((times[ipt] - t0[ipv] + 0.5 * p[ipv]) / p[ipv])
-            tc = times[ipt] - (t0[ipv] + epoch * p[ipv])
-            if not (bt1 <= tc <= bt4):
-                flux[ipv, :, ipt] = 1.0
-            else:
-                for isample in range(1, nsamples[0] + 1):
-                    time_offset = exptimes[0] * ((isample - 0.5) / nsamples[0] - 0.5)
-                    z = sep_c(tc + time_offset, xyc)
-                    ap0, kappa = ccia(1.0, kmean, z)
-                    dadk = 2.0*kmean*kappa
-                    for ipb in range(npb):
-                        if z <= 1.0 - kmax:
-                            for ipb in range(npb):
-                                iplanet = ldm_lookup(z / (1.0 + kmean), gc, n1, coef[ipb])
-                                flux[ipv, ipb, ipt] += (istar[ipv, ipb] - iplanet * ap0 * afac[ipb]) / istar[ipv, ipb]
-                        else:
-                            for ipb in range(npb):
-                                iplanet = ldm_lookup(z / (1.0 + kmean), gc, n1, coef[ipb])
-                                flux[ipv, ipb, ipt] += (istar[ipv, ipb] - iplanet * (
-                                            ap0 + (k[ipv, ipb] - kmean) * dadk)) / istar[ipv, ipb]
-                flux[ipv, :, ipt] /= nsamples[0]
+    flux = _tsmodel_parallel(times, k, t0, p, a, i, e, w, nsamples, exptimes, ldp, istar, pt0, pdt, rules, ng)
     set_num_threads(nthreads_current)
-    return flux
+    return flux

pytransit/models/roadrunner/tsmodel.py

@@ -28,7 +28,7 @@
 #  along with this program.  If not, see <https://www.gnu.org/licenses/>.
 from typing import Union, List, Optional
 
-from numpy import ndarray, isscalar, atleast_2d, atleast_1d, array, zeros, pi
+from numpy import ndarray, isscalar, atleast_2d, atleast_1d, array, zeros, pi, ascontiguousarray
 from pytransit.models.roadrunner.common import calculate_weights_3d
 from scipy.integrate import trapezoid
 
@@ -154,6 +154,7 @@ class TransmissionSpectroscopyModel(RoadRunnerModel):
                 for ipv in range(npv):
                     for ipb in range(npb):
                         istar[ipv, ipb] = 2 * pi * trapezoid(self._ldz * ldpi[ipv, ipb], self._ldz)
+        ldp = ascontiguousarray(ldp)
 
         if self.parallel:
             flux = tsmodel_parallel(self.time, k, t0, p, a, i, e, w, self.nsamples, self.exptimes,

{pytransit-2.9.2.dist-info → pytransit-2.10.0.dist-info}/METADATA

@@ -1,6 +1,6 @@
 Metadata-Version: 2.4
 Name: PyTransit
-Version: 2.9.2
+Version: 2.10.0
 Summary: Fast and painless exoplanet transit light curve modelling.
 Author-email: Hannu Parviainen <hpparvi@gmail.com>
 License: GPLv2
@@ -42,6 +42,7 @@ PyTransit
 =========
 
 [![Licence](http://img.shields.io/badge/license-GPLv2-blue.svg?style=flat)](http://www.gnu.org/licenses/gpl-2.0.html)
+[![Documentation Status](https://readthedocs.org/projects/pytransit/badge/?version=stable)](https://pytransit.readthedocs.io/en/stable/)
 [![MNRAS](https://img.shields.io/badge/MNRAS-10.1093%2Fmnras%2Fstv894-blue.svg)](http://mnras.oxfordjournals.org/content/450/3/3233)
 [![arXiv](http://img.shields.io/badge/arXiv-1504.07433-blue.svg?style=flat)](http://arxiv.org/abs/1504.07433)
 [![ASCL](https://img.shields.io/badge/ASCL-A1505.024-blue.svg?style=flat)](http://ascl.net/1505.024)

{pytransit-2.9.2.dist-info → pytransit-2.10.0.dist-info}/RECORD

@@ -1,4 +1,5 @@
-pytransit/__init__.py,sha256=wbx421r999MaoFLG2jQT9xxuLiO12CjD3rAqCnR85XQ,3406
+pytransit/__init__.py,sha256=IVzjo9sDVfDGfmXNobVblugXWyRZgeWbyTNivP4uNO8,3488
+pytransit/_version.py,sha256=Lxklg9HTweyVZKDc1n479t5P2z89i9gcwyKNcMCstwg,530
 pytransit/limb_darkening.py,sha256=OeBPhUWXGKe_QFWcuSMv0b1sR6HbCdrYtTlYhPX-L_c,3538
 pytransit/contamination/__init__.py,sha256=l2TQSaQSc5S0-6YVB_VXu1ddr4_H8GAw4ypNBV2lgNY,1603
 pytransit/contamination/contamination.py,sha256=-BvHXWTFz3FA7p25YA0ABE5OQhwwZFEsM4oIsgjP7nM,14891
@@ -11,7 +12,7 @@ pytransit/lpf/__init__.py,sha256=SxPyGKv6BYIAR-C1c7oKfEsyHhbsjNofiFNMEWBOOPA,126
 pytransit/lpf/cntlpf.py,sha256=CW88sAfLfeRyT1YWNQqiiG1tQhdFfLoyAa-k0dTEFs0,5915
 pytransit/lpf/eclipselpf.py,sha256=ry4aAijqCuJ2AkRma80KPQT6s0xeteSdyvpinlTK7BU,7924
 pytransit/lpf/logposteriorfunction.py,sha256=4J9cYc5taqCksO8ePDMNAlOrZIXck5k6ZeFEKmEXzPQ,18764
-pytransit/lpf/lpf.py,sha256=GJWOH8sN45s2ytoauk9xOXacZECN3G60sQZqLXcuoWs,26778
+pytransit/lpf/lpf.py,sha256=DtAJO64fK3tq-VZ_3qK1uCD92PSXeb1kvEoJ8a-EwHw,26779
 pytransit/lpf/mclpf.py,sha256=5hYRu5m5pyf9DKihhI4uKXQ-Wr6vCVI3ZCl6cTaHXgE,2271
 pytransit/lpf/mptlpf.py,sha256=nVqP1wwnKv8DJkrWwIEm5OPzA-LImlQLX4l3IwafewQ,6254
 pytransit/lpf/ocllpf.py,sha256=7U1hb_V7M1UNDNlQamv9sGmzuuE2SFX1q5J7DGQUKjE,3273
@@ -90,20 +91,20 @@ pytransit/models/opencl/ma_uniform.cl,sha256=p56RTe7y5W4d07ARvNRZc4CkM6JP0Sc5p3c
 pytransit/models/opencl/orbits.cl,sha256=MG1QVUXSgnubjUi1JWvz_9AJIXh6nUll1tA4EOfxnZA,4230
 pytransit/models/opencl/qpower2.cl,sha256=KDBItLYAyegPupO4YHjBUlzgVP_mr2st5aLtd-5bItQ,4498
 pytransit/models/roadrunner/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
-pytransit/models/roadrunner/common.py,sha256=7CPykGwD-zC1tCXQCoCQnCCzyCworu3K102hhRP8miI,23734
+pytransit/models/roadrunner/common.py,sha256=YdYymuXvpC1XbAlECJVyxympbi3VWLWq7MSpvSth2v4,26872
 pytransit/models/roadrunner/ecintersection.py,sha256=I-QX-j3cf2nyeyGO_m3thOBgtZlJ8i9sFxCBjqLb-os,20461
 pytransit/models/roadrunner/esmodel.py,sha256=L-9xI-Jo_qXyRFQ1F_IeyJsDS-CbdYhH92Al-FNh21E,5474
 pytransit/models/roadrunner/model_eclipse.py,sha256=UIDCQZPQT52tfad9XtJgJpzAPorly3NvRDWOJMCdI18,3125
 pytransit/models/roadrunner/model_ecspec.py,sha256=uao37pDV2a0_rnyqPuMjkEbj6lUG5Swiu1i_LY7N1wU,2651
 pytransit/models/roadrunner/model_full.py,sha256=o8VB7e1cQdnrd81zPmgG0xdDOhoDPUQ-rlab4xNyYu8,6849
 pytransit/models/roadrunner/model_simple.py,sha256=qEZiDxgXwvnfznzWYlmeLYBnlR8Bp4E1rnC-8t95L-8,5843
-pytransit/models/roadrunner/model_trspec.py,sha256=CruVo9LqDkcnwUp6_uzAE_n5sHxWJoTuRJypVqS6Rj0,8077
+pytransit/models/roadrunner/model_trspec.py,sha256=IOR6OY41vv8qTmSC-DGJ2AD5W1QbiWWoOXU2nyilKF4,6696
 pytransit/models/roadrunner/opmodel.py,sha256=v-Ptc-E2a4Vbe5d0qjjAVJT9rUTX9hb2YItC8Qm_w-U,13876
 pytransit/models/roadrunner/opmodel_full.py,sha256=baH78JCwkWJJE4Yjms2Jf64MVd9_nYPhdibezeJQxPI,12588
 pytransit/models/roadrunner/rrmodel.cl,sha256=MGBDT4bW8vHU0wSjOzF4GFK_YzZlcoMUFvwz_1B8X1Q,16782
 pytransit/models/roadrunner/rrmodel.py,sha256=0qcjayQ7S5EHRUM5KP3qqLAF-97RiV0lKrvPPIihTnM,16180
 pytransit/models/roadrunner/rrmodel_cl.py,sha256=SdnarHFI_yetl8xY1u7Rs6Tgh9p70RZpEYFymtDgvTw,27921
-pytransit/models/roadrunner/tsmodel.py,sha256=S1xlWJyDhNUnAL2tm-RqCe6hsadR-je4xgX18R5II6I,7461
+pytransit/models/roadrunner/tsmodel.py,sha256=FbbgZSD_T6vuuNXkzaokGMmDiIFc0GZ7HhZmk3W2sFQ,7517
 pytransit/orbits/__init__.py,sha256=Asb5OytXI9duU0MGF5ufV0QCD9m5Ku_S-pPOkqQIFfM,884
 pytransit/orbits/orbits.py,sha256=e8Gt3VL8dk9HXnsEDfE6gsjRtE7QyvNgCKYkdTFvk-w,4936
 pytransit/orbits/orbits_py.py,sha256=BqP89e5s7sXe-Q2-SIYFTD5lbJfrHBNRV5xmrNHC0o8,20423
@@ -136,8 +137,8 @@ pytransit/utils/io/__init__.py,sha256=Wh74NDb7IA_cXpvnQBtoAccZXYzm4oqvHJp8lHgWf7
 pytransit/utils/io/base.py,sha256=VH8LPu-Fzjv2kOftf91py_sRlgJBHMUtu0QyCI2kymc,15984
 pytransit/utils/io/lcdata.py,sha256=PIy3iKbuCj_ytk8IZ9bIFsrwNxLnGC2VtO336eg4tJI,32924
 pytransit/utils/io/rvdata.py,sha256=9hx0HjN5xmX6jvu73on6mLYZl50aCzuuz9RoKhJaaIg,7229
-pytransit-2.9.2.dist-info/licenses/LICENSE,sha256=s6DMeIC72jjJFdFX2MTU3lCz_4UCUQAcCMtxAqYBiv0,18045
-pytransit-2.9.2.dist-info/METADATA,sha256=O8qI3Uc0S0tDUqJJf0-Z-UN5oTvV2AdOeMTw-qLz6e4,5718
-pytransit-2.9.2.dist-info/WHEEL,sha256=YVMoNqKzERt-wjUZwJ33xBGAwnFl-4cqbYkTtWa4itE,91
-pytransit-2.9.2.dist-info/top_level.txt,sha256=Y4wC10N9yBBNslkDDVW6JTmTjMRD5GPgu992lWf6rwk,10
-pytransit-2.9.2.dist-info/RECORD,,
+pytransit-2.10.0.dist-info/licenses/LICENSE,sha256=s6DMeIC72jjJFdFX2MTU3lCz_4UCUQAcCMtxAqYBiv0,18045
+pytransit-2.10.0.dist-info/METADATA,sha256=Gnd204JSUeLfeShGWVYvOS8TjE-c8kO73Ruv5JH8A0s,5856
+pytransit-2.10.0.dist-info/WHEEL,sha256=YVMoNqKzERt-wjUZwJ33xBGAwnFl-4cqbYkTtWa4itE,91
+pytransit-2.10.0.dist-info/top_level.txt,sha256=Y4wC10N9yBBNslkDDVW6JTmTjMRD5GPgu992lWf6rwk,10
+pytransit-2.10.0.dist-info/RECORD,,