pdist2 — Compute pairwise distances between two sets of observations.

pdist2(X,Y) returns a size(X,1) by size(Y,1) matrix of distances between rows of X and rows of Y.

Syntax

D = pdist2(X, Y)
D = pdist2(X, Y, Distance, DistParameter, Name, Value)

Inputs

NameTypeRequiredDefaultDescription
XNumericArrayYesObservation matrix with observations in rows.
YNumericArrayYesSecond observation matrix with observations in rows.
optionsAnyVariadicDistance metric, metric parameter, or Smallest/Largest selection options.

Returns

NameTypeDescription
DNumericArrayPairwise distances.
DNumericArraySelected pairwise distances.
INumericArrayOne-based row indices from X for selected distances.

Returned values from pdist2 depend on how many outputs the caller requests.

Errors

IdentifierWhenMessage
RunMat:distance:InvalidArgumentInputs, dimensions, metrics, metric parameters, or selection options are malformed.distance helper: invalid argument
RunMat:distance:InternalRunMat cannot allocate or construct a distance output.distance helper: internal error

How pdist2 works

  • X and Y must be real numeric vectors or 2-D matrices with the same number of columns.
  • The default metric is Euclidean distance.
  • Supported metrics match pdist: "euclidean", "squaredeuclidean", "cityblock", "chebychev", "minkowski", "seuclidean", "mahalanobis", "cosine", "correlation", "hamming", "jaccard", and "spearman".
  • pdist2(X,Y,"minkowski",P) uses positive finite exponent P.
  • pdist2(X,Y,"seuclidean",S) accepts a standard-deviation vector and defaults to scale estimated from X. pdist2(X,Y,"mahalanobis",C) accepts a symmetric positive definite covariance matrix and defaults to covariance estimated from X.
  • "Smallest",K and "Largest",K return the K smallest or largest distances per row of Y. [D,I] = pdist2(...) also returns one-based row indices from X for the selected distances.

Examples

Compute distances between two sets

X = [0 0; 2 0];
Y = [1 0; 3 0];
D = pdist2(X,Y)

Expected output:

D is a 2-by-2 distance matrix.

Use squared Euclidean distance

D = pdist2(X,Y,"squaredeuclidean")

Expected output:

D contains squared distances.

Keep nearest distances per query row

[D,I] = pdist2(X,Y,"euclidean","Smallest",1)

Expected output:

D is a 1-by-size(Y,1) matrix of nearest distances; I contains row indices from X.

Using pdist2 with coding agents

Open a RunMat example with live inputs, then ask the agent to explain how pdist2 changes the result.

Run a small pdist2 example, explain the result, then change one input and compare the output.

FAQ

Does Smallest return indices?

Yes. Request [D,I] to receive one-based row indices from X, matching MATLAB's nearest-distance output form.

Are custom distance functions supported?

No. RunMat currently supports the standard named numeric distance metrics.

Open-source implementation

Unlike proprietary runtimes, every RunMat function is open-source. Read exactly how pdist2 is executed, line by line, in Rust.

About RunMat

RunMat is an open-source runtime that executes MATLAB-syntax code blazing on any GPU. It is licensed under the Apache 2.0 license.

  • RunMat automatically optimizes your math for GPU execution on Apple, Nvidia, and AMD hardware. No code changes needed. Simulations that took hours now take minutes.
  • Start running code in seconds. RunMat runs in the browser, on the desktop, or from the CLI. No license server, no IT ticket.

Getting started · Benchmarks · Pricing

Download RunMat

Download RunMat for full performance, or use RunMat in your browser for zero setup.