  We address the following question of partial desingularization preserving
normal crossings. Given an algebraic (or analytic) variety X in characteristic
zero, can we find a finite sequence of blowings-up preserving the
normal-crossings locus of X, after which the transform X' of X has only
singularities from an explicit finite list of minimal singularities, which we
define using the determinants of circulant matrices. In the case of surfaces,
for example, the pinch point or Whitney umbrella is the only singularity needed
in addition to normal crossings.
  We develop techniques for factorization (splitting) of a monic polynomial
with regular (or analytic) coefficients, satisfying a generic normal crossings
hypothesis, which we use together with resolution of singularities techniques
to find local circulant normal forms of singularities. These techniques in
their current state are enough for a positive answer to the question above, for
dim X up to 4, or in arbitrary dimension if we preserve normal crossings only
of order at most three. In these cases, minimal singularities have smooth
normalization.
