On the development and implementation of high-order flux reconstruction schemes for computational fluid dynamics