Scalable Matrix-free Solver for 3D Polarized Radiative Transfer in Stellar Atmospheres