  We formulate and analyze an adaptive algorithm for isogeometric analysis with
hierarchical B-splines for weakly-singular boundary integral equations. We
prove that the employed weighted-residual error estimator is reliable and
converges at optimal algebraic rate. Numerical experiments with isogeometric
boundary elements for the 3D Poisson problem confirm the theoretical results,
which also cover general elliptic systems like linear elasticity.
