FA-70251 / Map projection transforms / Open access

Mollweide forward with Newton solve: northing scale · case 01

The map is squashed vertically.

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

ROOT CAUSE

The northing omits the sqrt(2) factor.

VERIFIED REPAIR

At the northing scale step restore `py = R * math.sqrt(2) * math.sin(theta)`, leaving the rest of the model unchanged.

Unsuccessful approach: sqrt(2) is restored but the sine is taken of the latitude instead of the auxiliary angle.

Case contract

Input [lon, lat, lon0]; sphere R = 6371000, dlon wrapped to [-180, 180). Solve 2t + sin 2t = pi sin(phi) by Newton from t = phi, at most 50 steps, stopping when the residual is below 1e-12, with derivative 2 + 2 cos 2t. x = R (2 sqrt 2 / pi) dlon cos t, y = R sqrt 2 sin t, rounded to 2 decimals.

Why this case matters

Global thematic maps use Mollweide for equal area; an under-converged auxiliary angle skews high latitudes.

1 / The failure

Exit 1

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    lon, lat, lon0 = x
    R = 6371000.0
    phi = math.radians(lat)
    dl = math.radians(((lon - lon0 + 180.0) % 360.0) - 180.0)
    theta = phi
    for _ in range(50):
        f = 2 * theta + math.sin(2 * theta) - math.pi * math.sin(phi)
        if abs(f) < 1e-12:
            break
        theta -= f / (2 + 2 * math.cos(2 * theta))
    px = R * 2 * math.sqrt(2) / math.pi * dl * math.cos(theta)
    py = R * math.sin(theta)
    return [round(px, 2), round(py, 2)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('control #0', [101.361, 11.005, 160.0], [-5803523.67, 1355983.62]), ('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #7', [114.805, 18.464, -90.0], [-15037610.8, 2265254.5]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #10', [-134.836, -65.77, 160.0], [3720190.56, -7401319.17]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #13', [-49.492, 24.945, 0.0], [-4663467.36, 3043363.54]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #16', [120.053, 58.23, -90.0], [-10039242.4, 6698546.76]), ('control #17', [164.893, -9.572, 10.0], [15372850.89, -1180095.58]), ('control #18', [35.53, 12.622, 10.0], [2517518.51, 1554061.2]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])]]
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 fixture	Actual	Expected	Outcome
control #0	`[-5803523.67, 958825.21]`	`[-5803523.67, 1355983.62]`	Failed
control #1	`[6102963.12, 1390858.41]`	`[6102963.12, 1966970.83]`	Failed
control #2	`[-7888806.65, -4726817.65]`	`[-7888806.65, -6684729.63]`	Failed
control #3	`[-4994731.15, -4719032.68]`	`[-4994731.15, -6673720.01]`	Failed
control #4	`[3649809.48, 5758872.18]`	`[3649809.48, 8144275.14]`	Failed
control #5	`[10593963.16, 4828485.2]`	`[10593963.16, 6828509.25]`	Failed
control #6	`[-5572506.91, -5500523.53]`	`[-5572506.91, -7778914.97]`	Failed
control #22	`[0.0, 0.0]`	`[0.0, 0.0]`	Passed

SHA-256 / a34cd36a9103a4d548485ff8c110d87b401d416f223b49f2478944a27ea95215

2 / The unsuccessful fix

Exit 1

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    lon, lat, lon0 = x
    R = 6371000.0
    phi = math.radians(lat)
    dl = math.radians(((lon - lon0 + 180.0) % 360.0) - 180.0)
    theta = phi
    for _ in range(50):
        f = 2 * theta + math.sin(2 * theta) - math.pi * math.sin(phi)
        if abs(f) < 1e-12:
            break
        theta -= f / (2 + 2 * math.cos(2 * theta))
    px = R * 2 * math.sqrt(2) / math.pi * dl * math.cos(theta)
    py = R * math.sqrt(2) * math.sin(phi)
    return [round(px, 2), round(py, 2)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('control #0', [101.361, 11.005, 160.0], [-5803523.67, 1355983.62]), ('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #7', [114.805, 18.464, -90.0], [-15037610.8, 2265254.5]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #10', [-134.836, -65.77, 160.0], [3720190.56, -7401319.17]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #13', [-49.492, 24.945, 0.0], [-4663467.36, 3043363.54]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #16', [120.053, 58.23, -90.0], [-10039242.4, 6698546.76]), ('control #17', [164.893, -9.572, 10.0], [15372850.89, -1180095.58]), ('control #18', [35.53, 12.622, 10.0], [2517518.51, 1554061.2]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])]]
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 fixture	Actual	Expected	Outcome
control #0	`[-5803523.67, 1719952.2]`	`[-5803523.67, 1355983.62]`	Failed
control #1	`[6102963.12, 2484387.02]`	`[6102963.12, 1966970.83]`	Failed
control #2	`[-7888806.65, -7648199.01]`	`[-7888806.65, -6684729.63]`	Failed
control #3	`[-4994731.15, -7638790.83]`	`[-4994731.15, -6673720.01]`	Failed
control #4	`[3649809.48, 8692524.52]`	`[3649809.48, 8144275.14]`	Failed
control #5	`[10593963.16, 7769298.51]`	`[10593963.16, 6828509.25]`	Failed
control #6	`[-5572506.91, -8475053.5]`	`[-5572506.91, -7778914.97]`	Failed
control #22	`[0.0, 0.0]`	`[0.0, 0.0]`	Passed

SHA-256 / b2818b424d86db9767e8d2de5b5d3e9ae44ef06eba73f7a8dd0552a71050dbd0

3 / The verified repair

Exit 0

"""Failure Map reference implementation. Python standard library only."""
import json
import math
N = 1
observations = []
def solve(x):
    lon, lat, lon0 = x
    R = 6371000.0
    phi = math.radians(lat)
    dl = math.radians(((lon - lon0 + 180.0) % 360.0) - 180.0)
    theta = phi
    for _ in range(50):
        f = 2 * theta + math.sin(2 * theta) - math.pi * math.sin(phi)
        if abs(f) < 1e-12:
            break
        theta -= f / (2 + 2 * math.cos(2 * theta))
    px = R * 2 * math.sqrt(2) / math.pi * dl * math.cos(theta)
    py = R * math.sqrt(2) * math.sin(theta)
    return [round(px, 2), round(py, 2)]
def check(label, actual, expected):
    observations.append({"check": label, "actual": actual, "expected": expected, "passed": actual == expected})
fixtures = [[('control #0', [101.361, 11.005, 160.0], [-5803523.67, 1355983.62]), ('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #1', [-137.531, 16.006, 160.0], [6102963.12, 1966970.83]), ('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #6', [-110.318, -70.158, 0.0], [-5572506.91, -7778914.97]), ('control #7', [114.805, 18.464, -90.0], [-15037610.8, 2265254.5]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #2', [152.471, -58.088, -90.0], [-7888806.65, -6684729.63]), ('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #8', [-149.861, -30.081, 160.0], [4589310.14, -3649245.65]), ('control #9', [101.273, 37.68, 160.0], [-5084573.46, 4523399.69]), ('control #10', [-134.836, -65.77, 160.0], [3720190.56, -7401319.17]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #3', [-64.263, -57.975, 10.0], [-4994731.15, -6673720.01]), ('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #11', [148.892, -11.799, 0.0], [14710485.98, -1453302.05]), ('control #12', [-27.35, -61.095, 10.0], [-2367959.82, -6972924.49]), ('control #13', [-49.492, 24.945, 0.0], [-4663467.36, 3043363.54]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])], [('control #4', [95.241, 74.746, 10.0], [3649809.48, 8144275.14]), ('control #5', [172.208, 59.576, 10.0], [10593963.16, 6828509.25]), ('control #14', [165.148, 70.359, 10.0], [7787627.95, 7795585.92]), ('control #15', [-71.762, 5.799, 160.0], [12797412.04, 715745.6]), ('control #16', [120.053, 58.23, -90.0], [-10039242.4, 6698546.76]), ('control #17', [164.893, -9.572, 10.0], [15372850.89, -1180095.58]), ('control #18', [35.53, 12.622, 10.0], [2517518.51, 1554061.2]), ('control #22', [0.0, 0.0, 0.0], [0.0, 0.0])]]
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 fixture	Actual	Expected	Outcome
control #0	`[-5803523.67, 1355983.62]`	`[-5803523.67, 1355983.62]`	Passed
control #1	`[6102963.12, 1966970.83]`	`[6102963.12, 1966970.83]`	Passed
control #2	`[-7888806.65, -6684729.63]`	`[-7888806.65, -6684729.63]`	Passed
control #3	`[-4994731.15, -6673720.01]`	`[-4994731.15, -6673720.01]`	Passed
control #4	`[3649809.48, 8144275.14]`	`[3649809.48, 8144275.14]`	Passed
control #5	`[10593963.16, 6828509.25]`	`[10593963.16, 6828509.25]`	Passed
control #6	`[-5572506.91, -7778914.97]`	`[-5572506.91, -7778914.97]`	Passed
control #22	`[0.0, 0.0]`	`[0.0, 0.0]`	Passed

SHA-256 / 650a2d4376dadaf054d34dfc3a3046ff49c1f0ebc3e15693c13a26a498842c57

Verification & scope

Stipulated deterministic toy contract on a bounded input domain; results are rounded as stated and no conformance with any published standard or library is claimed. 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:19.073855+00:00.

Case digest / 9ee2cdf58031a05c8fda51aa6b3ae991941f2a00898e33849e0ef60980715a54