On the Development and Implementation of High-Order Flux Reconstruction Schemes for Computational Fluid Dynamics