  We introduce GYRE, a new open-source stellar oscillation code which solves
the adiabatic/non-adiabatic pulsation equations using a novel Magnus Multiple
Shooting (MMS) numerical scheme. The code has a global error scaling of up to
6th order in the grid spacing, and can therefore achieve high accuracy with few
grid points. It is moreover robust and efficiently makes use of multiple
processor cores and/or nodes. We present an example calculation using GYRE, and
discuss recent work to integrate GYRE into the asteroseismic optimization
module of the MESA stellar evolution code.
