  We study the density of states (DOS) in diffusive superconductors with
pointlike magnetic impurities of arbitrary strength described by the Poissonian
statistics. The mean-field theory predicts a nontrivial structure of the DOS
with the continuum of quasiparticle states and (possibly) the impurity band. In
this approximation, all the spectral edges are hard, marking distinct
boundaries between spectral regions of finite and zero DOS. Considering
instantons in the replica sigma-model technique, we calculate the average DOS
beyond the mean-field level and determine the smearing of the spectral edges
due interplay of fluctuations of potential and nonpotential disorder. The
latter, represented by inhomogeneity in the concentration of magnetic
impurities, affects the subgap DOS in two ways: via fluctuations of the
pair-breaking strength and via induced fluctuations of the order parameter. In
limiting cases, we reproduce previously reported results for the subgap DOS in
disordered superconductors with strong magnetic impurities.
