I have compared some distance calculations between [`geopy`][1] and [PostGIS `ST_Distance_Sphere`][2], and it seems that PostGIS `ST_Distance` uses a calculation that is similar to the great circle distance rather than Vincenty.

Geopy Vincenty:

    >>> from geopy.distance import great_circle, vincenty
    >>> p1=(32.083659,34.806480)
    >>> p2=(32.0843373,34.8065365)
    >>> vincenty_distance = vincenty(p1, p2).meters
    >>> print vincenty_distance
    75.4044140365

Geopy Great Circle:

    >>> great_circle_distance = great_circle(p1, p2).meters
    >>> print great_circle_distance
    75.6324224394

PostGIS `ST_Distance_Sphere`:

    => SELECT ST_Distance_Sphere(
              ST_SetSRID(ST_MakePoint(34.806480, 32.083659), 4326),
              ST_SetSRID(ST_MakePoint(34.8065365, 32.0843373), 4326));
     st_distance_sphere
    --------------------
           75.611223453
    (1 row)


**Does PostGIS have a Vincenty Distance calculation?**


  [1]: http://gis.stackexchange.com/questions/84885/whats-the-difference-between-vincenty-and-great-circle-distance-calculations
  [2]: http://postgis.refractions.net/docs/ST_Distance_Sphere.html