pianistpy 0.9.0__tar.gz

pianistpy-0.9.0/AUTHORS.md

@@ -0,0 +1,7 @@
+# Authors
+
+PIANISTpy is written and maintained by
+
+- Barry T. Chiang (Yale University, 2026), barrytchiang@gmail.com
+
+To cite PIANISTpy, cite the paper that describes it (see the citation file).

pianistpy-0.9.0/LICENSE

@@ -0,0 +1,29 @@
+BSD 3-Clause License
+
+Copyright (c) 2026, Barry T. Chiang
+All rights reserved.
+
+Redistribution and use in source and binary forms, with or without
+modification, are permitted provided that the following conditions are met:
+
+1. Redistributions of source code must retain the above copyright notice, this
+   list of conditions and the following disclaimer.
+
+2. Redistributions in binary form must reproduce the above copyright notice,
+   this list of conditions and the following disclaimer in the documentation
+   and/or other materials provided with the distribution.
+
+3. Neither the name of the copyright holder nor the names of its
+   contributors may be used to endorse or promote products derived from
+   this software without specific prior written permission.
+
+THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
+FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

pianistpy-0.9.0/PKG-INFO

@@ -0,0 +1,224 @@
+Metadata-Version: 2.4
+Name: pianistpy
+Version: 0.9.0
+Summary: Distribution functions, particles and observables of anisotropic spherical systems.
+Author-email: "Barry T. Chiang" <barrytchiang@gmail.com>
+License: BSD 3-Clause License
+        
+        Copyright (c) 2026, Barry T. Chiang
+        All rights reserved.
+        
+        Redistribution and use in source and binary forms, with or without
+        modification, are permitted provided that the following conditions are met:
+        
+        1. Redistributions of source code must retain the above copyright notice, this
+           list of conditions and the following disclaimer.
+        
+        2. Redistributions in binary form must reproduce the above copyright notice,
+           this list of conditions and the following disclaimer in the documentation
+           and/or other materials provided with the distribution.
+        
+        3. Neither the name of the copyright holder nor the names of its
+           contributors may be used to endorse or promote products derived from
+           this software without specific prior written permission.
+        
+        THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"
+        AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE
+        IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE
+        DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT HOLDER OR CONTRIBUTORS BE LIABLE
+        FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL
+        DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR
+        SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER
+        CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,
+        OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE USE
+        OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.
+        
+Project-URL: Homepage, https://pianistpy.com
+Project-URL: Documentation, https://pianistpy.com/docs/
+Project-URL: Source, https://github.com/barrytchiang/PIANISTpy
+Project-URL: Issues, https://github.com/barrytchiang/PIANISTpy/issues
+Requires-Python: >=3.10
+Description-Content-Type: text/markdown
+License-File: LICENSE
+License-File: AUTHORS.md
+Requires-Dist: numpy>=2.0.0
+Requires-Dist: scipy>=1.14.0
+Requires-Dist: h5py>=3.11.0
+Requires-Dist: astropy>=7.2.0
+Requires-Dist: mpmath>=0.17
+Requires-Dist: sympy>=1.2
+Requires-Dist: jax>=0.6.1
+Requires-Dist: diffrax>=0.7.0
+Requires-Dist: optimistix>=0.0.11
+Requires-Dist: filelock
+Dynamic: license-file
+
+<p align="center">
+  <img src="https://raw.githubusercontent.com/barrytchiang/PIANISTpy/main/docs/_static/logo.png" alt="PIANISTpy" width="320">
+</p>
+
+# PIANISTpy
+
+[![Tests](https://github.com/barrytchiang/PIANISTpy/actions/workflows/tests.yml/badge.svg)](https://github.com/barrytchiang/PIANISTpy/actions/workflows/tests.yml)
+[![Coverage](https://codecov.io/gh/barrytchiang/PIANISTpy/branch/main/graph/badge.svg)](https://codecov.io/gh/barrytchiang/PIANISTpy)
+[![Documentation](https://img.shields.io/badge/docs-pianistpy.com-00356B)](https://pianistpy.com)
+[![License](https://img.shields.io/badge/license-BSD--3--Clause-blue)](LICENSE)
+
+**P**article **I**nitial conditions with **ANIS**o**T**ropy: distribution functions, particle initial conditions
+and observables of spherical systems with anisotropic velocities, truncated self-consistently in phase space.
+
+PIANISTpy builds equilibrium models of spherical systems (dark matter halos and subhalos, star clusters, dwarf
+galaxies) with a velocity anisotropy you choose. The outer boundary of a model is set in phase space, so the
+truncated model is itself an equilibrium and needs no relaxation run before it is used.
+
+> Status: version 0.9.0, a pre-release. Version 1.0.0, released with the paper that describes the method,
+> adds the certified accuracy tiers on both backends, the remaining capabilities and the citation.
+
+## What it does
+
+- **Density profiles.** Any double power law in its three slopes, with the Navarro-Frenk-White, Hernquist and
+  Plummer profiles as named presets.
+- **Velocity anisotropy.** Isotropic, constant, Osipkov-Merritt, and the generalized Osipkov-Merritt (Cuddeford)
+  models, radially biased (Type I) or tangentially biased (Type II).
+- **Distribution functions.** The distribution function (DF) of each model, from the analytic inversion of its
+  density, returned whether or not it is physical, with a statement of where it is not.
+- **Truncation in phase space.** Method I truncates a model's own DF at an energy you give, or at the energy that
+  yields a target mass or a target radius. Method II gives a new anisotropy to a density that is already
+  truncated.
+- **Particles.** Positions and velocities drawn by inversion, with no rejection step, and reproducible particle
+  by particle.
+- **Observables.** Intrinsic and line-of-sight velocity dispersions, line-of-sight velocity distributions, and
+  their averages over apertures.
+- **Superposition.** Distribution functions inverted from one truncated density (Method II children of one
+  parent), combined with positive weights that sum to one into one equilibrium model.
+- **Inference.** The log-likelihood of a model against data, as a JAX function of the components' weights and
+  the scales of the units.
+
+## Installation
+
+```
+pip install pianistpy
+```
+
+PIANISTpy needs Python 3.10 or later, NumPy, SciPy, h5py, Astropy, mpmath, SymPy, JAX, Diffrax, Optimistix and filelock.
+The oldest versions the test suite passes with are NumPy 2.0.0, SciPy 1.14.0, h5py 3.11.0, Astropy 7.2.0,
+mpmath 0.17, SymPy 1.2, and together JAX 0.6.1, Diffrax 0.7.0 and Optimistix 0.0.11 (these three are bound as
+one set; older Diffrax or Optimistix releases pass only with a newer JAX); `pip` enforces them. It runs on Linux,
+macOS and Windows.
+
+## Quick start
+
+The example of the paper's first figure: an NFW halo with a constant velocity anisotropy of 0.25, truncated in
+phase space at the virial mass of a concentration-10 halo, presented as a Milky Way-mass halo.
+
+```python
+import pianistpy as pp
+import astropy.units as u
+
+# an untruncated NFW density profile with a constant velocity anisotropy of 0.25
+halo = pp.Component(pp.DoublePowerLaw.nfw(), pp.Anisotropy(beta_c=0.25))
+
+# the untruncated distribution function
+df = pp.compute_df(halo)
+
+# Method I: truncate it so that the truncated mass equals the virial mass of the
+# untruncated NFW halo of concentration 10 (the mass inside 10 scale radii)
+model = pp.truncate(df, pp.Truncation(mass_t=1.0, concentration=10))
+
+# physical units: a halo of 1e12 solar masses with R_vir = c * r_0 = 200 kpc
+model = model.to("physical", mass_0=1e12 * u.Msun, r_0=20 * u.kpc)
+
+# draw ten million equal-mass particles centered at the origin and write them to a file
+particles = pp.Sampler(model).sample(10_000_000, seed=1)
+particles.write("NFW_Const0p25_Mt1e12_PARIC.h5")
+
+# the line-of-sight velocity dispersion at three projected radii
+sigma = pp.sigma_los(model, [2, 20, 100] * u.kpc)
+```
+
+`pp.solve(halo, truncation)` runs the two steps at once.
+
+## Documentation
+
+The user guide, the notebook tutorials, the technical notes and the API reference are at
+[pianistpy.com](https://pianistpy.com).
+
+## Example notebooks
+
+Tutorials for each workflow, and notebooks that reproduce the paper's figures, are in
+[PIANISTpy-notebooks](https://github.com/barrytchiang/PIANISTpy-notebooks), tagged with the package version
+they were run with.
+
+## Results and files
+
+Every result (a distribution function, a model, a set of particles) is an immutable object that holds what it
+was computed from: the inputs, the numerical settings and the version of the code. Its `write` method stores
+it in a single HDF5 file, and `pp.read(path)` returns an equal object, whatever kind of result the file holds.
+Nothing is written unless you ask.
+
+A `Model` is always a valid equilibrium. A request for which none exists raises `UnphysicalModelError`, and the
+exception carries the distribution function so that you can inspect it.
+
+## Units
+
+Everything is computed and returned in internal units, the dimensionless normalization of Widrow (2000) in which
+the scale radius and the central potential of the untruncated profile are 1. A result's `to` method presents it
+in model units (G = M_0 = r_0 = 1) or, through Astropy, in physical units, without recomputing it, and what is
+computed from it is returned in that system.
+
+Three inputs of a truncation are given in the units that are natural for them, as their documentation states: the
+truncation energy `k_t` in internal units, a target radius `r_t` in units of the scale radius, and a target mass
+`mass_t` in units of M_0, together with the concentration that defines M_0.
+
+## Accuracy
+
+Every result carries a validation report, a short and fixed summary of the checks it passed, with its tier: a
+check within the accuracy goal (a relative error of 1e-6) is silent; above the goal, within the ceiling (1e-4) or
+within the hard limit (1e-2), the result is returned with a warning in the package's log and as an
+`AccuracyWarning` (Python's standard filter silences or escalates it), and the flag stays with the result and with
+everything computed from it; beyond the hard limit the result is refused.
+
+The numerical settings come in two tested tiers, both held to these checks by the test suite:
+
+```python
+pp.Settings()            # the default tier
+pp.Settings(quick=True)  # a faster tier for exploration
+```
+
+## Session switches
+
+None of these changes a number.
+
+```python
+pp.set_verbosity("info")   # "quiet" (the default), "info" or "debug"
+pp.set_timing(True)        # print the wall time of each stage
+pp.set_workers("half")     # "auto" (the default), "full", "half" or a number of cores
+pp.set_compilation_cache(False)  # the compilation cache is on by default; a path moves it
+pp.set_platform("cpu")     # "auto" (the default), "cpu" (the CPU-only mode) or "gpu"; before the first computation
+```
+
+The CPU-only mode gives the same numbers on every machine; the environment variable `PIANISTPY_PLATFORM=cpu`
+sets it before the package is imported. Each file records the backend that made it.
+
+Hardware notes. The workers' switch sets how many cores the package's parallel work uses; `"auto"` uses every
+allowed core inside a batch job and half of them anywhere else, and no value changes a number. A GPU accelerates the
+two batched table builds, the speed table (2.7 to 7 times an 8-core CPU on an RTX 3080 Ti) and the line-of-sight
+velocity table (1.0 to 1.4 times a 32-core node); the distribution function and the potential solve gain nothing from
+a GPU and run at their own speed on the CPU. The CPU-only install is complete and is the right one for a laptop. JAX's
+CUDA builds are for Linux: on Windows the GPU is reached through the Windows Subsystem for Linux with the Linux
+wheels, which pays only for large table builds; a native Windows Python runs the CPU-only mode. A GPU and a CPU give
+results that differ at rounding on the kernels; the accuracy goal holds on both, and the file's provenance says which
+made it. The compilation cache lives in the home folder by default (`~/.cache/pianistpy/jax`) and can fill a small
+one; `set_compilation_cache` with a path moves it.
+
+## Citing
+
+If you use PIANISTpy, cite the paper that describes it. The citation entry is added with the 1.0.0 release.
+
+## License
+
+BSD 3-Clause. See `LICENSE`.
+
+## Contributing
+
+See `CONTRIBUTING.md` and the code conventions in `STYLE_GUIDE.md`.

pianistpy-0.9.0/README.md

@@ -0,0 +1,169 @@
+<p align="center">
+  <img src="https://raw.githubusercontent.com/barrytchiang/PIANISTpy/main/docs/_static/logo.png" alt="PIANISTpy" width="320">
+</p>
+
+# PIANISTpy
+
+[![Tests](https://github.com/barrytchiang/PIANISTpy/actions/workflows/tests.yml/badge.svg)](https://github.com/barrytchiang/PIANISTpy/actions/workflows/tests.yml)
+[![Coverage](https://codecov.io/gh/barrytchiang/PIANISTpy/branch/main/graph/badge.svg)](https://codecov.io/gh/barrytchiang/PIANISTpy)
+[![Documentation](https://img.shields.io/badge/docs-pianistpy.com-00356B)](https://pianistpy.com)
+[![License](https://img.shields.io/badge/license-BSD--3--Clause-blue)](LICENSE)
+
+**P**article **I**nitial conditions with **ANIS**o**T**ropy: distribution functions, particle initial conditions
+and observables of spherical systems with anisotropic velocities, truncated self-consistently in phase space.
+
+PIANISTpy builds equilibrium models of spherical systems (dark matter halos and subhalos, star clusters, dwarf
+galaxies) with a velocity anisotropy you choose. The outer boundary of a model is set in phase space, so the
+truncated model is itself an equilibrium and needs no relaxation run before it is used.
+
+> Status: version 0.9.0, a pre-release. Version 1.0.0, released with the paper that describes the method,
+> adds the certified accuracy tiers on both backends, the remaining capabilities and the citation.
+
+## What it does
+
+- **Density profiles.** Any double power law in its three slopes, with the Navarro-Frenk-White, Hernquist and
+  Plummer profiles as named presets.
+- **Velocity anisotropy.** Isotropic, constant, Osipkov-Merritt, and the generalized Osipkov-Merritt (Cuddeford)
+  models, radially biased (Type I) or tangentially biased (Type II).
+- **Distribution functions.** The distribution function (DF) of each model, from the analytic inversion of its
+  density, returned whether or not it is physical, with a statement of where it is not.
+- **Truncation in phase space.** Method I truncates a model's own DF at an energy you give, or at the energy that
+  yields a target mass or a target radius. Method II gives a new anisotropy to a density that is already
+  truncated.
+- **Particles.** Positions and velocities drawn by inversion, with no rejection step, and reproducible particle
+  by particle.
+- **Observables.** Intrinsic and line-of-sight velocity dispersions, line-of-sight velocity distributions, and
+  their averages over apertures.
+- **Superposition.** Distribution functions inverted from one truncated density (Method II children of one
+  parent), combined with positive weights that sum to one into one equilibrium model.
+- **Inference.** The log-likelihood of a model against data, as a JAX function of the components' weights and
+  the scales of the units.
+
+## Installation
+
+```
+pip install pianistpy
+```
+
+PIANISTpy needs Python 3.10 or later, NumPy, SciPy, h5py, Astropy, mpmath, SymPy, JAX, Diffrax, Optimistix and filelock.
+The oldest versions the test suite passes with are NumPy 2.0.0, SciPy 1.14.0, h5py 3.11.0, Astropy 7.2.0,
+mpmath 0.17, SymPy 1.2, and together JAX 0.6.1, Diffrax 0.7.0 and Optimistix 0.0.11 (these three are bound as
+one set; older Diffrax or Optimistix releases pass only with a newer JAX); `pip` enforces them. It runs on Linux,
+macOS and Windows.
+
+## Quick start
+
+The example of the paper's first figure: an NFW halo with a constant velocity anisotropy of 0.25, truncated in
+phase space at the virial mass of a concentration-10 halo, presented as a Milky Way-mass halo.
+
+```python
+import pianistpy as pp
+import astropy.units as u
+
+# an untruncated NFW density profile with a constant velocity anisotropy of 0.25
+halo = pp.Component(pp.DoublePowerLaw.nfw(), pp.Anisotropy(beta_c=0.25))
+
+# the untruncated distribution function
+df = pp.compute_df(halo)
+
+# Method I: truncate it so that the truncated mass equals the virial mass of the
+# untruncated NFW halo of concentration 10 (the mass inside 10 scale radii)
+model = pp.truncate(df, pp.Truncation(mass_t=1.0, concentration=10))
+
+# physical units: a halo of 1e12 solar masses with R_vir = c * r_0 = 200 kpc
+model = model.to("physical", mass_0=1e12 * u.Msun, r_0=20 * u.kpc)
+
+# draw ten million equal-mass particles centered at the origin and write them to a file
+particles = pp.Sampler(model).sample(10_000_000, seed=1)
+particles.write("NFW_Const0p25_Mt1e12_PARIC.h5")
+
+# the line-of-sight velocity dispersion at three projected radii
+sigma = pp.sigma_los(model, [2, 20, 100] * u.kpc)
+```
+
+`pp.solve(halo, truncation)` runs the two steps at once.
+
+## Documentation
+
+The user guide, the notebook tutorials, the technical notes and the API reference are at
+[pianistpy.com](https://pianistpy.com).
+
+## Example notebooks
+
+Tutorials for each workflow, and notebooks that reproduce the paper's figures, are in
+[PIANISTpy-notebooks](https://github.com/barrytchiang/PIANISTpy-notebooks), tagged with the package version
+they were run with.
+
+## Results and files
+
+Every result (a distribution function, a model, a set of particles) is an immutable object that holds what it
+was computed from: the inputs, the numerical settings and the version of the code. Its `write` method stores
+it in a single HDF5 file, and `pp.read(path)` returns an equal object, whatever kind of result the file holds.
+Nothing is written unless you ask.
+
+A `Model` is always a valid equilibrium. A request for which none exists raises `UnphysicalModelError`, and the
+exception carries the distribution function so that you can inspect it.
+
+## Units
+
+Everything is computed and returned in internal units, the dimensionless normalization of Widrow (2000) in which
+the scale radius and the central potential of the untruncated profile are 1. A result's `to` method presents it
+in model units (G = M_0 = r_0 = 1) or, through Astropy, in physical units, without recomputing it, and what is
+computed from it is returned in that system.
+
+Three inputs of a truncation are given in the units that are natural for them, as their documentation states: the
+truncation energy `k_t` in internal units, a target radius `r_t` in units of the scale radius, and a target mass
+`mass_t` in units of M_0, together with the concentration that defines M_0.
+
+## Accuracy
+
+Every result carries a validation report, a short and fixed summary of the checks it passed, with its tier: a
+check within the accuracy goal (a relative error of 1e-6) is silent; above the goal, within the ceiling (1e-4) or
+within the hard limit (1e-2), the result is returned with a warning in the package's log and as an
+`AccuracyWarning` (Python's standard filter silences or escalates it), and the flag stays with the result and with
+everything computed from it; beyond the hard limit the result is refused.
+
+The numerical settings come in two tested tiers, both held to these checks by the test suite:
+
+```python
+pp.Settings()            # the default tier
+pp.Settings(quick=True)  # a faster tier for exploration
+```
+
+## Session switches
+
+None of these changes a number.
+
+```python
+pp.set_verbosity("info")   # "quiet" (the default), "info" or "debug"
+pp.set_timing(True)        # print the wall time of each stage
+pp.set_workers("half")     # "auto" (the default), "full", "half" or a number of cores
+pp.set_compilation_cache(False)  # the compilation cache is on by default; a path moves it
+pp.set_platform("cpu")     # "auto" (the default), "cpu" (the CPU-only mode) or "gpu"; before the first computation
+```
+
+The CPU-only mode gives the same numbers on every machine; the environment variable `PIANISTPY_PLATFORM=cpu`
+sets it before the package is imported. Each file records the backend that made it.
+
+Hardware notes. The workers' switch sets how many cores the package's parallel work uses; `"auto"` uses every
+allowed core inside a batch job and half of them anywhere else, and no value changes a number. A GPU accelerates the
+two batched table builds, the speed table (2.7 to 7 times an 8-core CPU on an RTX 3080 Ti) and the line-of-sight
+velocity table (1.0 to 1.4 times a 32-core node); the distribution function and the potential solve gain nothing from
+a GPU and run at their own speed on the CPU. The CPU-only install is complete and is the right one for a laptop. JAX's
+CUDA builds are for Linux: on Windows the GPU is reached through the Windows Subsystem for Linux with the Linux
+wheels, which pays only for large table builds; a native Windows Python runs the CPU-only mode. A GPU and a CPU give
+results that differ at rounding on the kernels; the accuracy goal holds on both, and the file's provenance says which
+made it. The compilation cache lives in the home folder by default (`~/.cache/pianistpy/jax`) and can fill a small
+one; `set_compilation_cache` with a path moves it.
+
+## Citing
+
+If you use PIANISTpy, cite the paper that describes it. The citation entry is added with the 1.0.0 release.
+
+## License
+
+BSD 3-Clause. See `LICENSE`.
+
+## Contributing
+
+See `CONTRIBUTING.md` and the code conventions in `STYLE_GUIDE.md`.

pianistpy-0.9.0/pianistpy/__init__.py

@@ -0,0 +1,77 @@
+"""PIANISTpy: equilibrium models of anisotropic spherical systems."""
+
+__all__ = [
+    "AccuracyError",
+    "AccuracyWarning",
+    "Anisotropy",
+    "Aperture",
+    "Component",
+    "DistributionFunction",
+    "DoublePowerLaw",
+    "InputError",
+    "LosTable",
+    "Model",
+    "PianistpyError",
+    "PianistpyWarning",
+    "ResultError",
+    "Sampler",
+    "Settings",
+    "Truncation",
+    "UnphysicalModelError",
+    "__version__",
+    "compute_df",
+    "log",
+    "platform",
+    "read",
+    "set_compilation_cache",
+    "set_timing",
+    "set_verbosity",
+    "set_platform",
+    "set_workers",
+    "sigma_los",
+    "sigma_los_aperture",
+    "solve",
+    "superpose",
+    "truncate",
+]
+
+# The reader of data files registers itself on import.
+from . import inference as inference
+from ._version import __version__
+from .anisotropy import Anisotropy
+from .df import (
+    Component,
+    DistributionFunction,
+    Model,
+    Truncation,
+    compute_df,
+    solve,
+    truncate,
+)
+from .exceptions import (
+    AccuracyError,
+    AccuracyWarning,
+    InputError,
+    PianistpyError,
+    PianistpyWarning,
+    ResultError,
+    UnphysicalModelError,
+)
+from .io import read
+from .observables import Aperture, LosTable, sigma_los, sigma_los_aperture
+from .particles import Sampler
+from .profiles import DoublePowerLaw
+from .runtime import (
+    configure_jax,
+    log,
+    platform,
+    set_compilation_cache,
+    set_platform,
+    set_timing,
+    set_verbosity,
+    set_workers,
+)
+from .settings import Settings
+from .superposition import superpose
+
+configure_jax()

pianistpy-0.9.0/pianistpy/_version.py

@@ -0,0 +1,5 @@
+"""Version of the package."""
+
+__all__ = ["__version__"]
+
+__version__ = "0.9.0"

pianistpy-0.9.0/pianistpy/anisotropy/__init__.py

@@ -0,0 +1,5 @@
+"""Velocity anisotropy models."""
+
+__all__ = ["Anisotropy"]
+
+from .core import Anisotropy