  By employing a novel generalization of the inverse scattering transform
method known as the unified transform or Fokas method, it can be shown that the
solution of certain physically significant boundary value problems for the
elliptic sine-Gordon equation, as well as for the elliptic version of the Ernst
equation, can be expressed in terms of the solution of appropriate $2 \times
2$-matrix Riemann--Hilbert (RH) problems. These RH problems are defined in
terms of certain functions, called spectral functions, which involve the given
boundary conditions, but also unknown boundary values. For arbitrary boundary
conditions, the determination of these unknown boundary values requires the
analysis of a nonlinear Fredholm integral equation. However, there exist
particular boundary conditions, called linearizable, for which it is possible
to bypass this nonlinear step and to characterize the spectral functions
directly in terms of the given boundary conditions. Here, we review the
implementation of this effective procedure for the following linearizable
boundary value problems: (a) the elliptic sine-Gordon equation in a semi-strip
with zero Dirichlet boundary values on the unbounded sides and with constant
Dirichlet boundary value on the bounded side; (b) the elliptic Ernst equation
with boundary conditions corresponding to a uniformly rotating disk of dust;
(c) the elliptic Ernst equation with boundary conditions corresponding to a
disk rotating uniformly around a central black hole; (d) the elliptic Ernst
equation with vanishing Neumann boundary values on a rotating disk.
