Vincenty's formulae¶

Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a).

Version: v1 · 2026-09-28 · History
Domain-specific #: 12796
Domain group: Natural Sciences
Origin domain: Geology & Earth Sciences
Subdomain: Geodesy → Geology & Earth Sciences

Core Idea¶

Vincenty's formulae is treated here as the recurring cross-domain formal modeling identity summarized by this source-grounded definition: Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a).

Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a). They are based on the assumption that the figure of the Earth is an oblate spheroid, and hence are more accurate than methods that assume a spherical Earth, such as great-circle distance. The first (direct) method computes the location of a point that is a given distance and azimuth (direction) from another point.

The second (inverse) method computes the geographical distance and azimuth between two given points. They have been widely used in geodesy because they are accurate to within 0.5 mm (0.020in) on the Earth ellipsoid. If the standard 2-argument arctangent atan2 function is used, then these values are usually handled correctly.

For Vincenty's formulae, the abstraction is narrower than the article's general subject matter: a positive case must preserve Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a). Retaining only the name, a familiar example, or a downstream effect is insufficient. The specialist roles and tests remain anchored in cross-domain formal modeling, which is why this identity is domain-specific rather than prime.

Structural Signature¶

Sig role-phrases:

Defining carrier — Between two nearly antipodal points, the iterative formula may fail to converge; this will occur when the first guess at λ as computed by the equation above is greater than π in absolute value.
Constitutive relation — Legendre showed that an ellipsoidal geodesic can be exactly mapped to a great circle on the auxiliary sphere by mapping the geographic latitude to reduced latitude and setting the azimuth of the great circle equal to that of the geodesic.
Operating condition — The longitude on the ellipsoid and the distance along the geodesic are then given in terms of the longitude on the sphere and the arc length along the great circle by simple integrals.
Recognition evidence — Bessel and Helmert gave rapidly converging series for these integrals, which allow the geodesic to be computed with arbitrary accuracy.
Admissible variation — The expressions were put in Horner (or nested) form, since this allows polynomials to be evaluated using only a single temporary register.
Characteristic consequence — This requires about 130 iterations to give a result accurate to 1 mm.
Failure boundary — An example of an incorrect result is provided by the NGS online utility, which returns a distance that is about 5 km too long.

What It Is Not¶

Not the whole field of cross-domain formal modeling. The node requires the specific identity stated by Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a).
Not an over-broad reading. This converges to the correct result 19944127.421 m after about 60 iterations; however, in other cases many thousands of iterations are required.
Not an over-broad reading. Finally, simple iterative techniques were used to solve the implicit equations in the direct and inverse methods; even though these are slow (and in the case of the inverse method it sometimes does not converge), they result in the least increase in code size.
Not an over-broad reading. Vincenty's goal was to express existing algorithms for geodesics on an ellipsoid in a form that minimized the program length (Vincenty 1975a).
Not automatically Haversine Formula. Retrieval proximity does not establish equivalence; the two identities must be compared by carrier, operation, and failure boundary.

Scope of Application¶

Vincenty's formulae applies literally inside cross-domain formal modeling wherever the source-defined carrier and relation can be established. Its documented habitats include:

Direct problem. If the standard 2-argument arctangent atan2 function is used, then these values are usually handled correctly.
Background. Finally, simple iterative techniques were used to solve the implicit equations in the direct and inverse methods; even though these are slow (and in the case of the inverse method it sometimes does not converge), they result in the least increase in code size.
Documented setting. Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a).
Background. The expressions were put in Horner (or nested) form, since this allows polynomials to be evaluated using only a single temporary register.
Nearly antipodal points. Depending on how the inverse method is implemented, the algorithm might return the correct result (19936288.579 m), an incorrect result, or an error indicator.
Nearly antipodal points. Vincenty suggested a method of accelerating the convergence in such cases (Rapp, 1993).

Outside cross-domain formal modeling, the name should be retained only when these same operational conditions survive; otherwise the comparison belongs to the broader parent Representation or should be marked as analogy.

Clarity¶

A clear use of Vincenty's formulae names the carrier, the operative relation, and the conditions under which the source treats the identity as present. The minimal definition is Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a). The strongest recognition evidence in the frozen account is: Bessel and Helmert gave rapidly converging series for these integrals, which allow the geodesic to be computed with arbitrary accuracy. A report should distinguish that evidence from a proxy, consequence, or common implementation. It should also state the qualification This converges to the correct result 19944127.421 m after about 60 iterations; however, in other cases many thousands of iterations are required. so that a reader can reproduce the classification rather than infer it from topical resemblance.

Manages Complexity¶

Vincenty's formulae compresses multiple cross-domain formal modeling details into a stable diagnostic relation. The source shows both the central mechanism—legendre showed that an ellipsoidal geodesic can be exactly mapped to a great circle on the auxiliary sphere by mapping the geographic latitude to reduced latitude and setting the azimuth of the great circle equal to that of the geodesic.—and the practical consequence—this requires about 130 iterations to give a result accurate to 1 mm. This compression makes cases comparable while leaving parameters, conventions, exceptions, and evidential quality explicit. It is lossy by design: local history and implementation details may be omitted only when they do not alter the defining relation.

Abstract Reasoning¶

Type the carrier. Identify the cross-domain formal modeling entities to which the claim applies.
State the relation. Use the source-grounded identity: Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a).
Check operation and conditions. The longitude on the ellipsoid and the distance along the geodesic are then given in terms of the longitude on the sphere and the arc length along the great circle by simple integrals.
Demand recognition evidence. Bessel and Helmert gave rapidly converging series for these integrals, which allow the geodesic to be computed with arbitrary accuracy.
Test variation. Change an implementation or setting while preserving the expressions were put in Horner (or nested) form, since this allows polynomials to be evaluated using only a single temporary register.
Run the collapse test. Remove the defining operation; if the label still seems equally apt, only a topic or correlate was retained.
Reduce cautiously. When the specialist conditions cannot be carried, route the residual comparison to Representation.

Knowledge Transfer¶

Within the home domain. Knowledge about Vincenty's formulae transfers literally when a new case preserves the same carrier type, relation, and recognition test. If the standard 2-argument arctangent atan2 function is used, then these values are usually handled correctly. Finally, simple iterative techniques were used to solve the implicit equations in the direct and inverse methods; even though these are slow (and in the case of the inverse method it sometimes does not converge), they result in the least increase in code size.

Beyond the home domain. No canonical parent is asserted for Vincenty's formulae. An outside case receives the specialist name only when the same typed roles and rejection conditions can be filled literally; otherwise the comparison remains an analogy pending later graph densification.

Examples¶

Canonical¶

Vincenty suggested a method of accelerating the convergence in such cases (Rapp, 1993). This case is canonical because it supplies a concrete carrier and lets the defining relation be checked rather than merely named.

Mapped back: carrier → the entities in the documented case; operation → Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a); recognition evidence → Bessel and Helmert gave rapidly converging series for these integrals, which allow the geodesic to be computed with arbitrary accuracy

Applied / In Practice¶

In an unpublished report, Vincenty (1975b) gave an alternative iterative scheme to handle such cases. The applied case shows how the identity is used under a second setting or qualification while keeping the same operative relation.

Mapped back: changed setting → Nearly antipodal points; invariant → Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a); boundary → the case exits the class when this converges to the correct result 19944127.421 m after about 60 iterations; however, in other cases many thousands of iterations are required

Structural Tensions¶

T1 — Stable identity versus admissible variation. This converges to the correct result 19944127.421 m after about 60 iterations; however, in other cases many thousands of iterations are required. The tension matters because emphasizing only one side either dissolves the identity or overstates what the evidence and domain conventions warrant.

Diagnostic: Which changes preserve the defining relation, and which replace it?

T2 — Recognition versus proxy. Finally, simple iterative techniques were used to solve the implicit equations in the direct and inverse methods; even though these are slow (and in the case of the inverse method it sometimes does not converge), they result in the least increase in code size. The tension matters because emphasizing only one side either dissolves the identity or overstates what the evidence and domain conventions warrant.

Diagnostic: Does the cited evidence establish the identity or only a correlated sign?

T3 — Definition versus implementation. Vincenty's goal was to express existing algorithms for geodesics on an ellipsoid in a form that minimized the program length (Vincenty 1975a). The tension matters because emphasizing only one side either dissolves the identity or overstates what the evidence and domain conventions warrant.

Diagnostic: Is the observed implementation constitutive, optional, or merely common?

T4 — Scope versus overextension. His unpublished report (1975b) mentions the use of a Wang 720 desk calculator, which had only a few kilobytes of memory. The tension matters because emphasizing only one side either dissolves the identity or overstates what the evidence and domain conventions warrant.

Diagnostic: Can every claimed application fill the same typed roles without metaphor?

T5 — Transfer versus domain accent. Between two nearly antipodal points, the iterative formula may fail to converge; this will occur when the first guess at λ as computed by the equation above is greater than π in absolute value. The tension matters because emphasizing only one side either dissolves the identity or overstates what the evidence and domain conventions warrant.

Diagnostic: Does the receiving case instantiate Vincenty's formulae literally, co-instantiate Representation, or only resemble it?

T6 — Autonomy versus reduction. Legendre showed that an ellipsoidal geodesic can be exactly mapped to a great circle on the auxiliary sphere by mapping the geographic latitude to reduced latitude and setting the azimuth of the great circle equal to that of the geodesic. The tension matters because emphasizing only one side either dissolves the identity or overstates what the evidence and domain conventions warrant.

Diagnostic: What does Vincenty's formulae distinguish that the broader parent Representation leaves together?

Structural–Framed Character¶

Vincenty's formulae is mixed or framed-leaning. Its structural side is the repeatable organization summarized by Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a). Its framed side is the cross-domain formal modeling vocabulary that fixes the carrier, evidence, exceptions, and admissible transformations.

Evaluative weight: the identity can be stated descriptively even when applications carry practical stakes. Human-practice dependence: the source-grounded carrier determines whether the relation exists independently or is constituted by a practice. Institutional origin: disciplinary conventions stabilize the name and test. Vocabulary portability: The longitude on the ellipsoid and the distance along the geodesic are then given in terms of the longitude on the sphere and the arc length along the great circle by simple integrals. Import versus recognition: literal transfer requires the same mechanism; shape alone is analogy.

Its portable skeleton is Representation. Its character: a recurring specialist identity whose thin organization can be abstracted, while its operational meaning remains domain-bound.

Structural Core vs. Domain Accent¶

What is skeletal. Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a). The stable skeleton is the typed relation expressed in that definition and the entry's recognition and collapse tests. The source identifies these operative conditions: Between two nearly antipodal points, the iterative formula may fail to converge; this will occur when the first guess at λ as computed by the equation above is greater than π in absolute value. Legendre showed that an ellipsoidal geodesic can be exactly mapped to a great circle on the auxiliary sphere by mapping the geographic latitude to reduced latitude and setting the azimuth of the great circle equal to that of the geodesic. It further constrains recognition and variation through: The longitude on the ellipsoid and the distance along the geodesic are then given in terms of the longitude on the sphere and the arc length along the great circle by simple integrals. Bessel and Helmert gave rapidly converging series for these integrals, which allow the geodesic to be computed with arbitrary accuracy.

What is domain-bound. cross-domain formal modeling supplies the operative entities, technical vocabulary, warrants, and exceptions that make Vincenty's formulae literal. Its documented scope includes the condition that If the standard 2-argument arctangent atan2 function is used, then these values are usually handled correctly. Another bounded application condition is that Finally, simple iterative techniques were used to solve the implicit equations in the direct and inverse methods; even though these are slow (and in the case of the inverse method it sometimes does not converge), they result in the least increase in code size. These are not decorative examples; they determine which carrier and evidence can fill the abstraction's roles.

Why no parent is asserted. Removing those specialist details does not currently yield one live catalog node that is a necessary genus for every instance. The entry is therefore approved as unparented rather than attached by topical resemblance. Its collapse evidence remains specific—The expressions were put in Horner (or nested) form, since this allows polynomials to be evaluated using only a single temporary register.—and future graph densification may discover a defensible relation only if it preserves that boundary.

This entry is a kind of Algorithm.

Approved unparented node. No current live node supplies a defensible necessary genus or structural prerequisite for Vincenty's formulae. The reviewed identity is: Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a). The accelerated suggestion was declined because topical or lexical similarity does not establish hierarchy; the node is admitted without a parent pending later graph densification.
Related reasoning operations. Evidence, representation, comparison, classification, transformation, or evaluation may participate in particular cases, but participation does not make any one of them a necessary parent of every instance.

Relationships to Other Abstractions¶

Current abstraction Vincenty's formulae Domain-specific

Parents (1) — more general patterns this builds on

Vincenty's formulae is a kind of Algorithm Prime

Vincenty's formulae is a domain-specific instance of algorithm under its frozen identity. The complete catalog already supplies this broader identity.

Hierarchy paths (2) — routes to 2 parentless roots

Vincenty's formulae → Algorithm → Function (Mapping)

Show alternative path (1)

Neighborhood in Abstraction Space¶

Vincenty's formulae sits in a moderately populated region (48^th percentile for distinctiveness): it has near-neighbors but no dense thicket of look-alikes.

Family — Geometric Figures & Constructions (32 abstractions)

Nearest neighbors

Albers Equal-Area Conic Projection — 0.88
Gnomonic Projection — 0.86
Oblate Spheroidal Coordinates — 0.86
Ellipse — 0.86
Laplace expansion (potential) — 0.86

Computed from structural-signature embeddings · 2026-10-08

Not to Be Confused With¶

Representation. The parent omits the specialist differentia. Tell: Can the case establish Vincenty's formulae are two related iterative methods used in geodesy to calculate the distance between two points on the surface of a spheroid, developed by Thaddeus Vincenty (1975a)?
Haversine Formula. A half-angle spherical-trigonometry relation that converts two latitude–longitude positions into their central angle and great-circle arc distance, with explicit radius, angle-unit, and numerical-boundary controls. Tell: Which entry's carrier, operation, and failure condition are satisfied?
Trilateration. Locate an unknown point by intersecting distance constraints from known reference points, using redundant ranges and an uncertainty model when real measurements do not meet at one exact solution. Tell: Which entry's carrier, operation, and failure condition are satisfied?
Orthogonal coordinates. A curvilinear coordinate system whose coordinate curves or hypersurfaces meet mutually at right angles, making the metric tensor diagonal in the coordinate basis. Tell: Which entry's carrier, operation, and failure condition are satisfied?
A measurement, proxy, or consequence. Those may provide evidence without being the identity. Tell: Would Vincenty's formulae remain present if the detector or downstream effect changed?
A metaphorical analogue. A similar shape outside cross-domain formal modeling lacks the specialist mechanism. Tell: Do the native roles transfer literally, or only the parent Representation?

References¶

Frozen Wikipedia discovery revision: https://en.wikipedia.org/wiki/Vincenty%27s_formulae (revision 1362852778).
Preserved source candidate: http://www.ngs.noaa.gov/TOOLS/Inv_Fwd/Inv_Fwd.html
Preserved source candidate: https://geographiclib.sourceforge.io/geodesic-papers/helmert80-en.html
Preserved source candidate: https://geographiclib.sourceforge.io/geod.html
Preserved source candidate: https://geographiclib.sourceforge.io/geod-addenda.html
Preserved source candidate: https://books.google.com/books?id=-d0EAAAAQAAJ&pg=PA130-IA4
Preserved source candidate: http://hdl.handle.net/1811/24409
Preserved source candidate: http://www.ngs.noaa.gov/PUBS_LIB/inverse.pdf
Preserved source candidate: https://geographiclib.sourceforge.io/geodesic-papers/vincenty75b.pdf

The frozen Wikipedia revision is discovery provenance. The retained source set was reviewed for identity, formal or operational relation, and scope. The encyclopedia's structural synthesis is bounded to those claims; a thin authority surface is recorded as a nonblocking source-strengthening repair rather than concealed.