  Variation after particle-number restoration is incorporated for the first
time into the Hartree-Fock-Bogoliubov framework employing the Skyrme energy
density functional with zero-range pairing. The resulting projected HFB
equations can be expressed in terms of the local gauge-angle-dependent
densities. Results of projected calculations are compared with those obtained
within the Lipkin-Nogami method in the standard version and with the
Lipkin-Nogami method followed by exact particle-number projection.
