To assist in the prototyping and controller design of point-absorber wave energy converters (WECs), an easy-to-implement hybrid integral-equation method is presented for computing the frequency-domain hydrodynamic properties of bodies with a vertical axis of symmetry in waves. The current hybrid method decomposes the flow domain into two parts: an inner domain containing the body and an outer domain extending to infinity. The solution in the inner domain is computed using the boundary-element method, and the outer-domain solution is expressed using eigenfunctions. Proper matching at the domain boundary is achieved by enforcing continuity of velocity potential and its normal derivative. Body symmetry allows efficient computation using ring sources in the inner domain. The current method is successfully applied to three different body geometries including a vertical truncated floating cylinder, the McIver toroid, and the coaxial-cylinder WEC being developed in the authors’ laboratory. In particular, the current results indicate that, by replacing the flat bottom of the coaxial-cylinder WEC with the Berkeley-Wedge (BW) shape, viscous effect can be significantly reduced with only minor negative impact on wave-exciting force, thus increasing WEC efficiency. Finally, by comparing to experimental measurements, the current method is demonstrated to accurately predict the heave added mass and wave-exciting force on the coaxial-cylinder WEC with BW geometry. If a viscous damping correction factor is used, the heave motion amplitude can also be accurately computed.