FA-17276 / Floating-point arithmetic / Open access
Complex exponential subtracts a binary exponent as though natural · case 01
Complex exponential subtracts a binary exponent as though natural.
ROOT CAUSE
Complex exponential subtracts a binary exponent as though natural. The faulty expression is f=x-k.
VERIFIED REPAIR
Apply the contract at this fault site using f=x-k*math.log(2).
Unsuccessful approach: The attempted local correction f=x/k if k else x still violates the explicit regression fixtures.
Case contract
Compute exp(x+iy) for finite components using binary exponent splitting. Fixtures keep each final component finite or zero, including cases where exp(x) alone overflows. Return two rendered components; exact trigonometric zeros are preserved. Finite results are rendered to eleven significant decimal digits; modeled domain violations and arithmetic errors are explicit strings.
Why this case matters
An offline floating representation model isolates a reproducible arithmetic fault.
1 / The failure
Exit 1"""Failure Map reference implementation. Python standard library only."""
import json
import math
import struct
def render(x):
if math.isnan(x): return 'nan'
if math.isinf(x): return '-infinity' if x<0 else '+infinity'
return format(x,'.11g')
N = 1
observations = []
def solve(x,y):
try:
k=math.floor(x/math.log(2))
f=x-k
scale=math.exp(f)
c=math.cos(y); s=math.sin(y)
real=math.ldexp(scale*c,k)
imag=math.ldexp(scale*s,k)
return [render(real),render(imag)]
except (ValueError, OverflowError, ZeroDivisionError, TypeError):
return "arithmetic-error"
def check(label, actual, expected):
observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
check('overflowing common factor', solve(710.0,math.pi/4), [render(math.exp(710-math.log(2))*(2*math.cos(math.pi/4))),render(math.exp(710-math.log(2))*(2*math.sin(math.pi/4)))])
check('negative real projection', solve(float(N),2.0), [render(math.exp(N)*math.cos(2)),render(math.exp(N)*math.sin(2))])
check('ordinary', solve(float(N),0.25), [render(math.exp(N)*math.cos(0.25)),render(math.exp(N)*math.sin(0.25))])
check('negative phase', solve(float(N),-0.25), [render(math.exp(N)*math.cos(-0.25)),render(math.exp(N)*math.sin(-0.25))])
check('negative real', solve(-float(N),0.25), [render(math.exp(-N)*math.cos(0.25)),render(math.exp(-N)*math.sin(0.25))])
check('negative zero phase', solve(float(N),-0.0), [render(math.exp(N)),"-0"])
check('zero', solve(0.0,0.0), ["1","0"])
check('tiny common factor', solve(-744.0,0.25), [render(math.exp(-744)*math.cos(0.25)),render(math.exp(-744)*math.sin(0.25))])
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 fixture | Actual | Expected | Outcome |
|---|---|---|---|
| overflowing common factor | ['5.4416830713e+171', '5.4416830713e+171'] | ['1.5796728483e+308', '1.5796728483e+308'] | Failed |
| negative real projection | ['-0.83229367309', '1.8185948537'] | ['-1.1312043838', '2.471726672'] | Failed |
| ordinary | ['1.9378248434', '0.49480791851'] | ['2.6337770293', '0.67251368673'] | Failed |
| negative phase | ['1.9378248434', '-0.49480791851'] | ['2.6337770293', '-0.67251368673'] | Failed |
| negative real | ['0.65844425733', '0.16812842168'] | ['0.35644296024', '0.091014830274'] | Failed |
| negative zero phase | ['2', '-0'] | ['2.7182818285', '-0'] | Failed |
| zero | ['1', '0'] | ['1', '0'] | Passed |
| tiny common factor | ['9.936837979e-181', '2.5372913004e-181'] | ['9.8813129168e-324', '0'] | Failed |
SHA-256 / 422b92bc28a797426fb049b47cbb89d5b928226feb51de61802e4e3e8979620e
2 / The unsuccessful fix
Exit 1"""Failure Map reference implementation. Python standard library only."""
import json
import math
import struct
def render(x):
if math.isnan(x): return 'nan'
if math.isinf(x): return '-infinity' if x<0 else '+infinity'
return format(x,'.11g')
N = 1
observations = []
def solve(x,y):
try:
k=math.floor(x/math.log(2))
f=x/k if k else x
scale=math.exp(f)
c=math.cos(y); s=math.sin(y)
real=math.ldexp(scale*c,k)
imag=math.ldexp(scale*s,k)
return [render(real),render(imag)]
except (ValueError, OverflowError, ZeroDivisionError, TypeError):
return "arithmetic-error"
def check(label, actual, expected):
observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
check('overflowing common factor', solve(710.0,math.pi/4), [render(math.exp(710-math.log(2))*(2*math.cos(math.pi/4))),render(math.exp(710-math.log(2))*(2*math.sin(math.pi/4)))])
check('negative real projection', solve(float(N),2.0), [render(math.exp(N)*math.cos(2)),render(math.exp(N)*math.sin(2))])
check('ordinary', solve(float(N),0.25), [render(math.exp(N)*math.cos(0.25)),render(math.exp(N)*math.sin(0.25))])
check('negative phase', solve(float(N),-0.25), [render(math.exp(N)*math.cos(-0.25)),render(math.exp(N)*math.sin(-0.25))])
check('negative real', solve(-float(N),0.25), [render(math.exp(-N)*math.cos(0.25)),render(math.exp(-N)*math.sin(0.25))])
check('negative zero phase', solve(float(N),-0.0), [render(math.exp(N)),"-0"])
check('zero', solve(0.0,0.0), ["1","0"])
check('tiny common factor', solve(-744.0,0.25), [render(math.exp(-744)*math.cos(0.25)),render(math.exp(-744)*math.sin(0.25))])
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 fixture | Actual | Expected | Outcome |
|---|---|---|---|
| overflowing common factor | arithmetic-error | ['1.5796728483e+308', '1.5796728483e+308'] | Failed |
| negative real projection | ['-2.2624087675', '4.943453344'] | ['-1.1312043838', '2.471726672'] | Failed |
| ordinary | ['5.2675540586', '1.3450273735'] | ['2.6337770293', '0.67251368673'] | Failed |
| negative phase | ['5.2675540586', '-1.3450273735'] | ['2.6337770293', '-0.67251368673'] | Failed |
| negative real | ['0.39936662978', '0.10197504252'] | ['0.35644296024', '0.091014830274'] | Failed |
| negative zero phase | ['5.4365636569', '-0'] | ['2.7182818285', '-0'] | Failed |
| zero | ['1', '0'] | ['1', '0'] | Passed |
| tiny common factor | ['9.8813129168e-324', '0'] | ['9.8813129168e-324', '0'] | Passed |
SHA-256 / 481c8f97e5ddf7749b010d9938fa5d5f7468e10ed3a64e67af992b9d156dc44b
3 / The verified repair
Exit 0"""Failure Map reference implementation. Python standard library only."""
import json
import math
import struct
def render(x):
if math.isnan(x): return 'nan'
if math.isinf(x): return '-infinity' if x<0 else '+infinity'
return format(x,'.11g')
N = 1
observations = []
def solve(x,y):
try:
k=math.floor(x/math.log(2))
f=x-k*math.log(2)
scale=math.exp(f)
c=math.cos(y); s=math.sin(y)
real=math.ldexp(scale*c,k)
imag=math.ldexp(scale*s,k)
return [render(real),render(imag)]
except (ValueError, OverflowError, ZeroDivisionError, TypeError):
return "arithmetic-error"
def check(label, actual, expected):
observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
check('overflowing common factor', solve(710.0,math.pi/4), [render(math.exp(710-math.log(2))*(2*math.cos(math.pi/4))),render(math.exp(710-math.log(2))*(2*math.sin(math.pi/4)))])
check('negative real projection', solve(float(N),2.0), [render(math.exp(N)*math.cos(2)),render(math.exp(N)*math.sin(2))])
check('ordinary', solve(float(N),0.25), [render(math.exp(N)*math.cos(0.25)),render(math.exp(N)*math.sin(0.25))])
check('negative phase', solve(float(N),-0.25), [render(math.exp(N)*math.cos(-0.25)),render(math.exp(N)*math.sin(-0.25))])
check('negative real', solve(-float(N),0.25), [render(math.exp(-N)*math.cos(0.25)),render(math.exp(-N)*math.sin(0.25))])
check('negative zero phase', solve(float(N),-0.0), [render(math.exp(N)),"-0"])
check('zero', solve(0.0,0.0), ["1","0"])
check('tiny common factor', solve(-744.0,0.25), [render(math.exp(-744)*math.cos(0.25)),render(math.exp(-744)*math.sin(0.25))])
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 fixture | Actual | Expected | Outcome |
|---|---|---|---|
| overflowing common factor | ['1.5796728483e+308', '1.5796728483e+308'] | ['1.5796728483e+308', '1.5796728483e+308'] | Passed |
| negative real projection | ['-1.1312043838', '2.471726672'] | ['-1.1312043838', '2.471726672'] | Passed |
| ordinary | ['2.6337770293', '0.67251368673'] | ['2.6337770293', '0.67251368673'] | Passed |
| negative phase | ['2.6337770293', '-0.67251368673'] | ['2.6337770293', '-0.67251368673'] | Passed |
| negative real | ['0.35644296024', '0.091014830274'] | ['0.35644296024', '0.091014830274'] | Passed |
| negative zero phase | ['2.7182818285', '-0'] | ['2.7182818285', '-0'] | Passed |
| zero | ['1', '0'] | ['1', '0'] | Passed |
| tiny common factor | ['9.8813129168e-324', '0'] | ['9.8813129168e-324', '0'] | Passed |
SHA-256 / 89198c78290ab9a83ddcd2c3e75088ee557eeb226ffc27747f970eb6f8bb402b
Verification & scope
Controlled binary64 or explicitly stipulated miniature format; no hardware exception flags or platform floating environment are modeled. 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:39:44.730375+00:00.
Case digest / 030b7fdce425ccc584e678eedbc793e2c8e0749bf662656b7b2327d204dc7132