We describe a new algorithm for computing Whitney stratifications of complex
projective varieties. The main ingredients are (a) an algebraic criterion, due
to L\^e and Teissier, which reformulates Whitney regularity in terms of
conormal spaces and maps, and (b) a new interpretation of this conormal
criterion via primary decomposition, which can be practically implemented on a
computer. We show that this algorithm improves upon the existing state of the
art by several orders of magnitude, even for relatively small input varieties.
En route, we introduce related algorithms for efficiently stratifying affine
varieties, flags on a given variety, and algebraic maps.Comment: There is an error in the published version of the article (Found
Comput Math, 2022) which has been fixed in this update. Section 3 is entirely
new, but the downstream results Sections 4-6 remain largely the same. We have
also updated the Runtimes and Complexity estimates in Section 7. The def. of
the integral closure of an ideal has also been correcte