function xy = perform_geodesic_embedding(X,options)

% perform_geodesic_embedding - compute embedding using distance to seed points.
%
%   xy = perform_geodesic_embedding(X,options);
%
%   You can provide:
%       options.landmarks: the list of the numbers of the 
%           extremal points used to compute the embedding.
%           The dimension of the embedding is m=(nbr.landmarks)/2
%       options.dim: the embedding dimension. The dim*2 needed lamdmarks
%           are then computed automatically using farthest points.
%
%   Copyright (c) 2006 Gabriel Peyr