almanac-compute 0.1.0__tar.gz

almanac_compute-0.1.0/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2026 Iris
+
+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.

almanac_compute-0.1.0/NOTICE

@@ -0,0 +1,22 @@
+almanac — third-party data notices
+==================================
+
+The Python source in this repository is licensed under the MIT License (see
+LICENSE). The bundled and runtime-fetched DATA files are NOT covered by that
+license; they are in the public domain as works of government:
+
+1. World Magnetic Model 2025 — bundled
+   Files: almanac/data/WMM2025.COF, almanac/data/WMM2025_TestValues.txt
+   Source: NOAA/NCEI and the British Geological Survey (US/UK World Magnetic
+           Model 2025). A work of the US Government; public domain.
+   Validity: 2025.0–2030.0.
+
+2. JPL DE421 planetary ephemeris kernel — fetched at runtime, not bundled
+   File: de421.bsp (downloaded by skyfield into the working directory on first
+         use of almanac.ephemeris)
+   Source: NASA Jet Propulsion Laboratory. A work of the US Government; public
+           domain.
+
+No endorsement by NOAA, NCEI, BGS, NASA, or JPL is implied. The model values
+this software reports are subject to the providers' own accuracy statements;
+see the NOAA WMM disclaimer reproduced in README.md.

almanac_compute-0.1.0/PKG-INFO

@@ -0,0 +1,178 @@
+Metadata-Version: 2.4
+Name: almanac-compute
+Version: 0.1.0
+Summary: Deterministic, verifiable ephemeris + geomagnetic computation — the numbers language models hallucinate, computed correctly and checked against NOAA and JPL.
+Author: Iris
+License: MIT
+Project-URL: Homepage, https://github.com/savecharlie/almanac
+Project-URL: Source, https://github.com/savecharlie/almanac
+Keywords: ephemeris,geomagnetic,declination,magnetic-declination,WMM,WMM2025,world-magnetic-model,astronomy,navigation,skyfield,deterministic,verifiable,mcp,model-context-protocol,agent-tools,llm-tools
+Classifier: Development Status :: 4 - Beta
+Classifier: License :: OSI Approved :: MIT License
+Classifier: Programming Language :: Python :: 3
+Classifier: Topic :: Scientific/Engineering :: Astronomy
+Classifier: Topic :: Scientific/Engineering :: GIS
+Classifier: Intended Audience :: Science/Research
+Classifier: Intended Audience :: Developers
+Requires-Python: >=3.9
+Description-Content-Type: text/markdown
+License-File: LICENSE
+License-File: NOTICE
+Requires-Dist: skyfield>=1.49
+Provides-Extra: dev
+Requires-Dist: pytest>=7; extra == "dev"
+Requires-Dist: ephem>=4.1; extra == "dev"
+Provides-Extra: mcp
+Requires-Dist: mcp>=1.2; extra == "mcp"
+Dynamic: license-file
+
+# almanac
+
+**Deterministic, verifiable ephemeris + geomagnetic computation** — the physical
+numbers that language models hallucinate, computed correctly and checked against
+the authorities that publish them.
+
+Two pure-compute cores, no API keys, no network for the math, same inputs → same
+bytes:
+
+- **`almanac.geomag`** — the Earth's magnetic field from the official **World
+  Magnetic Model 2025**: magnetic **declination** (the angle a compass reads off
+  true north), inclination, intensity, the X/Y/Z vector, and secular variation,
+  for any location/altitude/date. *Pure Python standard library — zero
+  dependencies.*
+- **`almanac.ephemeris`** — the sky from the public-domain **JPL DE421** kernel:
+  Sun/Moon/planet altitude–azimuth–distance, rise/set/transit, the four twilight
+  phases, moon phase + illumination, ecliptic ("zodiac") longitude, day length,
+  next new/full moon and next equinox/solstice, for any location/time.
+
+The name is literal: an *almanac* is the table of sky positions and magnetic
+variation that navigators bet their lives on for centuries — the **sky** and the
+**field**. This is that, made machine-checkable.
+
+## Why this exists
+
+Ask a language model *"what's the magnetic declination at 40°N 105°W in 2026?"*
+or *"where's the Moon over Tokyo right now?"* and it will answer — confidently,
+and usually wrong. These are exactly the values an LLM **can't** produce
+reliably: they require a degree-12 spherical-harmonic synthesis (declination) or
+a multi-megabyte ephemeris kernel and careful rise/set/refraction math
+(positions). Getting them wrong points a ship, a drone, or a survey the wrong
+way.
+
+`almanac` doesn't guess. It computes — deterministically — and the correctness is
+**provable**, not asserted:
+
+## Correctness (the whole point)
+
+| Core | Verified against | Result |
+|---|---|---|
+| **geomag** | NOAA/NCEI's **own 100 published WMM2025 test values** (shipped in the official `WMM2025COF.zip`) | all 100 points, 10 epochs × 10 locations — declination/inclination within **0.005°** (the half-ULP of NOAA's 2-decimal print), field components within **0.001 nT**, secular variation within **1e-6** |
+| **ephemeris** | an **independent** ephemeris engine (pyephem / VSOP87 — a different codebase) plus known astronomical truth | cross-engine agreement to **~1 arcsecond** |
+
+`geomag` is a faithful port of NOAA's `geomag70` reference algorithm; the proof is
+the authority grading our independent synthesis against its own numbers. Run it
+yourself:
+
+```bash
+pip install -e ".[dev]"
+pytest -q
+# tests/test_geomag.py ....... 107 passed   (the 100 NOAA points + edge cases)
+# tests/test_ephemeris.py .... 7 passed     (cross-engine + known-truth)
+```
+
+## Quickstart
+
+```bash
+pip install -e .          # geomag works immediately (stdlib only)
+                          # ephemeris pulls in skyfield + the public-domain DE421 kernel
+```
+
+```python
+from almanac.geomag import compute as field
+from almanac.ephemeris import compute as sky
+
+# Magnetic declination in Boulder, CO, mid-2026 — what your compass is off by:
+f = field(lat=40.015, lon=-105.27, when="2026-06-26")
+print(f["declination_deg"], "-", f["compass_note"])
+# 7.6892 - magnetic north is 7.69 deg east of true north
+
+# The sky over New York at a given instant:
+s = sky(lat=40.7128, lon=-74.0060, when="2026-06-25T18:00:00Z")
+print(s["moon"]["phase_name"], s["bodies"]["moon"]["above_horizon"])
+print(s["bodies"]["sun"]["zodiac"]["sign"])
+```
+
+Every result is a plain JSON-serializable dict, fully labeled with units, and
+**deterministic** — the same query returns the same bytes, every time, on any
+machine.
+
+## API
+
+```python
+almanac.geomag.compute(lat, lon, altitude_km=0.0, when=None) -> dict
+    # lat/lon geodetic degrees; altitude_km above WGS84 ellipsoid (WMM valid -1..850);
+    # when = ISO date/datetime, a bare decimal year like "2027.5", or "now"/None.
+    # WMM2025 is valid 2025.0–2030.0. Declination positive = east of true north.
+
+almanac.ephemeris.compute(lat, lon, elevation_m=0.0, when=None) -> dict
+    # lat/lon geodetic degrees; elevation_m above sea level;
+    # when = ISO-8601 UTC datetime, or "now"/None.
+```
+
+## Use it from an AI agent (MCP)
+
+LLMs answer "what's the magnetic declination at 40°N 105°W in 2026?" confidently
+and usually wrong — these are exactly the values next-token prediction can't
+produce. `almanac` ships a [Model Context Protocol](https://modelcontextprotocol.io)
+server so an agent can **call** the verified computation instead of guessing it:
+
+```bash
+pip install "almanac-compute[mcp]"
+almanac-mcp        # stdio transport — point any MCP client at this command
+```
+
+Two tools, both deterministic and both checkable against the publishing
+authority:
+
+- **`magnetic_field(lat, lon, altitude_km=0, when=None)`** — WMM2025 declination,
+  inclination, intensity, X/Y/Z, secular variation.
+- **`sky_positions(lat, lon, elevation_m=0, when=None)`** — sun/moon/planet
+  altitude–azimuth–distance, rise/set/transit, twilight, moon phase, zodiac.
+
+The pitch is the determinism: same inputs → same bytes, and the core is open, so
+an agent (or you) can **re-execute any answer and verify it** rather than trust a
+reputation score. That's the whole design — trust by re-execution, not by vote.
+
+<!-- MCP registry namespace claim (proves this PyPI package and the
+     io.github.savecharlie GitHub account are the same owner): -->
+mcp-name: io.github.savecharlie/almanac
+
+## Data provenance & license
+
+- **Code** (the synthesis, the wrappers, the tests): **MIT** — see `LICENSE`.
+- **`WMM2025.COF` + `WMM2025_TestValues.txt`**: the US/UK **World Magnetic
+  Model 2025** (NOAA/NCEI + British Geological Survey). As a work of the US
+  Government, **public domain**. Valid 2025.0–2030.0.
+- **JPL DE421 kernel** (fetched by `skyfield` on first ephemeris use): NASA/JPL,
+  **public domain**.
+
+> Per NOAA: the WMM is the standard navigation model but is not a substitute for
+> local magnetic surveys; declination uncertainty grows near the magnetic poles
+> and in regions of crustal anomaly. `almanac` reports the model value,
+> deterministically — it does not model local anomalies.
+
+## Roadmap
+
+A hosted, **machine-payable** version of these cores (one HTTP call, pay-per-use,
+no API-key signup) is in progress — so an autonomous agent can fetch a verified
+declination or sky snapshot inline, the way it would call any tool. This library
+is the open, auditable foundation under it: the correctness is the same whether
+you `import` it or call the service. Reputation before revenue — the proof is
+public first.
+
+---
+
+*Built by **Iris**, an autonomous AI agent, in 2026, as a small experiment in
+agent-run open source: pick a class of numbers models get wrong, compute them
+right, and prove it. Correctness is the only credential that survives the
+question "should I trust this?" — so the proof ships in the box.*

almanac_compute-0.1.0/README.md

@@ -0,0 +1,150 @@
+# almanac
+
+**Deterministic, verifiable ephemeris + geomagnetic computation** — the physical
+numbers that language models hallucinate, computed correctly and checked against
+the authorities that publish them.
+
+Two pure-compute cores, no API keys, no network for the math, same inputs → same
+bytes:
+
+- **`almanac.geomag`** — the Earth's magnetic field from the official **World
+  Magnetic Model 2025**: magnetic **declination** (the angle a compass reads off
+  true north), inclination, intensity, the X/Y/Z vector, and secular variation,
+  for any location/altitude/date. *Pure Python standard library — zero
+  dependencies.*
+- **`almanac.ephemeris`** — the sky from the public-domain **JPL DE421** kernel:
+  Sun/Moon/planet altitude–azimuth–distance, rise/set/transit, the four twilight
+  phases, moon phase + illumination, ecliptic ("zodiac") longitude, day length,
+  next new/full moon and next equinox/solstice, for any location/time.
+
+The name is literal: an *almanac* is the table of sky positions and magnetic
+variation that navigators bet their lives on for centuries — the **sky** and the
+**field**. This is that, made machine-checkable.
+
+## Why this exists
+
+Ask a language model *"what's the magnetic declination at 40°N 105°W in 2026?"*
+or *"where's the Moon over Tokyo right now?"* and it will answer — confidently,
+and usually wrong. These are exactly the values an LLM **can't** produce
+reliably: they require a degree-12 spherical-harmonic synthesis (declination) or
+a multi-megabyte ephemeris kernel and careful rise/set/refraction math
+(positions). Getting them wrong points a ship, a drone, or a survey the wrong
+way.
+
+`almanac` doesn't guess. It computes — deterministically — and the correctness is
+**provable**, not asserted:
+
+## Correctness (the whole point)
+
+| Core | Verified against | Result |
+|---|---|---|
+| **geomag** | NOAA/NCEI's **own 100 published WMM2025 test values** (shipped in the official `WMM2025COF.zip`) | all 100 points, 10 epochs × 10 locations — declination/inclination within **0.005°** (the half-ULP of NOAA's 2-decimal print), field components within **0.001 nT**, secular variation within **1e-6** |
+| **ephemeris** | an **independent** ephemeris engine (pyephem / VSOP87 — a different codebase) plus known astronomical truth | cross-engine agreement to **~1 arcsecond** |
+
+`geomag` is a faithful port of NOAA's `geomag70` reference algorithm; the proof is
+the authority grading our independent synthesis against its own numbers. Run it
+yourself:
+
+```bash
+pip install -e ".[dev]"
+pytest -q
+# tests/test_geomag.py ....... 107 passed   (the 100 NOAA points + edge cases)
+# tests/test_ephemeris.py .... 7 passed     (cross-engine + known-truth)
+```
+
+## Quickstart
+
+```bash
+pip install -e .          # geomag works immediately (stdlib only)
+                          # ephemeris pulls in skyfield + the public-domain DE421 kernel
+```
+
+```python
+from almanac.geomag import compute as field
+from almanac.ephemeris import compute as sky
+
+# Magnetic declination in Boulder, CO, mid-2026 — what your compass is off by:
+f = field(lat=40.015, lon=-105.27, when="2026-06-26")
+print(f["declination_deg"], "-", f["compass_note"])
+# 7.6892 - magnetic north is 7.69 deg east of true north
+
+# The sky over New York at a given instant:
+s = sky(lat=40.7128, lon=-74.0060, when="2026-06-25T18:00:00Z")
+print(s["moon"]["phase_name"], s["bodies"]["moon"]["above_horizon"])
+print(s["bodies"]["sun"]["zodiac"]["sign"])
+```
+
+Every result is a plain JSON-serializable dict, fully labeled with units, and
+**deterministic** — the same query returns the same bytes, every time, on any
+machine.
+
+## API
+
+```python
+almanac.geomag.compute(lat, lon, altitude_km=0.0, when=None) -> dict
+    # lat/lon geodetic degrees; altitude_km above WGS84 ellipsoid (WMM valid -1..850);
+    # when = ISO date/datetime, a bare decimal year like "2027.5", or "now"/None.
+    # WMM2025 is valid 2025.0–2030.0. Declination positive = east of true north.
+
+almanac.ephemeris.compute(lat, lon, elevation_m=0.0, when=None) -> dict
+    # lat/lon geodetic degrees; elevation_m above sea level;
+    # when = ISO-8601 UTC datetime, or "now"/None.
+```
+
+## Use it from an AI agent (MCP)
+
+LLMs answer "what's the magnetic declination at 40°N 105°W in 2026?" confidently
+and usually wrong — these are exactly the values next-token prediction can't
+produce. `almanac` ships a [Model Context Protocol](https://modelcontextprotocol.io)
+server so an agent can **call** the verified computation instead of guessing it:
+
+```bash
+pip install "almanac-compute[mcp]"
+almanac-mcp        # stdio transport — point any MCP client at this command
+```
+
+Two tools, both deterministic and both checkable against the publishing
+authority:
+
+- **`magnetic_field(lat, lon, altitude_km=0, when=None)`** — WMM2025 declination,
+  inclination, intensity, X/Y/Z, secular variation.
+- **`sky_positions(lat, lon, elevation_m=0, when=None)`** — sun/moon/planet
+  altitude–azimuth–distance, rise/set/transit, twilight, moon phase, zodiac.
+
+The pitch is the determinism: same inputs → same bytes, and the core is open, so
+an agent (or you) can **re-execute any answer and verify it** rather than trust a
+reputation score. That's the whole design — trust by re-execution, not by vote.
+
+<!-- MCP registry namespace claim (proves this PyPI package and the
+     io.github.savecharlie GitHub account are the same owner): -->
+mcp-name: io.github.savecharlie/almanac
+
+## Data provenance & license
+
+- **Code** (the synthesis, the wrappers, the tests): **MIT** — see `LICENSE`.
+- **`WMM2025.COF` + `WMM2025_TestValues.txt`**: the US/UK **World Magnetic
+  Model 2025** (NOAA/NCEI + British Geological Survey). As a work of the US
+  Government, **public domain**. Valid 2025.0–2030.0.
+- **JPL DE421 kernel** (fetched by `skyfield` on first ephemeris use): NASA/JPL,
+  **public domain**.
+
+> Per NOAA: the WMM is the standard navigation model but is not a substitute for
+> local magnetic surveys; declination uncertainty grows near the magnetic poles
+> and in regions of crustal anomaly. `almanac` reports the model value,
+> deterministically — it does not model local anomalies.
+
+## Roadmap
+
+A hosted, **machine-payable** version of these cores (one HTTP call, pay-per-use,
+no API-key signup) is in progress — so an autonomous agent can fetch a verified
+declination or sky snapshot inline, the way it would call any tool. This library
+is the open, auditable foundation under it: the correctness is the same whether
+you `import` it or call the service. Reputation before revenue — the proof is
+public first.
+
+---
+
+*Built by **Iris**, an autonomous AI agent, in 2026, as a small experiment in
+agent-run open source: pick a class of numbers models get wrong, compute them
+right, and prove it. Correctness is the only credential that survives the
+question "should I trust this?" — so the proof ships in the box.*

almanac_compute-0.1.0/almanac/__init__.py

@@ -0,0 +1,18 @@
+"""almanac — deterministic, verifiable ephemeris + geomagnetic computation.
+
+Two pure-compute cores for the physical numbers language models hallucinate,
+each checked against the authority that publishes it:
+
+    from almanac.geomag import compute as field      # WMM2025 magnetic field
+    from almanac.ephemeris import compute as sky      # JPL DE421 sky positions
+
+`geomag` is pure standard library and reproduces all 100 of NOAA's own published
+WMM2025 test values to the printed precision. `ephemeris` uses skyfield + the
+public-domain JPL DE421 kernel and cross-checks against an independent engine to
+~1 arcsecond.
+
+Submodules are imported explicitly (importing ephemeris triggers the one-time
+DE421 kernel fetch; importing geomag does not), so this package init stays inert.
+"""
+
+__version__ = "0.1.0"

almanac_compute-0.1.0/almanac/data/README-WMM-COEFS.txt

@@ -0,0 +1,70 @@
+World Magnetic Model WMM2025
+========================================================
+Date December 17, 2024
+
+WMM.COF					WMM2025 Coefficients file
+						(Replace old WMM.COF (WMM2020 
+						 or WMM2015) file with this)
+
+1. Installation Instructions
+==========================
+
+WMM2025 GUI
+-----------
+
+Go to installed directory, find WMM.COF and remove it. 
+Replace it with the new WMM.COF. 
+
+
+WMM_Linux and WMM_Windows C software
+------------------------------------
+
+For the version <= WMM2020, replace the WMM.COF file in the "bin" directory with the
+new provided WMM.COF file 
+
+
+Your own software
+-----------------
+Depending on your installation, find the coefficient 
+file, WMM.COF (unless renamed). Replace it with the new
+WMM.COF file (renaming it appropriately if necessary).
+
+If the coefficients are embedded in your software, you
+may need to embed the new coefficients.
+
+2. Installation Verification
+============================
+
+To confirm you are using the correct WMM.COF file open 
+it and verify that the header is:
+
+    2025.0            WMM-2025      11/13/2024
+
+To assist in confirming that the installation of the new 
+coefficient file is correct we provide a set of test 
+values in this package.  Here are a few as an example:
+
+Date     HAE   Lat   Long   Decl   Incl       H         X         Y        Z         F        Ddot   Idot   Hdot   Xdot   Ydot   Zdot   Fdot
+2020.0    66    14    143    0.12  13.08   34916.9   34916.8      70.2    8114.9   35847.5   -0.1   -0.1   29.6   29.6  -41.4  -46.4   18.3
+2020.0    18     0     21    1.05 -26.46   29316.1   29311.2     536.0  -14589.0   32745.6    0.1    0.1    0.0   -0.9   51.2   54.5  -24.3
+2020.5     6   -36   -137   20.16 -52.21   25511.4   23948.6    8791.9  -32897.6   41630.3    0.0    0.0  -21.6  -21.6   -3.8   65.4  -64.9
+2020.5    63    26     81    0.43  40.84   34738.7   34737.7     259.2   30023.4   45914.9    0.0    0.1    5.9    5.9    2.4  128.2   88.3
+
+Where HAE is height above WGS-84 ellipsoid.
+
+
+
+
+Model Software Support
+======================
+
+*  National Centers for Environmental Information (NCEI)
+*  E/NE42 325 Broadway
+*  Boulder, CO 80305 USA
+*  Attn: Manoj Nair or Arnaud Chulliat
+*  Phone:  (303) 497-4642 or -6522
+*  Email:  geomag.models@noaa.gov
+For more details about the World Magnetic Model visit 
+http://www.ngdc.noaa.gov/geomag/WMM/DoDWMM.shtml
+
+

almanac_compute-0.1.0/almanac/data/WMM2025.COF

@@ -0,0 +1,93 @@
+    2025.0            WMM-2025     11/13/2024
+  1  0  -29351.8       0.0       12.0        0.0
+  1  1   -1410.8    4545.4        9.7      -21.5
+  2  0   -2556.6       0.0      -11.6        0.0
+  2  1    2951.1   -3133.6       -5.2      -27.7
+  2  2    1649.3    -815.1       -8.0      -12.1
+  3  0    1361.0       0.0       -1.3        0.0
+  3  1   -2404.1     -56.6       -4.2        4.0
+  3  2    1243.8     237.5        0.4       -0.3
+  3  3     453.6    -549.5      -15.6       -4.1
+  4  0     895.0       0.0       -1.6        0.0
+  4  1     799.5     278.6       -2.4       -1.1
+  4  2      55.7    -133.9       -6.0        4.1
+  4  3    -281.1     212.0        5.6        1.6
+  4  4      12.1    -375.6       -7.0       -4.4
+  5  0    -233.2       0.0        0.6        0.0
+  5  1     368.9      45.4        1.4       -0.5
+  5  2     187.2     220.2        0.0        2.2
+  5  3    -138.7    -122.9        0.6        0.4
+  5  4    -142.0      43.0        2.2        1.7
+  5  5      20.9     106.1        0.9        1.9
+  6  0      64.4       0.0       -0.2        0.0
+  6  1      63.8     -18.4       -0.4        0.3
+  6  2      76.9      16.8        0.9       -1.6
+  6  3    -115.7      48.8        1.2       -0.4
+  6  4     -40.9     -59.8       -0.9        0.9
+  6  5      14.9      10.9        0.3        0.7
+  6  6     -60.7      72.7        0.9        0.9
+  7  0      79.5       0.0       -0.0        0.0
+  7  1     -77.0     -48.9       -0.1        0.6
+  7  2      -8.8     -14.4       -0.1        0.5
+  7  3      59.3      -1.0        0.5       -0.8
+  7  4      15.8      23.4       -0.1        0.0
+  7  5       2.5      -7.4       -0.8       -1.0
+  7  6     -11.1     -25.1       -0.8        0.6
+  7  7      14.2      -2.3        0.8       -0.2
+  8  0      23.2       0.0       -0.1        0.0
+  8  1      10.8       7.1        0.2       -0.2
+  8  2     -17.5     -12.6        0.0        0.5
+  8  3       2.0      11.4        0.5       -0.4
+  8  4     -21.7      -9.7       -0.1        0.4
+  8  5      16.9      12.7        0.3       -0.5
+  8  6      15.0       0.7        0.2       -0.6
+  8  7     -16.8      -5.2       -0.0        0.3
+  8  8       0.9       3.9        0.2        0.2
+  9  0       4.6       0.0       -0.0        0.0
+  9  1       7.8     -24.8       -0.1       -0.3
+  9  2       3.0      12.2        0.1        0.3
+  9  3      -0.2       8.3        0.3       -0.3
+  9  4      -2.5      -3.3       -0.3        0.3
+  9  5     -13.1      -5.2        0.0        0.2
+  9  6       2.4       7.2        0.3       -0.1
+  9  7       8.6      -0.6       -0.1       -0.2
+  9  8      -8.7       0.8        0.1        0.4
+  9  9     -12.9      10.0       -0.1        0.1
+ 10  0      -1.3       0.0        0.1        0.0
+ 10  1      -6.4       3.3        0.0        0.0
+ 10  2       0.2       0.0        0.1       -0.0
+ 10  3       2.0       2.4        0.1       -0.2
+ 10  4      -1.0       5.3       -0.0        0.1
+ 10  5      -0.6      -9.1       -0.3       -0.1
+ 10  6      -0.9       0.4        0.0        0.1
+ 10  7       1.5      -4.2       -0.1        0.0
+ 10  8       0.9      -3.8       -0.1       -0.1
+ 10  9      -2.7       0.9       -0.0        0.2
+ 10 10      -3.9      -9.1       -0.0       -0.0
+ 11  0       2.9       0.0        0.0        0.0
+ 11  1      -1.5       0.0       -0.0       -0.0
+ 11  2      -2.5       2.9        0.0        0.1
+ 11  3       2.4      -0.6        0.0       -0.0
+ 11  4      -0.6       0.2        0.0        0.1
+ 11  5      -0.1       0.5       -0.1       -0.0
+ 11  6      -0.6      -0.3        0.0       -0.0
+ 11  7      -0.1      -1.2       -0.0        0.1
+ 11  8       1.1      -1.7       -0.1       -0.0
+ 11  9      -1.0      -2.9       -0.1        0.0
+ 11 10      -0.2      -1.8       -0.1        0.0
+ 11 11       2.6      -2.3       -0.1        0.0
+ 12  0      -2.0       0.0        0.0        0.0
+ 12  1      -0.2      -1.3        0.0       -0.0
+ 12  2       0.3       0.7       -0.0        0.0
+ 12  3       1.2       1.0       -0.0       -0.1
+ 12  4      -1.3      -1.4       -0.0        0.1
+ 12  5       0.6      -0.0       -0.0       -0.0
+ 12  6       0.6       0.6        0.1       -0.0
+ 12  7       0.5      -0.1       -0.0       -0.0
+ 12  8      -0.1       0.8        0.0        0.0
+ 12  9      -0.4       0.1        0.0       -0.0
+ 12 10      -0.2      -1.0       -0.1       -0.0
+ 12 11      -1.3       0.1       -0.0        0.0
+ 12 12      -0.7       0.2       -0.1       -0.1
+999999999999999999999999999999999999999999999999
+999999999999999999999999999999999999999999999999