FA-16791 / Floating-point arithmetic / Open access
Cosine-pi maps odd integers to the even extremum · case 01
Cosine-pi maps odd integers to the even extremum.
ROOT CAUSE
Cosine-pi maps odd integers to the even extremum. The faulty expression is if r==1: return '1'.
VERIFIED REPAIR
Apply the contract at this fault site using if r==1: return '-1'.
Unsuccessful approach: The attempted local correction if r==1: return '0' still violates the explicit regression fixtures.
Case contract
Evaluate cos(pi*x) for finite x after modulo-two reduction. Half integers return positive zero; integers are exact; near half integers use a sine of the distance to the half integer. 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):
try:
r=abs(math.remainder(x,2.0))
if r==0: return '1'
if r==1: return '1'
if r==0.5: return '0'
if r<0.5:
result=math.sin(math.pi*(0.5-r))
else:
result=-math.sin(math.pi*(r-0.5))
return render(result)
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('even integer', solve(2.0*N), "1")
check('odd integer', solve(2.0*N+1), "-1")
check('negative integer', solve(-2.0*N-1), "-1")
check('half integer', solve(2.0*N+0.5), "0")
check('quarter', solve(0.25), render(math.sqrt(0.5)))
check('three quarters', solve(0.75), render(-math.sqrt(0.5)))
check('near half below', solve(0.5-N*2.0**-50), render(math.sin(math.pi*N*2.0**-50)))
check('near half above', solve(0.5+N*2.0**-50), render(-math.sin(math.pi*N*2.0**-50)))
check('negative half', solve(-0.5), "0")
check('negative quarter', solve(-0.25), render(math.sqrt(0.5)))
check('large integer', solve(2.0**52+2*N), "1")
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 |
|---|---|---|---|
| even integer | 1 | 1 | Passed |
| odd integer | 1 | -1 | Failed |
| negative integer | 1 | -1 | Failed |
| half integer | 0 | 0 | Passed |
| quarter | 0.70710678119 | 0.70710678119 | Passed |
| three quarters | -0.70710678119 | -0.70710678119 | Passed |
| near half below | 2.7902947984e-15 | 2.7902947984e-15 | Passed |
| near half above | -2.7902947984e-15 | -2.7902947984e-15 | Passed |
| negative half | 0 | 0 | Passed |
| negative quarter | 0.70710678119 | 0.70710678119 | Passed |
| large integer | 1 | 1 | Passed |
SHA-256 / b077f0b48591e58070bf03a9e50069c81bfa11b0dcf89c233bbc2a624671adba
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):
try:
r=abs(math.remainder(x,2.0))
if r==0: return '1'
if r==1: return '0'
if r==0.5: return '0'
if r<0.5:
result=math.sin(math.pi*(0.5-r))
else:
result=-math.sin(math.pi*(r-0.5))
return render(result)
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('even integer', solve(2.0*N), "1")
check('odd integer', solve(2.0*N+1), "-1")
check('negative integer', solve(-2.0*N-1), "-1")
check('half integer', solve(2.0*N+0.5), "0")
check('quarter', solve(0.25), render(math.sqrt(0.5)))
check('three quarters', solve(0.75), render(-math.sqrt(0.5)))
check('near half below', solve(0.5-N*2.0**-50), render(math.sin(math.pi*N*2.0**-50)))
check('near half above', solve(0.5+N*2.0**-50), render(-math.sin(math.pi*N*2.0**-50)))
check('negative half', solve(-0.5), "0")
check('negative quarter', solve(-0.25), render(math.sqrt(0.5)))
check('large integer', solve(2.0**52+2*N), "1")
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 |
|---|---|---|---|
| even integer | 1 | 1 | Passed |
| odd integer | 0 | -1 | Failed |
| negative integer | 0 | -1 | Failed |
| half integer | 0 | 0 | Passed |
| quarter | 0.70710678119 | 0.70710678119 | Passed |
| three quarters | -0.70710678119 | -0.70710678119 | Passed |
| near half below | 2.7902947984e-15 | 2.7902947984e-15 | Passed |
| near half above | -2.7902947984e-15 | -2.7902947984e-15 | Passed |
| negative half | 0 | 0 | Passed |
| negative quarter | 0.70710678119 | 0.70710678119 | Passed |
| large integer | 1 | 1 | Passed |
SHA-256 / bd3fcdc9d68f49193b3b9468de611343e90d3f2b22831cdb2d450c1916124c42
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):
try:
r=abs(math.remainder(x,2.0))
if r==0: return '1'
if r==1: return '-1'
if r==0.5: return '0'
if r<0.5:
result=math.sin(math.pi*(0.5-r))
else:
result=-math.sin(math.pi*(r-0.5))
return render(result)
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('even integer', solve(2.0*N), "1")
check('odd integer', solve(2.0*N+1), "-1")
check('negative integer', solve(-2.0*N-1), "-1")
check('half integer', solve(2.0*N+0.5), "0")
check('quarter', solve(0.25), render(math.sqrt(0.5)))
check('three quarters', solve(0.75), render(-math.sqrt(0.5)))
check('near half below', solve(0.5-N*2.0**-50), render(math.sin(math.pi*N*2.0**-50)))
check('near half above', solve(0.5+N*2.0**-50), render(-math.sin(math.pi*N*2.0**-50)))
check('negative half', solve(-0.5), "0")
check('negative quarter', solve(-0.25), render(math.sqrt(0.5)))
check('large integer', solve(2.0**52+2*N), "1")
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 |
|---|---|---|---|
| even integer | 1 | 1 | Passed |
| odd integer | -1 | -1 | Passed |
| negative integer | -1 | -1 | Passed |
| half integer | 0 | 0 | Passed |
| quarter | 0.70710678119 | 0.70710678119 | Passed |
| three quarters | -0.70710678119 | -0.70710678119 | Passed |
| near half below | 2.7902947984e-15 | 2.7902947984e-15 | Passed |
| near half above | -2.7902947984e-15 | -2.7902947984e-15 | Passed |
| negative half | 0 | 0 | Passed |
| negative quarter | 0.70710678119 | 0.70710678119 | Passed |
| large integer | 1 | 1 | Passed |
SHA-256 / 41416c11b6aa6861cfaa3b93f83ae6b9d6dd75bceb32e3074e0a4dba029dbad1
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:39.881604+00:00.
Case digest / 7ef6edf8bcd6636149e746600d07cfa4029edd040330499e3e15e2bdc1b1a5e2