photomancy 0.0.1__tar.gz

photomancy-0.0.1/.gitignore

@@ -0,0 +1,221 @@
+# Byte-compiled / optimized / DLL files
+__pycache__/
+*.py[codz]
+*$py.class
+
+# C extensions
+*.so
+
+# Distribution / packaging
+.Python
+build/
+develop-eggs/
+dist/
+downloads/
+eggs/
+.eggs/
+lib/
+lib64/
+parts/
+sdist/
+var/
+wheels/
+share/python-wheels/
+*.egg-info/
+.installed.cfg
+*.egg
+MANIFEST
+
+# PyInstaller
+#   Usually these files are written by a python script from a template
+#   before PyInstaller builds the exe, so as to inject date/other infos into it.
+*.manifest
+*.spec
+
+# Installer logs
+pip-log.txt
+pip-delete-this-directory.txt
+
+# Unit test / coverage reports
+htmlcov/
+.tox/
+.nox/
+.coverage
+.coverage.*
+.cache
+nosetests.xml
+coverage.xml
+*.cover
+*.py.cover
+.hypothesis/
+.pytest_cache/
+cover/
+
+# Translations
+*.mo
+*.pot
+
+# Django stuff:
+*.log
+local_settings.py
+db.sqlite3
+db.sqlite3-journal
+
+# Flask stuff:
+instance/
+.webassets-cache
+
+# Scrapy stuff:
+.scrapy
+
+# Sphinx documentation
+docs/_build/
+
+# PyBuilder
+.pybuilder/
+target/
+
+# Jupyter Notebook
+.ipynb_checkpoints
+
+# IPython
+profile_default/
+ipython_config.py
+
+# pyenv
+#   For a library or package, you might want to ignore these files since the code is
+#   intended to run in multiple environments; otherwise, check them in:
+# .python-version
+
+# pipenv
+#   According to pypa/pipenv#598, it is recommended to include Pipfile.lock in version control.
+#   However, in case of collaboration, if having platform-specific dependencies or dependencies
+#   having no cross-platform support, pipenv may install dependencies that don't work, or not
+#   install all needed dependencies.
+# Pipfile.lock
+
+# UV
+#   Similar to Pipfile.lock, it is generally recommended to include uv.lock in version control.
+#   This is especially recommended for binary packages to ensure reproducibility, and is more
+#   commonly ignored for libraries.
+# uv.lock
+
+# poetry
+#   Similar to Pipfile.lock, it is generally recommended to include poetry.lock in version control.
+#   This is especially recommended for binary packages to ensure reproducibility, and is more
+#   commonly ignored for libraries.
+#   https://python-poetry.org/docs/basic-usage/#commit-your-poetrylock-file-to-version-control
+# poetry.lock
+# poetry.toml
+
+# pdm
+#   Similar to Pipfile.lock, it is generally recommended to include pdm.lock in version control.
+#   pdm recommends including project-wide configuration in pdm.toml, but excluding .pdm-python.
+#   https://pdm-project.org/en/latest/usage/project/#working-with-version-control
+# pdm.lock
+# pdm.toml
+.pdm-python
+.pdm-build/
+
+# pixi
+#   Similar to Pipfile.lock, it is generally recommended to include pixi.lock in version control.
+# pixi.lock
+#   Pixi creates a virtual environment in the .pixi directory, just like venv module creates one
+#   in the .venv directory. It is recommended not to include this directory in version control.
+.pixi
+
+# PEP 582; used by e.g. github.com/David-OConnor/pyflow and github.com/pdm-project/pdm
+__pypackages__/
+
+# Celery stuff
+celerybeat-schedule
+celerybeat.pid
+
+# Redis
+*.rdb
+*.aof
+*.pid
+
+# RabbitMQ
+mnesia/
+rabbitmq/
+rabbitmq-data/
+
+# ActiveMQ
+activemq-data/
+
+# SageMath parsed files
+*.sage.py
+
+# Environments
+.env
+.envrc
+.venv
+env/
+venv/
+ENV/
+env.bak/
+venv.bak/
+
+# Spyder project settings
+.spyderproject
+.spyproject
+
+# Rope project settings
+.ropeproject
+
+# mkdocs documentation
+/site
+
+# mypy
+.mypy_cache/
+.dmypy.json
+dmypy.json
+
+# Pyre type checker
+.pyre/
+
+# pytype static type analyzer
+.pytype/
+
+# Cython debug symbols
+cython_debug/
+
+# PyCharm
+#   JetBrains specific template is maintained in a separate JetBrains.gitignore that can
+#   be found at https://github.com/github/gitignore/blob/main/Global/JetBrains.gitignore
+#   and can be added to the global gitignore or merged into this file.  For a more nuclear
+#   option (not recommended) you can uncomment the following to ignore the entire idea folder.
+# .idea/
+
+# Abstra
+#   Abstra is an AI-powered process automation framework.
+#   Ignore directories containing user credentials, local state, and settings.
+#   Learn more at https://abstra.io/docs
+.abstra/
+
+# Visual Studio Code
+#   Visual Studio Code specific template is maintained in a separate VisualStudioCode.gitignore 
+#   that can be found at https://github.com/github/gitignore/blob/main/Global/VisualStudioCode.gitignore
+#   and can be added to the global gitignore or merged into this file. However, if you prefer, 
+#   you could uncomment the following to ignore the entire vscode folder
+# .vscode/
+# Temporary file for partial code execution
+tempCodeRunnerFile.py
+
+# Ruff stuff:
+.ruff_cache/
+
+# PyPI configuration file
+.pypirc
+
+# Marimo
+marimo/_static/
+marimo/_lsp/
+__marimo__/
+
+# Streamlit
+.streamlit/secrets.toml
+
+# hatch-vcs generated version file
+src/photomancy/_version.py

photomancy-0.0.1/.pre-commit-config.yaml

@@ -0,0 +1,20 @@
+repos:
+  - repo: https://github.com/pre-commit/pre-commit-hooks
+    rev: "v6.0.0"
+    hooks:
+      - id: trailing-whitespace
+      - id: name-tests-test
+        args: [--pytest-test-first]
+      - id: end-of-file-fixer
+  - repo: https://github.com/astral-sh/ruff-pre-commit
+    rev: v0.15.8
+    hooks:
+      - id: ruff-check
+        args: [--fix]
+      - id: ruff-format
+  - repo: https://github.com/compilerla/conventional-pre-commit
+    rev: v4.4.0
+    hooks:
+      - id: conventional-pre-commit
+        stages: [commit-msg]
+        args: []

photomancy-0.0.1/.python-version

@@ -0,0 +1 @@
+3.12

photomancy-0.0.1/.readthedocs.yaml

@@ -0,0 +1,20 @@
+# Required
+version: 2
+
+# Set the OS, Python version and other tools you might need
+build:
+  os: ubuntu-22.04
+  tools:
+    python: "3.12"
+
+python:
+  install:
+    - method: pip
+      path: .
+      extra_requirements:
+        # Install with the [docs] flag for sphinx docs specific extensions
+        - docs
+
+# Build documentation in the "docs/" directory with Sphinx
+sphinx:
+  configuration: docs/conf.py

photomancy-0.0.1/CHANGELOG.md

@@ -0,0 +1,17 @@
+# Changelog
+
+## 0.0.1 (2026-06-18)
+
+
+### Features
+
+* add Backend/Posterior protocol and Laplace backend ([c18e80c](https://github.com/CoreySpohn/photomancy/commit/c18e80c1104e447dbcd2dcdd61346959db5ec662))
+* **core:** add scene-PyTree logdensity assembly ([77f80ed](https://github.com/CoreySpohn/photomancy/commit/77f80edea1cd82f5dbef604dfa92b5ae37759723))
+* scaffold photomancy and port orbit fitting from orbix ([d3fe243](https://github.com/CoreySpohn/photomancy/commit/d3fe2436403568c8928788ef43f9a14fe7fae130))
+
+
+### Miscellaneous Chores
+
+* release 0.0.1 ([25be93d](https://github.com/CoreySpohn/photomancy/commit/25be93d29830aab5f51b6482c1c6950acd2eee01))
+
+## Changelog

photomancy-0.0.1/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Corey Spohn
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

photomancy-0.0.1/PKG-INFO

@@ -0,0 +1,86 @@
+Metadata-Version: 2.4
+Name: photomancy
+Version: 0.0.1
+Summary: Divination from light: Bayesian inference and value-of-information for HWO direct imaging
+Project-URL: Homepage, https://github.com/CoreySpohn/photomancy
+Project-URL: Issues, https://github.com/CoreySpohn/photomancy/issues
+Author-email: Corey Spohn <corey.a.spohn@nasa.gov>
+License: MIT License
+        
+        Copyright (c) 2026 Corey Spohn
+        
+        Permission is hereby granted, free of charge, to any person obtaining a copy
+        of this software and associated documentation files (the "Software"), to deal
+        in the Software without restriction, including without limitation the rights
+        to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+        copies of the Software, and to permit persons to whom the Software is
+        furnished to do so, subject to the following conditions:
+        
+        The above copyright notice and this permission notice shall be included in all
+        copies or substantial portions of the Software.
+        
+        THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+        IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+        FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+        AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+        LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+        OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+        SOFTWARE.
+License-File: LICENSE
+Classifier: Development Status :: 2 - Pre-Alpha
+Classifier: Intended Audience :: Science/Research
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering :: Astronomy
+Requires-Python: >=3.11
+Requires-Dist: blackjax>=1.2
+Requires-Dist: equinox
+Requires-Dist: hwoutils
+Requires-Dist: jax>=0.8.1
+Requires-Dist: jaxlib>=0.8.1
+Requires-Dist: jaxtyping
+Requires-Dist: numpy
+Requires-Dist: numpyro>=0.15.0
+Requires-Dist: orbix
+Requires-Dist: skyscapes
+Provides-Extra: dev
+Requires-Dist: pre-commit; extra == 'dev'
+Provides-Extra: docs
+Requires-Dist: matplotlib; extra == 'docs'
+Requires-Dist: myst-nb; extra == 'docs'
+Requires-Dist: sphinx; extra == 'docs'
+Requires-Dist: sphinx-autoapi; extra == 'docs'
+Requires-Dist: sphinx-autodoc-typehints; extra == 'docs'
+Requires-Dist: sphinx-book-theme; extra == 'docs'
+Provides-Extra: imaging
+Requires-Dist: coronalyze; extra == 'imaging'
+Provides-Extra: npe
+Requires-Dist: flowjax; extra == 'npe'
+Provides-Extra: ns
+Requires-Dist: jaxns; extra == 'ns'
+Provides-Extra: test
+Requires-Dist: arviz; extra == 'test'
+Requires-Dist: hypothesis; extra == 'test'
+Requires-Dist: nox; extra == 'test'
+Requires-Dist: pytest; extra == 'test'
+Requires-Dist: pytest-cov; extra == 'test'
+Description-Content-Type: text/markdown
+
+# photomancy
+
+> Divination from light.
+
+A JAX-native Bayesian inference and value-of-information engine for the Habitable
+Worlds Observatory direct-imaging simulation suite. `orbix` builds the geometry,
+`skyscapes` builds the scene, and **photomancy** divines the scene back from the
+data: posteriors, evidence, and the next-best observation, over orbits, disks, and
+(later) atmospheres and images.
+
+The engine is forward-model agnostic. A fit is a `logdensity` over a partitioned
+scene PyTree, assembled from three plug-ins -- a forward model, a likelihood, and a
+prior -- and run through a uniform `Backend` (Laplace mixture, NUTS, adaptive
+tempered SMC, MCLMC, Pathfinder) that returns one `Posterior` exposing `.sample`,
+`.log_prob`, and `.evidence`.
+
+Status: early development. Orbit fitting is implemented; disk, atmosphere, and
+image-domain fitting are planned.

photomancy-0.0.1/README.md

@@ -0,0 +1,18 @@
+# photomancy
+
+> Divination from light.
+
+A JAX-native Bayesian inference and value-of-information engine for the Habitable
+Worlds Observatory direct-imaging simulation suite. `orbix` builds the geometry,
+`skyscapes` builds the scene, and **photomancy** divines the scene back from the
+data: posteriors, evidence, and the next-best observation, over orbits, disks, and
+(later) atmospheres and images.
+
+The engine is forward-model agnostic. A fit is a `logdensity` over a partitioned
+scene PyTree, assembled from three plug-ins -- a forward model, a likelihood, and a
+prior -- and run through a uniform `Backend` (Laplace mixture, NUTS, adaptive
+tempered SMC, MCLMC, Pathfinder) that returns one `Posterior` exposing `.sample`,
+`.log_prob`, and `.evidence`.
+
+Status: early development. Orbit fitting is implemented; disk, atmosphere, and
+image-domain fitting are planned.

photomancy-0.0.1/pyproject.toml

@@ -0,0 +1,74 @@
+[build-system]
+requires = ['hatchling', "hatch-fancy-pypi-readme", "hatch-vcs"]
+build-backend = "hatchling.build"
+
+[project]
+name = "photomancy"
+description = "Divination from light: Bayesian inference and value-of-information for HWO direct imaging"
+authors = [{ name = "Corey Spohn", email = "corey.a.spohn@nasa.gov" }]
+dynamic = ['readme', 'version']
+requires-python = ">=3.11"
+license = { file = "LICENSE" }
+classifiers = [
+    "Development Status :: 2 - Pre-Alpha",
+    "Intended Audience :: Science/Research",
+    "License :: OSI Approved :: MIT License",
+    "Programming Language :: Python :: 3",
+    "Topic :: Scientific/Engineering :: Astronomy",
+]
+dependencies = [
+    "blackjax>=1.2",
+    "equinox",
+    "hwoutils",
+    "jax>=0.8.1",
+    "jaxlib>=0.8.1",
+    "jaxtyping",
+    "numpy",
+    "numpyro>=0.15.0",
+    "orbix",
+    "skyscapes",
+]
+
+[project.optional-dependencies]
+dev = ["pre-commit"]
+docs = [
+    "sphinx",
+    "myst-nb",
+    "sphinx-book-theme",
+    "sphinx-autoapi",
+    "sphinx_autodoc_typehints",
+    "matplotlib",
+]
+test = ["nox", "pytest", "hypothesis", "pytest-cov", "arviz"]
+# Deferred extras (declared, not installed in the BlackJAX-only first cut):
+ns = ["jaxns"]          # nested-sampling evidence backend
+npe = ["flowjax"]       # amortized neural posterior estimation
+imaging = ["coronalyze"]  # differentiable detection/extraction likelihood helpers
+
+[project.urls]
+Homepage = "https://github.com/CoreySpohn/photomancy"
+Issues = "https://github.com/CoreySpohn/photomancy/issues"
+
+[tool.hatch.version]
+source = "vcs"
+
+[tool.hatch.build.hooks.vcs]
+version-file = "src/photomancy/_version.py"
+
+[tool.hatch.metadata.hooks.fancy-pypi-readme]
+content-type = "text/markdown"
+
+[[tool.hatch.metadata.hooks.fancy-pypi-readme.fragments]]
+path = "README.md"
+
+[tool.hatch.build.targets.wheel]
+packages = ["src/photomancy"]
+
+[tool.hatch.build.targets.sdist]
+exclude = ["/scripts", "/docs", "/tests", "/.github"]
+
+[tool.ruff.lint]
+select = ["B", "D", "E", "F", "I", "UP", "RUF"]
+
+[tool.ruff.lint.pydocstyle]
+convention = "google"

photomancy-0.0.1/src/photomancy/__init__.py

@@ -0,0 +1,19 @@
+"""photomancy: divination from light.
+
+A JAX-native Bayesian inference and value-of-information engine for fitting
+astrophysical scenes (orbits, disks, and later atmospheres and images) to
+direct-imaging data. orbix builds the geometry, skyscapes builds the scene, and
+photomancy divines the scene back from the data: posteriors, evidence, and the
+next-best observation.
+
+The engine is forward-model agnostic. A fit is a ``logdensity`` over a partitioned
+scene PyTree, assembled from a forward model, a likelihood, and a prior, and run
+through a uniform ``Backend`` that returns one ``Posterior``.
+"""
+
+try:
+    from photomancy._version import __version__
+except ImportError:  # pragma: no cover - generated by hatch-vcs at build time
+    __version__ = "0.0.0"
+
+__all__ = ["__version__"]

photomancy-0.0.1/src/photomancy/_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 = '0.0.1'
+__version_tuple__ = version_tuple = (0, 0, 1)
+
+__commit_id__ = commit_id = None

photomancy-0.0.1/src/photomancy/backends/__init__.py

@@ -0,0 +1,11 @@
+"""photomancy inference backends.
+
+Each backend wraps one inference method (Laplace today; NUTS, adaptive-tempered
+SMC, MCLMC, Pathfinder to come) behind a uniform ``run(logdensity, init, key)``
+that returns a Posterior.
+"""
+
+from photomancy.backends.base import AbstractBackend
+from photomancy.backends.laplace import LaplaceBackend
+
+__all__ = ["AbstractBackend", "LaplaceBackend"]

photomancy-0.0.1/src/photomancy/backends/base.py

@@ -0,0 +1,25 @@
+"""The Backend protocol: how inference is run on a logdensity."""
+
+from abc import abstractmethod
+
+import equinox as eqx
+
+
+class AbstractBackend(eqx.Module):
+    """Runs inference on a flat logdensity and returns a Posterior.
+
+    Backends are ``eqx.Module`` config objects (hyperparameters as fields) with a
+    pure ``run``. A backend sees only the flat ``logdensity`` -- never the scene or
+    the forward model -- which is what keeps the engine forward-model agnostic.
+    """
+
+    @abstractmethod
+    def run(self, logdensity, init, key=None):
+        """Run inference on ``logdensity`` from ``init``, returning a Posterior.
+
+        Args:
+            logdensity: ``z -> scalar`` log-density over the flat parameter position.
+            init: Initial flat position. Shape ``(d,)``.
+            key: PRNG key (unused by deterministic backends such as Laplace).
+        """
+        raise NotImplementedError

photomancy-0.0.1/src/photomancy/backends/laplace.py

@@ -0,0 +1,31 @@
+"""Laplace backend: MAP + Hessian -> GaussianPosterior."""
+
+import jax
+import jax.numpy as jnp
+from jax.scipy.optimize import minimize
+
+from photomancy.backends.base import AbstractBackend
+from photomancy.posterior import GaussianPosterior
+
+
+class LaplaceBackend(AbstractBackend):
+    """MAP + Laplace approximation: a Gaussian posterior around the mode.
+
+    Optimizes the logdensity to its MAP, takes the Hessian there for the
+    covariance, and returns a ``GaussianPosterior`` carrying the Laplace
+    log-evidence. Exact for Gaussian targets.
+    """
+
+    def run(self, logdensity, init, key=None):
+        """Fit the Laplace approximation and return a GaussianPosterior."""
+
+        def neg(z):
+            return -logdensity(z)
+
+        z_map = minimize(neg, init, method="BFGS").x
+        hess = jax.hessian(neg)(z_map)
+        cov = jnp.linalg.inv(hess)
+        d = z_map.shape[0]
+        _, logdet_cov = jnp.linalg.slogdet(cov)
+        log_z = logdensity(z_map) + 0.5 * d * jnp.log(2.0 * jnp.pi) + 0.5 * logdet_cov
+        return GaussianPosterior(mean=z_map, cov=cov, evidence=log_z)

photomancy-0.0.1/src/photomancy/core/__init__.py

@@ -0,0 +1,10 @@
+"""photomancy.core: the domain-agnostic inference engine.
+
+Assembles a logdensity over a partitioned scene PyTree from plug-in forward
+models, likelihoods, and priors. The partition / ravel boundary helpers and the
+Backend protocol + unified Posterior build on top of this.
+"""
+
+from photomancy.core.model import build_logdensity, build_scene_logdensity
+
+__all__ = ["build_logdensity", "build_scene_logdensity"]