  New families of composition methods with processing of order 4 and 6 are
presented and analyzed. They are specifically designed to be used for the
numerical integration of differential equations whose vector field is separated
into three or more parts which are explicitly solvable. The new schemes are
shown to be more efficient than previous state-of-the-art splitting methods.
