  We show how to compute the two-loop sunset integrals at finite volume, for
non-degenerate masses and non-zero momentum. We present results for all
integrals that appear in the Chiral Perturbation Therory ($\chi$PT) calculation
of the pseudoscalar meson masses and decay constants at NNLO, including the
case of Partially Quenched $\chi$PT. We also provide numerical implementations
of the finite-volume sunset integrals, and review the results for one-loop
integrals at finite volume.
