  We present a number of new contributions to the topic of constructing
efficient higher-order splitting methods for the numerical integration of
evolution equations. Particular schemes are constructed via setup and solution
of polynomial systems for the splitting coefficients. To this end we use and
modify a recent approach for generating these systems for a large class of
splittings. In particular, various types of pairs of schemes intended for use
in adaptive integrators are constructed.
