  In recent years, an increasing amount of attention is being paid to the
gravitational few-body problem and its applications to astrophysical scenarios.
Among the main reasons for this renewed interest there is large number of newly
discovered exoplanets and the detection of gravitational waves. Here, we
present two numerical codes to model three- and few-body systems, called
TSUNAMI and OKINAMI. The TSUNAMI code is a direct few-body code with
algorithmic regularization, tidal forces and post-Newtonian corrections.
OKINAMI is a secular, double-averaged code for stable hierarchical triples. We
describe the main methods implemented in our codes, and review our recent
results and applications to gravitational-wave astronomy, planetary science and
statistical escape theories.
