Here are described four solvers for time-harmonic electromagnetic fields in checkerboard patterns. A pattern is built by four squares with constant permittivity, ε 1 or ε 2 . It is enclosed by conducting walls or is a unit cell of a periodic structure. The field is represented in two ways: by B, the transverse component of the magnetic induction, and by A, the magnetic vector potential in Lorenz gauge. B and A satisfy Helmholtz equations in each square as well as transmission and boundary conditions (BCs). These governing equations yield eigensolutions B and A, which are found to be H 1 at worst. Variational versions of the governing equations are introduced. The weak formulations for B are standard, while those for A are new. They imply that the derivative transmission and BCs are satisfied weakly on interfaces between regions with different permittivity. Eigenpairs are computed approximately by spectral element methods. They yield mutually consistent eigenpairs. However, only about half of the eigenpairs (ω 2 , A) correspond to eigenpairs (ω 2 , B). For each set of BCs, the first few eigenfrequencies ω are given by tables, and some of the eigenfunctions are presented by contour plots.