FAILURE MAP
← Case archive

FA-69636 / Orbital propagation / Open access

Stumpff functions for universal variables: Series branch is used far from zero · case 01

Moderate z values lose accuracy in the truncated series.

Verified by executionVariant 1 · 7 checks per implementationDownload source bundle ↓JSON ↗

ROOT CAUSE

The closed forms are used only for |z|>1.

VERIFIED REPAIR

Switch to closed forms outside |z|<=1e-6.

Unsuccessful approach: A 0.1 threshold still leaves truncation error in the reported digits.

Case contract

Input z. For z>1e-6: C=(1-cos s)/z, S=(s-sin s)/s^3 with s=sqrt z; for z<-1e-6: C=(cosh s-1)/(-z), S=(sinh s-s)/s^3 with s=sqrt(-z); otherwise series C=1/2-z/24+z^2/720, S=1/6-z/120+z^2/5040. Return [C,S] rounded to 10.

Why this case matters

Orbit determination and mission planning chain many small conversions; one wrong branch or unit silently moves a spacecraft by kilometres.

1 / The failure

Exit 1
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    z=x
    if z>1.0:
        s=math.sqrt(z)
        c=(1-math.cos(s))/z; S=(s-math.sin(s))/s**3
    elif z<-1.0:
        s=math.sqrt(-z)
        c=(math.cosh(s)-1)/(-z); S=(math.sinh(s)-s)/s**3
    else:
        c=1/2-z/24+z*z/720; S=1/6-z/120+z*z/5040
    return [round(c,10),round(S,10)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 30.0', 30.0, [0.0102526963, 0.0377242152]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658]), ('stumpff functions for universal variables 0.0', 0.0, [0.5, 0.1666666667]), ('stumpff functions for universal variables -1e-07', -1e-07, [0.5000000042, 0.1666666675])], [('stumpff functions for universal variables -0.5', -0.5, [0.521183673, 0.1708832825]), ('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658]), ('stumpff functions for universal variables 0.0', 0.0, [0.5, 0.1666666667]), ('stumpff functions for universal variables -1e-07', -1e-07, [0.5000000042, 0.1666666675]), ('stumpff functions for universal variables -4.0', -4.0, [0.6905489228, 0.203357551]), ('stumpff functions for universal variables -25.0', -25.0, [2.928397941, 0.5536256846])], [('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables -0.5', -0.5, [0.521183673, 0.1708832825]), ('stumpff functions for universal variables -4.0', -4.0, [0.6905489228, 0.203357551]), ('stumpff functions for universal variables -25.0', -25.0, [2.928397941, 0.5536256846]), ('stumpff functions for universal variables 9.8696', 9.8696, [0.2026424576, 0.1013212062]), ('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665])], [('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 9.8696', 9.8696, [0.2026424576, 0.1013212062]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665]), ('stumpff functions for universal variables -3e-06', -3e-06, [0.500000125, 0.1666666917]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603])], [('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665]), ('stumpff functions for universal variables -3e-06', -3e-06, [0.500000125, 0.1666666917]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 30.0', 30.0, [0.0102526963, 0.0377242152]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658])]]
for label, args, expected in fixtures[N-1]:
    check(label, solve(args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
stumpff functions for universal variables 0.5[0.4795138889, 0.1625496032][0.4795108058, 0.1625492603]Failed
stumpff functions for universal variables 0.05[0.4979201389, 0.166250496][0.4979201358, 0.1662504957]Failed
stumpff functions for universal variables 4.0[0.3540367091, 0.1363378216][0.3540367091, 0.1363378216]Passed
stumpff functions for universal variables 30.0[0.0102526963, 0.0377242152][0.0102526963, 0.0377242152]Passed
stumpff functions for universal variables 1e-07[0.4999999958, 0.1666666658][0.4999999958, 0.1666666658]Passed
stumpff functions for universal variables 0.0[0.5, 0.1666666667][0.5, 0.1666666667]Passed
stumpff functions for universal variables -1e-07[0.5000000042, 0.1666666675][0.5000000042, 0.1666666675]Passed

SHA-256 / e33c114b9cf740c58230dfcb3d40d0231562352840bd7b22d5bcf2470f2c8ef1

2 / The unsuccessful fix

Exit 1
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    z=x
    if z>0.1:
        s=math.sqrt(z)
        c=(1-math.cos(s))/z; S=(s-math.sin(s))/s**3
    elif z<-0.1:
        s=math.sqrt(-z)
        c=(math.cosh(s)-1)/(-z); S=(math.sinh(s)-s)/s**3
    else:
        c=1/2-z/24+z*z/720; S=1/6-z/120+z*z/5040
    return [round(c,10),round(S,10)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 30.0', 30.0, [0.0102526963, 0.0377242152]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658]), ('stumpff functions for universal variables 0.0', 0.0, [0.5, 0.1666666667]), ('stumpff functions for universal variables -1e-07', -1e-07, [0.5000000042, 0.1666666675])], [('stumpff functions for universal variables -0.5', -0.5, [0.521183673, 0.1708832825]), ('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658]), ('stumpff functions for universal variables 0.0', 0.0, [0.5, 0.1666666667]), ('stumpff functions for universal variables -1e-07', -1e-07, [0.5000000042, 0.1666666675]), ('stumpff functions for universal variables -4.0', -4.0, [0.6905489228, 0.203357551]), ('stumpff functions for universal variables -25.0', -25.0, [2.928397941, 0.5536256846])], [('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables -0.5', -0.5, [0.521183673, 0.1708832825]), ('stumpff functions for universal variables -4.0', -4.0, [0.6905489228, 0.203357551]), ('stumpff functions for universal variables -25.0', -25.0, [2.928397941, 0.5536256846]), ('stumpff functions for universal variables 9.8696', 9.8696, [0.2026424576, 0.1013212062]), ('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665])], [('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 9.8696', 9.8696, [0.2026424576, 0.1013212062]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665]), ('stumpff functions for universal variables -3e-06', -3e-06, [0.500000125, 0.1666666917]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603])], [('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665]), ('stumpff functions for universal variables -3e-06', -3e-06, [0.500000125, 0.1666666917]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 30.0', 30.0, [0.0102526963, 0.0377242152]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658])]]
for label, args, expected in fixtures[N-1]:
    check(label, solve(args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
stumpff functions for universal variables 0.5[0.4795108058, 0.1625492603][0.4795108058, 0.1625492603]Passed
stumpff functions for universal variables 0.05[0.4979201389, 0.166250496][0.4979201358, 0.1662504957]Failed
stumpff functions for universal variables 4.0[0.3540367091, 0.1363378216][0.3540367091, 0.1363378216]Passed
stumpff functions for universal variables 30.0[0.0102526963, 0.0377242152][0.0102526963, 0.0377242152]Passed
stumpff functions for universal variables 1e-07[0.4999999958, 0.1666666658][0.4999999958, 0.1666666658]Passed
stumpff functions for universal variables 0.0[0.5, 0.1666666667][0.5, 0.1666666667]Passed
stumpff functions for universal variables -1e-07[0.5000000042, 0.1666666675][0.5000000042, 0.1666666675]Passed

SHA-256 / 4b90f14f875f95ac8c1b79b10729208cbcb0ff767303af95fd721782895bca42

3 / The verified repair

Exit 0
"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    z=x
    if z>1e-6:
        s=math.sqrt(z)
        c=(1-math.cos(s))/z; S=(s-math.sin(s))/s**3
    elif z<-1e-6:
        s=math.sqrt(-z)
        c=(math.cosh(s)-1)/(-z); S=(math.sinh(s)-s)/s**3
    else:
        c=1/2-z/24+z*z/720; S=1/6-z/120+z*z/5040
    return [round(c,10),round(S,10)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 30.0', 30.0, [0.0102526963, 0.0377242152]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658]), ('stumpff functions for universal variables 0.0', 0.0, [0.5, 0.1666666667]), ('stumpff functions for universal variables -1e-07', -1e-07, [0.5000000042, 0.1666666675])], [('stumpff functions for universal variables -0.5', -0.5, [0.521183673, 0.1708832825]), ('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658]), ('stumpff functions for universal variables 0.0', 0.0, [0.5, 0.1666666667]), ('stumpff functions for universal variables -1e-07', -1e-07, [0.5000000042, 0.1666666675]), ('stumpff functions for universal variables -4.0', -4.0, [0.6905489228, 0.203357551]), ('stumpff functions for universal variables -25.0', -25.0, [2.928397941, 0.5536256846])], [('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables -0.5', -0.5, [0.521183673, 0.1708832825]), ('stumpff functions for universal variables -4.0', -4.0, [0.6905489228, 0.203357551]), ('stumpff functions for universal variables -25.0', -25.0, [2.928397941, 0.5536256846]), ('stumpff functions for universal variables 9.8696', 9.8696, [0.2026424576, 0.1013212062]), ('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665])], [('stumpff functions for universal variables -0.05', -0.05, [0.5020868087, 0.1670838297]), ('stumpff functions for universal variables 9.8696', 9.8696, [0.2026424576, 0.1013212062]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665]), ('stumpff functions for universal variables -3e-06', -3e-06, [0.500000125, 0.1666666917]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603])], [('stumpff functions for universal variables 0.5', 0.5, [0.4795108058, 0.1625492603]), ('stumpff functions for universal variables 0.05', 0.05, [0.4979201358, 0.1662504957]), ('stumpff functions for universal variables 2e-06', 2e-06, [0.4999999167, 0.16666665]), ('stumpff functions for universal variables -3e-06', -3e-06, [0.500000125, 0.1666666917]), ('stumpff functions for universal variables 4.0', 4.0, [0.3540367091, 0.1363378216]), ('stumpff functions for universal variables 30.0', 30.0, [0.0102526963, 0.0377242152]), ('stumpff functions for universal variables 1e-07', 1e-07, [0.4999999958, 0.1666666658])]]
for label, args, expected in fixtures[N-1]:
    check(label, solve(args), expected)
print(json.dumps({"observations": observations, "passed": all(x["passed"] for x in observations)}, ensure_ascii=False))
raise SystemExit(0 if all(x["passed"] for x in observations) else 1)
Boundary fixtureActualExpectedOutcome
stumpff functions for universal variables 0.5[0.4795108058, 0.1625492603][0.4795108058, 0.1625492603]Passed
stumpff functions for universal variables 0.05[0.4979201358, 0.1662504957][0.4979201358, 0.1662504957]Passed
stumpff functions for universal variables 4.0[0.3540367091, 0.1363378216][0.3540367091, 0.1363378216]Passed
stumpff functions for universal variables 30.0[0.0102526963, 0.0377242152][0.0102526963, 0.0377242152]Passed
stumpff functions for universal variables 1e-07[0.4999999958, 0.1666666658][0.4999999958, 0.1666666658]Passed
stumpff functions for universal variables 0.0[0.5, 0.1666666667][0.5, 0.1666666667]Passed
stumpff functions for universal variables -1e-07[0.5000000042, 0.1666666675][0.5000000042, 0.1666666675]Passed

SHA-256 / 30688794f626f12a813838cda158e3a82fdf456d674956e4938e23cb0e6c1190

Verification & scope

A deterministic toy two-body model with stipulated constants and conventions; not flight dynamics software or a validated SGP4 implementation. This reproducer isolates one failure mechanism. Results cover the supplied fixtures. Variants within a family share a test contract and should remain grouped when constructing evaluation splits. Related mechanisms with a shared evaluation_group must also remain together; these controlled models are not independent production incidents.

Observations recorded using Python 3.12.14 at 2026-09-29T14:48:13.223183+00:00.

Case digest / 04d4faec96077f8f101a4e85ac7833759f336025e268c927ef8a72ea251e48e7