We develop an efficient and robust numerical scheme to compute multi-fronts in one-dimensional Real Ginzburg-Landau equations that range from well-separated to strongly interacting and colliding. The scheme is based on the global centre-manifold reduction where one considers an initial sum of fronts plus a remainder function (not necessarily small) and applying a suitable projection based on the neutral Eigenmodes of each front. Such a scheme efficiently captures the weakly interacting tails of the fronts. Furthermore, as the fronts become strongly interacting, we show how they may be added to the remainder function to accurately compute through collisions. We then present results of our numerical scheme applied to various real Ginzburg Landau equations where we observe colliding fronts, travelling fronts and fronts converging to bound states. Finally, we discuss how this numerical scheme can be extended to general PDE systems and other multi-localised structures.