Scalable approximation and solvers for ionic electrodiffusion in cellular geometries