  Orbital and spin ordering phenomena in strongly correlated systems are
commonly studied using the local-density approximation + dynamical mean-field
theory approach. Typically, however, such simulations are restricted to
simplified models (density-density Coulomb interactions, high symmetry
couplings and few-band models). In this work we implement an efficient general
hybridization-expansion continuous-time quantum Monte Carlo impurity solver
(Krylov approach) which allows us to investigate orbital and spin ordering in a
more realistic setting, including interactions that are often neglected (e.g.,
spin-flip and pair-hopping terms), enlarged basis sets (full d versus eg),
low-symmetry distortions, and reaching the very low-temperature (experimental)
regime. We use this solver to study ordering phenomena in a selection of
exemplary low-symmetry transition-metal oxides: LaMnO3 and rare-earth
manganites as well as the perovskites CaVO3 and YTiO3. We show that spin-flip
and pair hopping terms do not affect the Kugel-Khomskii orbital-order melting
transition in rare-earth manganites, or the suppression of orbital fluctuations
driven by crystal field and Coulomb repulsion. For the Mott insulator YTiO3 we
find a ferromagnetic transition temperature 50 K, in remarkably good agreement
with experiments. For LaMnO3 we show that the classical t2g-spin approximation,
commonly adopted for studying manganites, yields indeed an occupied eg orbital
in very good agreement with that obtained for the full d 5-orbital Hubbard
model, while the spin-spin e_g-t_{2g} correlation function calculated from the
full d model is 0.74, very close to the value expected for aligned eg and t2g
spins; the eg spectral function matrix is also well reproduced. Finally, we
show that the t2g screening reduces the eg-eg Coulomb repulsion by about 10%
