Automatic Generation of Finite Difference Field Solvers for Toroidally Confined Plasmas. Commun. Comput. Phys. 2020;28(3):1176-1188. doi:10.4208/cicp.OA-2017-0151