We review the current status of general relativistic studies for coalescences of black hole--neutron star binaries. First, high-precision computations of black hole--neutron star binaries in quasiequilibrium circular orbits are summarized, focusing on the quasiequilibrium sequences and the mass-shedding limit. Next, the current status of numerical-relativity simulations for the merger of black hole--neutron star binaries is described. We summarize our understanding for the merger process, tidal disruption and its criterion, properties of the merger remnant and ejected material, gravitational waveforms, and gravitational-wave spectra. We also discuss expected electromagnetic counterparts to black hole--neutron star coalescences.