FA-70232 / Map projection transforms / Member archive

Mollweide forward with Newton solve: newton derivative · case 02

Mid-latitude points jitter between two wrong positions depending on iteration parity.

Member previewVariant 2 · 3 implementations · 8 checks per implementation

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.

One recorded failure

Sample boundary fixture

This sample comes from the broken implementation of a controlled reproducer.

Boundary fixture	Actual	Expected	Outcome
control #1	`[5974515.86, 2662514.4]`	`[6102963.12, 1966970.83]`	Failed

MEMBER ARCHIVE

The complete case is available to members.

This record includes three runnable implementations, regression fixtures, execution results, and source hashes.

Member access is invitation-based. Sign in with your invited account to inspect the sources.