The major challenges for a fully polarized radiative transfer driven approach to Zeeman-Doppler imaging are still the enormous computational requirements. In every cycle of the iterative interplay between the forward process (spectral synthesis) and the inverse process (derivative based optimization) the Stokes profile synthesis requires several thousand evaluations of the polarized radiative transfer equation for a given stellar surface model. To cope with these computational demands and to allow for the incorporation of a full Stokes profile synthesis into Doppler- and Zeeman-Doppler imaging applications as well as into large scale solar Stokes profile