We apply high-order mixed finite element discretization techniques and their associated preconditioned iterative solvers to the Variable Eddington Factor (VEF) equations in two spatial dimensions. The VEF discretizations are coupled a Discontinuous Galerkin (DG) of discrete ordinates transport equation form effective linear algorithms that compatible with (curved) meshes. This combination is mo...